diff --git a/Mesh/directions3D.cpp b/Mesh/directions3D.cpp
index ea27ca3179b8e6227f8a0e41006cd43ca9322082..5c17379fa9679e5490a58dd671bfac18967690bd 100644
--- a/Mesh/directions3D.cpp
+++ b/Mesh/directions3D.cpp
@@ -508,14 +508,6 @@ void Frame_field::initFace(GFace* gf)
   // compute cumulative cross-data "vertices x elements" for the whole contour of gf
   std::list<GEdge*> edges = gf->edges();
   vertex_to_elements.clear();
-  std::set<GEdge*> mySet;
-  std::list<GEdge*>::iterator it = edges.begin();
-  while(it != edges.end()){
-    mySet.insert(*it);
-    ++it;
-  }
-  edges.clear();
-  edges.insert(edges.begin(), mySet.begin(), mySet.end());
   for( std::list<GEdge*>::const_iterator it=edges.begin(); it!=edges.end(); it++){
     build_vertex_to_elements(*it,false);
   }
diff --git a/Mesh/meshGFaceBamg.cpp b/Mesh/meshGFaceBamg.cpp
index 250a0a53e08c528bd659dd396a080c93e29b2b11..9803db2da32546864c56af4e9f65c1a09c23e427 100644
--- a/Mesh/meshGFaceBamg.cpp
+++ b/Mesh/meshGFaceBamg.cpp
@@ -70,21 +70,11 @@ static void computeMeshMetricsForBamg(GFace *gf, int numV,
 }
 
 
-void meshGFaceBamg(GFace *gf){
-
-   std::list<GEdge*> edges = gf->edges();
-
-  //Replace edges by their compounds
-   std::set<GEdge*> mySet;
-   std::list<GEdge*>::iterator it = edges.begin();
-   while(it != edges.end()){
-     mySet.insert(*it);
-     ++it;
-   }
-   edges.clear();
-   edges.insert(edges.begin(), mySet.begin(), mySet.end());
-   std::set<MVertex*> bcVertex;
-   for (std::list<GEdge*>::iterator it = edges.begin(); it != edges.end(); it++){
+void meshGFaceBamg(GFace *gf)
+{
+  std::list<GEdge*> edges = gf->edges();
+  std::set<MVertex*> bcVertex;
+  for (std::list<GEdge*>::iterator it = edges.begin(); it != edges.end(); it++){
     for (unsigned int i = 0; i < (*it)->lines.size(); i++){
       bcVertex.insert((*it)->lines[i]->getVertex(0));
       bcVertex.insert((*it)->lines[i]->getVertex(1));