diff --git a/Geo/GRegion.cpp b/Geo/GRegion.cpp
index 021b14259ad43726e89e5a2a981aff85e5112b11..19a4916cf1ec47f80827d1130ddc46ffa0135d7e 100644
--- a/Geo/GRegion.cpp
+++ b/Geo/GRegion.cpp
@@ -310,7 +310,7 @@ double GRegion::computeSolidProperties (std::vector<double> cg,
   double surface = 0;
   cg[0] = cg[1] = cg[2] = 0.0;
   for ( ; it != l_faces.end(); ++it,++itdir){    
-    printf("face %d dir %d %d elements\n",(*it)->tag(),*itdir,(int)(*it)->triangles.size());
+    //printf("face %d dir %d %d elements\n",(*it)->tag(),*itdir,(int)(*it)->triangles.size());
     for (unsigned int i = 0; i < (*it)->triangles.size(); ++i){
       MTriangle *e = (*it)->triangles[i];
       int npt;
diff --git a/benchmarks/stl/skull.geo b/benchmarks/stl/skull.geo
index 8863d81c0f7cb049fa18a56a0e569ca2d07179d1..636cb452055c0d56fc0edf0bf5c09e6f1cfccc7c 100644
--- a/benchmarks/stl/skull.geo
+++ b/benchmarks/stl/skull.geo
@@ -1,5 +1,6 @@
 Mesh.RemeshParametrization=1; //(0) harmonic (1) conformal 
 Mesh.RemeshAlgorithm=1; //(0) nosplit (1) automatic (2) split metis
+
 Mesh.CharacteristicLengthFactor=0.1;
 Mesh.Algorithm3D = 4; //Frontal (4) Delaunay(1)