From f97de58c7e55127e270338b5c8b7af43fd2d0f80 Mon Sep 17 00:00:00 2001
From: Bruno <Bruno.Levy@inria.fr>
Date: Mon, 23 Sep 2024 22:47:32 +0200
Subject: [PATCH] Snap coordinates to floating points before tetrahedralizing
 if pre-processing and non-zero tolerance are set.

---
 src/lib/geogram/mesh/mesh_tetrahedralize.cpp | 10 ++++++++++
 1 file changed, 10 insertions(+)

diff --git a/src/lib/geogram/mesh/mesh_tetrahedralize.cpp b/src/lib/geogram/mesh/mesh_tetrahedralize.cpp
index 6096cb7a682e..97ee29361f12 100644
--- 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);