Enhance _bez_csx4.py and _bez_ccx4.py with AABB pre-filters and parameter tolerances for boundary and isolated intersection detection. Refactor boolean2d.py for endpoint snapping accuracy. Add _nurbs_join.py to support NURBS curve joining operations.

Author: sth-v

examples/ssx/bez_ssx6_baseline.py

@@ -157,6 +157,7 @@ def max_residual(branches: list, S1: np.ndarray, S2: np.ndarray, rational: bool)
 
         for i in range(len(stuv)):
             s, t, u, v = stuv[i]
+
             p1 = eval_surface(S1, s, t, rational=rational)
             p2 = eval_surface(S2, u, v, rational=rational)
             d = float(np.linalg.norm(p1 - p2))