diff --git a/Mesh/meshGRegion.cpp b/Mesh/meshGRegion.cpp
index d04b788cfa2152d5a2c6dc9399b3a4ba69318e43..a8e5902dc5556b30ebf8e3245035fd72a7b84d3a 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 cabedbfd06f5b47f1557b79547278828b09f9373..0536ceb23e8d408d101081bc06356f0f92fc5e38 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);