Commit 059a27d7 by Christophe Geuzaine

Merge branch 'Fix_invalid_cavities_delaunay_3d' into 'master'

Fix invalid cavities delaunay 3d

See merge request !62
parents e1b338cd 7cfaaa3e
Pipeline #931 passed with stage
in 58 minutes 54 seconds
......@@ -980,7 +980,7 @@ Tet* walk (Tet *t, Vert *v, int maxx, double &totSearch, int thread)
(double*)f.V[1],
(double*)f.V[2],
(double*)v);
if (val >=-1.e-09) count++;
if (val >=0.0) count++;
if (val < _min){
NEIGH = iNeigh;
_min = val;
......@@ -1340,7 +1340,7 @@ void delaunayTrgl (const unsigned int numThreads,
}
}
if (invalidCavities[0])Msg::Warning("%d invalid cavities",invalidCavities[0]);
if (invalidCavities[0])Msg::Error("%d invalid cavities",invalidCavities[0]);
#if defined(_VERBOSE)
printf("average searches per point %12.5E\n",totSearchGlob/Npts);
......
......@@ -158,6 +158,7 @@ bool tetgenmesh::reconstructmesh(void *p)
std::vector<MTetrahedron*> tets;
delaunayMeshIn3D(_vertices, tets, false);
if (Msg::GetErrorCount()) return false;
Msg::Debug("Points have been tetrahedralized");
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment