prevent negative areas, due to delaunay triangles at edge next to one another

Jammy2211 · Jammy2211 · commit 12d7ee74064b · 2026-02-11T17:06:30.000Z
diff --git a/autoarray/structures/mesh/delaunay_2d.py b/autoarray/structures/mesh/delaunay_2d.py
@@ -59,7 +59,9 @@ def scipy_delaunay(points_np, query_points_np, use_voronoi_areas, areas_factor):
             xp=np,
         )
 
+    areas = np.maximum(areas, 0.0)
     split_point_areas = areas_factor * np.sqrt(areas)
+    split_point_areas = np.where(np.isfinite(split_point_areas), split_point_areas, 0.0)
 
     # ---------- Compute split cross points for Split regularization ----------
     split_points = split_points_from(

Original file line number	Diff line number	Diff line change
`@@ -59,7 +59,9 @@ def scipy_delaunay(points_np, query_points_np, use_voronoi_areas, areas_factor):`
`59`	`59`	`xp=np,`
`60`	`60`	`)`
`61`	`61`
	`62`	`+ areas = np.maximum(areas, 0.0)`
`62`	`63`	`split_point_areas = areas_factor * np.sqrt(areas)`
	`64`	`+ split_point_areas = np.where(np.isfinite(split_point_areas), split_point_areas, 0.0)`
`63`	`65`
`64`	`66`	`# ---------- Compute split cross points for Split regularization ----------`
`65`	`67`	`split_points = split_points_from(`