Use same tolerance in "Tet* walk (Tet *t, Vert *v, int maxx, double &totSearch,…

Use same tolerance in "Tet* walk (Tet *t, Vert *v, int maxx, double &totSearch, int thread)" as in "static Tet* tetContainsV (Vert *v, cavityContainer &cavity)" as now the initial tet from the walk function may be some flat tets away. When building up the cavity, flat tets between the initial tet and the tet containing the vertex may not include the vertex in their sphere, resulting in that they won't end up in the cavity and so the function "tetContainsV" will not find any tet. What will result in invalid cavities, unable to recover the boundary and finally, no mesh.

Use same tolerance in "Tet* walk (Tet t, Vert v, int maxx, double &totSearch,…
7cfaaa3e · Wendy Merks-Swolfs · 0ce61dd3 · 7cfaaa3e
Commit 7cfaaa3e authored Nov 30, 2017 by Wendy Merks-Swolfs
Hide whitespace changes
Inline Side-by-side

Showing with 1 addition and 1 deletion

delaunay3d.cpp Mesh/delaunay3d.cpp +1 -1

No files found.
--- a/Mesh/delaunay3d.cpp
+++ b/Mesh/delaunay3d.cpp
@@ -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;