Snap coordinates to floating points before tetrahedralizing if

pre-processing and non-zero tolerance are set.
BrunoLevy · Sep 23, 2024 · f97de58 · f97de58
1 parent 32418e8
commit f97de58
Showing 1 changed file with 10 additions and 0 deletions.
diff --git a/src/lib/geogram/mesh/mesh_tetrahedralize.cpp b/src/lib/geogram/mesh/mesh_tetrahedralize.cpp
@@ -95,6 +95,16 @@ namespace GEO {
             }
 
             if(preprocess) {
+
+		// Snap to floating-point coords
+		if(epsilon > 0) {
+		    index_t N = M.vertices.nb() * M.vertices.dimension();
+		    double* coords = M.vertices.point_ptr(0);
+		    for(index_t i=0; i<N; ++i) {
+			coords[i] = double(float(coords[i]));
+		    }
+		}
+
                 mesh_repair(M, MESH_REPAIR_DEFAULT, epsilon);
                 fill_holes(M, max_hole_area);
                 MeshSurfaceIntersection intersection(M);