reset surface element linked list when deleting surface elements

libsrc/meshing/meshclass.cpp

@@ -137,6 +137,16 @@ namespace netgen
   }
 
 
+  void Mesh :: ClearSurfaceElements()
+  { 
+    surfelements.SetSize(0); 
+    for (int i = 0; i < facedecoding.Size(); i++)
+      facedecoding[i].firstelement = -1;
+
+    timestamp = NextTimeStamp();
+  }
+
+
 
   PointIndex Mesh :: AddPoint (const Point3d & p, int layer)
   { 

libsrc/meshing/meshclass.hpp

@@ -191,11 +191,7 @@ public:
   void DeleteMesh();
   
   ///
-  void ClearSurfaceElements()
-  { 
-    surfelements.SetSize(0); 
-    timestamp = NextTimeStamp();
-  }
+  void ClearSurfaceElements();
 
   ///
   void ClearVolumeElements()

libsrc/visualization/vsmesh.cpp

@@ -136,8 +136,6 @@ namespace netgen
 
     glPolygonMode (GL_FRONT_AND_BACK, GL_FILL);
 
-
-
     // glPolygonOffset (1,10);
     glPolygonOffset (2,2);
     glEnable (GL_POLYGON_OFFSET_FILL);
@@ -210,9 +208,6 @@ namespace netgen
     glDisable (GL_LINE_SMOOTH);
 
 
-
-
-
     if (vispar.drawoutline)
       {
 	glPolygonOffset (1, 1);
@@ -897,8 +892,6 @@ namespace netgen
     static int timer = NgProfiler::CreateTimer ("Mesh::BuildFilledList");
     NgProfiler::RegionTimer reg (timer);
 
-    // cout << "buildilled" << endl;
-
 #ifdef PARALLELGL
     cout << "buildfillelist, id = " << id << ", nse = " << mesh -> GetNSE() << endl;
 
@@ -1010,7 +1003,6 @@ namespace netgen
 	const FaceDescriptor & fd = mesh->GetFaceDescriptor(faceindex);
 	mesh->GetSurfaceElementsOfFace (faceindex, seia);
 
-	
 	// Philippose - 06/07/2009
 	// Modified the colour system to integrate the face colours into 
 	// the mesh data structure, rather than limit it to the OCC geometry 
@@ -1392,7 +1384,6 @@ namespace netgen
     linetimestamp = NextTimeStamp();
 
 #ifdef STLGEOM
-    xxx
     STLGeometry * stlgeometry = dynamic_cast<STLGeometry*> (ng_geometry);
     bool checkvicinity = (stlgeometry != NULL) && stldoctor.showvicinity;
 #endif