geometry: added floating point tolerance back into interior sectioning

devitocodes · Jun 29, 2023 · 0635e0a · 0635e0a
1 parent 19c885b
commit 0635e0a
Showing 1 changed file with 3 additions and 1 deletion.
diff --git a/schism/geometry/geometry.py b/schism/geometry/geometry.py
@@ -236,8 +236,10 @@ def _get_interior_mask(self):
             _, i = kdtree.query(query, workers=-1)
 
             # Check nearest boundary point not in box
+            # Errs on the side of caution (cutoff slightly larger
+            # according to floating point precision)
             interior = np.any(np.abs(query - bp[i])
-                              > cutoff*np.array(spacing)[np.newaxis, :],
+                              > (cutoff+_feps)*np.array(spacing)[np.newaxis, :],
                               axis=1)
 
             # Fill interior mask away from boundary