- Standardize data types in numeric operations across Bézier and NURBS methods for consistency and precision.

- Add `examples/ssx/nurbs_nurbs_intersection_10.py` to demonstrate NURBS intersection use cases.
- Refactor `ndinterval.c` to fix method names and refine interval operations.
- Introduce `_deflate.py` for advanced intersection processing with improved derivative handling and multi-dimensional Bernstein patch methods.
- Optimize memory usage by caching frequently used computations in `fdm.py` and `aabb.py`.

Author: sth-v

mmcore/geom/implicit/sdf.py

@@ -790,4 +790,4 @@ def value_grad_at(
     res=sdf_intersection_deepest_point(t1, t2, t1.bounds(),t2.bounds(),batch_size=256)
     print(res)
 
-    print(time.time()-s)
+    print(time.time()-s)

mmcore/numeric/intersection/_deflate.py

@@ -1359,94 +1359,5 @@ def trace_one_direction(sign):
     print("γ rows:", gamma_rows)
     print("points:", len(curve_3d))
     print("first/last:", curve_3d[0], curve_3d[-1])
-    s1 = bern_to_nurbs_bezier(ts1,rational=False)
-    s2 = bern_to_nurbs_bezier(ts2,rational=False)
-    try:
-
-        from dataclasses import dataclass, field
-
-
-        @dataclass
-        class PointMaterial:
-            color: tuple[float, float, float, float] = field(default=(1.0, 1.0, 1.0, 1.0))
-            size: int = 8
-
-
-        @dataclass
-        class ControlNetMaterial:
-            color: tuple[float, float, float, float] = field(default=(1.0, 1.0, 1.0, 1.0))
-            control_point_material: PointMaterial = field(default_factory=PointMaterial)
-
-
-        @dataclass
-        class CurveMaterial:
-            color: tuple[float, float, float, float] = field(default=(1.0, 1.0, 1.0, 1.0))
-            show_control_net: bool = field(default=False)
-            control_net_material: ControlNetMaterial = field(default_factory=ControlNetMaterial)
-
-
-        @dataclass
-        class WiresMaterial(CurveMaterial):
-
-            u_count: int = field(default=0)
-            v_count: int = field(default=0)
-
-
-        @dataclass
-        class SurfaceMaterial:
-
-            color: tuple[float, float, float, float] = field(default=(0.5, 0.5, 0.9, 0.05))
-            wires_material: WiresMaterial = field(
-                default_factory=lambda: WiresMaterial((1.0, 1.0, 1.0, 1.0), show_control_net=False, u_count=0,
-                                                      v_count=0)
-            )
-            show_control_net: bool = field(default=False)
-            control_net_material: ControlNetMaterial = field(default_factory=ControlNetMaterial)
-
-
-        # Defaults
-        ssx_point_material = PointMaterial(color=(0.0, 1.0, 0.5, 1.0), size=12)
-        ssx_branch_material: CurveMaterial = CurveMaterial((0.0, 1.0, 0.5, 1.0),
-                                                           show_control_net=True,
-                                                           control_net_material=ControlNetMaterial((0.0, 1.0, 0.5, 0.7),
-
-                                                                                                   control_point_material=PointMaterial(
-                                                                                                       (0.0, 1.0, 0.5,
-                                                                                                        0.4), size=8)
-                                                                                                   )
-                                                           )
-
-        surface_material: SurfaceMaterial = SurfaceMaterial()
-
-        from mmcore.extras.renderer.renderer3d import Viewer,OrbitCamera
-        from mmcore.geom.bvh.lbvh import AABB
-        bb = AABB.from_points(s1.control_points.reshape(-1, 3)).merge(AABB.from_points(s2.control_points.reshape(-1, 3)))
-        viewer = Viewer(camera=OrbitCamera(target=bb.centroid(), distance=np.linalg.norm(bb.diag()) * 2, near=1.0))
-
-        viewer.add_nurbs_surface(s1, color=surface_material.wires_material.color, surface_color=surface_material.color,
-                                 u_count=surface_material.wires_material.u_count,
-                                 v_count=surface_material.wires_material.v_count)
-        viewer.add_nurbs_surface(
-            s2, color=surface_material.wires_material.color, surface_color=surface_material.color,
-            u_count=surface_material.wires_material.u_count, v_count=surface_material.wires_material.v_count
-
-        )
-        from mmcore.geom._nurbs_interp import interpolate_curve
-        from mmcore.geom._nurbs_eval import NURBSCurveTuple
-
-        control_points,kv=interpolate_curve(curve_3d,min(3,len(curve_3d)-1))
-        branch=NURBSCurveTuple(order=min(3,len(curve_3d)-1)+1,knot=np.array(kv),control_points=np.array(control_points),weights= np.ones(len(control_points)) )
-        viewer.add_nurbs_curve(branch, color=ssx_branch_material.color)
-        if ssx_branch_material.show_control_net:
-                for p in branch.control_points:
-                    viewer.add_point3d(p,
-                                       color=ssx_branch_material.control_net_material.control_point_material.color,
-                                       size_px=ssx_branch_material.control_net_material.control_point_material.size)
-        viewer.run()
-    except Exception as e:
-        traceback.print_exc()
-
 
-    print(s1)
-    print(s2)
     print(np.array(curve_3d).tolist())