remove dead code

Mesh/BoundaryLayers.cpp

@@ -32,7 +32,7 @@ static double GetAveEdgeLength(std::vector<MVertex*> &elem_verts)
   if(!size)
     return 0.0;
   for( int i = 0; i < size-1; i++ )
-    ave += elem_verts[i]->distance(elem_verts[i+1]);    
+    ave += elem_verts[i]->distance(elem_verts[i+1]);
   ave += elem_verts[0]->distance(elem_verts[size-1]);
   ave /= size;
   return ave;
@@ -47,7 +47,7 @@ static void addExtrudeNormals(std::vector<T*> &elements, int invert,
     Msg::Error("Boundary layer index should be 0 or 1");
     return;
   }
- 
+
   if(octree && !gouraud){ // get extrusion direction from post-processing view
     // Trevor Strickler modified this section heavily
     std::set<MVertex*> verts;
@@ -117,10 +117,10 @@ static void addExtrudeNormals(std::vector<T*> &elements, int invert,
 typedef std::set<std::pair<bool, std::pair<int, int> > > infoset;
 
 // Trevor Strickler Modified this function
-//skipScaleCalcMap maps an entity tag to a flag telling whether to skip the 
+//skipScaleCalcMap maps an entity tag to a flag telling whether to skip the
 // scale calc when extruding only that entity.  The flag is false when an extrusion
 // is not scaleLast when in a boundary layer that has at least one scaleLast region.
-// Effectively, this makes the vertices on the boundary between a scaled and not 
+// Effectively, this makes the vertices on the boundary between a scaled and not
 // scaled region 'average' between being scaled and not scaled.
 template<class T>
 static void addExtrudeNormals(std::set<T*> &entities,
@@ -245,7 +245,7 @@ static void checkDepends(GModel *m, GFace *f, std::set<GFace*> &dep)
     dep.insert(from);
     checkDepends(m, from, dep);
   }
-  
+
   // Added by Trevor Strickler for compound face extrusion
   if( f->geomType() == GEntity::CompoundSurface ){
     std::list<GFace*> compounds = ((GFaceCompound*)(f))->getCompounds();
@@ -259,18 +259,19 @@ static void checkDepends(GModel *m, GFace *f, std::set<GFace*> &dep)
       checkDepends(m, *itgf, dep);
     }
   }
-  
+
 }
 
 // Trevor Strickler
-static unsigned int FixErasedExtrScaleFlags(GModel *m, std::map<int, bool> &faceSkipScaleCalc, std::map<int, bool> &edgeSkipScaleCalc)
+static unsigned int FixErasedExtrScaleFlags(GModel *m, std::map<int, bool> &faceSkipScaleCalc,
+                                            std::map<int, bool> &edgeSkipScaleCalc)
 {
   unsigned int num_changed = 0;
   std::set<GRegion *, GEntityLessThan>::iterator itreg;
   // fix all extruded faces bordering ScaleLast regions
   for( itreg = m->firstRegion(); itreg != m->lastRegion(); itreg++ ){
     ExtrudeParams *r_ep = (*itreg)->meshAttributes.extrude;
-    if(!r_ep || !r_ep->mesh.ExtrudeMesh || r_ep->geo.Mode != EXTRUDED_ENTITY 
+    if(!r_ep || !r_ep->mesh.ExtrudeMesh || r_ep->geo.Mode != EXTRUDED_ENTITY
         || !r_ep->mesh.ScaleLast )
       continue;
     std::list<GFace *> reg_faces = (*itreg)->faces();
@@ -308,7 +309,7 @@ static unsigned int FixErasedExtrScaleFlags(GModel *m, std::map<int, bool> &face
       }
     }
   }
-  
+
   return num_changed;
 }
 
@@ -320,7 +321,7 @@ int Mesh2DWithBoundaryLayers(GModel *m)
   std::map<int, bool> faceSkipScaleCalc, edgeSkipScaleCalc; // Trevor Strickler
   ExtrudeParams::calcLayerScaleFactor[0] = 0; // Trevor Strickler
   ExtrudeParams::calcLayerScaleFactor[1] = 0; // Trevor Strickler
-  
+
   // 2D boundary layers
   for(GModel::eiter it = m->firstEdge(); it != m->lastEdge(); it++){
     GEdge *ge = *it;
@@ -467,7 +468,7 @@ int Mesh2DWithBoundaryLayers(GModel *m)
           vdest = ge->getEndVertex();
         }
         GPoint p = vsrc->point();
-	
+
         ep->Extrude(ep->mesh.NbLayer - 1, ep->mesh.NbElmLayer[ep->mesh.NbLayer - 1],
                     p.x(), p.y(), p.z());
         vdest->setPosition(p);
@@ -517,105 +518,3 @@ int Mesh2DWithBoundaryLayers(GModel *m)
   return 1;
 }
 
-// give it a try : add one quad layer on the
-  /*
-void addOneLayerOnContour(GFace *gf, GVertex *gv){
-, int nbLayers, double hplus, double factor){
-  // for each vertex
-  std::map<MVertex*,std::vector<MVertex*> >layers;
-  std::vector<MQuadrangle*> newQuads;
-  std::vector<MTriangle*> newTris;
-
-  std::list<GEdgeLoop>::iterator it = gf->edgeLoops.begin();
-  for (; it != gf->edgeLoops.end(); ++it){
-    bool found = false;
-    std::list<GEdge*> ed;
-    for (GEdgeLoop::iter it2 = it->begin(); it2 != it->end(); ++it2){
-      if (it2->ge->getBeginVertex() == gv || it2->ge->getEndVertex() == gv) {
-	found = true;
-      }
-      ed.push_back(it2->ge);
-    }
-    // we found an edge loop with the GVertex that was specified
-    if (found){
-      // compute model vertices that will produce fans
-      for (GEdgeLoop::iter it2 = it->begin(); it2 != it->end(); ++it2){
-	GEdgeLoop::iter it3 = it2; ++it3;
-	GVertex *gv = it2->getEndVertex();
-	GEdgeSigned *before,*after = *it2;
-	if (it3 == it->end()){
-	  before = *(it->begin());
-	}
-	else{
-	  before = *it2;
-	}
-      }
-
-      for (std::list<GEdge*>::iterator it = ed.begin(); it != ed.end(); ++it){
-	GEdge *ge = *it;
-	for (int i=0;i<ge->lines.size();i++){
-	  SPoint2 p[2];
-	  reparamMeshEdgeOnFace ( ge->lines[i]->getVertex(0), ge->lines[i]->getVertex(1),gf,p[0],p[1]);
-	  MVertex *vd[2];
-	  for (int j=0;j<2;j++){
-	    MVertex *v = ge->lines[i]->getVertex(j);
-	    std::map<MVertex*,MVertex*>::iterator itv = duplicates.find(v);
-	    if (itv == duplicates.end()){
-	      vd[j] = new MFaceVertex(v->x(),v->y(),v->z(),gf,p[j].x(),p[j].y());
-	      duplicates[v] = vd[j];
-	      gf->mesh_vertices.push_back(vd[j]);
-	    }
-	    else
-	      vd[j] = itv->second;
-	  }
-	  newQuads.push_back(new MQuadrangle(ge->lines[i]->getVertex(0), ge->lines[i]->getVertex(1),vd[1],vd[0]));
-	}
-      }
-      for (int i=0;i<gf->quadrangles.size();i++){
-	MQuadrangle *q = gf->quadrangles[i];
-	MVertex *vs[4];
-	for (int j=0;j<4;j++){
-	  MVertex *v = q->getVertex(j);
-	  std::map<MVertex*,MVertex*>::iterator itv = duplicates.find(v);
-	  if (itv == duplicates.end()){
-	    vs[j] = v;
-	  }
-	  else{
-	    vs[j] = itv->second;
-	  }
-	}
-	newQuads.push_back(new MQuadrangle(vs[0],vs[1],vs[2],vs[3]));
-	delete q;
-      }
-      for (int i=0;i<gf->triangles.size();i++){
-	MTriangle *t = gf->triangles[i];
-	MVertex *vs[3];
-	for (int j=0;j<3;j++){
-	  MVertex *v = t->getVertex(j);
-	  std::map<MVertex*,MVertex*>::iterator itv = duplicates.find(v);
-	  if (itv == duplicates.end()){
-	    vs[j] = v;
-	  }
-	  else{
-	    vs[j] = itv->second;
-	  }
-	}
-	newTris.push_back(new MTriangle(vs[0],vs[1],vs[2]));
-	delete t;
-      }
-
-      gf->triangles = newTris;
-      gf->quadrangles = newQuads;
-    }
-  }
-}
-*/
-
-void addBoundaryLayers(GFace *gf)
-{
-  //if (backgroundMesh::current()){
-  //}
-  // first compute the cross field if it is not computed yet
-  // start from a selection of edges and create points in the boundary layer
-  // connect everybody with delaunay 
-}