only one boundarytio remesh and refill holes with rbf

Geo/GModel.cpp

@@ -1508,7 +1508,7 @@ void GModel::createTopologyFromRegions(std::vector<discreteRegion*> &discRegions
   Msg::Debug("Done creating topology from regions");
 }
 
-void GModel::createTopologyFromFaces(std::vector<discreteFace*> &discFaces)
+void GModel::createTopologyFromFaces(std::vector<discreteFace*> &discFaces, int onlyOneBoundary)
 {
   Msg::Debug("Creating topology from faces...");
 
@@ -1597,6 +1597,11 @@ void GModel::createTopologyFromFaces(std::vector<discreteFace*> &discFaces)
     std::vector<std::vector<MEdge> > boundaries;
     int nbBounds = connectedSurfaceBoundaries(myEdges, boundaries);   
 
+    //EMI RBF fix
+    if (onlyOneBoundary){
+      nbBounds = 1;
+    }
+
     // create new discrete edges
     for (int ib = 0; ib < nbBounds; ib++){
       int numE = getMaxElementaryNumber(1) + 1;

Geo/GModel.h

@@ -351,7 +351,7 @@ class GModel
   // create topology from mesh
   void createTopologyFromMesh();
   void createTopologyFromRegions(std::vector<discreteRegion*> &discRegions);
-  void createTopologyFromFaces(std::vector<discreteFace*> &pFaces);
+  void createTopologyFromFaces(std::vector<discreteFace*> &pFaces, int onlyOneBoundary=0);
   void makeDiscreteRegionsSimplyConnected();
   void makeDiscreteFacesSimplyConnected();