@itwin/core-geometry 4.0.0-dev.99 → 4.1.0-dev.0

package/lib/cjs/geometry3d/CoincidentGeometryOps.js

@@ -5,6 +5,10 @@
 *--------------------------------------------------------------------------------------------*/
 Object.defineProperty(exports, "__esModule", { value: true });
 exports.CoincidentGeometryQuery = void 0;
+/** @packageDocumentation
+ * @module CartesianGeometry
+ */
+const core_bentley_1 = require("@itwin/core-bentley");
 const CurveLocationDetail_1 = require("../curve/CurveLocationDetail");
 const Geometry_1 = require("../Geometry");
 const AngleSweep_1 = require("./AngleSweep");
@@ -62,71 +66,71 @@ class CoincidentGeometryQuery {
      * @param pointB1 end point of segment B
      */
     coincidentSegmentRangeXY(pointA0, pointA1, pointB0, pointB1, restrictToBounds = true) {
-        const detailAOnB = this.projectPointToSegmentXY(pointA0, pointB0, pointB1);
-        if (pointA0.distanceXY(detailAOnB.point) > this._tolerance)
+        const detailA0OnB = this.projectPointToSegmentXY(pointA0, pointB0, pointB1);
+        if (pointA0.distanceXY(detailA0OnB.point) > this._tolerance)
             return undefined;
         const detailA1OnB = this.projectPointToSegmentXY(pointA1, pointB0, pointB1);
         if (pointA1.distanceXY(detailA1OnB.point) > this._tolerance)
             return undefined;
-        const detailBOnA = this.projectPointToSegmentXY(pointB0, pointA0, pointA1);
-        if (pointB0.distanceXY(detailBOnA.point) > this._tolerance)
+        const detailB0OnA = this.projectPointToSegmentXY(pointB0, pointA0, pointA1);
+        if (pointB0.distanceXY(detailB0OnA.point) > this._tolerance)
             return undefined;
         const detailB1OnA = this.projectPointToSegmentXY(pointB1, pointA0, pointA1);
         if (pointB1.distanceXY(detailB1OnA.point) > this._tolerance)
             return undefined;
-        detailAOnB.fraction1 = detailA1OnB.fraction;
-        detailAOnB.point1 = detailA1OnB.point; // capture -- detailB0OnA is not reused.
-        detailBOnA.fraction1 = detailB1OnA.fraction;
-        detailBOnA.point1 = detailB1OnA.point;
+        detailA0OnB.fraction1 = detailA1OnB.fraction;
+        detailA0OnB.point1 = detailA1OnB.point; // capture -- detailA1OnB is not reused.
+        detailB0OnA.fraction1 = detailB1OnA.fraction;
+        detailB0OnA.point1 = detailB1OnA.point;
         if (!restrictToBounds)
-            return CurveLocationDetail_1.CurveLocationDetailPair.createCapture(detailBOnA, detailAOnB);
-        const segment = Segment1d_1.Segment1d.create(detailBOnA.fraction, detailBOnA.fraction1);
+            return CurveLocationDetail_1.CurveLocationDetailPair.createCapture(detailB0OnA, detailA0OnB);
+        const segment = Segment1d_1.Segment1d.create(detailB0OnA.fraction, detailB0OnA.fraction1);
         if (segment.clampDirectedTo01()) {
-            segment.reverseIfNeededForDeltaSign(1.0);
             const f0 = segment.x0;
             const f1 = segment.x1;
-            const h0 = detailBOnA.inverseInterpolateFraction(f0);
-            const h1 = detailBOnA.inverseInterpolateFraction(f1);
+            const h0 = detailB0OnA.inverseInterpolateFraction(f0);
+            const h1 = detailB0OnA.inverseInterpolateFraction(f1);
             // recompute fractions and points..
-            CoincidentGeometryQuery.assignDetailInterpolatedFractionsAndPoints(detailBOnA, f0, f1, pointA0, pointA1);
-            CoincidentGeometryQuery.assignDetailInterpolatedFractionsAndPoints(detailAOnB, h0, h1, pointB0, pointB1);
-            return CurveLocationDetail_1.CurveLocationDetailPair.createCapture(detailBOnA, detailAOnB);
+            CoincidentGeometryQuery.assignDetailInterpolatedFractionsAndPoints(detailB0OnA, f0, f1, pointA0, pointA1, f0 > f1);
+            CoincidentGeometryQuery.assignDetailInterpolatedFractionsAndPoints(detailA0OnB, h0, h1, pointB0, pointB1, h0 > h1);
+            return CurveLocationDetail_1.CurveLocationDetailPair.createCapture(detailB0OnA, detailA0OnB);
         }
         else {
             if (segment.signedDelta() < 0.0) {
-                if (detailBOnA.point.isAlmostEqual(pointA0)) {
-                    detailBOnA.collapseToStart();
-                    detailAOnB.collapseToStart();
-                    return CurveLocationDetail_1.CurveLocationDetailPair.createCapture(detailBOnA, detailAOnB);
+                if (detailB0OnA.point.isAlmostEqual(pointA0, this.tolerance)) {
+                    detailB0OnA.collapseToStart();
+                    detailA0OnB.collapseToStart();
+                    return CurveLocationDetail_1.CurveLocationDetailPair.createCapture(detailB0OnA, detailA0OnB);
                 }
-                if (detailBOnA.point1.isAlmostEqual(pointA1)) {
-                    detailBOnA.collapseToEnd();
-                    detailAOnB.collapseToEnd();
-                    return CurveLocationDetail_1.CurveLocationDetailPair.createCapture(detailBOnA, detailAOnB);
+                if (detailB0OnA.point1.isAlmostEqual(pointA1, this.tolerance)) {
+                    detailB0OnA.collapseToEnd();
+                    detailA0OnB.collapseToEnd();
+                    return CurveLocationDetail_1.CurveLocationDetailPair.createCapture(detailB0OnA, detailA0OnB);
                 }
             }
             else {
-                if (detailBOnA.point.isAlmostEqual(pointA1)) {
-                    detailBOnA.collapseToStart();
-                    detailAOnB.collapseToEnd();
-                    return CurveLocationDetail_1.CurveLocationDetailPair.createCapture(detailBOnA, detailAOnB);
+                if (detailB0OnA.point.isAlmostEqual(pointA1, this.tolerance)) {
+                    detailB0OnA.collapseToStart();
+                    detailA0OnB.collapseToEnd();
+                    return CurveLocationDetail_1.CurveLocationDetailPair.createCapture(detailB0OnA, detailA0OnB);
                 }
-                if (detailBOnA.point1.isAlmostEqual(pointA0)) {
-                    detailBOnA.collapseToEnd();
-                    detailAOnB.collapseToStart();
-                    return CurveLocationDetail_1.CurveLocationDetailPair.createCapture(detailBOnA, detailAOnB);
+                if (detailB0OnA.point1.isAlmostEqual(pointA0, this.tolerance)) {
+                    detailB0OnA.collapseToEnd();
+                    detailA0OnB.collapseToStart();
+                    return CurveLocationDetail_1.CurveLocationDetailPair.createCapture(detailB0OnA, detailA0OnB);
                 }
             }
         }
         return undefined;
     }
     /**
-     * Create a CurveLocationDetailPair from . . .
+     * Create a CurveLocationDetailPair for a coincident interval of two overlapping curves
      * @param cpA curveA
-     * @param cpB curve B
-     * @param fractionsOnA fractions of an active section of curveA
-     * @param fractionB0 fraction of an original containing B interval
-     * @param fractionB1 end fraction of an original containing B interval
+     * @param cpB curveB
+     * @param fractionsOnA coincident interval of curveB in fraction space of curveA
+     * @param fractionB0 curveB start in fraction space of curveA
+     * @param fractionB1 curveB end in fraction space of curveA
+     * @param reverse whether curveB and curveA have opposite direction
      */
     createDetailPair(cpA, cpB, fractionsOnA, fractionB0, fractionB1, reverse) {
         const deltaB = fractionB1 - fractionB0;
@@ -135,10 +139,8 @@ class CoincidentGeometryQuery {
         if (g0 !== undefined && g1 !== undefined) {
             const detailA = CurveLocationDetail_1.CurveLocationDetail.createCurveEvaluatedFractionFraction(cpA, fractionsOnA.x0, fractionsOnA.x1);
             const detailB = CurveLocationDetail_1.CurveLocationDetail.createCurveEvaluatedFractionFraction(cpB, g0, g1);
-            if (reverse) {
+            if (reverse)
                 detailA.swapFractionsAndPoints();
-                detailB.swapFractionsAndPoints();
-            }
             return CurveLocationDetail_1.CurveLocationDetailPair.createCapture(detailA, detailB);
         }
         return undefined;
@@ -156,11 +158,11 @@ class CoincidentGeometryQuery {
      * @param arcA
      * @param arcB
      * @param _restrictToBounds
-     * @return 0, 1, or 2 overlap intervals.
+     * @return 0, 1, or 2 overlap points/intervals
      */
     coincidentArcIntersectionXY(arcA, arcB, _restrictToBounds = true) {
         let result;
-        if (arcA.center.isAlmostEqual(arcB.center)) {
+        if (arcA.center.isAlmostEqual(arcB.center, this.tolerance)) {
             const matrixBToA = arcA.matrixRef.multiplyMatrixInverseMatrix(arcB.matrixRef);
             if (matrixBToA) {
                 const ux = matrixBToA.at(0, 0);
@@ -174,25 +176,49 @@ class CoincidentGeometryQuery {
                 if (Geometry_1.Geometry.isAlmostEqualNumber(ru, 1.0)
                     && Geometry_1.Geometry.isAlmostEqualNumber(rv, 1.0)
                     && Geometry_1.Geometry.isAlmostEqualNumber(0, dot)) {
-                    const alphaB0Radians = Math.atan2(uy, ux); // angular position of arcB 0 point in A sweep
-                    const sweepDirection = cross > 0 ? 1.0 : -1.0; // 1 if arcB's parameter space sweeps forward, -1 if reverse
-                    const betaStartRadians = alphaB0Radians + sweepDirection * arcB.sweep.startRadians;
-                    const betaEndRadians = alphaB0Radians + sweepDirection * arcB.sweep.endRadians;
+                    const alphaB0Radians = Math.atan2(uy, ux); // angular position of arcB 0 point in arcA sweep
+                    const sweepDirection = cross > 0 ? 1.0 : -1.0; // 1 if arcB parameter space sweeps in same direction as arcA, -1 if opposite
+                    const betaStartRadians = alphaB0Radians + sweepDirection * arcB.sweep.startRadians; // arcB start in arcA parameter space
+                    const betaEndRadians = alphaB0Radians + sweepDirection * arcB.sweep.endRadians; // arcB end in arcA parameter space
                     const fractionSpacesReversed = (sweepDirection * arcA.sweep.sweepRadians * arcB.sweep.sweepRadians) < 0;
                     const sweepB = AngleSweep_1.AngleSweep.createStartEndRadians(betaStartRadians, betaEndRadians);
                     const sweepA = arcA.sweep;
                     const fractionPeriodA = sweepA.fractionPeriod();
-                    const fractionB0 = sweepA.radiansToPositivePeriodicFraction(sweepB.startRadians);
-                    const fractionSweep = sweepB.sweepRadians / sweepA.sweepRadians;
-                    const fractionB1 = fractionB0 + fractionSweep;
+                    const fractionB0 = sweepA.radiansToPositivePeriodicFraction(sweepB.startRadians); // arcB start in arcA fraction space
+                    (0, core_bentley_1.assert)(fractionB0 >= 0.0);
+                    const fractionSweep = sweepB.sweepRadians / sweepA.sweepRadians; // arcB sweep in arcA fraction space
+                    const fractionB1 = fractionB0 + fractionSweep; // arcB end in arcA fraction space
                     const fractionSweepB = Segment1d_1.Segment1d.create(fractionB0, fractionB1);
-                    if (fractionSweepB.clampDirectedTo01())
-                        result = this.appendDetailPair(result, this.createDetailPair(arcA, arcB, fractionSweepB, fractionB0, fractionB1, fractionSpacesReversed));
-                    if (fractionB1 > fractionPeriodA) {
-                        const fractionSweepBWrap = Segment1d_1.Segment1d.create(fractionB0 - fractionPeriodA, fractionB1 - fractionPeriodA);
-                        if (fractionSweepBWrap.clampDirectedTo01())
-                            result = this.appendDetailPair(result, this.createDetailPair(arcA, arcB, fractionSweepBWrap, fractionB0, fractionB1, fractionSpacesReversed));
-                    }
+                    /** lambda to add a coincident interval or isolated intersection, given inputs in arcA fraction space
+                     * @param arcBInArcAFractionSpace span of arcB in arcA fraction space. On return, clamped to [0,1] if nontrivial.
+                     * @param testStartOfArcA if no nontrivial coincident interval was found, look for an isolated intersection at the start (true) or end (false) of arcA
+                     * @returns whether a detail pair was appended to result
+                     */
+                    const appendCoincidentIntersection = (arcBInArcAFractionSpace, testStartOfArcA) => {
+                        const size = result ? result.length : 0;
+                        const arcBStart = arcBInArcAFractionSpace.x0;
+                        const arcBEnd = arcBInArcAFractionSpace.x1;
+                        if (arcBInArcAFractionSpace.clampDirectedTo01() && !Geometry_1.Geometry.isSmallRelative(arcBInArcAFractionSpace.absoluteDelta())) {
+                            result = this.appendDetailPair(result, this.createDetailPair(arcA, arcB, arcBInArcAFractionSpace, arcBStart, arcBEnd, fractionSpacesReversed));
+                        }
+                        else { // test isolated intersection
+                            const testStartOfArcB = fractionSpacesReversed ? testStartOfArcA : !testStartOfArcA;
+                            const arcAPt = this._point0 = testStartOfArcA ? arcA.startPoint(this._point0) : arcA.endPoint(this._point0);
+                            const arcBPt = this._point1 = testStartOfArcB ? arcB.startPoint(this._point1) : arcB.endPoint(this._point1);
+                            if (arcAPt.isAlmostEqual(arcBPt, this.tolerance)) {
+                                const detailA = CurveLocationDetail_1.CurveLocationDetail.createCurveFractionPoint(arcA, testStartOfArcA ? 0 : 1, arcAPt);
+                                const detailB = CurveLocationDetail_1.CurveLocationDetail.createCurveFractionPoint(arcB, testStartOfArcB ? 0 : 1, arcBPt);
+                                result = this.appendDetailPair(result, CurveLocationDetail_1.CurveLocationDetailPair.createCapture(detailA, detailB));
+                            }
+                        }
+                        return result !== undefined && result.length > size;
+                    };
+                    appendCoincidentIntersection(fractionSweepB, false); // compute overlap in strict interior, or at end of arcA
+                    // check overlap at start of arcA with a periodic shift of fractionSweepB
+                    if (fractionB1 >= fractionPeriodA)
+                        appendCoincidentIntersection(Segment1d_1.Segment1d.create(fractionB0 - fractionPeriodA, fractionB1 - fractionPeriodA), true);
+                    else if (fractionB0 === 0.0)
+                        appendCoincidentIntersection(Segment1d_1.Segment1d.create(fractionB0 + fractionPeriodA, fractionB1 + fractionPeriodA), true);
                 }
             }
         }

package/lib/cjs/geometry3d/CoincidentGeometryOps.js.map

@@ -1 +1 @@
-{"version":3,"file":"CoincidentGeometryOps.js","sourceRoot":"","sources":["../../../src/geometry3d/CoincidentGeometryOps.ts"],"names":[],"mappings":";AAAA;;;+FAG+F;;;AAM/F,sEAA4F;AAE5F,0CAAuC;AACvC,6CAA0C;AAC1C,uDAAsD;AACtD,2CAAwC;AAExC;;;;;GAKG;AACH,MAAa,uBAAuB;IAIlC,IAAW,SAAS;QAClB,OAAO,IAAI,CAAC,UAAU,CAAC;IACzB,CAAC;IACD,YAAoB,YAAoB,mBAAQ,CAAC,mBAAmB;QAClE,IAAI,CAAC,UAAU,GAAG,SAAS,CAAC;IAC9B,CAAC;IACM,MAAM,CAAC,MAAM,CAAC,YAAoB,mBAAQ,CAAC,mBAAmB;QACnE,OAAO,IAAI,uBAAuB,CAAC,SAAS,CAAC,CAAC;IAChD,CAAC;IACD;;;OAGG;IACI,MAAM,CAAC,0CAA0C,CAAC,MAA2B,EAAE,EAAU,EAAE,EAAU,EAC1G,MAAe,EAAE,MAAe,EAAE,OAAgB,KAAK;QACvD,IAAI,IAAI,EAAE;YACR,MAAM,CAAC,QAAQ,GAAG,EAAE,CAAC;YACrB,MAAM,CAAC,SAAS,GAAG,EAAE,CAAC;SACvB;aAAM;YACL,MAAM,CAAC,QAAQ,GAAG,EAAE,CAAC;YACrB,MAAM,CAAC,SAAS,GAAG,EAAE,CAAC;SACvB;QACD,MAAM,CAAC,KAAK,GAAG,MAAM,CAAC,WAAW,CAAC,MAAM,CAAC,QAAQ,EAAE,MAAM,EAAE,MAAM,CAAC,KAAK,CAAC,CAAC;QACzE,MAAM,CAAC,MAAM,GAAG,MAAM,CAAC,WAAW,CAAC,MAAM,CAAC,SAAS,EAAE,MAAM,EAAE,MAAM,CAAC,MAAM,CAAC,CAAC;IAE9E,CAAC;IAED;;OAEG;IACI,uBAAuB,CAAC,UAAmB,EAAE,MAAe,EAAE,MAAe;QAClF,IAAI,CAAC,QAAQ,GAAG,0BAAQ,CAAC,cAAc,CAAC,MAAM,EAAE,MAAM,EAAE,IAAI,CAAC,QAAQ,CAAC,CAAC;QACvE,IAAI,CAAC,QAAQ,GAAG,0BAAQ,CAAC,cAAc,CAAC,MAAM,EAAE,UAAU,EAAE,IAAI,CAAC,QAAQ,CAAC,CAAC;QAC3E,MAAM,KAAK,GAAG,IAAI,CAAC,QAAQ,CAAC,YAAY,CAAC,IAAI,CAAC,QAAQ,CAAC,CAAC;QACxD,MAAM,KAAK,GAAG,IAAI,CAAC,QAAQ,CAAC,YAAY,CAAC,IAAI,CAAC,QAAQ,CAAC,CAAC;QACxD,MAAM,QAAQ,GAAG,mBAAQ,CAAC,kBAAkB,CAAC,KAAK,EAAE,KAAK,EAAE,GAAG,CAAC,CAAC;QAChE,OAAO,yCAAmB,CAAC,wBAAwB,CAAC,SAAS,EAAE,QAAQ,EACrE,MAAM,CAAC,WAAW,CAAC,QAAQ,EAAE,MAAM,CAAC,CAAC,CAAC;IAC1C,CAAC;IACD;;;;;;;OAOG;IACI,wBAAwB,CAAC,OAAgB,EAAE,OAAgB,EAAE,OAAgB,EAAE,OAAgB,EAAE,mBAA4B,IAAI;QACtI,MAAM,UAAU,GAAG,IAAI,CAAC,uBAAuB,CAAC,OAAO,EAAE,OAAO,EAAE,OAAO,CAAC,CAAC;QAC3E,IAAI,OAAO,CAAC,UAAU,CAAC,UAAU,CAAC,KAAK,CAAC,GAAG,IAAI,CAAC,UAAU;YACxD,OAAO,SAAS,CAAC;QACnB,MAAM,WAAW,GAAG,IAAI,CAAC,uBAAuB,CAAC,OAAO,EAAE,OAAO,EAAE,OAAO,CAAC,CAAC;QAC5E,IAAI,OAAO,CAAC,UAAU,CAAC,WAAW,CAAC,KAAK,CAAC,GAAG,IAAI,CAAC,UAAU;YACzD,OAAO,SAAS,CAAC;QAEnB,MAAM,UAAU,GAAG,IAAI,CAAC,uBAAuB,CAAC,OAAO,EAAE,OAAO,EAAE,OAAO,CAAC,CAAC;QAC3E,IAAI,OAAO,CAAC,UAAU,CAAC,UAAU,CAAC,KAAK,CAAC,GAAG,IAAI,CAAC,UAAU;YACxD,OAAO,SAAS,CAAC;QACnB,MAAM,WAAW,GAAG,IAAI,CAAC,uBAAuB,CAAC,OAAO,EAAE,OAAO,EAAE,OAAO,CAAC,CAAC;QAC5E,IAAI,OAAO,CAAC,UAAU,CAAC,WAAW,CAAC,KAAK,CAAC,GAAG,IAAI,CAAC,UAAU;YACzD,OAAO,SAAS,CAAC;QAEnB,UAAU,CAAC,SAAS,GAAG,WAAW,CAAC,QAAQ,CAAC;QAC5C,UAAU,CAAC,MAAM,GAAG,WAAW,CAAC,KAAK,CAAC,CAAE,wCAAwC;QAChF,UAAU,CAAC,SAAS,GAAG,WAAW,CAAC,QAAQ,CAAC;QAC5C,UAAU,CAAC,MAAM,GAAG,WAAW,CAAC,KAAK,CAAC;QACtC,IAAI,CAAC,gBAAgB;YACnB,OAAO,6CAAuB,CAAC,aAAa,CAAC,UAAU,EAAE,UAAU,CAAC,CAAC;QAEvE,MAAM,OAAO,GAAG,qBAAS,CAAC,MAAM,CAAC,UAAU,CAAC,QAAQ,EAAE,UAAU,CAAC,SAAS,CAAC,CAAC;QAC5E,IAAI,OAAO,CAAC,iBAAiB,EAAE,EAAE;YAC/B,OAAO,CAAC,2BAA2B,CAAC,GAAG,CAAC,CAAC;YACzC,MAAM,EAAE,GAAG,OAAO,CAAC,EAAE,CAAC;YACtB,MAAM,EAAE,GAAG,OAAO,CAAC,EAAE,CAAC;YACtB,MAAM,EAAE,GAAG,UAAU,CAAC,0BAA0B,CAAC,EAAE,CAAC,CAAC;YACrD,MAAM,EAAE,GAAG,UAAU,CAAC,0BAA0B,CAAC,EAAE,CAAC,CAAC;YACrD,mCAAmC;YACnC,uBAAuB,CAAC,0CAA0C,CAAC,UAAU,EAAE,EAAE,EAAE,EAAE,EAAE,OAAO,EAAE,OAAO,CAAC,CAAC;YACzG,uBAAuB,CAAC,0CAA0C,CAAC,UAAU,EAAE,EAAE,EAAE,EAAE,EAAE,OAAO,EAAE,OAAO,CAAC,CAAC;YACzG,OAAO,6CAAuB,CAAC,aAAa,CAAC,UAAU,EAAE,UAAU,CAAC,CAAC;SACtE;aAAM;YACL,IAAI,OAAO,CAAC,WAAW,EAAE,GAAG,GAAG,EAAE;gBAC/B,IAAI,UAAU,CAAC,KAAK,CAAC,aAAa,CAAC,OAAO,CAAC,EAAE;oBAC3C,UAAU,CAAC,eAAe,EAAE,CAAC;oBAC7B,UAAU,CAAC,eAAe,EAAE,CAAC;oBAC7B,OAAO,6CAAuB,CAAC,aAAa,CAAC,UAAU,EAAE,UAAU,CAAC,CAAC;iBACtE;gBACD,IAAI,UAAU,CAAC,MAAM,CAAC,aAAa,CAAC,OAAO,CAAC,EAAE;oBAC5C,UAAU,CAAC,aAAa,EAAE,CAAC;oBAC3B,UAAU,CAAC,aAAa,EAAE,CAAC;oBAC3B,OAAO,6CAAuB,CAAC,aAAa,CAAC,UAAU,EAAE,UAAU,CAAC,CAAC;iBACtE;aACF;iBAAM;gBACL,IAAI,UAAU,CAAC,KAAK,CAAC,aAAa,CAAC,OAAO,CAAC,EAAE;oBAC3C,UAAU,CAAC,eAAe,EAAE,CAAC;oBAC7B,UAAU,CAAC,aAAa,EAAE,CAAC;oBAC3B,OAAO,6CAAuB,CAAC,aAAa,CAAC,UAAU,EAAE,UAAU,CAAC,CAAC;iBACtE;gBACD,IAAI,UAAU,CAAC,MAAM,CAAC,aAAa,CAAC,OAAO,CAAC,EAAE;oBAC5C,UAAU,CAAC,aAAa,EAAE,CAAC;oBAC3B,UAAU,CAAC,eAAe,EAAE,CAAC;oBAC7B,OAAO,6CAAuB,CAAC,aAAa,CAAC,UAAU,EAAE,UAAU,CAAC,CAAC;iBACtE;aACF;SACF;QACD,OAAO,SAAS,CAAC;IACnB,CAAC;IACD;;;;;;;OAOG;IACK,gBAAgB,CAAC,GAAmB,EAAE,GAAmB,EAAE,YAAuB,EAAE,UAAkB,EAAE,UAAkB,EAAE,OAAgB;QAClJ,MAAM,MAAM,GAAG,UAAU,GAAG,UAAU,CAAC;QACvC,MAAM,EAAE,GAAG,mBAAQ,CAAC,yBAAyB,CAAC,YAAY,CAAC,EAAE,GAAG,UAAU,EAAE,MAAM,CAAC,CAAC;QACpF,MAAM,EAAE,GAAG,mBAAQ,CAAC,yBAAyB,CAAC,YAAY,CAAC,EAAE,GAAG,UAAU,EAAE,MAAM,CAAC,CAAC;QACpF,IAAI,EAAE,KAAK,SAAS,IAAI,EAAE,KAAK,SAAS,EAAE;YACxC,MAAM,OAAO,GAAG,yCAAmB,CAAC,oCAAoC,CAAC,GAAG,EAAE,YAAY,CAAC,EAAE,EAAE,YAAY,CAAC,EAAE,CAAC,CAAC;YAChH,MAAM,OAAO,GAAG,yCAAmB,CAAC,oCAAoC,CAAC,GAAG,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC;YACtF,IAAI,OAAO,EAAE;gBACX,OAAO,CAAC,sBAAsB,EAAE,CAAC;gBACjC,OAAO,CAAC,sBAAsB,EAAE,CAAC;aAClC;YACD,OAAO,6CAAuB,CAAC,aAAa,CAAC,OAAO,EAAE,OAAO,CAAC,CAAC;SAChE;QACD,OAAO,SAAS,CAAC;IACnB,CAAC;IACO,gBAAgB,CAAC,MAA6C,EAAE,IAAyC;QAC/G,IAAI,IAAI,KAAK,SAAS;YACpB,OAAO,MAAM,CAAC;QAChB,IAAI,MAAM,KAAK,SAAS;YACtB,OAAO,CAAC,IAAI,CAAC,CAAC;QAChB,MAAM,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC;QAClB,OAAO,MAAM,CAAC;IAChB,CAAC;IACD;;;;;;OAMG;IACI,2BAA2B,CAAC,IAAW,EAAE,IAAW,EAAE,oBAA6B,IAAI;QAC5F,IAAI,MAA6C,CAAC;QAClD,IAAI,IAAI,CAAC,MAAM,CAAC,aAAa,CAAC,IAAI,CAAC,MAAM,CAAC,EAAE;YAC1C,MAAM,UAAU,GAAG,IAAI,CAAC,SAAS,CAAC,2BAA2B,CAAC,IAAI,CAAC,SAAS,CAAC,CAAC;YAC9E,IAAI,UAAU,EAAE;gBACd,MAAM,EAAE,GAAG,UAAU,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;gBAAC,MAAM,EAAE,GAAG,UAAU,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;gBAC/D,MAAM,EAAE,GAAG,UAAU,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;gBAAC,MAAM,EAAE,GAAG,UAAU,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;gBAC/D,MAAM,EAAE,GAAG,mBAAQ,CAAC,YAAY,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC;gBACzC,MAAM,EAAE,GAAG,mBAAQ,CAAC,YAAY,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC;gBACzC,MAAM,GAAG,GAAG,mBAAQ,CAAC,cAAc,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC;gBACpD,MAAM,KAAK,GAAG,mBAAQ,CAAC,gBAAgB,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC;gBACxD,IAAI,mBAAQ,CAAC,mBAAmB,CAAC,EAAE,EAAE,GAAG,CAAC;uBACpC,mBAAQ,CAAC,mBAAmB,CAAC,EAAE,EAAE,GAAG,CAAC;uBACrC,mBAAQ,CAAC,mBAAmB,CAAC,CAAC,EAAE,GAAG,CAAC,EAAE;oBACzC,MAAM,cAAc,GAAG,IAAI,CAAC,KAAK,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC,CAAO,8CAA8C;oBAC/F,MAAM,cAAc,GAAG,KAAK,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAG,4DAA4D;oBAC7G,MAAM,gBAAgB,GAAG,cAAc,GAAG,cAAc,GAAG,IAAI,CAAC,KAAK,CAAC,YAAY,CAAC;oBACnF,MAAM,cAAc,GAAG,cAAc,GAAG,cAAc,GAAG,IAAI,CAAC,KAAK,CAAC,UAAU,CAAC;oBAC/E,MAAM,sBAAsB,GAAG,CAAC,cAAc,GAAG,IAAI,CAAC,KAAK,CAAC,YAAY,GAAG,IAAI,CAAC,KAAK,CAAC,YAAY,CAAC,GAAG,CAAC,CAAC;oBACxG,MAAM,MAAM,GAAG,uBAAU,CAAC,qBAAqB,CAAC,gBAAgB,EAAE,cAAc,CAAC,CAAC;oBAClF,MAAM,MAAM,GAAG,IAAI,CAAC,KAAK,CAAC;oBAC1B,MAAM,eAAe,GAAG,MAAM,CAAC,cAAc,EAAE,CAAC;oBAChD,MAAM,UAAU,GAAG,MAAM,CAAC,iCAAiC,CAAC,MAAM,CAAC,YAAY,CAAC,CAAC;oBACjF,MAAM,aAAa,GAAG,MAAM,CAAC,YAAY,GAAG,MAAM,CAAC,YAAY,CAAC;oBAChE,MAAM,UAAU,GAAG,UAAU,GAAG,aAAa,CAAC;oBAC9C,MAAM,cAAc,GAAG,qBAAS,CAAC,MAAM,CAAC,UAAU,EAAE,UAAU,CAAC,CAAC;oBAChE,IAAI,cAAc,CAAC,iBAAiB,EAAE;wBACpC,MAAM,GAAG,IAAI,CAAC,gBAAgB,CAAC,MAAM,EAAE,IAAI,CAAC,gBAAgB,CAAC,IAAI,EAAE,IAAI,EAAE,cAAc,EAAE,UAAU,EAAE,UAAU,EAAE,sBAAsB,CAAC,CAAC,CAAC;oBAC5I,IAAI,UAAU,GAAG,eAAe,EAAE;wBAChC,MAAM,kBAAkB,GAAG,qBAAS,CAAC,MAAM,CAAC,UAAU,GAAG,eAAe,EAAE,UAAU,GAAG,eAAe,CAAC,CAAC;wBACxG,IAAI,kBAAkB,CAAC,iBAAiB,EAAE;4BACxC,MAAM,GAAG,IAAI,CAAC,gBAAgB,CAAC,MAAM,EAAE,IAAI,CAAC,gBAAgB,CAAC,IAAI,EAAE,IAAI,EAAE,kBAAkB,EAAE,UAAU,EAAE,UAAU,EAAE,sBAAsB,CAAC,CAAC,CAAC;qBACjJ;iBACF;aACF;SACF;QACD,OAAO,MAAM,CAAC;IAChB,CAAC;CACF;AA3LD,0DA2LC","sourcesContent":["/*---------------------------------------------------------------------------------------------\r\n* Copyright (c) Bentley Systems, Incorporated. All rights reserved.\r\n* See LICENSE.md in the project root for license terms and full copyright notice.\r\n*--------------------------------------------------------------------------------------------*/\r\n\r\n/** @packageDocumentation\r\n * @module CartesianGeometry\r\n */\r\nimport { Arc3d } from \"../curve/Arc3d\";\r\nimport { CurveLocationDetail, CurveLocationDetailPair } from \"../curve/CurveLocationDetail\";\r\nimport { CurvePrimitive } from \"../curve/CurvePrimitive\";\r\nimport { Geometry } from \"../Geometry\";\r\nimport { AngleSweep } from \"./AngleSweep\";\r\nimport { Point3d, Vector3d } from \"./Point3dVector3d\";\r\nimport { Segment1d } from \"./Segment1d\";\r\n\r\n/**\r\n * `CoincidentGeometryQuery` has methods useful in testing for overlapping geometry.\r\n * * Each instance carries tolerance information that can be reused over extended call sequences.\r\n * * These methods are expected to be called internally by CurveCurve intersection methods.\r\n * @internal\r\n */\r\nexport class CoincidentGeometryQuery {\r\n  private _vectorU?: Vector3d;\r\n  private _vectorV?: Vector3d;\r\n  private _tolerance: number;\r\n  public get tolerance(): number {\r\n    return this._tolerance;\r\n  }\r\n  private constructor(tolerance: number = Geometry.smallMetricDistance) {\r\n    this._tolerance = tolerance;\r\n  }\r\n  public static create(tolerance: number = Geometry.smallMetricDistance): CoincidentGeometryQuery {\r\n    return new CoincidentGeometryQuery(tolerance);\r\n  }\r\n  /**\r\n   * * Assign both the fraction and fraction1 values in the detail, possibly swapped.\r\n   * * reevaluate the points as simple interpolation between given points.\r\n   */\r\n  public static assignDetailInterpolatedFractionsAndPoints(detail: CurveLocationDetail, f0: number, f1: number,\r\n    pointA: Point3d, pointB: Point3d, swap: boolean = false) {\r\n    if (swap) {\r\n      detail.fraction = f1;\r\n      detail.fraction1 = f0;\r\n    } else {\r\n      detail.fraction = f0;\r\n      detail.fraction1 = f1;\r\n    }\r\n    detail.point = pointA.interpolate(detail.fraction, pointB, detail.point);\r\n    detail.point1 = pointA.interpolate(detail.fraction1, pointB, detail.point1);\r\n\r\n  }\r\n\r\n  /** Return a curve location detail with projection of a `spacePoint` to the line segment with `pointA` and `pointB`\r\n   *\r\n   */\r\n  public projectPointToSegmentXY(spacePoint: Point3d, pointA: Point3d, pointB: Point3d): CurveLocationDetail {\r\n    this._vectorU = Vector3d.createStartEnd(pointA, pointB, this._vectorU);\r\n    this._vectorV = Vector3d.createStartEnd(pointA, spacePoint, this._vectorV);\r\n    const uDotU = this._vectorU.dotProductXY(this._vectorU);\r\n    const uDotV = this._vectorU.dotProductXY(this._vectorV);\r\n    const fraction = Geometry.safeDivideFraction(uDotV, uDotU, 0.0);\r\n    return CurveLocationDetail.createCurveFractionPoint(undefined, fraction,\r\n      pointA.interpolate(fraction, pointB));\r\n  }\r\n  /**\r\n   * * project `pointA0` and `pointA1` onto the segment with `pointB0` and `pointB1`\r\n   * * In the returned detail pair, the `detailA` has fractions along segmentA, and `detailB` has fractions along segment B\r\n   * @param pointA0 start point of segment A\r\n   * @param pointA1 end point of segment A\r\n   * @param pointB0 start point of segment B\r\n   * @param pointB1 end point of segment B\r\n   */\r\n  public coincidentSegmentRangeXY(pointA0: Point3d, pointA1: Point3d, pointB0: Point3d, pointB1: Point3d, restrictToBounds: boolean = true): CurveLocationDetailPair | undefined {\r\n    const detailAOnB = this.projectPointToSegmentXY(pointA0, pointB0, pointB1);\r\n    if (pointA0.distanceXY(detailAOnB.point) > this._tolerance)\r\n      return undefined;\r\n    const detailA1OnB = this.projectPointToSegmentXY(pointA1, pointB0, pointB1);\r\n    if (pointA1.distanceXY(detailA1OnB.point) > this._tolerance)\r\n      return undefined;\r\n\r\n    const detailBOnA = this.projectPointToSegmentXY(pointB0, pointA0, pointA1);\r\n    if (pointB0.distanceXY(detailBOnA.point) > this._tolerance)\r\n      return undefined;\r\n    const detailB1OnA = this.projectPointToSegmentXY(pointB1, pointA0, pointA1);\r\n    if (pointB1.distanceXY(detailB1OnA.point) > this._tolerance)\r\n      return undefined;\r\n\r\n    detailAOnB.fraction1 = detailA1OnB.fraction;\r\n    detailAOnB.point1 = detailA1OnB.point;  // capture -- detailB0OnA is not reused.\r\n    detailBOnA.fraction1 = detailB1OnA.fraction;\r\n    detailBOnA.point1 = detailB1OnA.point;\r\n    if (!restrictToBounds)\r\n      return CurveLocationDetailPair.createCapture(detailBOnA, detailAOnB);\r\n\r\n    const segment = Segment1d.create(detailBOnA.fraction, detailBOnA.fraction1);\r\n    if (segment.clampDirectedTo01()) {\r\n      segment.reverseIfNeededForDeltaSign(1.0);\r\n      const f0 = segment.x0;\r\n      const f1 = segment.x1;\r\n      const h0 = detailBOnA.inverseInterpolateFraction(f0);\r\n      const h1 = detailBOnA.inverseInterpolateFraction(f1);\r\n      // recompute fractions and points..\r\n      CoincidentGeometryQuery.assignDetailInterpolatedFractionsAndPoints(detailBOnA, f0, f1, pointA0, pointA1);\r\n      CoincidentGeometryQuery.assignDetailInterpolatedFractionsAndPoints(detailAOnB, h0, h1, pointB0, pointB1);\r\n      return CurveLocationDetailPair.createCapture(detailBOnA, detailAOnB);\r\n    } else {\r\n      if (segment.signedDelta() < 0.0) {\r\n        if (detailBOnA.point.isAlmostEqual(pointA0)) {\r\n          detailBOnA.collapseToStart();\r\n          detailAOnB.collapseToStart();\r\n          return CurveLocationDetailPair.createCapture(detailBOnA, detailAOnB);\r\n        }\r\n        if (detailBOnA.point1.isAlmostEqual(pointA1)) {\r\n          detailBOnA.collapseToEnd();\r\n          detailAOnB.collapseToEnd();\r\n          return CurveLocationDetailPair.createCapture(detailBOnA, detailAOnB);\r\n        }\r\n      } else {\r\n        if (detailBOnA.point.isAlmostEqual(pointA1)) {\r\n          detailBOnA.collapseToStart();\r\n          detailAOnB.collapseToEnd();\r\n          return CurveLocationDetailPair.createCapture(detailBOnA, detailAOnB);\r\n        }\r\n        if (detailBOnA.point1.isAlmostEqual(pointA0)) {\r\n          detailBOnA.collapseToEnd();\r\n          detailAOnB.collapseToStart();\r\n          return CurveLocationDetailPair.createCapture(detailBOnA, detailAOnB);\r\n        }\r\n      }\r\n    }\r\n    return undefined;\r\n  }\r\n  /**\r\n   * Create a CurveLocationDetailPair from . . .\r\n   * @param cpA curveA\r\n   * @param cpB curve B\r\n   * @param fractionsOnA fractions of an active section of curveA\r\n   * @param fractionB0 fraction of an original containing B interval\r\n   * @param fractionB1 end fraction of an original containing B interval\r\n   */\r\n  private createDetailPair(cpA: CurvePrimitive, cpB: CurvePrimitive, fractionsOnA: Segment1d, fractionB0: number, fractionB1: number, reverse: boolean): CurveLocationDetailPair | undefined {\r\n    const deltaB = fractionB1 - fractionB0;\r\n    const g0 = Geometry.conditionalDivideFraction(fractionsOnA.x0 - fractionB0, deltaB);\r\n    const g1 = Geometry.conditionalDivideFraction(fractionsOnA.x1 - fractionB0, deltaB);\r\n    if (g0 !== undefined && g1 !== undefined) {\r\n      const detailA = CurveLocationDetail.createCurveEvaluatedFractionFraction(cpA, fractionsOnA.x0, fractionsOnA.x1);\r\n      const detailB = CurveLocationDetail.createCurveEvaluatedFractionFraction(cpB, g0, g1);\r\n      if (reverse) {\r\n        detailA.swapFractionsAndPoints();\r\n        detailB.swapFractionsAndPoints();\r\n      }\r\n      return CurveLocationDetailPair.createCapture(detailA, detailB);\r\n    }\r\n    return undefined;\r\n  }\r\n  private appendDetailPair(result: CurveLocationDetailPair[] | undefined, pair: CurveLocationDetailPair | undefined): CurveLocationDetailPair[] | undefined {\r\n    if (pair === undefined)\r\n      return result;\r\n    if (result === undefined)\r\n      return [pair];\r\n    result.push(pair);\r\n    return result;\r\n  }\r\n  /**\r\n   * Test if 2 arcs have coinciding portions.\r\n   * @param arcA\r\n   * @param arcB\r\n   * @param _restrictToBounds\r\n   * @return 0, 1, or 2 overlap intervals.\r\n   */\r\n  public coincidentArcIntersectionXY(arcA: Arc3d, arcB: Arc3d, _restrictToBounds: boolean = true): CurveLocationDetailPair[] | undefined {\r\n    let result: CurveLocationDetailPair[] | undefined;\r\n    if (arcA.center.isAlmostEqual(arcB.center)) {\r\n      const matrixBToA = arcA.matrixRef.multiplyMatrixInverseMatrix(arcB.matrixRef);\r\n      if (matrixBToA) {\r\n        const ux = matrixBToA.at(0, 0); const uy = matrixBToA.at(1, 0);\r\n        const vx = matrixBToA.at(0, 1); const vy = matrixBToA.at(1, 1);\r\n        const ru = Geometry.hypotenuseXY(ux, uy);\r\n        const rv = Geometry.hypotenuseXY(vx, vy);\r\n        const dot = Geometry.dotProductXYXY(ux, uy, vx, vy);\r\n        const cross = Geometry.crossProductXYXY(ux, uy, vx, vy);\r\n        if (Geometry.isAlmostEqualNumber(ru, 1.0)\r\n          && Geometry.isAlmostEqualNumber(rv, 1.0)\r\n          && Geometry.isAlmostEqualNumber(0, dot)) {\r\n          const alphaB0Radians = Math.atan2(uy, ux);       // angular position of arcB 0 point in A sweep\r\n          const sweepDirection = cross > 0 ? 1.0 : -1.0;   // 1 if arcB's parameter space sweeps forward, -1 if reverse\r\n          const betaStartRadians = alphaB0Radians + sweepDirection * arcB.sweep.startRadians;\r\n          const betaEndRadians = alphaB0Radians + sweepDirection * arcB.sweep.endRadians;\r\n          const fractionSpacesReversed = (sweepDirection * arcA.sweep.sweepRadians * arcB.sweep.sweepRadians) < 0;\r\n          const sweepB = AngleSweep.createStartEndRadians(betaStartRadians, betaEndRadians);\r\n          const sweepA = arcA.sweep;\r\n          const fractionPeriodA = sweepA.fractionPeriod();\r\n          const fractionB0 = sweepA.radiansToPositivePeriodicFraction(sweepB.startRadians);\r\n          const fractionSweep = sweepB.sweepRadians / sweepA.sweepRadians;\r\n          const fractionB1 = fractionB0 + fractionSweep;\r\n          const fractionSweepB = Segment1d.create(fractionB0, fractionB1);\r\n          if (fractionSweepB.clampDirectedTo01())\r\n            result = this.appendDetailPair(result, this.createDetailPair(arcA, arcB, fractionSweepB, fractionB0, fractionB1, fractionSpacesReversed));\r\n          if (fractionB1 > fractionPeriodA) {\r\n            const fractionSweepBWrap = Segment1d.create(fractionB0 - fractionPeriodA, fractionB1 - fractionPeriodA);\r\n            if (fractionSweepBWrap.clampDirectedTo01())\r\n              result = this.appendDetailPair(result, this.createDetailPair(arcA, arcB, fractionSweepBWrap, fractionB0, fractionB1, fractionSpacesReversed));\r\n          }\r\n        }\r\n      }\r\n    }\r\n    return result;\r\n  }\r\n}\r\n"]}
+{"version":3,"file":"CoincidentGeometryOps.js","sourceRoot":"","sources":["../../../src/geometry3d/CoincidentGeometryOps.ts"],"names":[],"mappings":";AAAA;;;+FAG+F;;;AAE/F;;GAEG;AACH,sDAA6C;AAE7C,sEAA4F;AAE5F,0CAAuC;AACvC,6CAA0C;AAC1C,uDAAsD;AACtD,2CAAwC;AAExC;;;;;GAKG;AACH,MAAa,uBAAuB;IAMlC,IAAW,SAAS;QAClB,OAAO,IAAI,CAAC,UAAU,CAAC;IACzB,CAAC;IACD,YAAoB,YAAoB,mBAAQ,CAAC,mBAAmB;QAClE,IAAI,CAAC,UAAU,GAAG,SAAS,CAAC;IAC9B,CAAC;IACM,MAAM,CAAC,MAAM,CAAC,YAAoB,mBAAQ,CAAC,mBAAmB;QACnE,OAAO,IAAI,uBAAuB,CAAC,SAAS,CAAC,CAAC;IAChD,CAAC;IACD;;;OAGG;IACI,MAAM,CAAC,0CAA0C,CAAC,MAA2B,EAAE,EAAU,EAAE,EAAU,EAC1G,MAAe,EAAE,MAAe,EAAE,OAAgB,KAAK;QACvD,IAAI,IAAI,EAAE;YACR,MAAM,CAAC,QAAQ,GAAG,EAAE,CAAC;YACrB,MAAM,CAAC,SAAS,GAAG,EAAE,CAAC;SACvB;aAAM;YACL,MAAM,CAAC,QAAQ,GAAG,EAAE,CAAC;YACrB,MAAM,CAAC,SAAS,GAAG,EAAE,CAAC;SACvB;QACD,MAAM,CAAC,KAAK,GAAG,MAAM,CAAC,WAAW,CAAC,MAAM,CAAC,QAAQ,EAAE,MAAM,EAAE,MAAM,CAAC,KAAK,CAAC,CAAC;QACzE,MAAM,CAAC,MAAM,GAAG,MAAM,CAAC,WAAW,CAAC,MAAM,CAAC,SAAS,EAAE,MAAM,EAAE,MAAM,CAAC,MAAM,CAAC,CAAC;IAE9E,CAAC;IAED;;OAEG;IACI,uBAAuB,CAAC,UAAmB,EAAE,MAAe,EAAE,MAAe;QAClF,IAAI,CAAC,QAAQ,GAAG,0BAAQ,CAAC,cAAc,CAAC,MAAM,EAAE,MAAM,EAAE,IAAI,CAAC,QAAQ,CAAC,CAAC;QACvE,IAAI,CAAC,QAAQ,GAAG,0BAAQ,CAAC,cAAc,CAAC,MAAM,EAAE,UAAU,EAAE,IAAI,CAAC,QAAQ,CAAC,CAAC;QAC3E,MAAM,KAAK,GAAG,IAAI,CAAC,QAAQ,CAAC,YAAY,CAAC,IAAI,CAAC,QAAQ,CAAC,CAAC;QACxD,MAAM,KAAK,GAAG,IAAI,CAAC,QAAQ,CAAC,YAAY,CAAC,IAAI,CAAC,QAAQ,CAAC,CAAC;QACxD,MAAM,QAAQ,GAAG,mBAAQ,CAAC,kBAAkB,CAAC,KAAK,EAAE,KAAK,EAAE,GAAG,CAAC,CAAC;QAChE,OAAO,yCAAmB,CAAC,wBAAwB,CAAC,SAAS,EAAE,QAAQ,EACrE,MAAM,CAAC,WAAW,CAAC,QAAQ,EAAE,MAAM,CAAC,CAAC,CAAC;IAC1C,CAAC;IACD;;;;;;;OAOG;IACI,wBAAwB,CAAC,OAAgB,EAAE,OAAgB,EAAE,OAAgB,EAAE,OAAgB,EAAE,mBAA4B,IAAI;QACtI,MAAM,WAAW,GAAG,IAAI,CAAC,uBAAuB,CAAC,OAAO,EAAE,OAAO,EAAE,OAAO,CAAC,CAAC;QAC5E,IAAI,OAAO,CAAC,UAAU,CAAC,WAAW,CAAC,KAAK,CAAC,GAAG,IAAI,CAAC,UAAU;YACzD,OAAO,SAAS,CAAC;QACnB,MAAM,WAAW,GAAG,IAAI,CAAC,uBAAuB,CAAC,OAAO,EAAE,OAAO,EAAE,OAAO,CAAC,CAAC;QAC5E,IAAI,OAAO,CAAC,UAAU,CAAC,WAAW,CAAC,KAAK,CAAC,GAAG,IAAI,CAAC,UAAU;YACzD,OAAO,SAAS,CAAC;QAEnB,MAAM,WAAW,GAAG,IAAI,CAAC,uBAAuB,CAAC,OAAO,EAAE,OAAO,EAAE,OAAO,CAAC,CAAC;QAC5E,IAAI,OAAO,CAAC,UAAU,CAAC,WAAW,CAAC,KAAK,CAAC,GAAG,IAAI,CAAC,UAAU;YACzD,OAAO,SAAS,CAAC;QACnB,MAAM,WAAW,GAAG,IAAI,CAAC,uBAAuB,CAAC,OAAO,EAAE,OAAO,EAAE,OAAO,CAAC,CAAC;QAC5E,IAAI,OAAO,CAAC,UAAU,CAAC,WAAW,CAAC,KAAK,CAAC,GAAG,IAAI,CAAC,UAAU;YACzD,OAAO,SAAS,CAAC;QAEnB,WAAW,CAAC,SAAS,GAAG,WAAW,CAAC,QAAQ,CAAC;QAC7C,WAAW,CAAC,MAAM,GAAG,WAAW,CAAC,KAAK,CAAC,CAAE,wCAAwC;QACjF,WAAW,CAAC,SAAS,GAAG,WAAW,CAAC,QAAQ,CAAC;QAC7C,WAAW,CAAC,MAAM,GAAG,WAAW,CAAC,KAAK,CAAC;QACvC,IAAI,CAAC,gBAAgB;YACnB,OAAO,6CAAuB,CAAC,aAAa,CAAC,WAAW,EAAE,WAAW,CAAC,CAAC;QAEzE,MAAM,OAAO,GAAG,qBAAS,CAAC,MAAM,CAAC,WAAW,CAAC,QAAQ,EAAE,WAAW,CAAC,SAAS,CAAC,CAAC;QAC9E,IAAI,OAAO,CAAC,iBAAiB,EAAE,EAAE;YAC/B,MAAM,EAAE,GAAG,OAAO,CAAC,EAAE,CAAC;YACtB,MAAM,EAAE,GAAG,OAAO,CAAC,EAAE,CAAC;YACtB,MAAM,EAAE,GAAG,WAAW,CAAC,0BAA0B,CAAC,EAAE,CAAC,CAAC;YACtD,MAAM,EAAE,GAAG,WAAW,CAAC,0BAA0B,CAAC,EAAE,CAAC,CAAC;YACtD,mCAAmC;YACnC,uBAAuB,CAAC,0CAA0C,CAAC,WAAW,EAAE,EAAE,EAAE,EAAE,EAAE,OAAO,EAAE,OAAO,EAAE,EAAE,GAAG,EAAE,CAAC,CAAC;YACnH,uBAAuB,CAAC,0CAA0C,CAAC,WAAW,EAAE,EAAE,EAAE,EAAE,EAAE,OAAO,EAAE,OAAO,EAAE,EAAE,GAAG,EAAE,CAAC,CAAC;YACnH,OAAO,6CAAuB,CAAC,aAAa,CAAC,WAAW,EAAE,WAAW,CAAC,CAAC;SACxE;aAAM;YACL,IAAI,OAAO,CAAC,WAAW,EAAE,GAAG,GAAG,EAAE;gBAC/B,IAAI,WAAW,CAAC,KAAK,CAAC,aAAa,CAAC,OAAO,EAAE,IAAI,CAAC,SAAS,CAAC,EAAE;oBAC5D,WAAW,CAAC,eAAe,EAAE,CAAC;oBAC9B,WAAW,CAAC,eAAe,EAAE,CAAC;oBAC9B,OAAO,6CAAuB,CAAC,aAAa,CAAC,WAAW,EAAE,WAAW,CAAC,CAAC;iBACxE;gBACD,IAAI,WAAW,CAAC,MAAM,CAAC,aAAa,CAAC,OAAO,EAAE,IAAI,CAAC,SAAS,CAAC,EAAE;oBAC7D,WAAW,CAAC,aAAa,EAAE,CAAC;oBAC5B,WAAW,CAAC,aAAa,EAAE,CAAC;oBAC5B,OAAO,6CAAuB,CAAC,aAAa,CAAC,WAAW,EAAE,WAAW,CAAC,CAAC;iBACxE;aACF;iBAAM;gBACL,IAAI,WAAW,CAAC,KAAK,CAAC,aAAa,CAAC,OAAO,EAAE,IAAI,CAAC,SAAS,CAAC,EAAE;oBAC5D,WAAW,CAAC,eAAe,EAAE,CAAC;oBAC9B,WAAW,CAAC,aAAa,EAAE,CAAC;oBAC5B,OAAO,6CAAuB,CAAC,aAAa,CAAC,WAAW,EAAE,WAAW,CAAC,CAAC;iBACxE;gBACD,IAAI,WAAW,CAAC,MAAM,CAAC,aAAa,CAAC,OAAO,EAAE,IAAI,CAAC,SAAS,CAAC,EAAE;oBAC7D,WAAW,CAAC,aAAa,EAAE,CAAC;oBAC5B,WAAW,CAAC,eAAe,EAAE,CAAC;oBAC9B,OAAO,6CAAuB,CAAC,aAAa,CAAC,WAAW,EAAE,WAAW,CAAC,CAAC;iBACxE;aACF;SACF;QACD,OAAO,SAAS,CAAC;IACnB,CAAC;IACD;;;;;;;;OAQG;IACK,gBAAgB,CAAC,GAAmB,EAAE,GAAmB,EAAE,YAAuB,EAAE,UAAkB,EAAE,UAAkB,EAAE,OAAgB;QAClJ,MAAM,MAAM,GAAG,UAAU,GAAG,UAAU,CAAC;QACvC,MAAM,EAAE,GAAG,mBAAQ,CAAC,yBAAyB,CAAC,YAAY,CAAC,EAAE,GAAG,UAAU,EAAE,MAAM,CAAC,CAAC;QACpF,MAAM,EAAE,GAAG,mBAAQ,CAAC,yBAAyB,CAAC,YAAY,CAAC,EAAE,GAAG,UAAU,EAAE,MAAM,CAAC,CAAC;QACpF,IAAI,EAAE,KAAK,SAAS,IAAI,EAAE,KAAK,SAAS,EAAE;YACxC,MAAM,OAAO,GAAG,yCAAmB,CAAC,oCAAoC,CAAC,GAAG,EAAE,YAAY,CAAC,EAAE,EAAE,YAAY,CAAC,EAAE,CAAC,CAAC;YAChH,MAAM,OAAO,GAAG,yCAAmB,CAAC,oCAAoC,CAAC,GAAG,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC;YACtF,IAAI,OAAO;gBACT,OAAO,CAAC,sBAAsB,EAAE,CAAC;YACnC,OAAO,6CAAuB,CAAC,aAAa,CAAC,OAAO,EAAE,OAAO,CAAC,CAAC;SAChE;QACD,OAAO,SAAS,CAAC;IACnB,CAAC;IACO,gBAAgB,CAAC,MAA6C,EAAE,IAAyC;QAC/G,IAAI,IAAI,KAAK,SAAS;YACpB,OAAO,MAAM,CAAC;QAChB,IAAI,MAAM,KAAK,SAAS;YACtB,OAAO,CAAC,IAAI,CAAC,CAAC;QAChB,MAAM,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC;QAClB,OAAO,MAAM,CAAC;IAChB,CAAC;IACD;;;;;;OAMG;IACI,2BAA2B,CAAC,IAAW,EAAE,IAAW,EAAE,oBAA6B,IAAI;QAC5F,IAAI,MAA6C,CAAC;QAClD,IAAI,IAAI,CAAC,MAAM,CAAC,aAAa,CAAC,IAAI,CAAC,MAAM,EAAE,IAAI,CAAC,SAAS,CAAC,EAAE;YAC1D,MAAM,UAAU,GAAG,IAAI,CAAC,SAAS,CAAC,2BAA2B,CAAC,IAAI,CAAC,SAAS,CAAC,CAAC;YAC9E,IAAI,UAAU,EAAE;gBACd,MAAM,EAAE,GAAG,UAAU,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;gBAAC,MAAM,EAAE,GAAG,UAAU,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;gBAC/D,MAAM,EAAE,GAAG,UAAU,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;gBAAC,MAAM,EAAE,GAAG,UAAU,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;gBAC/D,MAAM,EAAE,GAAG,mBAAQ,CAAC,YAAY,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC;gBACzC,MAAM,EAAE,GAAG,mBAAQ,CAAC,YAAY,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC;gBACzC,MAAM,GAAG,GAAG,mBAAQ,CAAC,cAAc,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC;gBACpD,MAAM,KAAK,GAAG,mBAAQ,CAAC,gBAAgB,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC;gBACxD,IAAI,mBAAQ,CAAC,mBAAmB,CAAC,EAAE,EAAE,GAAG,CAAC;uBACpC,mBAAQ,CAAC,mBAAmB,CAAC,EAAE,EAAE,GAAG,CAAC;uBACrC,mBAAQ,CAAC,mBAAmB,CAAC,CAAC,EAAE,GAAG,CAAC,EAAE;oBACzC,MAAM,cAAc,GAAG,IAAI,CAAC,KAAK,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC,CAAQ,iDAAiD;oBACnG,MAAM,cAAc,GAAG,KAAK,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAI,6EAA6E;oBAC/H,MAAM,gBAAgB,GAAG,cAAc,GAAG,cAAc,GAAG,IAAI,CAAC,KAAK,CAAC,YAAY,CAAC,CAAG,qCAAqC;oBAC3H,MAAM,cAAc,GAAG,cAAc,GAAG,cAAc,GAAG,IAAI,CAAC,KAAK,CAAC,UAAU,CAAC,CAAO,mCAAmC;oBACzH,MAAM,sBAAsB,GAAG,CAAC,cAAc,GAAG,IAAI,CAAC,KAAK,CAAC,YAAY,GAAG,IAAI,CAAC,KAAK,CAAC,YAAY,CAAC,GAAG,CAAC,CAAC;oBACxG,MAAM,MAAM,GAAG,uBAAU,CAAC,qBAAqB,CAAC,gBAAgB,EAAE,cAAc,CAAC,CAAC;oBAClF,MAAM,MAAM,GAAG,IAAI,CAAC,KAAK,CAAC;oBAC1B,MAAM,eAAe,GAAG,MAAM,CAAC,cAAc,EAAE,CAAC;oBAChD,MAAM,UAAU,GAAG,MAAM,CAAC,iCAAiC,CAAC,MAAM,CAAC,YAAY,CAAC,CAAC,CAAG,oCAAoC;oBACxH,IAAA,qBAAM,EAAC,UAAU,IAAI,GAAG,CAAC,CAAC;oBAC1B,MAAM,aAAa,GAAG,MAAM,CAAC,YAAY,GAAG,MAAM,CAAC,YAAY,CAAC,CAAoB,oCAAoC;oBACxH,MAAM,UAAU,GAAG,UAAU,GAAG,aAAa,CAAC,CAAsC,kCAAkC;oBACtH,MAAM,cAAc,GAAG,qBAAS,CAAC,MAAM,CAAC,UAAU,EAAE,UAAU,CAAC,CAAC;oBAEhE;;;;uBAIG;oBACH,MAAM,4BAA4B,GAAG,CAAC,uBAAkC,EAAE,eAAwB,EAAW,EAAE;wBAC7G,MAAM,IAAI,GAAG,MAAM,CAAC,CAAC,CAAC,MAAM,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC;wBACxC,MAAM,SAAS,GAAG,uBAAuB,CAAC,EAAE,CAAC;wBAC7C,MAAM,OAAO,GAAG,uBAAuB,CAAC,EAAE,CAAC;wBAC3C,IAAI,uBAAuB,CAAC,iBAAiB,EAAE,IAAI,CAAC,mBAAQ,CAAC,eAAe,CAAC,uBAAuB,CAAC,aAAa,EAAE,CAAC,EAAE;4BACrH,MAAM,GAAG,IAAI,CAAC,gBAAgB,CAAC,MAAM,EAAE,IAAI,CAAC,gBAAgB,CAAC,IAAI,EAAE,IAAI,EAAE,uBAAuB,EAAE,SAAS,EAAE,OAAO,EAAE,sBAAsB,CAAC,CAAC,CAAC;yBAChJ;6BAAM,EAAG,6BAA6B;4BACrC,MAAM,eAAe,GAAG,sBAAsB,CAAC,CAAC,CAAC,eAAe,CAAC,CAAC,CAAC,CAAC,eAAe,CAAC;4BACpF,MAAM,MAAM,GAAG,IAAI,CAAC,OAAO,GAAG,eAAe,CAAC,CAAC,CAAC,IAAI,CAAC,UAAU,CAAC,IAAI,CAAC,OAAO,CAAC,CAAA,CAAC,CAAC,IAAI,CAAC,QAAQ,CAAC,IAAI,CAAC,OAAO,CAAC,CAAC;4BAC3G,MAAM,MAAM,GAAG,IAAI,CAAC,OAAO,GAAG,eAAe,CAAC,CAAC,CAAC,IAAI,CAAC,UAAU,CAAC,IAAI,CAAC,OAAO,CAAC,CAAA,CAAC,CAAC,IAAI,CAAC,QAAQ,CAAC,IAAI,CAAC,OAAO,CAAC,CAAC;4BAC3G,IAAI,MAAM,CAAC,aAAa,CAAC,MAAM,EAAE,IAAI,CAAC,SAAS,CAAC,EAAE;gCAChD,MAAM,OAAO,GAAG,yCAAmB,CAAC,wBAAwB,CAAC,IAAI,EAAE,eAAe,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC;gCACpG,MAAM,OAAO,GAAG,yCAAmB,CAAC,wBAAwB,CAAC,IAAI,EAAE,eAAe,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC;gCACpG,MAAM,GAAG,IAAI,CAAC,gBAAgB,CAAC,MAAM,EAAE,6CAAuB,CAAC,aAAa,CAAC,OAAO,EAAE,OAAO,CAAC,CAAC,CAAC;6BACjG;yBACF;wBACD,OAAO,MAAM,KAAK,SAAS,IAAI,MAAM,CAAC,MAAM,GAAG,IAAI,CAAC;oBACtD,CAAC,CAAC;oBAEF,4BAA4B,CAAC,cAAc,EAAE,KAAK,CAAC,CAAC,CAAE,wDAAwD;oBAE9G,yEAAyE;oBACzE,IAAI,UAAU,IAAI,eAAe;wBAC/B,4BAA4B,CAAC,qBAAS,CAAC,MAAM,CAAC,UAAU,GAAG,eAAe,EAAE,UAAU,GAAG,eAAe,CAAC,EAAE,IAAI,CAAC,CAAC;yBAC9G,IAAI,UAAU,KAAK,GAAG;wBACzB,4BAA4B,CAAC,qBAAS,CAAC,MAAM,CAAC,UAAU,GAAG,eAAe,EAAE,UAAU,GAAG,eAAe,CAAC,EAAE,IAAI,CAAC,CAAC;iBACpH;aACF;SACF;QACD,OAAO,MAAM,CAAC;IAChB,CAAC;CACF;AArND,0DAqNC","sourcesContent":["/*---------------------------------------------------------------------------------------------\r\n* Copyright (c) Bentley Systems, Incorporated. All rights reserved.\r\n* See LICENSE.md in the project root for license terms and full copyright notice.\r\n*--------------------------------------------------------------------------------------------*/\r\n\r\n/** @packageDocumentation\r\n * @module CartesianGeometry\r\n */\r\nimport { assert } from \"@itwin/core-bentley\";\r\nimport { Arc3d } from \"../curve/Arc3d\";\r\nimport { CurveLocationDetail, CurveLocationDetailPair } from \"../curve/CurveLocationDetail\";\r\nimport { CurvePrimitive } from \"../curve/CurvePrimitive\";\r\nimport { Geometry } from \"../Geometry\";\r\nimport { AngleSweep } from \"./AngleSweep\";\r\nimport { Point3d, Vector3d } from \"./Point3dVector3d\";\r\nimport { Segment1d } from \"./Segment1d\";\r\n\r\n/**\r\n * `CoincidentGeometryQuery` has methods useful in testing for overlapping geometry.\r\n * * Each instance carries tolerance information that can be reused over extended call sequences.\r\n * * These methods are expected to be called internally by CurveCurve intersection methods.\r\n * @internal\r\n */\r\nexport class CoincidentGeometryQuery {\r\n  private _vectorU?: Vector3d;\r\n  private _vectorV?: Vector3d;\r\n  private _point0?: Point3d;\r\n  private _point1?: Point3d;\r\n  private _tolerance: number;\r\n  public get tolerance(): number {\r\n    return this._tolerance;\r\n  }\r\n  private constructor(tolerance: number = Geometry.smallMetricDistance) {\r\n    this._tolerance = tolerance;\r\n  }\r\n  public static create(tolerance: number = Geometry.smallMetricDistance): CoincidentGeometryQuery {\r\n    return new CoincidentGeometryQuery(tolerance);\r\n  }\r\n  /**\r\n   * * Assign both the fraction and fraction1 values in the detail, possibly swapped.\r\n   * * reevaluate the points as simple interpolation between given points.\r\n   */\r\n  public static assignDetailInterpolatedFractionsAndPoints(detail: CurveLocationDetail, f0: number, f1: number,\r\n    pointA: Point3d, pointB: Point3d, swap: boolean = false) {\r\n    if (swap) {\r\n      detail.fraction = f1;\r\n      detail.fraction1 = f0;\r\n    } else {\r\n      detail.fraction = f0;\r\n      detail.fraction1 = f1;\r\n    }\r\n    detail.point = pointA.interpolate(detail.fraction, pointB, detail.point);\r\n    detail.point1 = pointA.interpolate(detail.fraction1, pointB, detail.point1);\r\n\r\n  }\r\n\r\n  /** Return a curve location detail with projection of a `spacePoint` to the line segment with `pointA` and `pointB`\r\n   *\r\n   */\r\n  public projectPointToSegmentXY(spacePoint: Point3d, pointA: Point3d, pointB: Point3d): CurveLocationDetail {\r\n    this._vectorU = Vector3d.createStartEnd(pointA, pointB, this._vectorU);\r\n    this._vectorV = Vector3d.createStartEnd(pointA, spacePoint, this._vectorV);\r\n    const uDotU = this._vectorU.dotProductXY(this._vectorU);\r\n    const uDotV = this._vectorU.dotProductXY(this._vectorV);\r\n    const fraction = Geometry.safeDivideFraction(uDotV, uDotU, 0.0);\r\n    return CurveLocationDetail.createCurveFractionPoint(undefined, fraction,\r\n      pointA.interpolate(fraction, pointB));\r\n  }\r\n  /**\r\n   * * project `pointA0` and `pointA1` onto the segment with `pointB0` and `pointB1`\r\n   * * In the returned detail pair, the `detailA` has fractions along segmentA, and `detailB` has fractions along segment B\r\n   * @param pointA0 start point of segment A\r\n   * @param pointA1 end point of segment A\r\n   * @param pointB0 start point of segment B\r\n   * @param pointB1 end point of segment B\r\n   */\r\n  public coincidentSegmentRangeXY(pointA0: Point3d, pointA1: Point3d, pointB0: Point3d, pointB1: Point3d, restrictToBounds: boolean = true): CurveLocationDetailPair | undefined {\r\n    const detailA0OnB = this.projectPointToSegmentXY(pointA0, pointB0, pointB1);\r\n    if (pointA0.distanceXY(detailA0OnB.point) > this._tolerance)\r\n      return undefined;\r\n    const detailA1OnB = this.projectPointToSegmentXY(pointA1, pointB0, pointB1);\r\n    if (pointA1.distanceXY(detailA1OnB.point) > this._tolerance)\r\n      return undefined;\r\n\r\n    const detailB0OnA = this.projectPointToSegmentXY(pointB0, pointA0, pointA1);\r\n    if (pointB0.distanceXY(detailB0OnA.point) > this._tolerance)\r\n      return undefined;\r\n    const detailB1OnA = this.projectPointToSegmentXY(pointB1, pointA0, pointA1);\r\n    if (pointB1.distanceXY(detailB1OnA.point) > this._tolerance)\r\n      return undefined;\r\n\r\n    detailA0OnB.fraction1 = detailA1OnB.fraction;\r\n    detailA0OnB.point1 = detailA1OnB.point;  // capture -- detailA1OnB is not reused.\r\n    detailB0OnA.fraction1 = detailB1OnA.fraction;\r\n    detailB0OnA.point1 = detailB1OnA.point;\r\n    if (!restrictToBounds)\r\n      return CurveLocationDetailPair.createCapture(detailB0OnA, detailA0OnB);\r\n\r\n    const segment = Segment1d.create(detailB0OnA.fraction, detailB0OnA.fraction1);\r\n    if (segment.clampDirectedTo01()) {\r\n      const f0 = segment.x0;\r\n      const f1 = segment.x1;\r\n      const h0 = detailB0OnA.inverseInterpolateFraction(f0);\r\n      const h1 = detailB0OnA.inverseInterpolateFraction(f1);\r\n      // recompute fractions and points..\r\n      CoincidentGeometryQuery.assignDetailInterpolatedFractionsAndPoints(detailB0OnA, f0, f1, pointA0, pointA1, f0 > f1);\r\n      CoincidentGeometryQuery.assignDetailInterpolatedFractionsAndPoints(detailA0OnB, h0, h1, pointB0, pointB1, h0 > h1);\r\n      return CurveLocationDetailPair.createCapture(detailB0OnA, detailA0OnB);\r\n    } else {\r\n      if (segment.signedDelta() < 0.0) {\r\n        if (detailB0OnA.point.isAlmostEqual(pointA0, this.tolerance)) {\r\n          detailB0OnA.collapseToStart();\r\n          detailA0OnB.collapseToStart();\r\n          return CurveLocationDetailPair.createCapture(detailB0OnA, detailA0OnB);\r\n        }\r\n        if (detailB0OnA.point1.isAlmostEqual(pointA1, this.tolerance)) {\r\n          detailB0OnA.collapseToEnd();\r\n          detailA0OnB.collapseToEnd();\r\n          return CurveLocationDetailPair.createCapture(detailB0OnA, detailA0OnB);\r\n        }\r\n      } else {\r\n        if (detailB0OnA.point.isAlmostEqual(pointA1, this.tolerance)) {\r\n          detailB0OnA.collapseToStart();\r\n          detailA0OnB.collapseToEnd();\r\n          return CurveLocationDetailPair.createCapture(detailB0OnA, detailA0OnB);\r\n        }\r\n        if (detailB0OnA.point1.isAlmostEqual(pointA0, this.tolerance)) {\r\n          detailB0OnA.collapseToEnd();\r\n          detailA0OnB.collapseToStart();\r\n          return CurveLocationDetailPair.createCapture(detailB0OnA, detailA0OnB);\r\n        }\r\n      }\r\n    }\r\n    return undefined;\r\n  }\r\n  /**\r\n   * Create a CurveLocationDetailPair for a coincident interval of two overlapping curves\r\n   * @param cpA curveA\r\n   * @param cpB curveB\r\n   * @param fractionsOnA coincident interval of curveB in fraction space of curveA\r\n   * @param fractionB0 curveB start in fraction space of curveA\r\n   * @param fractionB1 curveB end in fraction space of curveA\r\n   * @param reverse whether curveB and curveA have opposite direction\r\n   */\r\n  private createDetailPair(cpA: CurvePrimitive, cpB: CurvePrimitive, fractionsOnA: Segment1d, fractionB0: number, fractionB1: number, reverse: boolean): CurveLocationDetailPair | undefined {\r\n    const deltaB = fractionB1 - fractionB0;\r\n    const g0 = Geometry.conditionalDivideFraction(fractionsOnA.x0 - fractionB0, deltaB);\r\n    const g1 = Geometry.conditionalDivideFraction(fractionsOnA.x1 - fractionB0, deltaB);\r\n    if (g0 !== undefined && g1 !== undefined) {\r\n      const detailA = CurveLocationDetail.createCurveEvaluatedFractionFraction(cpA, fractionsOnA.x0, fractionsOnA.x1);\r\n      const detailB = CurveLocationDetail.createCurveEvaluatedFractionFraction(cpB, g0, g1);\r\n      if (reverse)\r\n        detailA.swapFractionsAndPoints();\r\n      return CurveLocationDetailPair.createCapture(detailA, detailB);\r\n    }\r\n    return undefined;\r\n  }\r\n  private appendDetailPair(result: CurveLocationDetailPair[] | undefined, pair: CurveLocationDetailPair | undefined): CurveLocationDetailPair[] | undefined {\r\n    if (pair === undefined)\r\n      return result;\r\n    if (result === undefined)\r\n      return [pair];\r\n    result.push(pair);\r\n    return result;\r\n  }\r\n  /**\r\n   * Test if 2 arcs have coinciding portions.\r\n   * @param arcA\r\n   * @param arcB\r\n   * @param _restrictToBounds\r\n   * @return 0, 1, or 2 overlap points/intervals\r\n   */\r\n  public coincidentArcIntersectionXY(arcA: Arc3d, arcB: Arc3d, _restrictToBounds: boolean = true): CurveLocationDetailPair[] | undefined {\r\n    let result: CurveLocationDetailPair[] | undefined;\r\n    if (arcA.center.isAlmostEqual(arcB.center, this.tolerance)) {\r\n      const matrixBToA = arcA.matrixRef.multiplyMatrixInverseMatrix(arcB.matrixRef);\r\n      if (matrixBToA) {\r\n        const ux = matrixBToA.at(0, 0); const uy = matrixBToA.at(1, 0);\r\n        const vx = matrixBToA.at(0, 1); const vy = matrixBToA.at(1, 1);\r\n        const ru = Geometry.hypotenuseXY(ux, uy);\r\n        const rv = Geometry.hypotenuseXY(vx, vy);\r\n        const dot = Geometry.dotProductXYXY(ux, uy, vx, vy);\r\n        const cross = Geometry.crossProductXYXY(ux, uy, vx, vy);\r\n        if (Geometry.isAlmostEqualNumber(ru, 1.0)\r\n          && Geometry.isAlmostEqualNumber(rv, 1.0)\r\n          && Geometry.isAlmostEqualNumber(0, dot)) {\r\n          const alphaB0Radians = Math.atan2(uy, ux);        // angular position of arcB 0 point in arcA sweep\r\n          const sweepDirection = cross > 0 ? 1.0 : -1.0;    // 1 if arcB parameter space sweeps in same direction as arcA, -1 if opposite\r\n          const betaStartRadians = alphaB0Radians + sweepDirection * arcB.sweep.startRadians;   // arcB start in arcA parameter space\r\n          const betaEndRadians = alphaB0Radians + sweepDirection * arcB.sweep.endRadians;       // arcB end in arcA parameter space\r\n          const fractionSpacesReversed = (sweepDirection * arcA.sweep.sweepRadians * arcB.sweep.sweepRadians) < 0;\r\n          const sweepB = AngleSweep.createStartEndRadians(betaStartRadians, betaEndRadians);\r\n          const sweepA = arcA.sweep;\r\n          const fractionPeriodA = sweepA.fractionPeriod();\r\n          const fractionB0 = sweepA.radiansToPositivePeriodicFraction(sweepB.startRadians);   // arcB start in arcA fraction space\r\n          assert(fractionB0 >= 0.0);\r\n          const fractionSweep = sweepB.sweepRadians / sweepA.sweepRadians;                    // arcB sweep in arcA fraction space\r\n          const fractionB1 = fractionB0 + fractionSweep;                                      // arcB end in arcA fraction space\r\n          const fractionSweepB = Segment1d.create(fractionB0, fractionB1);\r\n\r\n          /** lambda to add a coincident interval or isolated intersection, given inputs in arcA fraction space\r\n           * @param arcBInArcAFractionSpace span of arcB in arcA fraction space. On return, clamped to [0,1] if nontrivial.\r\n           * @param testStartOfArcA if no nontrivial coincident interval was found, look for an isolated intersection at the start (true) or end (false) of arcA\r\n           * @returns whether a detail pair was appended to result\r\n           */\r\n          const appendCoincidentIntersection = (arcBInArcAFractionSpace: Segment1d, testStartOfArcA: boolean): boolean => {\r\n            const size = result ? result.length : 0;\r\n            const arcBStart = arcBInArcAFractionSpace.x0;\r\n            const arcBEnd = arcBInArcAFractionSpace.x1;\r\n            if (arcBInArcAFractionSpace.clampDirectedTo01() && !Geometry.isSmallRelative(arcBInArcAFractionSpace.absoluteDelta())) {\r\n              result = this.appendDetailPair(result, this.createDetailPair(arcA, arcB, arcBInArcAFractionSpace, arcBStart, arcBEnd, fractionSpacesReversed));\r\n            } else {  // test isolated intersection\r\n              const testStartOfArcB = fractionSpacesReversed ? testStartOfArcA : !testStartOfArcA;\r\n              const arcAPt = this._point0 = testStartOfArcA ? arcA.startPoint(this._point0): arcA.endPoint(this._point0);\r\n              const arcBPt = this._point1 = testStartOfArcB ? arcB.startPoint(this._point1): arcB.endPoint(this._point1);\r\n              if (arcAPt.isAlmostEqual(arcBPt, this.tolerance)) {\r\n                const detailA = CurveLocationDetail.createCurveFractionPoint(arcA, testStartOfArcA ? 0 : 1, arcAPt);\r\n                const detailB = CurveLocationDetail.createCurveFractionPoint(arcB, testStartOfArcB ? 0 : 1, arcBPt);\r\n                result = this.appendDetailPair(result, CurveLocationDetailPair.createCapture(detailA, detailB));\r\n              }\r\n            }\r\n            return result !== undefined && result.length > size;\r\n          };\r\n\r\n          appendCoincidentIntersection(fractionSweepB, false);  // compute overlap in strict interior, or at end of arcA\r\n\r\n          // check overlap at start of arcA with a periodic shift of fractionSweepB\r\n          if (fractionB1 >= fractionPeriodA)\r\n            appendCoincidentIntersection(Segment1d.create(fractionB0 - fractionPeriodA, fractionB1 - fractionPeriodA), true);\r\n          else if (fractionB0 === 0.0)\r\n            appendCoincidentIntersection(Segment1d.create(fractionB0 + fractionPeriodA, fractionB1 + fractionPeriodA), true);\r\n        }\r\n      }\r\n    }\r\n    return result;\r\n  }\r\n}\r\n"]}

package/lib/cjs/geometry3d/Matrix3d.js

@@ -2387,7 +2387,7 @@ class Matrix3d {
      * almost independent and matrix is invertible).
      */
     conditionNumber() {
-        const determinant = this.determinant();
+        const determinant = Math.abs(this.determinant());
         const columnMagnitudeSum = Geometry_1.Geometry.hypotenuseXYZ(this.coffs[0], this.coffs[3], this.coffs[6])
             + Geometry_1.Geometry.hypotenuseXYZ(this.coffs[1], this.coffs[4], this.coffs[7])
             + Geometry_1.Geometry.hypotenuseXYZ(this.coffs[2], this.coffs[5], this.coffs[8]);