pp

Mesh/meshGRegionBoundaryRecovery.cpp

@@ -130,7 +130,7 @@ bool tetgenmesh::reconstructmesh(void *p)
 
   Msg::Debug("Points have been tetrahedralized");
 
-  { //transfernodes();
+  {
     point pointloop;
     REAL x, y, z;
     int i;
@@ -147,7 +147,8 @@ bool tetgenmesh::reconstructmesh(void *p)
 	xmin = xmax = x;
 	ymin = ymax = y;
 	zmin = zmax = z;
-      } else {
+      }
+      else {
 	xmin = (x < xmin) ? x : xmin;
 	xmax = (x > xmax) ? x : xmax;
 	ymin = (y < ymin) ? y : ymin;
@@ -163,7 +164,7 @@ bool tetgenmesh::reconstructmesh(void *p)
     z = zmax - zmin;
     longest = sqrt(x * x + y * y + z * z);
     if (longest == 0.0) {
-      Msg::Warning("Error:  The point set is trivial.\n");
+      Msg::Warning("The point set is trivial");
       return true;
     }
 
@@ -171,7 +172,7 @@ bool tetgenmesh::reconstructmesh(void *p)
     if (minedgelength == 0.0) {
       minedgelength = longest * b->epsilon;
     }
-  } // transfernodes();
+  }
 
   point *idx2verlist;
 
@@ -182,7 +183,7 @@ bool tetgenmesh::reconstructmesh(void *p)
     idx2verlist[0] = dummypoint; // Let 0th-entry be dummypoint.
   }
 
-  if (1) {
+  {
     // Index the vertices.
     for (unsigned int i = 0; i < _vertices.size(); i++){
       _vertices[i]->setIndex(i);
@@ -282,7 +283,8 @@ bool tetgenmesh::reconstructmesh(void *p)
 		    bondflag++;
 		  }
 		}
-	      } else {
+	      }
+              else {
 		bondflag++;
 	      }
 	      enextself(checktet);
@@ -293,7 +295,8 @@ bool tetgenmesh::reconstructmesh(void *p)
 	      // All three faces at d in 'checktet' have been connected.
 	      // It can be removed from the link.
 	      prevchktet.tet[8 + prevchktet.ver] = tptr;
-	    } else {
+	    }
+            else {
 	      // Bakup the previous tet in the link.
 	      prevchktet = checktet;
 	    }
@@ -495,11 +498,10 @@ bool tetgenmesh::reconstructmesh(void *p)
       outmesh2medit("dump2");
     }
 
-  } // meshsurface()
+  }
 
   delete [] idx2verlist;
 
-  ////////////////////////////////////////////////////////
   // Boundary recovery.
 
   clock_t t;
@@ -526,14 +528,9 @@ bool tetgenmesh::reconstructmesh(void *p)
 
   Msg::Debug("Statistics:\n");
   Msg::Debug("  Input points: %ld", _vertices.size());
-  //if (b->refine) {
-  //  printf("  Input tetrahedra: %d\n", in->numberoftetrahedra);
-  //}
   if (b->plc) {
     Msg::Debug("  Input facets: %ld", f_list.size());
     Msg::Debug("  Input segments: %ld", e_list.size());
-    //printf("  Input holes: %d\n", in->numberofholes);
-    //printf("  Input regions: %d\n", in->numberofregions);
   }
 
   tetnumber = tetrahedrons->items - hullsize;
@@ -541,7 +538,8 @@ bool tetgenmesh::reconstructmesh(void *p)
 
   if (b->weighted) { // -w option
     Msg::Debug(" Mesh points: %ld", points->items - nonregularcount);
-  } else {
+  }
+  else {
     Msg::Debug(" Mesh points: %ld", points->items);
   }
   Msg::Debug("  Mesh tetrahedra: %ld", tetnumber);
@@ -569,7 +567,8 @@ bool tetgenmesh::reconstructmesh(void *p)
     if (st_segref_count > 0l) {
       Msg::Debug("  Steiner points on segments:  %ld", st_segref_count);
     }
-  } else {
+  }
+  else {
     Msg::Debug("  Convex hull faces: %ld", hullsize);
     if (meshhulledges > 0l) {
       Msg::Debug("  Convex hull edges: %ld", meshhulledges);
@@ -583,10 +582,8 @@ bool tetgenmesh::reconstructmesh(void *p)
   if (0) {
     outmesh2medit("dump");
   }
-  ////////////////////////////////////////////////////////
 
   {
-  ////////////////////////////////////////////////////////
     // Write mesh into to GRegion.
 
     Msg::Info("--> Write to GRegion...");