diff --git a/Mesh/BackgroundMesh.cpp b/Mesh/BackgroundMesh.cpp
index 79d9e2142a17dba9076655d7a0d1672078727379..22f997224501bb58eb027b6a1cbbfed5fb8946a8 100644
--- a/Mesh/BackgroundMesh.cpp
+++ b/Mesh/BackgroundMesh.cpp
@@ -111,7 +111,6 @@ static double max_edge_curvature(const GVertex *gv)
 
 static double max_surf_curvature(const GEdge *ge, double u)
 {
-  //printf("max surf \n");
   double val = 0;
   std::list<GFace *> faces = ge->faces();
   std::list<GFace *>::iterator it = faces.begin();
diff --git a/Mesh/meshGEdge.cpp b/Mesh/meshGEdge.cpp
index 700a7d7f7f1e79c41da2dd9e5b70f998c3b61d90..ae1c95c165f24cd96c9d4f11dd0a3256dca23184 100644
--- a/Mesh/meshGEdge.cpp
+++ b/Mesh/meshGEdge.cpp
@@ -22,7 +22,6 @@ typedef struct {
 
 static double F_Lc(GEdge *ge, double t)
 {
-  //printf("FLC \n");
   GPoint p = ge->point(t);
   double lc_here;
 
@@ -30,8 +29,6 @@ static double F_Lc(GEdge *ge, double t)
   double t_begin = bounds.low();
   double t_end = bounds.high();
 
-  //printf("FLC going to BGM ! \n");
- 
   if(t == t_begin)
     lc_here = BGM_MeshSize(ge->getBeginVertex(), t, 0, p.x(), p.y(), p.z());
   else if(t == t_end)
@@ -320,8 +317,6 @@ void meshGEdge::operator() (GEdge *ge)
                       CTX::instance()->mesh.lcIntegrationPrecision);
     N = std::max(ge->minimumMeshSegments() + 1, (int)(a + 1.));
   }
-  //printf("going to exit !!");
-  //exit(1);
 
   // force odd number of points for if blossom is used for
   // recombination