@itwin/core-geometry 5.4.0-dev.6 → 5.5.0-dev.1

package/lib/cjs/curve/internalContexts/CurveCurveIntersectXY.js

@@ -25,11 +25,13 @@ const Arc3d_1 = require("../Arc3d");
 const CurveChainWithDistanceIndex_1 = require("../CurveChainWithDistanceIndex");
 const CurveCollection_1 = require("../CurveCollection");
 const CurveLocationDetail_1 = require("../CurveLocationDetail");
+const CurvePrimitive_1 = require("../CurvePrimitive");
 const LineSegment3d_1 = require("../LineSegment3d");
 const LineString3d_1 = require("../LineString3d");
 const Loop_1 = require("../Loop");
 const Path_1 = require("../Path");
 const ProxyCurve_1 = require("../ProxyCurve");
+const TransitionSpiral3d_1 = require("../spiral/TransitionSpiral3d");
 // cspell:word XYRR
 /**
  * Handler class for XY intersections between _geometryB and another geometry.
@@ -739,7 +741,7 @@ class CurveCurveIntersectXY extends GeometryHandler_1.RecurseToCurvesGeometryHan
         return undefined;
     }
     /** Detail computation for linestring intersecting linestring. */
-    computeLineStringLineString(lsA, lsB, reversed) {
+    computeLineStringLineString(lsA, extendA0, extendA1, lsB, extendB0, extendB1, reversed) {
         const pointA0 = CurveCurveIntersectXY._workPointAA0;
         const pointA1 = CurveCurveIntersectXY._workPointAA1;
         const pointB0 = CurveCurveIntersectXY._workPointBB0;
@@ -763,7 +765,7 @@ class CurveCurveIntersectXY extends GeometryHandler_1.RecurseToCurvesGeometryHan
                 for (let ib = 1; ib < numB; ib++, pointB0.setFrom(pointB1), fB0 = fB1) {
                     lsB.pointAt(ib, pointB1);
                     fB1 = ib * dfB;
-                    this.dispatchSegmentSegment(lsA, ia === 1 && this._extendA0, pointA0, fA0, pointA1, fA1, (ia + 1) === numA && this._extendA1, lsB, ib === 1 && this._extendB0, pointB0, fB0, pointB1, fB1, (ib + 1) === numB && this._extendB1, reversed);
+                    this.dispatchSegmentSegment(lsA, ia === 1 && extendA0, pointA0, fA0, pointA1, fA1, (ia + 1) === numA && extendA1, lsB, ib === 1 && extendB0, pointB0, fB0, pointB1, fB1, (ib + 1) === numB && extendB1, reversed);
                 }
             }
         }
@@ -862,6 +864,14 @@ class CurveCurveIntersectXY extends GeometryHandler_1.RecurseToCurvesGeometryHan
         else if (this._geometryB instanceof BSplineCurve_1.BSplineCurve3d) {
             this.dispatchSegmentBsplineCurve(segmentA, this._extendA0, segmentA.point0Ref, 0.0, segmentA.point1Ref, 1.0, this._extendA1, this._geometryB, this._extendB0, this._extendB1, false);
         }
+        else if (this._geometryB instanceof TransitionSpiral3d_1.TransitionSpiral3d) {
+            const spiralApproximation = LineString3d_1.LineString3d.create();
+            this._geometryB.emitStrokes(spiralApproximation);
+            const numPreviousResults = this._results.length;
+            this.computeSegmentLineString(segmentA, this._extendA0, this._extendA1, spiralApproximation, false, false, false);
+            const numberOfNewResults = this._results.length - numPreviousResults;
+            this.refineSpiralResultsByNewton(segmentA, this._geometryB, numberOfNewResults);
+        }
         else if (this._geometryB instanceof CurveCollection_1.CurveCollection) {
             this.dispatchCurveCollection(segmentA, this.handleLineSegment3d.bind(this));
         }
@@ -874,7 +884,7 @@ class CurveCurveIntersectXY extends GeometryHandler_1.RecurseToCurvesGeometryHan
     handleLineString3d(lsA) {
         if (this._geometryB instanceof LineString3d_1.LineString3d) {
             const lsB = this._geometryB;
-            this.computeLineStringLineString(lsA, lsB, false);
+            this.computeLineStringLineString(lsA, this._extendA0, this._extendA1, lsB, this._extendB0, this._extendB1, false);
         }
         else if (this._geometryB instanceof LineSegment3d_1.LineSegment3d) {
             this.computeSegmentLineString(this._geometryB, this._extendB0, this._extendB1, lsA, this._extendA0, this._extendA1, true);
@@ -885,6 +895,14 @@ class CurveCurveIntersectXY extends GeometryHandler_1.RecurseToCurvesGeometryHan
         else if (this._geometryB instanceof BSplineCurve_1.BSplineCurve3d) {
             this.dispatchLineStringBSplineCurve(lsA, this._extendA0, this._extendA1, this._geometryB, this._extendB0, this._extendB1, false);
         }
+        else if (this._geometryB instanceof TransitionSpiral3d_1.TransitionSpiral3d) {
+            const spiralApproximation = LineString3d_1.LineString3d.create();
+            this._geometryB.emitStrokes(spiralApproximation);
+            const numPreviousResults = this._results.length;
+            this.computeLineStringLineString(lsA, this._extendA0, this._extendA1, spiralApproximation, false, false, false);
+            const numberOfNewResults = this._results.length - numPreviousResults;
+            this.refineSpiralResultsByNewton(lsA, this._geometryB, numberOfNewResults);
+        }
         else if (this._geometryB instanceof CurveCollection_1.CurveCollection) {
             this.dispatchCurveCollection(lsA, this.handleLineString3d.bind(this));
         }
@@ -907,6 +925,14 @@ class CurveCurveIntersectXY extends GeometryHandler_1.RecurseToCurvesGeometryHan
         else if (this._geometryB instanceof BSplineCurve_1.BSplineCurve3d) {
             this.dispatchArcBsplineCurve3d(arc0, this._extendA0, this._extendA1, this._geometryB, this._extendB0, this._extendB1, false);
         }
+        else if (this._geometryB instanceof TransitionSpiral3d_1.TransitionSpiral3d) {
+            const spiralApproximation = LineString3d_1.LineString3d.create();
+            this._geometryB.emitStrokes(spiralApproximation);
+            const numPreviousResults = this._results.length;
+            this.computeArcLineString(arc0, this._extendA0, this._extendA1, spiralApproximation, false, false, false);
+            const numberOfNewResults = this._results.length - numPreviousResults;
+            this.refineSpiralResultsByNewton(arc0, this._geometryB, numberOfNewResults);
+        }
         else if (this._geometryB instanceof CurveCollection_1.CurveCollection) {
             this.dispatchCurveCollection(arc0, this.handleArc3d.bind(this));
         }
@@ -929,6 +955,14 @@ class CurveCurveIntersectXY extends GeometryHandler_1.RecurseToCurvesGeometryHan
         else if (this._geometryB instanceof BSplineCurve_1.BSplineCurve3dBase) {
             this.dispatchBSplineCurve3dBSplineCurve3d(curve, this._geometryB, false);
         }
+        else if (this._geometryB instanceof TransitionSpiral3d_1.TransitionSpiral3d) {
+            const spiralApproximation = LineString3d_1.LineString3d.create();
+            this._geometryB.emitStrokes(spiralApproximation);
+            const numPreviousResults = this._results.length;
+            this.dispatchLineStringBSplineCurve(spiralApproximation, false, false, curve, this._extendA0, this._extendA1, true);
+            const numberOfNewResults = this._results.length - numPreviousResults;
+            this.refineSpiralResultsByNewton(curve, this._geometryB, numberOfNewResults);
+        }
         else if (this._geometryB instanceof CurveCollection_1.CurveCollection) {
             this.dispatchCurveCollection(curve, this.handleBSplineCurve3d.bind(this));
         }
@@ -937,6 +971,47 @@ class CurveCurveIntersectXY extends GeometryHandler_1.RecurseToCurvesGeometryHan
         }
         return undefined;
     }
+    /**
+     * Assuming the tail of `this._results` holds the given number of approximate solutions to the curve-spiral
+     * xy-intersection problem, replace each with its Newton refinement, unless it doesn't converge, in which case
+     * it is removed.
+     * @param spiral The transition spiral.
+     * @param otherCurve The other curve primitive.
+     * @param numberOfNewResults The number of results in the tail of `this._results` to be refined.
+     * @param reversed Whether the spiral is geometryA (true) or geometryB (false).
+     */
+    refineSpiralResultsByNewton(otherCurve, spiral, numberOfNewResults, reversed = false) {
+        const resultsToBeRefined = this._results.slice(this._results.length - numberOfNewResults);
+        this._results.length -= numberOfNewResults; // keep already refined results
+        for (const detail of resultsToBeRefined) {
+            let spiralFraction = reversed ? detail.detailA.fraction : detail.detailB.fraction;
+            let otherFraction = reversed ? detail.detailB.fraction : detail.detailA.fraction;
+            const xyMatchingFunction = new Newton_1.CurveCurveIntersectionXYRRToRRD(spiral, otherCurve);
+            const newtonSearcher = new Newton_1.Newton2dUnboundedWithDerivative(xyMatchingFunction);
+            newtonSearcher.setUV(spiralFraction, otherFraction);
+            if (newtonSearcher.runIterations()) {
+                spiralFraction = newtonSearcher.getU();
+                otherFraction = newtonSearcher.getV();
+                if (this.acceptFraction(false, spiralFraction, false) && this.acceptFraction(false, otherFraction, false))
+                    this.recordPointWithLocalFractions(otherFraction, otherCurve, 0, 1, spiralFraction, spiral, 0, 1, reversed);
+            }
+        }
+    }
+    /** Double dispatch handler for strongly typed spiral curve. */
+    handleTransitionSpiral(spiral) {
+        if (this._geometryB instanceof CurvePrimitive_1.CurvePrimitive) { // this also handles CurveChainWithDistanceIndex
+            const spiralApproximation = LineString3d_1.LineString3d.create();
+            spiral.emitStrokes(spiralApproximation);
+            const numPreviousResults = this._results.length;
+            this.handleLineString3d(spiralApproximation);
+            const numberOfNewResults = this._results.length - numPreviousResults;
+            this.refineSpiralResultsByNewton(this._geometryB, spiral, numberOfNewResults, true);
+        }
+        else if (this._geometryB instanceof CurveCollection_1.CurveCollection) {
+            this.dispatchCurveCollection(spiral, this.handleTransitionSpiral.bind(this));
+        }
+        return undefined;
+    }
     /** Double dispatch handler for strongly typed CurveChainWithDistanceIndex. */
     handleCurveChainWithDistanceIndex(chain) {
         super.handleCurveChainWithDistanceIndex(chain);