From 37790676a1f2281551250b7a4463c42699cc6741 Mon Sep 17 00:00:00 2001
From: Jean-Francois Remacle <jean-francois.remacle@uclouvain.be>
Date: Tue, 27 Sep 2016 19:02:09 +0000
Subject: [PATCH] debug hybrid meshes

---
 Mesh/meshGRegion.cpp               | 5 +++--
 Mesh/meshGRegionRelocateVertex.cpp | 3 ++-
 2 files changed, 5 insertions(+), 3 deletions(-)

diff --git a/Mesh/meshGRegion.cpp b/Mesh/meshGRegion.cpp
index d04b788cfa..a8e5902dc5 100644
--- a/Mesh/meshGRegion.cpp
+++ b/Mesh/meshGRegion.cpp
@@ -596,7 +596,6 @@ void MeshDelaunayVolumeTetgen(std::vector<GRegion*> &regions)
       return;
     }
     TransferTetgenMesh(gr, in, out, numberedV);
-    sqr.relocateVertices(gr,3);
   }
 
 
@@ -610,6 +609,8 @@ void MeshDelaunayVolumeTetgen(std::vector<GRegion*> &regions)
 
   // restore the initial set of faces
   gr->set(faces);
+  insertVerticesInRegion(gr,0,true);
+  sqr.relocateVertices(gr,3);
 
   
   // now do insertion of points
@@ -624,7 +625,7 @@ void MeshDelaunayVolumeTetgen(std::vector<GRegion*> &regions)
     int nbvertices_filler = (old_algo_hexa()) ?
       Filler::get_nbr_new_vertices() : Filler3D::get_nbr_new_vertices();
     if(!nbvertices_filler && !LpSmoother::get_nbr_interior_vertices()){
-      insertVerticesInRegion(gr,2000000000,true);
+      insertVerticesInRegion(gr,2000000000,false);
     }
   }
   // crete an initial mesh 
diff --git a/Mesh/meshGRegionRelocateVertex.cpp b/Mesh/meshGRegionRelocateVertex.cpp
index cabedbfd06..0536ceb23e 100644
--- a/Mesh/meshGRegionRelocateVertex.cpp
+++ b/Mesh/meshGRegionRelocateVertex.cpp
@@ -71,7 +71,8 @@ static int Stopping_Rule(double x0, double x1, double tol)
 
 double Maximize_Quality_Golden_Section( MVertex *ver, 
                                         double xTarget, double yTarget, double zTarget,
-                                        const std::vector<MElement*> &lt , double tol, double &q)
+                                        const std::vector<MElement*> &lt ,
+					double tol, double &q)
 {
   
   static const double lambda = 0.5 * (sqrt5 - 1.0);
-- 
GitLab