From f03117d819679ff46bebe3b23cb8595c54385da0 Mon Sep 17 00:00:00 2001
From: Emilie Marchandise <emilie.marchandise@uclouvain.be>
Date: Mon, 26 Mar 2012 14:53:54 +0000
Subject: [PATCH] fill holes with planar faces (create only triangles)
---
Geo/GFaceCompound.cpp | 3 +++
1 file changed, 3 insertions(+)
diff --git a/Geo/GFaceCompound.cpp b/Geo/GFaceCompound.cpp
index 51544b1814..a2fbfc862e 100644
--- a/Geo/GFaceCompound.cpp
+++ b/Geo/GFaceCompound.cpp
@@ -535,10 +535,13 @@ void GFaceCompound::fillNeumannBCS_Plane() const
GFace *newFace = GModel::current()->addPlanarFace(myEdgeLoops);
fillFaces.push_back(newFace);
int meshingAlgo = CTX::instance()->mesh.algo2d;
+ int recombine = CTX::instance()->mesh.recombineAll;
opt_mesh_algo2d(0, GMSH_SET, 1.0); //mesh adapt
+ opt_mesh_recombine_all(0, GMSH_SET, 0.0); //no recombination
meshGFace mgf;
mgf(newFace);
opt_mesh_algo2d(0, GMSH_SET, meshingAlgo);
+ opt_mesh_recombine_all(0, GMSH_SET, recombine);
for(unsigned int i = 0; i < newFace->triangles.size(); ++i){
loopfillTris.push_back(newFace->triangles[i]);
fillNodes.insert(newFace->triangles[i]->getVertex(0));
--
GitLab