From ca11b159564adead511de9b4ee523955b11fa7a2 Mon Sep 17 00:00:00 2001
From: Christophe Geuzaine <cgeuzaine@ulg.ac.be>
Date: Mon, 17 Nov 2008 16:31:51 +0000
Subject: [PATCH] fix bug with boundary layers introduced recently

---
 Mesh/HighOrder.cpp | 7 +++++++
 1 file changed, 7 insertions(+)

diff --git a/Mesh/HighOrder.cpp b/Mesh/HighOrder.cpp
index 661f4048ae..a3ec772d57 100644
--- a/Mesh/HighOrder.cpp
+++ b/Mesh/HighOrder.cpp
@@ -148,6 +148,7 @@ bool reparamOnFace(MVertex *v1, MVertex *v2, GFace *gf, SPoint2 &param1, SPoint2
   return false;
 }
 
+// FIXME: this should NOT be in HighOrder.cpp...
 bool reparamOnFace(MVertex *v, GFace *gf, SPoint2 &param)
 {
 
@@ -157,6 +158,12 @@ bool reparamOnFace(MVertex *v, GFace *gf, SPoint2 &param)
     return true;
   }
 
+  if(v->onWhat()->geomType() == GEntity::DiscreteCurve || 	 
+     v->onWhat()->geomType() == GEntity::BoundaryLayerCurve){ 	 
+    param = gf->parFromPoint(SPoint3(v->x(), v->y(), v->z()));
+    return true;
+  }
+
   if(v->onWhat()->dim() == 0){
     GVertex *gv = (GVertex*)v->onWhat();
 
-- 
GitLab