revert to bary

Author: Jammy2211

autoarray/structures/mesh/delaunay_2d.py

@@ -40,24 +40,27 @@ def scipy_delaunay(points_np, query_points_np, use_voronoi_areas, areas_factor):
 
     # ---------- Voronoi or Barycentric Areas used to weight split points ----------
 
-    # if use_voronoi_areas:
-    #
-    #     areas = voronoi_areas_numpy(
-    #         points,
-    #     )
-    #
-    #     max_area = np.percentile(areas, 90.0)
-    #
-    #     areas[areas == -1] = max_area
-    #     areas[areas > max_area] = max_area
-    #
-    # else:
+    if use_voronoi_areas:
 
-    areas = barycentric_dual_area_from(
-        points,
-        simplices,
-        xp=np,
-    )
+        try:
+            areas = voronoi_areas_numpy(points)
+        except Exception as e:
+            # Qhull precision problems -> fallback
+            print(f"[pure_callback] Voronoi failed ({type(e).__name__}); falling back to barycentric areas.")
+            areas = barycentric_dual_area_from(points, simplices, xp=np)
+
+        max_area = np.percentile(areas, 90.0)
+
+        areas[areas == -1] = max_area
+        areas[areas > max_area] = max_area
+
+    else:
+
+        areas = barycentric_dual_area_from(
+            points,
+            simplices,
+            xp=np,
+        )
 
     split_point_areas = areas_factor * np.sqrt(areas)