From 97d2ef35c62bd1a467e67423fc98f3af5fb1081a Mon Sep 17 00:00:00 2001
From: Bruno Levy <Bruno.Levy@inria.fr>
Date: Mon, 30 Oct 2023 10:22:38 +0100
Subject: [PATCH] mesh_tetrahedralize: automatically triangulate if surface was
 not triangulated.

---
 src/lib/geogram/basic/memory.h                     | 9 +++++++++
 src/lib/geogram/mesh/mesh_surface_intersection.cpp | 2 +-
 src/lib/geogram/mesh/mesh_tetrahedralize.cpp       | 6 +++---
 3 files changed, 13 insertions(+), 4 deletions(-)

diff --git a/src/lib/geogram/basic/memory.h b/src/lib/geogram/basic/memory.h
index 28cdb163de0f..fb128bf141ae 100644
--- a/src/lib/geogram/basic/memory.h
+++ b/src/lib/geogram/basic/memory.h
@@ -65,6 +65,15 @@
 
 #endif
 
+// Stack size depending on OS:
+// Linux: 10 Mb
+// Windows: 1 Mb
+// Mac OSX: 523 Kb
+// GEO_HAS_BIG_STACK is defined under Linux
+// and lets some of the functions that
+// manipulate exact precision numbers
+// allocate temporaries on the stack.
+
 #ifdef GEO_OS_LINUX
 #define GEO_HAS_BIG_STACK
 #endif
diff --git a/src/lib/geogram/mesh/mesh_surface_intersection.cpp b/src/lib/geogram/mesh/mesh_surface_intersection.cpp
index d2210f700785..36f708241ec4 100644
--- a/src/lib/geogram/mesh/mesh_surface_intersection.cpp
+++ b/src/lib/geogram/mesh/mesh_surface_intersection.cpp
@@ -492,7 +492,7 @@ namespace GEO {
                 for(index_t i=b; i<e; ++i) {
                     const IsectInfo& II = intersections[i];
                     
-                    // Each IsectIfo is either an individual vertex
+                    // Each IsectInfo is either an individual vertex
                     // or a segment with two vertices.
                     // Each vertex is represented combinatorially.
                     // The MeshInTriangle knows how to compute the
diff --git a/src/lib/geogram/mesh/mesh_tetrahedralize.cpp b/src/lib/geogram/mesh/mesh_tetrahedralize.cpp
index 7187220d1a89..9c7e0088f08f 100644
--- a/src/lib/geogram/mesh/mesh_tetrahedralize.cpp
+++ b/src/lib/geogram/mesh/mesh_tetrahedralize.cpp
@@ -73,10 +73,10 @@ namespace GEO {
         }
 
         if(!M.facets.are_simplices()) {
-            Logger::err("TetMeshing")
-                << "Mesh is not triangulated"
+            Logger::warn("TetMeshing")
+                << "Mesh is not triangulated (triangulating it)"
                 << std::endl;
-            return false;
+            tessellate_facets(M,3);
         }
 
         // in percent of bbox diagonal