Skip to content
Snippets Groups Projects
Commit c75a2522 authored by Christophe Geuzaine's avatar Christophe Geuzaine
Browse files

new tetgen boundary recovery code

parent f46816e3
No related branches found
No related tags found
No related merge requests found
......@@ -947,12 +947,7 @@ StringXNumber GeometryOptions_Number[] = {
StringXNumber MeshOptions_Number[] = {
{ F|O, "Algorithm" , opt_mesh_algo2d , ALGO_2D_AUTO ,
"2D mesh algorithm (1=MeshAdapt, 2=Automatic, 5=Delaunay, 6=Frontal, 7=bamg, 8=delquad)" },
{ F|O, "Algorithm3D" , opt_mesh_algo3d ,
#if defined(HAVE_TETGEN)
ALGO_3D_DELAUNAY ,
#else
ALGO_3D_FRONTAL ,
#endif
{ F|O, "Algorithm3D" , opt_mesh_algo3d , ALGO_3D_DELAUNAY ,
"3D mesh algorithm (1=Delaunay, 4=Frontal, 5=Frontal Delaunay, 6=Frontal Hex, 7=MMG3D, 9=R-tree)" },
{ F|O, "AngleSmoothNormals" , opt_mesh_angle_smooth_normals , 30.0 ,
"Threshold angle below which normals are not smoothed" },
......
......@@ -592,9 +592,12 @@ void MeshDelaunayVolumeTetgen(std::vector<GRegion*> &regions)
}
// uncomment this to test the new code
//#define NEW_CODE
#if !defined(HAVE_TETGEN)
#define NEW_CODE
#endif
static void MeshDelaunayVolumeNewCode(std::vector<GRegion*> &regions) {
static void MeshDelaunayVolumeNewCode(std::vector<GRegion*> &regions)
{
GRegion *gr = regions[0];
std::list<GFace*> faces = gr->faces();
std::set<GFace*> allFacesSet;
......@@ -610,9 +613,7 @@ static void MeshDelaunayVolumeNewCode(std::vector<GRegion*> &regions) {
gr->set(allFaces);
try{
meshGRegionBoundaryRecovery *init = new meshGRegionBoundaryRecovery();
init->reconstructmesh(gr);
delete init;
meshGRegionBoundaryRecovery(gr);
}
catch(int err){
Msg::Error("Could not recover boundary: error %d", err);
......
This diff is collapsed.
This diff is collapsed.
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment