From 6b672fecce50ae15ebd6c765007ba94cf5730f8b Mon Sep 17 00:00:00 2001
From: Amaury Johnan <amjohnen@gmail.com>
Date: Mon, 7 Apr 2014 18:51:34 +0000
Subject: [PATCH] skip computation of Jacobian for HO pyramids until Bezier
subdivision is implemented.
---
Plugin/AnalyseCurvedMesh.cpp | 9 +++++++++
1 file changed, 9 insertions(+)
diff --git a/Plugin/AnalyseCurvedMesh.cpp b/Plugin/AnalyseCurvedMesh.cpp
index f4b7c2e2d9..f95a255c09 100644
--- a/Plugin/AnalyseCurvedMesh.cpp
+++ b/Plugin/AnalyseCurvedMesh.cpp
@@ -475,6 +475,14 @@ void GMSH_AnalyseCurvedMeshPlugin::computeMinMax(MElement *const*el, int numEl,
_min_ratioJ = 1.7e308;
for (int k = 0; k < numEl; ++k) {
+ if (el[k]->getType() == TYPE_PYR && el[k]->getPolynomialOrder() > 1) {
+ static int i = 0;
+ if (++i == 1) {
+ Msg::Error("High-order pyramids skipped (subdivision not implemented).");
+ Msg::Error("Subdivision will come soon :) (before may 2014)");
+ }
+ continue;
+ }
#ifdef _ANALYSECURVEDMESH_BLAS_
jacBez.setAsProxy(jacBezB, k);
@@ -573,6 +581,7 @@ void GMSH_AnalyseCurvedMeshPlugin::computeMinMax(MElement *const*el, int numEl,
}
while (pqMax.size() > 0) {
+ break;
bj = pqMax.top();
pqMax.pop();
delete bj;
--
GitLab