From 1b974678867825e5fd0046ae610691cbacd7ed38 Mon Sep 17 00:00:00 2001
From: Jonathan Lambrechts <jonathan.lambrechts@uclouvain.be>
Date: Wed, 1 Apr 2015 16:19:40 +0000
Subject: [PATCH] dg : fix(?)  bug in commit 20542 : fix gLevelsetDistMesh with
 angle weighted normal, (please check this commit)

---
 Geo/gmshLevelset.cpp | 18 ++++++++++--------
 1 file changed, 10 insertions(+), 8 deletions(-)

diff --git a/Geo/gmshLevelset.cpp b/Geo/gmshLevelset.cpp
index 78beae10b6..05926f882a 100644
--- a/Geo/gmshLevelset.cpp
+++ b/Geo/gmshLevelset.cpp
@@ -979,14 +979,16 @@ double gLevelsetDistMesh::operator()(double x, double y, double z) const
       Msg::Error("Cannot compute a distance to an entity of dimension %d\n",
                  (*it)->getDim());
     }
-    if(fabs(distance) == fabs(minDistance)){
-      closestElements.push_back(*it);
-    }
-    if(fabs(distance) < fabs(minDistance)){
-      closestPoint = closePt;
-      minDistance = distance;
-      closestElements.clear();
-      closestElements.push_back(*it);
+    if ((*it)->getDim() == 2) {
+      if(fabs(distance) == fabs(minDistance)){
+        closestElements.push_back(*it);
+      }
+      if(fabs(distance) < fabs(minDistance)){
+        closestPoint = closePt;
+        minDistance = distance;
+        closestElements.clear();
+        closestElements.push_back(*it);
+      }
     }
   }
   if(closestElements.size() > 1){
-- 
GitLab