From 727b02f7c823fcc12a630277ca1350d11fab8e59 Mon Sep 17 00:00:00 2001
From: Thomas Toulorge <thomas.toulorge@mines-paristech.fr>
Date: Mon, 7 Dec 2015 19:50:18 +0000
Subject: [PATCH] Modified fast curving to act also on concave boundaries
---
contrib/HighOrderMeshOptimizer/OptHomFastCurving.cpp | 2 +-
1 file changed, 1 insertion(+), 1 deletion(-)
diff --git a/contrib/HighOrderMeshOptimizer/OptHomFastCurving.cpp b/contrib/HighOrderMeshOptimizer/OptHomFastCurving.cpp
index 2a56ba6eb1..58b9df9bc2 100644
--- a/contrib/HighOrderMeshOptimizer/OptHomFastCurving.cpp
+++ b/contrib/HighOrderMeshOptimizer/OptHomFastCurving.cpp
@@ -410,7 +410,7 @@ inline void insertIfCurved(MElement *el, std::list<MElement*> &bndEl)
SPoint3 xyzS(0.,0.,0.);
for (int iSF = 0; iSF < nV1; ++iSF) xyzS += xyz1[iSF]*f[iSF]; // Compute location of node in straight element
const SVector3 vec(xyzS, el->getVertex(iV)->point());
- const double normalDisp = dot(vec, normal); // Normal component of displacement
+ const double normalDisp = fabs(dot(vec, normal)); // Normal component of displacement
if (normalDisp > normalDispCurved) {
bndEl.push_back(el);
break;
--
GitLab