if directions==0 bug fix

spasmann · Jan 15, 2024 · 5ace684 · 5ace684
1 parent f6e5e9f
commit 5ace684
Show file tree

Hide file tree

Showing 2 changed files with 3 additions and 1 deletion.
diff --git a/mcdc/constant.py b/mcdc/constant.py
@@ -34,7 +34,7 @@
 INF = 1e10
 PI = math.acos(-1.0)
 SHIFT = 1e-10  # To ensure lattice, surface, and mesh crossings
-PREC = 1.0 + 1e-8  # Precision factor to determine if a distance is smaller
+PREC = 1.0 + 1e-5  # Precision factor to determine if a distance is smaller
 BANKMAX = 100  # Default maximum active bank
 
 # Domain Decomp mesh crossing flags

diff --git a/mcdc/kernel.py b/mcdc/kernel.py
@@ -28,6 +28,8 @@ def domain_crossing(P, mcdc):
         mesh = mcdc["technique"]["domain_mesh"]
         # Determine which dimension is crossed
         x, y, z, t, directions = mesh_crossing_evaluate(P, mesh)
+        if len(directions) == 0:
+            return
         flag = directions[0]
 
         if len(directions) > 1: