Use better bump for safe floor evaluation

libgeos · Sep 10, 2023 · 30dc8a1 · 30dc8a1
1 parent da6215c
commit 30dc8a1
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/capi/geos_ts_c.cpp b/capi/geos_ts_c.cpp
@@ -2983,7 +2983,7 @@ extern "C" {
                 //-- check for an integral scale
                 double scale = 1.0 / gridSize;
                 //-- add a small "bump" to ensure flooring works
-                double scaleInt = std::floor(scale + 0.1);
+                double scaleInt = std::floor(scale + 2 * GRIDSIZE_INTEGER_TOLERANCE);
                 //-- if scale factor is essentially integral, use the exact integer value
                 if (std::abs(scale - scaleInt) < GRIDSIZE_INTEGER_TOLERANCE) {
                     scale = scaleInt;