fix: domain error at acos [1,-1]

aicostsqo-io · Mar 14, 2024 · 04a1b62 · 04a1b62
1 parent 7a5f895
commit 04a1b62
Showing 1 changed file with 8 additions and 1 deletion.
diff --git a/marble-sculp/disc.py b/marble-sculp/disc.py
@@ -90,7 +90,14 @@ def baecher(
                 -math.sin(echis),
             ]
             fisher_dip_dev_angle = math.acos(
-                (fisher_constant + math.log(1 - np.random.uniform())) / fisher_constant
+                max(
+                    min(
+                        (fisher_constant + math.log(1 - np.random.uniform()))
+                        / fisher_constant,
+                        1,
+                    ),
+                    -1,
+                )
             )
             echis = (math.pi / 2) - (dip - fisher_dip_dev_angle)
             pole_dip_rotated = [