split points?

Author: Jammy2211

autoarray/structures/mesh/delaunay_2d.py

@@ -46,28 +46,6 @@ def scipy_delaunay(points_np, query_points_np, use_voronoi_areas, areas_factor):
             points,
         )
 
-        # --- DEBUG / SAFETY CHECK: Voronoi areas ---
-        if not np.isfinite(areas).all():
-            n_nan = np.isnan(areas).sum()
-            n_inf = np.isinf(areas).sum()
-            print(
-                f"[pure_callback] Voronoi areas NON-FINITE: "
-                f"n_nan={n_nan} n_inf={n_inf} "
-                f"min={np.nanmin(areas)} max={np.nanmax(areas)}"
-            )
-
-            # Save everything needed to reproduce offline
-            np.savez(
-                "callback_bad_voronoi_areas.npz",
-                points=points,
-                areas=areas,
-            )
-
-            raise FloatingPointError(
-                "voronoi_areas_numpy produced NaN/inf; "
-                "saved callback_bad_voronoi_areas.npz"
-            )
-
         max_area = np.percentile(areas, 90.0)
 
         areas[areas == -1] = max_area
@@ -89,6 +67,28 @@ def scipy_delaunay(points_np, query_points_np, use_voronoi_areas, areas_factor):
         area_weights=split_point_areas,
     )
 
+    # --- DEBUG / SAFETY CHECK: Voronoi areas ---
+    if not np.isfinite(split_points).all():
+        n_nan = np.isnan(split_points).sum()
+        n_inf = np.isinf(split_points).sum()
+        print(
+            f"[pure_callback] Voronoi split_points NON-FINITE: "
+            f"n_nan={n_nan} n_inf={n_inf} "
+            f"min={np.nanmin(split_points)} max={np.nanmax(split_points)}"
+        )
+
+        # Save everything needed to reproduce offline
+        np.savez(
+            "callback_bad_voronoi_split_points.npz",
+            points=points,
+            split_points=split_points,
+        )
+
+        raise FloatingPointError(
+            "voronoi_split_points_numpy produced NaN/inf; "
+            "saved callback_bad_voronoi_split_points.npz"
+        )
+
     # ---------- find_simplex for split cross points ----------
     split_points_idx = tri.find_simplex(split_points)