@itwin/core-geometry 4.3.0-dev.8 → 4.3.0

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

@@ -15,6 +15,7 @@ const GeometryHandler_1 = require("../../geometry3d/GeometryHandler");
 const GrowableFloat64Array_1 = require("../../geometry3d/GrowableFloat64Array");
 const Point3dVector3d_1 = require("../../geometry3d/Point3dVector3d");
 const Range_1 = require("../../geometry3d/Range");
+const Newton_1 = require("../../numerics/Newton");
 const Polynomials_1 = require("../../numerics/Polynomials");
 const Arc3d_1 = require("../Arc3d");
 const CurveChainWithDistanceIndex_1 = require("../CurveChainWithDistanceIndex");
@@ -22,7 +23,7 @@ const CurveCollection_1 = require("../CurveCollection");
 const CurveLocationDetail_1 = require("../CurveLocationDetail");
 const LineSegment3d_1 = require("../LineSegment3d");
 const LineString3d_1 = require("../LineString3d");
-// cspell:word XYRR
+// cspell:word XYRR currentdFdX
 /**
  * Handler class for XY close approach between _geometryB and another geometry.
  * * Approach means the XY distance (z is ignored) between _geometryB and another geometry.
@@ -357,10 +358,44 @@ class CurveCurveCloseApproachXY extends GeometryHandler_1.RecurseToCurvesGeometr
             this.testAndRecordProjection(cpB, fB1, pointB1, cpA, fA0, fA1, !reversed);
         }
     }
+    /**
+     * Return XY closest approach between a curve primitive and a point.
+     * Currently, this function only supports Arc3d and LineSegment.
+     * Note that this function doesn't handle endpoints.
+     */
+    getPointCurveClosestApproachXYNewton(curveP, pointQ) {
+        if (!(curveP instanceof Arc3d_1.Arc3d) && !(curveP instanceof LineSegment3d_1.LineSegment3d)) {
+            (0, core_bentley_1.assert)(!!"getPointCurveClosestApproachXYNewton only supports Arc3d and LineSegment");
+            return undefined;
+        }
+        const seeds = [0.2, 0.4, 0.6, 0.8]; // HEURISTIC: arcs have up to 4 perpendiculars; lines have only 1
+        const newtonEvaluator = new Newton_1.CurvePointCloseApproachXYRtoRD(curveP, pointQ);
+        const newtonSearcher = new Newton_1.Newton1dUnbounded(newtonEvaluator);
+        let minCloseApproachLength = Geometry_1.Geometry.largeCoordinateResult;
+        let minCurvePFraction;
+        let minPointP;
+        for (const seed of seeds) {
+            newtonSearcher.setX(seed);
+            if (newtonSearcher.runIterations()) {
+                const curvePFraction = newtonSearcher.getX();
+                if (this.acceptFraction(curvePFraction)) {
+                    const pointP = curveP.fractionToPoint(curvePFraction);
+                    const closeApproachLength = pointP.distanceSquaredXY(pointQ);
+                    if (closeApproachLength < minCloseApproachLength) {
+                        minCloseApproachLength = closeApproachLength;
+                        minCurvePFraction = curvePFraction;
+                        minPointP = pointP;
+                    }
+                }
+            }
+        }
+        if (minCurvePFraction && minPointP)
+            return CurveLocationDetail_1.CurveLocationDetail.createCurveFractionPoint(curveP, minCurvePFraction, minPointP);
+        return undefined;
+    }
     /** Find the closest approach between pointA and cpB. Add the approach if it's within fB0 and fB1. */
     testAndRecordProjection(cpA, fA, pointA, cpB, fB0, fB1, reversed) {
-        // NO NO NO -- this is 3D closest point --- need 2d !!
-        const detail = cpB.closestPoint(pointA, false);
+        const detail = this.getPointCurveClosestApproachXYNewton(cpB, pointA);
         if (detail) {
             const fB = Geometry_1.Geometry.restrictToInterval(detail.fraction, fB0, fB1);
             if (fB === detail.fraction) { // if fraction is within fB0 and fB1
@@ -389,13 +424,13 @@ class CurveCurveCloseApproachXY extends GeometryHandler_1.RecurseToCurvesGeometr
         this.computeSegmentSegment3D(cpA, pointA0, fractionA0, pointA1, fractionA1, cpB, pointB0, fractionB0, pointB1, fractionB1, reversed);
     }
     /**
-     * Low level dispatch of segment with arc.
+     * Low level dispatch of line segment with arc.
      * Find close approaches within maxDistance between a line segments (pointA0, pointA1) and an arc.
      * To consider:
      * 1) intersection between arc and segment.
-     * 2) arc endpoints to segment endpoints or arc endpoints projection to the segment.
-     * 3) line parallel to arc tangent.
-     * @param cpA curve A (line segment or line string)
+     * 2) endpoints to endpoints or endpoints projection to the other curve.
+     * 3) arc tangent parallel to line segment (or line string).
+     * @param cpA curve A (line segment or line string; if it is a line string, then the fractions must specify a segment)
      * @param pointA0 start point of the segment
      * @param fractionA0 fraction of the start of the segment
      * @param pointA1 end point of the segment
@@ -404,7 +439,7 @@ class CurveCurveCloseApproachXY extends GeometryHandler_1.RecurseToCurvesGeometr
      * @param reversed whether to reverse the details in the pair (e.g., so that detailB refers to geometryB).
      */
     dispatchSegmentArc(cpA, pointA0, fractionA0, pointA1, fractionA1, arc, reversed) {
-        // 1) intersection between arc and segment
+        // 1) intersection between arc and line segment (or string).
         // Suppose:
         // Arc: X = C + cU + sV where c = cos(theta) and s = sin(theta)
         // Line: contains points A0 and A1
@@ -437,9 +472,9 @@ class CurveCurveCloseApproachXY extends GeometryHandler_1.RecurseToCurvesGeometr
         }
         if (intersectionFound)
             return;
-        // 2) endpoints to endpoints or endpoints projection to the other curve
-        this.testAndRecordFractionalPairApproach(cpA, fractionA0, fractionA1, true, arc, 0, 1, false, reversed);
-        // 3) line parallel to arc tangent.
+        // 2) endpoints to endpoints or endpoints projection to the other curve.
+        this.testAndRecordFractionalPairApproach(cpA, fractionA0, fractionA1, true, arc, 0, 1, true, reversed);
+        // 3) arc tangent parallel to line segment (or string).
         // If line does not intersect the arc, then the closest (and/or the furthest) point on arc to the line is a
         // point where the tangent line on arc at that point is parallel to the line.
         const dotUT = data.vector0.crossProductStartEndXY(pointA0, pointA1);
@@ -455,63 +490,40 @@ class CurveCurveCloseApproachXY extends GeometryHandler_1.RecurseToCurvesGeometr
             }
         }
     }
-    /** Low level dispatch of circular arc with circular arc. radiusA must be larger than or equal to radiusB. */
-    dispatchCircularCircularOrdered(cpA, radiusA, cpB, radiusB, reversed) {
-        const c = cpA.center.distance(cpB.center);
-        const e = this._maxDistanceToAccept !== undefined ? this._maxDistanceToAccept : Geometry_1.Geometry.smallMetricDistance;
-        if (c > radiusA + radiusB + e) // distance between circles is more than max distance
-            return;
-        // TODO: 1) intersection between arcs
-        // 2) endpoints to endpoints
-        this.testAndRecordFractionalPairApproach(cpA, 0, 1, false, cpB, 0, 1, false, reversed);
-        // 3) line from one arc to another (perpendicular to arc tangents along center-center line)
-        if (!Geometry_1.Geometry.isSmallMetricDistance(c)) {
-            const vectorAB = Point3dVector3d_1.Vector3d.createStartEnd(cpA.center, cpB.center);
-            vectorAB.scaleInPlace(1.0 / c);
-            for (const rA of [-radiusA, radiusA]) {
-                for (const rB of [-radiusB, radiusB]) {
-                    const tangentDistance = c - rA + rB;
-                    if (tangentDistance < e) {
-                        const detailA = this.resolveDirectionToArcXYFraction(cpA, vectorAB, rA);
-                        if (detailA) {
-                            const detailB = this.resolveDirectionToArcXYFraction(cpB, vectorAB, rB);
-                            if (detailB)
-                                this.captureDetailPair(detailA, detailB, reversed);
-                        }
+    /** Solve Newton for 2 arcs and the given newtonEvaluator. */
+    solveArcArcNewton(curveP, curveQ, reversed, newtonEvaluator) {
+        const seedsU = [0.2, 0.4, 0.6, 0.8]; // HEURISTIC: 2 arcs have up to 4 perpendiculars/intersections
+        const seedsV = [0.2, 0.4, 0.6, 0.8];
+        const newtonSearcher = new Newton_1.Newton2dUnboundedWithDerivative(newtonEvaluator);
+        for (const seedU of seedsU) {
+            for (const seedV of seedsV) {
+                newtonSearcher.setUV(seedU, seedV);
+                if (newtonSearcher.runIterations()) {
+                    const curvePFraction = newtonSearcher.getU();
+                    const curveQFraction = newtonSearcher.getV();
+                    if (this.acceptFraction(curvePFraction) && this.acceptFraction(curveQFraction)) {
+                        this.recordPointWithLocalFractions(curvePFraction, curveP, 0, 1, curveQFraction, curveQ, 0, 1, reversed);
                     }
                 }
             }
         }
     }
-    /** Find the fractional point (if any) on the circular `arc` in the direction of `radialVector`. */
-    resolveDirectionToArcXYFraction(arc, radialVector, scale) {
-        // The scale ultimately only affects the direction --- easiest way to use it is two multiplies.
-        const c = scale * arc.matrixRef.columnDotXYZ(0, radialVector.x, radialVector.y, 0);
-        const s = scale * arc.matrixRef.columnDotXYZ(1, radialVector.x, radialVector.y, 0);
-        const radians = Math.atan2(s, c);
-        const fraction = arc.sweep.radiansToPositivePeriodicFraction(radians, 0);
-        if (fraction < 1.0)
-            return CurveLocationDetail_1.CurveLocationDetail.createCurveEvaluatedFraction(arc, fraction);
-        return undefined;
+    /** Find and store perpendicular line between 2 arcs. */
+    findPerpLineXYArcArcNewton(curveP, curveQ, reversed) {
+        const newtonEvaluator = new Newton_1.CurveCurveCloseApproachXYRRtoRRD(curveP, curveQ);
+        this.solveArcArcNewton(curveP, curveQ, reversed, newtonEvaluator);
     }
-    /** Low level dispatch of arc with arc. Only circular arcs are supported. */
+    /** Low level dispatch of arc with Arc3d. */
     dispatchArcArc(cpA, cpB, reversed) {
         const rangeA = cpA.range();
         const rangeB = cpB.range();
         rangeA.expandInPlace(this._maxDistanceToAccept);
         if (!rangeB.intersectsRangeXY(rangeA))
             return;
-        if (this._circularArcB) {
-            const radiusB = this._circularRadiusB;
-            const radiusA = cpA.circularRadiusXY();
-            if (radiusA !== undefined) {
-                if (radiusA >= radiusB)
-                    this.dispatchCircularCircularOrdered(cpA, radiusA, cpB, radiusB, reversed);
-                else
-                    this.dispatchCircularCircularOrdered(cpB, radiusB, cpA, radiusA, !reversed);
-                return;
-            }
-        }
+        // 1) endpoints to endpoints or endpoints projection to the other curve
+        this.testAndRecordFractionalPairApproach(cpA, 0, 1, true, cpB, 0, 1, true, reversed);
+        // 2) perpendicular line between 2 arcs (includes intersections)
+        this.findPerpLineXYArcArcNewton(cpA, cpB, reversed);
     }
     /** Low level dispatch of arc with (beziers of) a bspline curve */
     dispatchArcBsplineCurve3d(cpA, cpB, reversed) {

package/lib/cjs/curve/internalContexts/CurveCurveCloseApproachXY.js.map

@@ -1 +1 @@
-{"version":3,"file":"CurveCurveCloseApproachXY.js","sourceRoot":"","sources":["../../../../src/curve/internalContexts/CurveCurveCloseApproachXY.ts"],"names":[],"mappings":";AAAA;;;+FAG+F;;;AAE/F;;GAEG;AAEH,sDAA6C;AAC7C,6DAAgF;AAEhF,6CAA0C;AAC1C,sEAAkF;AAClF,gFAA6E;AAC7E,sEAAqE;AACrE,kDAAiD;AACjD,4DAAwE;AACxE,oCAAiC;AACjC,gFAA6E;AAC7E,wDAAqD;AACrD,gEAAyG;AAGzG,oDAAiD;AACjD,kDAA+C;AAE/C,mBAAmB;AAEnB;;;;;;;;;;;;;GAaG;AACH,MAAa,yBAA0B,SAAQ,gDAA8B;IAInE,YAAY,CAAC,SAA+B;QAClD,IAAI,CAAC,UAAU,GAAG,SAAS,CAAC;QAC5B,IAAI,CAAC,aAAa,GAAG,SAAS,CAAC;QAC/B,IAAI,CAAC,gBAAgB,GAAG,SAAS,CAAC;QAClC,IAAI,SAAS,YAAY,aAAK,EAAE;YAC9B,MAAM,CAAC,GAAG,SAAS,CAAC,gBAAgB,EAAE,CAAC;YACvC,IAAI,CAAC,KAAK,SAAS,EAAE;gBACnB,IAAI,CAAC,gBAAgB,GAAG,CAAC,CAAC;gBAC1B,IAAI,CAAC,aAAa,GAAG,SAAS,CAAC;aAChC;SACF;IACH,CAAC;IAoBD;;;OAGG;IACH,YAAmB,SAA+B;QAChD,KAAK,EAAE,CAAC;QACR,IAAI,CAAC,YAAY,CAAC,SAAS,CAAC,CAAC;QAC7B,IAAI,CAAC,mBAAmB,GAAG,mBAAQ,CAAC,0BAA0B,CAAC;QAC/D,IAAI,CAAC,QAAQ,GAAG,EAAE,CAAC;IACrB,CAAC;IACD,6EAA6E;IAC7E,IAAW,mBAAmB,CAAC,KAAyB;QACtD,IAAI,CAAC,oBAAoB,GAAG,KAAK,CAAC;QAClC,IAAI,IAAI,CAAC,oBAAoB,KAAK,SAAS,IAAI,IAAI,CAAC,oBAAoB,GAAG,CAAC;YAC1E,IAAI,CAAC,mBAAmB,GAAG,IAAI,CAAC,oBAAoB,GAAG,IAAI,CAAC,oBAAoB,CAAC;IACrF,CAAC;IACD,gFAAgF;IAChF,IAAW,mBAAmB;QAC5B,OAAO,IAAI,CAAC,oBAAoB,CAAC;IACnC,CAAC;IACD,mEAAmE;IACnE,IAAW,gBAAgB;QACzB,OAAO,IAAI,CAAC,oBAAoB,KAAK,SAAS,IAAI,IAAI,CAAC,oBAAoB,GAAG,CAAC,CAAC;IAClF,CAAC;IACD,iHAAiH;IAC1G,aAAa,CAAC,SAAmB;QACtC,IAAI,CAAC,YAAY,CAAC,SAAS,CAAC,CAAC;IAC/B,CAAC;IACD,8DAA8D;IACtD,cAAc,CAAC,QAAgB,EAAE,cAAsB,OAAO;QACpE,IAAI,QAAQ,GAAG,CAAC,WAAW;YACzB,OAAO,KAAK,CAAC;QACf,IAAI,QAAQ,GAAG,GAAG,GAAG,WAAW;YAC9B,OAAO,KAAK,CAAC;QACf,OAAO,IAAI,CAAC;IACd,CAAC;IACD;;;OAGG;IACI,iBAAiB,CAAC,eAAwB,KAAK;QACpD,MAAM,MAAM,GAAG,IAAI,CAAC,QAAQ,CAAC;QAC7B,IAAI,YAAY;YACd,IAAI,CAAC,QAAQ,GAAG,EAAE,CAAC;QACrB,OAAO,MAAM,CAAC;IAChB,CAAC;IACD;;;OAGG;IACK,8BAA8B,CACpC,GAAmB,EAAE,EAAU,EAAE,MAAe,EAAE,GAAmB,EAAE,EAAU,EAAE,MAAe,EAAE,QAAiB;QAErH,MAAM,EAAE,GAAG,MAAM,CAAC,iBAAiB,CAAC,MAAM,CAAC,CAAC;QAC5C,IAAI,EAAE,GAAG,IAAI,CAAC,mBAAmB,EAAE;YACjC,MAAM,OAAO,GAAG,yCAAmB,CAAC,wBAAwB,CAAC,GAAG,EAAE,EAAE,EAAE,MAAM,CAAC,CAAC;YAC9E,MAAM,OAAO,GAAG,yCAAmB,CAAC,wBAAwB,CAAC,GAAG,EAAE,EAAE,EAAE,MAAM,CAAC,CAAC;YAC9E,MAAM,IAAI,GAAG,6CAAuB,CAAC,aAAa,CAAC,OAAO,EAAE,OAAO,CAAC,CAAC;YACrE,IAAI,QAAQ;gBACV,IAAI,CAAC,WAAW,EAAE,CAAC;YACrB,IAAI,CAAC,QAAQ,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC;SAC1B;IACH,CAAC;IACD;;;;;;;;;;;;OAYG;IACK,6BAA6B,CACnC,cAAsB,EACtB,GAAmB,EACnB,UAAkB,EAClB,UAAkB,EAClB,cAAsB,EACtB,GAAmB,EACnB,UAAkB,EAClB,UAAkB,EAClB,QAAiB,EACjB,eAAqD;QAErD,IAAI,eAAe,EAAE,eAAe,CAAC;QACrC,IAAI,gBAAgB,EAAE,gBAAgB,CAAC;QACvC,MAAM,UAAU,GAAG,eAAe,KAAK,SAAS;YAC9C,eAAe,CAAC,OAAO,CAAC,YAAY;YACpC,eAAe,CAAC,OAAO,CAAC,YAAY,CAAC;QACvC,IAAI,UAAU,EAAE;YACd,eAAe,GAAG,mBAAQ,CAAC,WAAW,CAAC,UAAU,EAAE,eAAe,CAAC,OAAO,CAAC,QAAQ,EAAE,UAAU,CAAC,CAAC;YACjG,eAAe,GAAG,mBAAQ,CAAC,WAAW,CAAC,UAAU,EAAE,eAAe,CAAC,OAAO,CAAC,QAAQ,EAAE,UAAU,CAAC,CAAC;YACjG,gBAAgB,GAAG,mBAAQ,CAAC,WAAW,CAAC,UAAU,EAAE,eAAe,CAAC,OAAO,CAAC,SAAU,EAAE,UAAU,CAAC,CAAC;YACpG,gBAAgB,GAAG,mBAAQ,CAAC,WAAW,CAAC,UAAU,EAAE,eAAe,CAAC,OAAO,CAAC,SAAU,EAAE,UAAU,CAAC,CAAC;SACrG;aAAM;YACL,eAAe,GAAG,gBAAgB,GAAG,mBAAQ,CAAC,WAAW,CAAC,UAAU,EAAE,cAAc,EAAE,UAAU,CAAC,CAAC;YAClG,eAAe,GAAG,gBAAgB,GAAG,mBAAQ,CAAC,WAAW,CAAC,UAAU,EAAE,cAAc,EAAE,UAAU,CAAC,CAAC;SACnG;QACD,2CAA2C;QAC3C,MAAM,WAAW,GAAG,IAAI,CAAC,QAAQ,CAAC,MAAM,CAAC;QACzC,IAAI,WAAW,GAAG,CAAC,IAAI,CAAC,UAAU,EAAE;YAClC,MAAM,UAAU,GAAG,IAAI,CAAC,QAAQ,CAAC,WAAW,GAAG,CAAC,CAAC,CAAC,OAAO,CAAC;YAC1D,MAAM,UAAU,GAAG,IAAI,CAAC,QAAQ,CAAC,WAAW,GAAG,CAAC,CAAC,CAAC,OAAO,CAAC;YAC1D,IAAI,QAAQ,EAAE;gBACZ,IAAI,UAAU,CAAC,sBAAsB,CAAC,EAAC,KAAK,EAAE,GAAG,EAAE,QAAQ,EAAE,eAAe,EAAC,CAAC;oBAC5E,UAAU,CAAC,sBAAsB,CAAC,EAAC,KAAK,EAAE,GAAG,EAAE,QAAQ,EAAE,eAAe,EAAC,CAAC;oBAC1E,OAAO;aACV;iBAAM;gBACL,IAAI,UAAU,CAAC,sBAAsB,CAAC,EAAC,KAAK,EAAE,GAAG,EAAE,QAAQ,EAAE,eAAe,EAAC,CAAC;oBAC5E,UAAU,CAAC,sBAAsB,CAAC,EAAC,KAAK,EAAE,GAAG,EAAE,QAAQ,EAAE,eAAe,EAAC,CAAC;oBAC1E,OAAO;aACV;SACF;QACD,MAAM,OAAO,GAAG,yCAAmB,CAAC,wBAAwB,CAC1D,GAAG,EAAE,eAAe,EAAE,GAAG,CAAC,eAAe,CAAC,eAAe,CAAC,CAC3D,CAAC;QACF,MAAM,OAAO,GAAG,yCAAmB,CAAC,wBAAwB,CAC1D,GAAG,EAAE,eAAe,EAAE,GAAG,CAAC,eAAe,CAAC,eAAe,CAAC,CAC3D,CAAC;QACF,IAAI,UAAU,EAAE;YACd,OAAO,CAAC,sBAAsB,CAAC,gBAAgB,EAAE,GAAG,CAAC,eAAe,CAAC,gBAAgB,CAAC,CAAC,CAAC;YACxF,OAAO,CAAC,sBAAsB,CAAC,gBAAgB,EAAE,GAAG,CAAC,eAAe,CAAC,gBAAgB,CAAC,CAAC,CAAC;SACzF;aAAM;YACL,MAAM,EAAE,GAAG,OAAO,CAAC,KAAK,CAAC,iBAAiB,CAAC,OAAO,CAAC,KAAK,CAAC,CAAC;YAC1D,IAAI,EAAE,GAAG,IAAI,CAAC,mBAAmB;gBAC/B,OAAO;YACT,OAAO,CAAC,eAAe,CAAC,uCAAiB,CAAC,QAAQ,CAAC,CAAC;YACpD,OAAO,CAAC,eAAe,CAAC,uCAAiB,CAAC,QAAQ,CAAC,CAAC;SACrD;QACD,IAAI,QAAQ,EAAE;YACZ,IAAI,CAAC,QAAQ,CAAC,IAAI,CAAC,IAAI,6CAAuB,CAAC,OAAO,EAAE,OAAO,CAAC,CAAC,CAAC;SACnE;aAAM;YACL,IAAI,CAAC,QAAQ,CAAC,IAAI,CAAC,IAAI,6CAAuB,CAAC,OAAO,EAAE,OAAO,CAAC,CAAC,CAAC;SACnE;IACH,CAAC;IACD;;;;;;;;;;;;OAYG;IACK,6BAA6B,CACnC,IAA6B,EAC7B,GAAmB,EACnB,UAAkB,EAClB,UAAkB,EAClB,GAAmB,EACnB,UAAkB,EAClB,UAAkB,EAClB,QAAiB;QAEjB,MAAM,eAAe,GAAG,mBAAQ,CAAC,WAAW,CAAC,UAAU,EAAE,IAAI,CAAC,OAAO,CAAC,QAAQ,EAAE,UAAU,CAAC,CAAC;QAC5F,MAAM,eAAe,GAAG,mBAAQ,CAAC,WAAW,CAAC,UAAU,EAAE,IAAI,CAAC,OAAO,CAAC,QAAQ,EAAE,UAAU,CAAC,CAAC;QAC5F,uCAAuC;QACvC,MAAM,WAAW,GAAG,IAAI,CAAC,QAAQ,CAAC,MAAM,CAAC;QACzC,IAAI,WAAW,GAAG,CAAC,EAAE;YACnB,MAAM,UAAU,GAAG,IAAI,CAAC,QAAQ,CAAC,WAAW,GAAG,CAAC,CAAC,CAAC,OAAO,CAAC;YAC1D,MAAM,UAAU,GAAG,IAAI,CAAC,QAAQ,CAAC,WAAW,GAAG,CAAC,CAAC,CAAC,OAAO,CAAC;YAC1D,IAAI,QAAQ,EAAE;gBACZ,IAAI,UAAU,CAAC,sBAAsB,CAAC,EAAC,KAAK,EAAE,GAAG,EAAE,QAAQ,EAAE,eAAe,EAAC,CAAC;oBAC5E,UAAU,CAAC,sBAAsB,CAAC,EAAC,KAAK,EAAE,GAAG,EAAE,QAAQ,EAAE,eAAe,EAAC,CAAC;oBAC1E,OAAO;aACV;iBAAM;gBACL,IAAI,UAAU,CAAC,sBAAsB,CAAC,EAAC,KAAK,EAAE,GAAG,EAAE,QAAQ,EAAE,eAAe,EAAC,CAAC;oBAC5E,UAAU,CAAC,sBAAsB,CAAC,EAAC,KAAK,EAAE,GAAG,EAAE,QAAQ,EAAE,eAAe,EAAC,CAAC;oBAC1E,OAAO;aACV;SACF;QACD,oCAAoC;QACpC,yCAAmB,CAAC,4BAA4B,CAAC,GAAG,EAAE,eAAe,EAAE,IAAI,CAAC,OAAO,CAAC,CAAC;QACrF,yCAAmB,CAAC,4BAA4B,CAAC,GAAG,EAAE,eAAe,EAAE,IAAI,CAAC,OAAO,CAAC,CAAC;QACrF,IAAI,CAAC,OAAO,CAAC,eAAe,CAAC,uCAAiB,CAAC,QAAQ,CAAC,CAAC;QACzD,IAAI,CAAC,OAAO,CAAC,eAAe,CAAC,uCAAiB,CAAC,QAAQ,CAAC,CAAC;QACzD,IAAI,QAAQ;YACV,IAAI,CAAC,WAAW,EAAE,CAAC;QACrB,IAAI,CAAC,QAAQ,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC;IAC3B,CAAC;IACD;;;;;;OAMG;IACI,WAAW,CAChB,GAAmB,EAAE,GAAmB,EAAE,KAA4C,EAAE,QAAiB;QAEzG,IAAI,KAAK,KAAK,SAAS,EAAE;YACvB,KAAK,MAAM,CAAC,IAAI,KAAK,EAAE;gBACrB,IAAI,CAAC,6BAA6B,CAChC,CAAC,CAAC,OAAO,CAAC,QAAQ,EAAE,GAAG,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,OAAO,CAAC,QAAQ,EAAE,GAAG,EAAE,CAAC,EAAE,CAAC,EAAE,QAAQ,EAAE,CAAC,CAC1E,CAAC;aACH;SACF;IACH,CAAC;IACD;;;;;OAKG;IACI,iBAAiB,CACtB,OAAwC,EAAE,OAAwC,EAAE,QAAiB;QAErG,IAAI,OAAO,IAAI,OAAO,EAAE;YACtB,IAAI,QAAQ,EAAE;gBACZ,IAAI,CAAC,QAAQ,CAAC,IAAI,CAAC,6CAAuB,CAAC,aAAa,CAAC,OAAO,EAAE,OAAO,CAAC,CAAC,CAAC;aAC7E;iBAAM;gBACL,IAAI,CAAC,QAAQ,CAAC,IAAI,CAAC,6CAAuB,CAAC,aAAa,CAAC,OAAO,EAAE,OAAO,CAAC,CAAC,CAAC;aAC7E;SACF;IACH,CAAC;IACO,MAAM,CAAC,4BAA4B,CACzC,SAAiB,EACjB,MAAe,EACf,OAAgB,EAChB,OAAgB,EAChB,SAAiB,EACjB,kBAA0B,EAC1B,eAAwC;QAExC,IAAI,OAAO,GAAG,KAAK,CAAC;QACpB,IAAI,SAAS,GAAG,CAAC;YACf,SAAS,GAAG,CAAC,CAAC;aACX,IAAI,SAAS,GAAG,CAAC;YACpB,SAAS,GAAG,CAAC,CAAC;QAChB,IAAI,CAAC,WAAW,GAAG,OAAO,CAAC,WAAW,CAAC,SAAS,EAAE,OAAO,EAAE,IAAI,CAAC,WAAW,CAAC,CAAC;QAC7E,MAAM,eAAe,GAAG,IAAI,CAAC,WAAW,CAAC,iBAAiB,CAAC,MAAM,CAAC,CAAC;QACnE,IAAI,eAAe,IAAI,IAAI,CAAC,GAAG,CAAC,kBAAkB,EAAE,eAAe,CAAC,OAAO,CAAC,CAAC,CAAC,EAAE;YAC9E,eAAe,CAAC,OAAO,CAAC,KAAK,CAAC,SAAS,EAAE,MAAM,EAAE,SAAS,EAAE,eAAe,CAAC,CAAC;YAC7E,eAAe,CAAC,OAAO,CAAC,KAAK,CAAC,SAAS,EAAE,IAAI,CAAC,WAAW,EAAE,SAAS,EAAE,eAAe,CAAC,CAAC;YACvF,OAAO,GAAG,IAAI,CAAC;SAChB;QACD,OAAO,OAAO,CAAC;IACjB,CAAC;IACD;;;;;;;;;;OAUG;IACK,MAAM,CAAC,6BAA6B,CAC1C,EAAW,EACX,EAAW,EACX,EAAW,EACX,EAAW,EACX,kBAA0B;QAE1B,MAAM,EAAE,GAAG,EAAE,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC;QACvB,MAAM,EAAE,GAAG,EAAE,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC;QACvB,MAAM,EAAE,GAAG,EAAE,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC;QACvB,MAAM,EAAE,GAAG,EAAE,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC;QACvB,MAAM,IAAI,GAAG,EAAE,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC;QACzB,MAAM,IAAI,GAAG,EAAE,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC;QACzB,MAAM,IAAI,GAAG,EAAE,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC;QACzB,MAAM,IAAI,GAAG,EAAE,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC;QACzB,MAAM,IAAI,GAAG,EAAE,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC;QACzB,MAAM,IAAI,GAAG,EAAE,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC;QACzB,MAAM,IAAI,GAAG,mBAAQ,CAAC,gBAAgB,CAAC,EAAE,EAAE,EAAE,EAAE,IAAI,EAAE,IAAI,CAAC,CAAC;QAC3D,MAAM,IAAI,GAAG,mBAAQ,CAAC,gBAAgB,CAAC,EAAE,EAAE,EAAE,EAAE,IAAI,EAAE,IAAI,CAAC,CAAC;QAC3D,MAAM,IAAI,GAAG,CAAC,mBAAQ,CAAC,gBAAgB,CAAC,EAAE,EAAE,EAAE,EAAE,IAAI,EAAE,IAAI,CAAC,CAAC;QAC5D,MAAM,IAAI,GAAG,CAAC,mBAAQ,CAAC,gBAAgB,CAAC,EAAE,EAAE,EAAE,EAAE,IAAI,EAAE,IAAI,CAAC,CAAC;QAC5D,IAAI,IAAI,GAAG,IAAI,GAAG,GAAG,IAAI,IAAI,GAAG,IAAI,GAAG,GAAG,EAAE,EAAE,mDAAmD;YAC/F,MAAM,SAAS,GAAG,CAAC,IAAI,GAAG,CAAC,IAAI,GAAG,IAAI,CAAC,CAAC;YACxC,MAAM,SAAS,GAAG,CAAC,IAAI,GAAG,CAAC,IAAI,GAAG,IAAI,CAAC,CAAC;YACxC,OAAO,6CAAuB,CAAC,aAAa,CAC1C,yCAAmB,CAAC,wBAAwB,CAAC,SAAS,EAAE,SAAS,EAAE,EAAE,CAAC,WAAW,CAAC,SAAS,EAAE,EAAE,CAAC,CAAC,EACjG,yCAAmB,CAAC,wBAAwB,CAAC,SAAS,EAAE,SAAS,EAAE,EAAE,CAAC,WAAW,CAAC,SAAS,EAAE,EAAE,CAAC,CAAC,CAClG,CAAC;SACH;QACD,4FAA4F;QAC5F,MAAM,eAAe,GAAG,IAAI,6CAAuB,EAAE,CAAC;QACtD,eAAe,CAAC,OAAO,CAAC,CAAC,GAAG,CAAC,GAAG,kBAAkB,CAAC,CAAC,0CAA0C;QAC9F,IAAI,QAAQ,GAAG,KAAK,CAAC;QACrB,MAAM,EAAE,GAAG,mBAAQ,CAAC,mBAAmB,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC;QAChD,IAAI,IAAI,GAAG,IAAI,IAAI,kBAAkB,GAAG,EAAE,EAAE,EAAE,2BAA2B;YACvE,MAAM,SAAS,GAAG,mBAAQ,CAAC,cAAc,CAAC,EAAE,EAAE,EAAE,EAAE,IAAI,EAAE,IAAI,CAAC,GAAG,EAAE,CAAC;YACnE,IAAI,IAAI,CAAC,4BAA4B,CAAC,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,SAAS,EAAE,kBAAkB,EAAE,eAAe,CAAC;gBAClG,QAAQ,GAAG,IAAI,CAAC;SACnB;QACD,IAAI,IAAI,GAAG,IAAI,IAAI,kBAAkB,GAAG,EAAE,EAAE,EAAE,2BAA2B;YACvE,MAAM,SAAS,GAAG,mBAAQ,CAAC,cAAc,CAAC,EAAE,EAAE,EAAE,EAAE,IAAI,EAAE,IAAI,CAAC,GAAG,EAAE,CAAC;YACnE,IAAI,IAAI,CAAC,4BAA4B,CAAC,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,SAAS,EAAE,kBAAkB,EAAE,eAAe,CAAC;gBAClG,QAAQ,GAAG,IAAI,CAAC;SACnB;QACD,MAAM,EAAE,GAAG,mBAAQ,CAAC,mBAAmB,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC;QAChD,IAAI,IAAI,GAAG,IAAI,IAAI,kBAAkB,GAAG,EAAE,EAAE,EAAE,2BAA2B;YACvE,MAAM,SAAS,GAAG,CAAC,mBAAQ,CAAC,cAAc,CAAC,EAAE,EAAE,EAAE,EAAE,IAAI,EAAE,IAAI,CAAC,GAAG,EAAE,CAAC;YACpE,IAAI,IAAI,CAAC,4BAA4B,CAAC,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,SAAS,EAAE,kBAAkB,EAAE,eAAe,CAAC;gBAClG,QAAQ,GAAG,KAAK,CAAC;SACpB;QACD,IAAI,IAAI,GAAG,IAAI,IAAI,kBAAkB,GAAG,EAAE,EAAE,EAAE,2BAA2B;YACvE,MAAM,SAAS,GAAG,CAAC,mBAAQ,CAAC,cAAc,CAAC,EAAE,EAAE,EAAE,EAAE,IAAI,EAAE,IAAI,CAAC,GAAG,EAAE,CAAC;YACpE,IAAI,IAAI,CAAC,4BAA4B,CAAC,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,SAAS,EAAE,kBAAkB,EAAE,eAAe,CAAC;gBAClG,QAAQ,GAAG,KAAK,CAAC;SACpB;QACD,IAAI,eAAe,CAAC,OAAO,CAAC,CAAC,GAAG,kBAAkB;YAChD,OAAO,SAAS,CAAC;QACnB,IAAI,QAAQ;YACV,eAAe,CAAC,WAAW,EAAE,CAAC;QAChC,OAAO,eAAe,CAAC;IACzB,CAAC;IACD;;;;;;;;;;;;;OAaG;IACK,mCAAmC,CACzC,GAAmB,EACnB,GAAW,EACX,GAAW,EACX,iBAA0B,EAC1B,GAAmB,EACnB,GAAW,EACX,GAAW,EACX,iBAA0B,EAC1B,QAAiB;QAEjB,MAAM,OAAO,GAAG,GAAG,CAAC,eAAe,CAAC,GAAG,CAAC,CAAC;QACzC,MAAM,OAAO,GAAG,GAAG,CAAC,eAAe,CAAC,GAAG,CAAC,CAAC;QACzC,MAAM,OAAO,GAAG,GAAG,CAAC,eAAe,CAAC,GAAG,CAAC,CAAC;QACzC,MAAM,OAAO,GAAG,GAAG,CAAC,eAAe,CAAC,GAAG,CAAC,CAAC;QACzC,IAAI,CAAC,8BAA8B,CAAC,GAAG,EAAE,GAAG,EAAE,OAAO,EAAE,GAAG,EAAE,GAAG,EAAE,OAAO,EAAE,QAAQ,CAAC,CAAC;QACpF,IAAI,CAAC,8BAA8B,CAAC,GAAG,EAAE,GAAG,EAAE,OAAO,EAAE,GAAG,EAAE,GAAG,EAAE,OAAO,EAAE,QAAQ,CAAC,CAAC;QACpF,IAAI,CAAC,8BAA8B,CAAC,GAAG,EAAE,GAAG,EAAE,OAAO,EAAE,GAAG,EAAE,GAAG,EAAE,OAAO,EAAE,QAAQ,CAAC,CAAC;QACpF,IAAI,CAAC,8BAA8B,CAAC,GAAG,EAAE,GAAG,EAAE,OAAO,EAAE,GAAG,EAAE,GAAG,EAAE,OAAO,EAAE,QAAQ,CAAC,CAAC;QACpF,IAAI,iBAAiB,EAAE;YACrB,IAAI,CAAC,uBAAuB,CAAC,GAAG,EAAE,GAAG,EAAE,OAAO,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,QAAQ,CAAC,CAAC;YACzE,IAAI,CAAC,uBAAuB,CAAC,GAAG,EAAE,GAAG,EAAE,OAAO,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,QAAQ,CAAC,CAAC;SAC1E;QACD,IAAI,iBAAiB,EAAE;YACrB,IAAI,CAAC,uBAAuB,CAAC,GAAG,EAAE,GAAG,EAAE,OAAO,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC,QAAQ,CAAC,CAAC;YAC1E,IAAI,CAAC,uBAAuB,CAAC,GAAG,EAAE,GAAG,EAAE,OAAO,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC,QAAQ,CAAC,CAAC;SAC3E;IACH,CAAC;IACD,qGAAqG;IAC7F,uBAAuB,CAC7B,GAAmB,EAAE,EAAU,EAAE,MAAe,EAAE,GAAmB,EAAE,GAAW,EAAE,GAAW,EAAE,QAAiB;QAElH,sDAAsD;QACtD,MAAM,MAAM,GAAG,GAAG,CAAC,YAAY,CAAC,MAAM,EAAE,KAAK,CAAC,CAAC;QAC/C,IAAI,MAAM,EAAE;YACV,MAAM,EAAE,GAAG,mBAAQ,CAAC,kBAAkB,CAAC,MAAM,CAAC,QAAQ,EAAE,GAAG,EAAE,GAAG,CAAC,CAAC;YAClE,IAAI,EAAE,KAAK,MAAM,CAAC,QAAQ,EAAE,EAAE,oCAAoC;gBAChE,IAAI,CAAC,8BAA8B,CAAC,GAAG,EAAE,EAAE,EAAE,MAAM,EAAE,GAAG,EAAE,MAAM,CAAC,QAAQ,EAAE,MAAM,CAAC,KAAK,EAAE,QAAQ,CAAC,CAAC;aACpG;SACF;IACH,CAAC;IACD;;;;;OAKG;IACK,uBAAuB,CAC7B,GAAmB,EACnB,OAAgB,EAChB,UAAkB,EAClB,OAAgB,EAChB,UAAkB,EAClB,GAAmB,EACnB,OAAgB,EAChB,UAAkB,EAClB,OAAgB,EAChB,UAAkB,EAClB,QAAiB;QAEjB,kDAAkD;QAClD,MAAM,QAAQ,GAAG,yBAAyB,CAAC,6BAA6B,CACtE,OAAO,EAAE,OAAO,EAAE,OAAO,EAAE,OAAO,EAAE,IAAI,CAAC,mBAAmB,CAC7D,CAAC;QACF,uFAAuF;QACvF,IAAI,QAAQ,EAAE;YACZ,QAAQ,CAAC,OAAO,CAAC,QAAQ,CAAC,GAAG,CAAC,CAAC;YAC/B,QAAQ,CAAC,OAAO,CAAC,QAAQ,CAAC,GAAG,CAAC,CAAC;YAC/B,IAAI,CAAC,6BAA6B,CAAC,QAAQ,EAAE,GAAG,EAAE,UAAU,EAAE,UAAU,EAAE,GAAG,EAAE,UAAU,EAAE,UAAU,EAAE,QAAQ,CAAC,CAAC;SAClH;IACH,CAAC;IACD,kDAAkD;IAC1C,sBAAsB,CAC5B,GAAmB,EACnB,OAAgB,EAChB,UAAkB,EAClB,OAAgB,EAChB,UAAkB,EAClB,GAAmB,EACnB,OAAgB,EAChB,UAAkB,EAClB,OAAgB,EAChB,UAAkB,EAClB,QAAiB;QAEjB,IAAI,CAAC,uBAAuB,CAC1B,GAAG,EAAE,OAAO,EAAE,UAAU,EAAE,OAAO,EAAE,UAAU,EAC7C,GAAG,EAAE,OAAO,EAAE,UAAU,EAAE,OAAO,EAAE,UAAU,EAC7C,QAAQ,CACT,CAAC;IACJ,CAAC;IACD;;;;;;;;;;;;;;OAcG;IACK,kBAAkB,CACxB,GAAmB,EACnB,OAAgB,EAChB,UAAkB,EAClB,OAAgB,EAChB,UAAkB,EAClB,GAAU,EACV,QAAiB;QAEjB,0CAA0C;QAC1C,WAAW;QACX,+DAA+D;QAC/D,kCAAkC;QAClC,wGAAwG;QACxG,gFAAgF;QAChF,mBAAmB;QACnB,mBAAmB;QACnB,wBAAwB;QACxB,IAAI,iBAAiB,GAAG,KAAK,CAAC;QAC9B,MAAM,IAAI,GAAG,GAAG,CAAC,oBAAoB,EAAE,CAAC;QACxC,MAAM,YAAY,GAAG,OAAO,CAAC;QAC7B,MAAM,YAAY,GAAG,OAAO,CAAC;QAC7B,MAAM,KAAK,GAAG,mBAAQ,CAAC,gBAAgB,CAAC,YAAY,EAAE,CAAC,EAAE,YAAY,EAAE,CAAC,EAAE,IAAI,CAAC,MAAM,EAAE,CAAC,CAAC,CAAC,CAAC,kBAAkB;QAC7G,MAAM,IAAI,GAAG,mBAAQ,CAAC,gBAAgB,CAAC,YAAY,EAAE,CAAC,EAAE,YAAY,EAAE,CAAC,EAAE,IAAI,CAAC,OAAO,EAAE,CAAC,CAAC,CAAC,CAAC,iBAAiB;QAC5G,MAAM,KAAK,GAAG,mBAAQ,CAAC,gBAAgB,CAAC,YAAY,EAAE,CAAC,EAAE,YAAY,EAAE,CAAC,EAAE,IAAI,CAAC,QAAQ,EAAE,CAAC,CAAC,CAAC,CAAC,iBAAiB;QAC9G,MAAM,OAAO,GAAG,IAAI,2CAAoB,CAAC,CAAC,CAAC,CAAC;QAC5C,MAAM,KAAK,GAAG,IAAI,2CAAoB,CAAC,CAAC,CAAC,CAAC;QAC1C,MAAM,OAAO,GAAG,IAAI,2CAAoB,CAAC,CAAC,CAAC,CAAC;QAC5C,MAAM,QAAQ,GAAG,2BAAa,CAAC,yCAAyC,CAAE,qBAAqB;QAC7F,KAAK,EAAE,IAAI,EAAE,KAAK,EAAE,OAAO,EAAE,KAAK,EAAE,OAAO,CAC5C,CAAC;QACF,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,QAAQ,EAAE,CAAC,EAAE,EAAE;YACjC,MAAM,QAAQ,GAAG,IAAI,CAAC,MAAM,CAAC,WAAW,CACtC,IAAI,CAAC,OAAO,EAAE,OAAO,CAAC,gBAAgB,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,QAAQ,EAAE,KAAK,CAAC,gBAAgB,CAAC,CAAC,CAAC,CACpF,CAAC;YACF,MAAM,WAAW,GAAG,IAAI,CAAC,KAAK,CAAC,+BAA+B,CAAC,OAAO,CAAC,gBAAgB,CAAC,CAAC,CAAC,CAAC,CAAC;YAC5F,MAAM,YAAY,GAAG,yBAAW,CAAC,oCAAoC,CAAC,YAAY,EAAE,YAAY,EAAE,QAAQ,CAAC,CAAC;YAC5G,2FAA2F;YAC3F,IAAI,YAAY,KAAK,SAAS,IAAI,IAAI,CAAC,cAAc,CAAC,YAAY,CAAC,IAAI,IAAI,CAAC,cAAc,CAAC,WAAW,CAAC,EAAE;gBACvG,IAAI,CAAC,6BAA6B,CAChC,YAAY,EAAE,GAAG,EAAE,UAAU,EAAE,UAAU,EAAE,WAAW,EAAE,GAAG,EAAE,CAAC,EAAE,CAAC,EAAE,QAAQ,CAC5E,CAAC;gBACF,iBAAiB,GAAG,IAAI,CAAC;aAC1B;SACF;QACD,IAAI,iBAAiB;YACnB,OAAO;QACT,uEAAuE;QACvE,IAAI,CAAC,mCAAmC,CAAC,GAAG,EAAE,UAAU,EAAE,UAAU,EAAE,IAAI,EAAE,GAAG,EAAE,CAAC,EAAE,CAAC,EAAE,KAAK,EAAE,QAAQ,CAAC,CAAC;QACxG,mCAAmC;QACnC,2GAA2G;QAC3G,6EAA6E;QAC7E,MAAM,KAAK,GAAG,IAAI,CAAC,OAAO,CAAC,sBAAsB,CAAC,OAAO,EAAE,OAAO,CAAC,CAAC;QACpE,MAAM,KAAK,GAAG,IAAI,CAAC,QAAQ,CAAC,sBAAsB,CAAC,OAAO,EAAE,OAAO,CAAC,CAAC;QACrE,MAAM,eAAe,GAAG,IAAI,CAAC,KAAK,CAAC,KAAK,EAAE,KAAK,CAAC,CAAC;QACjD,KAAK,MAAM,QAAQ,IAAI,CAAC,eAAe,EAAE,eAAe,GAAG,IAAI,CAAC,EAAE,CAAC,EAAE;YACnE,MAAM,QAAQ,GAAG,IAAI,CAAC,MAAM,CAAC,WAAW,CAAC,IAAI,CAAC,OAAO,EAAE,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,EAAE,IAAI,CAAC,QAAQ,EAAE,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,CAAC,CAAC;YAC9G,MAAM,WAAW,GAAG,IAAI,CAAC,KAAK,CAAC,+BAA+B,CAAC,QAAQ,CAAC,CAAC;YACzE,MAAM,YAAY,GAAG,yBAAW,CAAC,oCAAoC,CAAC,YAAY,EAAE,YAAY,EAAE,QAAQ,CAAC,CAAC;YAC5G,2FAA2F;YAC3F,IAAI,YAAY,KAAK,SAAS,IAAI,IAAI,CAAC,cAAc,CAAC,YAAY,CAAC,IAAI,IAAI,CAAC,cAAc,CAAC,WAAW,CAAC,EAAE;gBACvG,IAAI,CAAC,6BAA6B,CAChC,YAAY,EAAE,GAAG,EAAE,UAAU,EAAE,UAAU,EAAE,WAAW,EAAE,GAAG,EAAE,CAAC,EAAE,CAAC,EAAE,QAAQ,CAC5E,CAAC;aACH;SACF;IACH,CAAC;IACD,6GAA6G;IACrG,+BAA+B,CACrC,GAAU,EAAE,OAAe,EAAE,GAAU,EAAE,OAAe,EAAE,QAAiB;QAE3E,MAAM,CAAC,GAAG,GAAG,CAAC,MAAM,CAAC,QAAQ,CAAC,GAAG,CAAC,MAAM,CAAC,CAAC;QAC1C,MAAM,CAAC,GAAG,IAAI,CAAC,oBAAoB,KAAK,SAAS,CAAC,CAAC,CAAC,IAAI,CAAC,oBAAoB,CAAC,CAAC,CAAC,mBAAQ,CAAC,mBAAmB,CAAC;QAC7G,IAAI,CAAC,GAAG,OAAO,GAAG,OAAO,GAAG,CAAC,EAAE,qDAAqD;YAClF,OAAO;QACT,qCAAqC;QACrC,4BAA4B;QAC5B,IAAI,CAAC,mCAAmC,CAAC,GAAG,EAAE,CAAC,EAAE,CAAC,EAAE,KAAK,EAAE,GAAG,EAAE,CAAC,EAAE,CAAC,EAAE,KAAK,EAAE,QAAQ,CAAC,CAAC;QACvF,2FAA2F;QAC3F,IAAI,CAAC,mBAAQ,CAAC,qBAAqB,CAAC,CAAC,CAAC,EAAE;YACtC,MAAM,QAAQ,GAAG,0BAAQ,CAAC,cAAc,CAAC,GAAG,CAAC,MAAM,EAAE,GAAG,CAAC,MAAM,CAAC,CAAC;YACjE,QAAQ,CAAC,YAAY,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC;YAC/B,KAAK,MAAM,EAAE,IAAI,CAAC,CAAC,OAAO,EAAE,OAAO,CAAC,EAAE;gBACpC,KAAK,MAAM,EAAE,IAAI,CAAC,CAAC,OAAO,EAAE,OAAO,CAAC,EAAE;oBACpC,MAAM,eAAe,GAAG,CAAC,GAAG,EAAE,GAAG,EAAE,CAAC;oBACpC,IAAI,eAAe,GAAG,CAAC,EAAE;wBACvB,MAAM,OAAO,GAAG,IAAI,CAAC,+BAA+B,CAAC,GAAG,EAAE,QAAQ,EAAE,EAAE,CAAC,CAAC;wBACxE,IAAI,OAAO,EAAE;4BACX,MAAM,OAAO,GAAG,IAAI,CAAC,+BAA+B,CAAC,GAAG,EAAE,QAAQ,EAAE,EAAE,CAAC,CAAC;4BACxE,IAAI,OAAO;gCACT,IAAI,CAAC,iBAAiB,CAAC,OAAO,EAAE,OAAO,EAAE,QAAQ,CAAC,CAAC;yBACtD;qBACF;iBACF;aACF;SACF;IACH,CAAC;IACD,mGAAmG;IAC3F,+BAA+B,CACrC,GAAU,EAAE,YAAsB,EAAE,KAAa;QAEjD,+FAA+F;QAC/F,MAAM,CAAC,GAAG,KAAK,GAAG,GAAG,CAAC,SAAS,CAAC,YAAY,CAAC,CAAC,EAAE,YAAY,CAAC,CAAC,EAAE,YAAY,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;QACnF,MAAM,CAAC,GAAG,KAAK,GAAG,GAAG,CAAC,SAAS,CAAC,YAAY,CAAC,CAAC,EAAE,YAAY,CAAC,CAAC,EAAE,YAAY,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;QACnF,MAAM,OAAO,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;QACjC,MAAM,QAAQ,GAAG,GAAG,CAAC,KAAK,CAAC,iCAAiC,CAAC,OAAO,EAAE,CAAC,CAAC,CAAC;QACzE,IAAI,QAAQ,GAAG,GAAG;YAChB,OAAO,yCAAmB,CAAC,4BAA4B,CAAC,GAAG,EAAE,QAAQ,CAAC,CAAC;QACzE,OAAO,SAAS,CAAC;IACnB,CAAC;IACD,4EAA4E;IACpE,cAAc,CAAC,GAAU,EAAE,GAAU,EAAE,QAAiB;QAC9D,MAAM,MAAM,GAAG,GAAG,CAAC,KAAK,EAAE,CAAC;QAC3B,MAAM,MAAM,GAAG,GAAG,CAAC,KAAK,EAAE,CAAC;QAC3B,MAAM,CAAC,aAAa,CAAC,IAAI,CAAC,oBAAqB,CAAC,CAAC;QACjD,IAAI,CAAC,MAAM,CAAC,iBAAiB,CAAC,MAAM,CAAC;YACnC,OAAO;QACT,IAAI,IAAI,CAAC,aAAa,EAAE;YACtB,MAAM,OAAO,GAAG,IAAI,CAAC,gBAAiB,CAAC;YACvC,MAAM,OAAO,GAAG,GAAG,CAAC,gBAAgB,EAAE,CAAC;YACvC,IAAI,OAAO,KAAK,SAAS,EAAE;gBACzB,IAAI,OAAO,IAAI,OAAO;oBACpB,IAAI,CAAC,+BAA+B,CAAC,GAAG,EAAE,OAAO,EAAE,GAAG,EAAE,OAAO,EAAE,QAAQ,CAAC,CAAC;;oBAE3E,IAAI,CAAC,+BAA+B,CAAC,GAAG,EAAE,OAAO,EAAE,GAAG,EAAE,OAAO,EAAE,CAAC,QAAQ,CAAC,CAAC;gBAC9E,OAAO;aACR;SACF;IACH,CAAC;IACD,kEAAkE;IAC1D,yBAAyB,CAAC,GAAU,EAAE,GAAmB,EAAE,QAAiB;QAClF,MAAM,EAAE,GAAG,2BAAY,CAAC,MAAM,EAAE,CAAC;QACjC,GAAG,CAAC,WAAW,CAAC,EAAE,CAAC,CAAC;QACpB,IAAI,CAAC,oBAAoB,CAAC,GAAG,EAAE,EAAE,EAAE,QAAQ,CAAC,CAAC;IAC/C,CAAC;IACD,2FAA2F;IACnF,oCAAoC,CAC1C,OAA2B,EAAE,OAA2B,EAAE,QAAiB;QAE3E,MAAM,GAAG,GAAG,2BAAY,CAAC,MAAM,EAAE,CAAC;QAClC,OAAO,CAAC,WAAW,CAAC,GAAG,CAAC,CAAC;QACzB,MAAM,GAAG,GAAG,2BAAY,CAAC,MAAM,EAAE,CAAC;QAClC,OAAO,CAAC,WAAW,CAAC,GAAG,CAAC,CAAC;QACzB,IAAI,CAAC,2BAA2B,CAAC,GAAG,EAAE,GAAG,EAAE,QAAQ,CAAC,CAAC;IACvD,CAAC;IACD,yEAAyE;IAClE,8BAA8B,CAAC,GAAiB,EAAE,MAAsB,EAAE,QAAiB;QAChG,MAAM,GAAG,GAAG,2BAAY,CAAC,MAAM,EAAE,CAAC;QAClC,MAAM,CAAC,WAAW,CAAC,GAAG,CAAC,CAAC;QACxB,IAAI,CAAC,2BAA2B,CAAC,GAAG,EAAE,GAAG,EAAE,QAAQ,CAAC,CAAC;IACvD,CAAC;IACD,sEAAsE;IAC/D,2BAA2B,CAAC,IAAmB,EAAE,MAAsB,EAAE,QAAiB;QAC/F,MAAM,GAAG,GAAG,2BAAY,CAAC,MAAM,EAAE,CAAC;QAClC,MAAM,CAAC,WAAW,CAAC,GAAG,CAAC,CAAC;QACxB,IAAI,CAAC,wBAAwB,CAAC,IAAI,EAAE,GAAG,EAAE,QAAQ,CAAC,CAAC;IACrD,CAAC;IACD,6DAA6D;IACtD,wBAAwB,CAAC,IAAmB,EAAE,GAAiB,EAAE,QAAiB;QACvF,MAAM,IAAI,GAAG,GAAG,CAAC,SAAS,EAAE,CAAC;QAC7B,MAAM,UAAU,GAAG,mBAAQ,CAAC,kBAAkB,CAAC,CAAC,EAAE,IAAI,GAAG,CAAC,EAAE,CAAC,CAAC,CAAC;QAC/D,MAAM,OAAO,GAAG,IAAI,CAAC,SAAS,CAAC;QAC/B,MAAM,OAAO,GAAG,IAAI,CAAC,SAAS,CAAC;QAC/B,MAAM,OAAO,GAAG,yBAAyB,CAAC,aAAa,CAAC;QACxD,MAAM,OAAO,GAAG,yBAAyB,CAAC,aAAa,CAAC;QACxD,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,GAAG,CAAC,EAAE,EAAE,CAAC,EAAE;YACjC,MAAM,GAAG,GAAG,CAAC,GAAG,UAAU,CAAC,CAAC,8BAA8B;YAC1D,MAAM,GAAG,GAAG,CAAC,CAAC,GAAG,CAAC,KAAK,IAAI,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,UAAU,CAAC,CAAE,qCAAqC;YACrG,GAAG,CAAC,YAAY,CAAC,+BAA+B,CAAC,CAAC,EAAE,OAAO,CAAC,CAAC;YAC7D,GAAG,CAAC,YAAY,CAAC,+BAA+B,CAAC,CAAC,GAAG,CAAC,EAAE,OAAO,CAAC,CAAC;YACjE,IAAI,CAAC,sBAAsB,CAAC,IAAI,EAAE,OAAO,EAAE,GAAG,EAAE,OAAO,EAAE,GAAG,EAAE,GAAG,EAAE,OAAO,EAAE,GAAG,EAAE,OAAO,EAAE,GAAG,EAAE,QAAQ,CAAC,CAAC;SAC1G;IACH,CAAC;IACD,yDAAyD;IAClD,oBAAoB,CAAC,IAAW,EAAE,GAAiB,EAAE,QAAiB;QAC3E,MAAM,MAAM,GAAG,IAAI,CAAC,KAAK,EAAE,CAAC;QAC5B,MAAM,MAAM,GAAG,GAAG,CAAC,KAAK,EAAE,CAAC;QAC3B,MAAM,CAAC,aAAa,CAAC,IAAI,CAAC,oBAAqB,CAAC,CAAC;QACjD,IAAI,CAAC,MAAM,CAAC,iBAAiB,CAAC,MAAM,CAAC;YACnC,OAAO;QACT,MAAM,OAAO,GAAG,yBAAyB,CAAC,aAAa,CAAC;QACxD,MAAM,OAAO,GAAG,yBAAyB,CAAC,aAAa,CAAC;QACxD,MAAM,IAAI,GAAG,GAAG,CAAC,SAAS,EAAE,CAAC;QAC7B,IAAI,IAAI,GAAG,CAAC,EAAE;YACZ,MAAM,GAAG,GAAG,GAAG,GAAG,CAAC,IAAI,GAAG,CAAC,CAAC,CAAC;YAC7B,IAAI,GAAG,CAAC;YACR,IAAI,GAAG,CAAC;YACR,GAAG,GAAG,GAAG,CAAC;YACV,GAAG,CAAC,OAAO,CAAC,CAAC,EAAE,OAAO,CAAC,CAAC;YACxB,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,IAAI,EAAE,EAAE,EAAE,EAAE,OAAO,CAAC,OAAO,CAAC,OAAO,CAAC,EAAE,GAAG,GAAG,GAAG,EAAE;gBACrE,GAAG,CAAC,OAAO,CAAC,EAAE,EAAE,OAAO,CAAC,CAAC;gBACzB,GAAG,GAAG,EAAE,GAAG,GAAG,CAAC;gBACf,IAAI,CAAC,kBAAkB,CAAC,GAAG,EAAE,OAAO,EAAE,GAAG,EAAE,OAAO,EAAE,GAAG,EAAE,IAAI,EAAE,CAAC,QAAQ,CAAC,CAAC;aAC3E;SACF;QACD,OAAO,SAAS,CAAC;IACnB,CAAC;IACD,8CAA8C;IACtC,uBAAuB,CAAC,KAAe,EAAE,YAAiC;QAChF,MAAM,KAAK,GAAG,IAAI,CAAC,UAAU,CAAC,CAAE,OAAO;QACvC,IAAI,CAAC,KAAK,IAAI,CAAC,KAAK,CAAC,QAAQ,IAAI,CAAC,CAAC,KAAK,YAAY,iCAAe,CAAC;YAClE,OAAO;QACT,KAAK,MAAM,KAAK,IAAI,KAAK,CAAC,QAAQ,EAAE;YAClC,IAAI,CAAC,aAAa,CAAC,KAAK,CAAC,CAAC;YAC1B,YAAY,CAAC,KAAK,CAAC,CAAC;SACrB;QACD,IAAI,CAAC,UAAU,GAAG,KAAK,CAAC,CAAE,UAAU;IACtC,CAAC;IACD,oFAAoF;IAC5E,mCAAmC,CAAC,KAAe,EAAE,YAAiC;QAC5F,IAAI,CAAC,IAAI,CAAC,UAAU,IAAI,CAAC,CAAC,IAAI,CAAC,UAAU,YAAY,yDAA2B,CAAC;YAC/E,OAAO;QACT,IAAI,KAAK,YAAY,yDAA2B,EAAE;YAChD,IAAA,qBAAM,EAAC,CAAC,CAAC,uDAAuD,CAAC,CAAC;YAClE,OAAO;SACR;QACD,MAAM,MAAM,GAAG,IAAI,CAAC,QAAQ,CAAC,MAAM,CAAC;QACpC,MAAM,KAAK,GAAG,IAAI,CAAC,UAAU,CAAC,CAAE,OAAO;QACvC,KAAK,MAAM,KAAK,IAAI,KAAK,CAAC,IAAI,CAAC,QAAQ,EAAE;YACvC,IAAI,CAAC,aAAa,CAAC,KAAK,CAAC,CAAC;YAC1B,YAAY,CAAC,KAAK,CAAC,CAAC;SACrB;QACD,IAAI,CAAC,aAAa,CAAC,KAAK,CAAC,CAAC,CAAE,UAAU;QACtC,IAAI,CAAC,QAAQ,GAAG,yDAA2B,CAAC,+BAA+B,CAAC,IAAI,CAAC,QAAQ,EAAE,MAAM,EAAE,SAAS,EAAE,KAAK,EAAE,IAAI,CAAC,CAAC;IAC7H,CAAC;IACD,0DAA0D;IAC1C,mBAAmB,CAAC,QAAuB;QACzD,IAAI,IAAI,CAAC,UAAU,YAAY,6BAAa,EAAE;YAC5C,MAAM,QAAQ,GAAG,IAAI,CAAC,UAAU,CAAC;YACjC,IAAI,CAAC,sBAAsB,CACzB,QAAQ,EAAE,QAAQ,CAAC,SAAS,EAAE,GAAG,EAAE,QAAQ,CAAC,SAAS,EAAE,GAAG,EAC1D,QAAQ,EAAE,QAAQ,CAAC,SAAS,EAAE,GAAG,EAAE,QAAQ,CAAC,SAAS,EAAE,GAAG,EAC1D,KAAK,CACN,CAAC;SACH;aAAM,IAAI,IAAI,CAAC,UAAU,YAAY,2BAAY,EAAE;YAClD,IAAI,CAAC,wBAAwB,CAAC,QAAQ,EAAE,IAAI,CAAC,UAAU,EAAE,KAAK,CAAC,CAAC;SACjE;aAAM,IAAI,IAAI,CAAC,UAAU,YAAY,aAAK,EAAE;YAC3C,IAAI,CAAC,kBAAkB,CAAC,QAAQ,EAAE,QAAQ,CAAC,SAAS,EAAE,GAAG,EAAE,QAAQ,CAAC,SAAS,EAAE,GAAG,EAAE,IAAI,CAAC,UAAU,EAAE,KAAK,CAAC,CAAC;SAC7G;aAAM,IAAI,IAAI,CAAC,UAAU,YAAY,6BAAc,EAAE;YACpD,IAAI,CAAC,2BAA2B,CAAC,QAAQ,EAAE,IAAI,CAAC,UAAU,EAAE,KAAK,CAAC,CAAC;SACpE;aAAM,IAAI,IAAI,CAAC,UAAU,YAAY,iCAAe,EAAE;YACrD,IAAI,CAAC,uBAAuB,CAAC,QAAQ,EAAE,IAAI,CAAC,mBAAmB,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC;SAC7E;aAAM,IAAI,IAAI,CAAC,UAAU,YAAY,yDAA2B,EAAE;YACjE,IAAI,CAAC,mCAAmC,CAAC,QAAQ,EAAE,IAAI,CAAC,mBAAmB,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC;SACzF;QACD,OAAO,SAAS,CAAC;IACnB,CAAC;IACD;;;;;;;;;;;;;;;OAeG;IACK,wBAAwB,CAAC,CAAS,EAAE,CAAS,EAAE,KAAc;QACnE,IAAI,MAAM,GAAG,CAAC,CAAC;QACf,IAAI,CAAC,GAAG,KAAK,CAAC,GAAG,CAAC,CAAC;YACjB,MAAM,GAAG,IAAI,CAAC;aACX,IAAI,CAAC,GAAG,KAAK,CAAC,IAAI,CAAC,CAAC;YACvB,MAAM,GAAG,IAAI,CAAC;QAChB,wBAAwB;QACxB,IAAI,CAAC,GAAG,KAAK,CAAC,GAAG,CAAC,CAAC;YACjB,MAAM,IAAI,IAAI,CAAC;aACZ,IAAI,CAAC,GAAG,KAAK,CAAC,IAAI,CAAC,CAAC;YACvB,MAAM,IAAI,IAAI,CAAC;QACjB,OAAO,MAAM,CAAC;IAChB,CAAC;IACD,0DAA0D;IAClD,2BAA2B,CAAC,GAAiB,EAAE,GAAiB,EAAE,QAAiB;QACzF,MAAM,MAAM,GAAG,GAAG,CAAC,KAAK,EAAE,CAAC;QAC3B,MAAM,MAAM,GAAG,GAAG,CAAC,KAAK,EAAE,CAAC;QAC3B,MAAM,CAAC,aAAa,CAAC,IAAI,CAAC,oBAAqB,CAAC,CAAC;QACjD,IAAI,CAAC,MAAM,CAAC,iBAAiB,CAAC,MAAM,CAAC;YACnC,OAAO;QACT,IAAI,KAAa,CAAC;QAClB,IAAI,KAAa,CAAC;QAClB,MAAM,OAAO,GAAG,eAAO,CAAC,UAAU,EAAE,CAAC;QACrC,MAAM,OAAO,GAAG,yBAAyB,CAAC,aAAa,CAAC;QACxD,MAAM,OAAO,GAAG,yBAAyB,CAAC,aAAa,CAAC;QACxD,MAAM,OAAO,GAAG,yBAAyB,CAAC,aAAa,CAAC;QACxD,MAAM,OAAO,GAAG,yBAAyB,CAAC,aAAa,CAAC;QACxD,MAAM,IAAI,GAAG,GAAG,CAAC,SAAS,EAAE,CAAC;QAC7B,MAAM,IAAI,GAAG,GAAG,CAAC,SAAS,EAAE,CAAC;QAC7B,IAAI,IAAI,GAAG,CAAC,IAAI,IAAI,GAAG,CAAC,EAAE;YACxB,MAAM,GAAG,GAAG,GAAG,GAAG,CAAC,IAAI,GAAG,CAAC,CAAC,CAAC;YAC7B,MAAM,GAAG,GAAG,GAAG,GAAG,CAAC,IAAI,GAAG,CAAC,CAAC,CAAC;YAC7B,IAAI,GAAG,GAAG,GAAG,CAAC;YACd,IAAI,GAAG,EAAE,GAAG,EAAE,GAAG,CAAC;YAClB,GAAG,CAAC,OAAO,CAAC,CAAC,EAAE,OAAO,CAAC,CAAC;YACxB,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,IAAI,EAAE,EAAE,EAAE,EAAE,OAAO,CAAC,OAAO,CAAC,OAAO,CAAC,EAAE,GAAG,GAAG,GAAG,EAAE;gBACrE,GAAG,GAAG,EAAE,GAAG,GAAG,CAAC;gBACf,GAAG,GAAG,GAAG,CAAC;gBACV,GAAG,CAAC,OAAO,CAAC,EAAE,EAAE,OAAO,CAAC,CAAC;gBACzB,kEAAkE;gBAClE,OAAO,CAAC,OAAO,EAAE,CAAC;gBAClB,OAAO,CAAC,WAAW,CAAC,OAAO,CAAC,CAAC;gBAC7B,OAAO,CAAC,WAAW,CAAC,OAAO,CAAC,CAAC;gBAC7B,OAAO,CAAC,aAAa,CAAC,IAAI,CAAC,oBAAqB,CAAC,CAAC;gBAClD,IAAI,OAAO,CAAC,iBAAiB,CAAC,MAAM,CAAC,EAAE;oBACrC,GAAG,CAAC,OAAO,CAAC,CAAC,EAAE,OAAO,CAAC,CAAC;oBACxB,KAAK,GAAG,IAAI,CAAC,wBAAwB,CAAC,OAAO,CAAC,CAAC,EAAE,OAAO,CAAC,CAAC,EAAE,OAAO,CAAC,CAAC;oBACrE,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,IAAI,EAAE,EAAE,EAAE,EAAE,OAAO,CAAC,OAAO,CAAC,OAAO,CAAC,EAAE,GAAG,GAAG,GAAG,EAAE,KAAK,GAAG,KAAK,EAAE;wBACpF,GAAG,CAAC,OAAO,CAAC,EAAE,EAAE,OAAO,CAAC,CAAC;wBACzB,KAAK,GAAG,IAAI,CAAC,wBAAwB,CAAC,OAAO,CAAC,CAAC,EAAE,OAAO,CAAC,CAAC,EAAE,OAAO,CAAC,CAAC;wBACrE,GAAG,GAAG,EAAE,GAAG,GAAG,CAAC;wBACf,sFAAsF;wBACtF,oFAAoF;wBACpF,IAAI,CAAC,KAAK,GAAG,KAAK,CAAC,KAAK,CAAC;4BACvB,IAAI,CAAC,sBAAsB,CAAC,GAAG,EAAE,OAAO,EAAE,GAAG,EAAE,OAAO,EAAE,GAAG,EAAE,GAAG,EAAE,OAAO,EAAE,GAAG,EAAE,OAAO,EAAE,GAAG,EAAE,QAAQ,CAAC,CAAC;qBAC3G;iBACF;aACF;SACF;IACH,CAAC;IACD,6DAA6D;IAC7C,kBAAkB,CAAC,GAAiB;QAClD,IAAI,IAAI,CAAC,UAAU,YAAY,2BAAY,EAAE;YAC3C,MAAM,GAAG,GAAG,IAAI,CAAC,UAAU,CAAC;YAC5B,IAAI,CAAC,2BAA2B,CAAC,GAAG,EAAE,GAAG,EAAE,KAAK,CAAC,CAAC;SACnD;aAAM,IAAI,IAAI,CAAC,UAAU,YAAY,6BAAa,EAAE;YACnD,IAAI,CAAC,wBAAwB,CAAC,IAAI,CAAC,UAAU,EAAE,GAAG,EAAE,IAAI,CAAC,CAAC;SAC3D;aAAM,IAAI,IAAI,CAAC,UAAU,YAAY,aAAK,EAAE;YAC3C,IAAI,CAAC,oBAAoB,CAAC,IAAI,CAAC,UAAU,EAAE,GAAG,EAAE,IAAI,CAAC,CAAC;SACvD;aAAM,IAAI,IAAI,CAAC,UAAU,YAAY,6BAAc,EAAE;YACpD,IAAI,CAAC,8BAA8B,CAAC,GAAG,EAAE,IAAI,CAAC,UAAU,EAAE,KAAK,CAAC,CAAC;SAClE;aAAM,IAAI,IAAI,CAAC,UAAU,YAAY,iCAAe,EAAE;YACrD,IAAI,CAAC,uBAAuB,CAAC,GAAG,EAAE,IAAI,CAAC,kBAAkB,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC;SACvE;aAAM,IAAI,IAAI,CAAC,UAAU,YAAY,yDAA2B,EAAE;YACjE,IAAI,CAAC,mCAAmC,CAAC,GAAG,EAAE,IAAI,CAAC,kBAAkB,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC;SACnF;QACD,OAAO,SAAS,CAAC;IACnB,CAAC;IACD,sDAAsD;IACtC,WAAW,CAAC,IAAW;QACrC,IAAI,IAAI,CAAC,UAAU,YAAY,6BAAa,EAAE;YAC5C,IAAI,CAAC,kBAAkB,CACrB,IAAI,CAAC,UAAU,EAAE,IAAI,CAAC,UAAU,CAAC,SAAS,EAAE,GAAG,EAAE,IAAI,CAAC,UAAU,CAAC,SAAS,EAAE,GAAG,EAAE,IAAI,EAAE,IAAI,CAC5F,CAAC;SACH;aAAM,IAAI,IAAI,CAAC,UAAU,YAAY,2BAAY,EAAE;YAClD,IAAI,CAAC,oBAAoB,CAAC,IAAI,EAAE,IAAI,CAAC,UAAU,EAAE,KAAK,CAAC,CAAC;SACzD;aAAM,IAAI,IAAI,CAAC,UAAU,YAAY,aAAK,EAAE;YAC3C,IAAI,CAAC,cAAc,CAAC,IAAI,EAAE,IAAI,CAAC,UAAU,EAAE,KAAK,CAAC,CAAC;SACnD;aAAM,IAAI,IAAI,CAAC,UAAU,YAAY,6BAAc,EAAE;YACpD,IAAI,CAAC,yBAAyB,CAAC,IAAI,EAAE,IAAI,CAAC,UAAU,EAAE,KAAK,CAAC,CAAC;SAC9D;aAAM,IAAI,IAAI,CAAC,UAAU,YAAY,iCAAe,EAAE;YACrD,IAAI,CAAC,uBAAuB,CAAC,IAAI,EAAE,IAAI,CAAC,WAAW,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC;SACjE;aAAM,IAAI,IAAI,CAAC,UAAU,YAAY,yDAA2B,EAAE;YACjE,IAAI,CAAC,mCAAmC,CAAC,IAAI,EAAE,IAAI,CAAC,WAAW,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC;SAC7E;QACD,OAAO,SAAS,CAAC;IACnB,CAAC;IACD,gEAAgE;IAChD,oBAAoB,CAAC,KAAqB;QACxD,IAAI,IAAI,CAAC,UAAU,YAAY,6BAAa,EAAE;YAC5C,IAAI,CAAC,2BAA2B,CAAC,IAAI,CAAC,UAAU,EAAE,KAAK,EAAE,IAAI,CAAC,CAAC;SAChE;aAAM,IAAI,IAAI,CAAC,UAAU,YAAY,2BAAY,EAAE;YAClD,IAAI,CAAC,8BAA8B,CAAC,IAAI,CAAC,UAAU,EAAE,KAAK,EAAE,IAAI,CAAC,CAAC;SACnE;aAAM,IAAI,IAAI,CAAC,UAAU,YAAY,aAAK,EAAE;YAC3C,IAAI,CAAC,yBAAyB,CAAC,IAAI,CAAC,UAAU,EAAE,KAAK,EAAE,IAAI,CAAC,CAAC;SAC9D;aAAM,IAAI,IAAI,CAAC,UAAU,YAAY,iCAAkB,EAAE;YACxD,IAAI,CAAC,oCAAoC,CAAC,KAAK,EAAE,IAAI,CAAC,UAAU,EAAE,KAAK,CAAC,CAAC;SAC1E;aAAM,IAAI,IAAI,CAAC,UAAU,YAAY,iCAAe,EAAE;YACrD,IAAI,CAAC,uBAAuB,CAAC,KAAK,EAAE,IAAI,CAAC,oBAAoB,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC;SAC3E;aAAM,IAAI,IAAI,CAAC,UAAU,YAAY,yDAA2B,EAAE;YACjE,IAAI,CAAC,mCAAmC,CAAC,KAAK,EAAE,IAAI,CAAC,oBAAoB,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC;SACvF;QACD,OAAO,SAAS,CAAC;IACnB,CAAC;IACD,8EAA8E;IAC9D,iCAAiC,CAAC,KAAkC;QAClF,KAAK,CAAC,iCAAiC,CAAC,KAAK,CAAC,CAAC;QAC/C,kIAAkI;QAClI,IAAI,CAAC,QAAQ,GAAG,yDAA2B,CAAC,+BAA+B,CAAC,IAAI,CAAC,QAAQ,EAAE,CAAC,EAAE,KAAK,EAAE,SAAS,EAAE,IAAI,CAAC,CAAC;IACxH,CAAC;IACD,8EAA8E;IAC9D,qBAAqB,CAAC,MAAuB;QAC3D;;;;;;;;;;;;UAYE;QACF,OAAO,SAAS,CAAC;IACnB,CAAC;;AA30Bc,uCAAa,GAAG,yBAAO,CAAC,MAAM,EAAE,CAAC;AACjC,uCAAa,GAAG,yBAAO,CAAC,MAAM,EAAE,CAAC;AACjC,uCAAa,GAAG,yBAAO,CAAC,MAAM,EAAE,CAAC;AACjC,uCAAa,GAAG,yBAAO,CAAC,MAAM,EAAE,CAAC;AACjC,qCAAW,GAAG,yBAAO,CAAC,MAAM,EAAE,CAAC;AAjCnC,8DAAyB","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 Curve\r\n */\r\n\r\nimport { assert } from \"@itwin/core-bentley\";\r\nimport { BSplineCurve3d, BSplineCurve3dBase } from \"../../bspline/BSplineCurve\";\r\nimport { BSplineCurve3dH } from \"../../bspline/BSplineCurve3dH\";\r\nimport { Geometry } from \"../../Geometry\";\r\nimport { RecurseToCurvesGeometryHandler } from \"../../geometry3d/GeometryHandler\";\r\nimport { GrowableFloat64Array } from \"../../geometry3d/GrowableFloat64Array\";\r\nimport { Point3d, Vector3d } from \"../../geometry3d/Point3dVector3d\";\r\nimport { Range3d } from \"../../geometry3d/Range\";\r\nimport { AnalyticRoots, SmallSystem } from \"../../numerics/Polynomials\";\r\nimport { Arc3d } from \"../Arc3d\";\r\nimport { CurveChainWithDistanceIndex } from \"../CurveChainWithDistanceIndex\";\r\nimport { CurveCollection } from \"../CurveCollection\";\r\nimport { CurveIntervalRole, CurveLocationDetail, CurveLocationDetailPair } from \"../CurveLocationDetail\";\r\nimport { CurvePrimitive } from \"../CurvePrimitive\";\r\nimport { AnyCurve } from \"../CurveTypes\";\r\nimport { LineSegment3d } from \"../LineSegment3d\";\r\nimport { LineString3d } from \"../LineString3d\";\r\n\r\n// cspell:word XYRR\r\n\r\n/**\r\n * Handler class for XY close approach between _geometryB and another geometry.\r\n * * Approach means the XY distance (z is ignored) between _geometryB and another geometry.\r\n * * Closest approach is a measure of the proximity of one curve to another. It's the length of the shortest line\r\n * segment perpendicular to both curves; if the curves intersect, the closest approach is zero. In the context of\r\n * this class, z-coordinates are ignored, so the closest approach is as seen in the top view. If you have coplanar\r\n * input curves and want to find closest approach in their plane, rotate them first into a plane parallel to the\r\n * xy-plane, then afterward, rotate the results back as required.\r\n * * Close approach can also be from a curve endpoint perpendicular to another curve or from a curve endpoint to\r\n * another curve endpoint.\r\n * * Instances are initialized and called from CurveCurve.\r\n * * geometryB is saved for later reference.\r\n * @internal\r\n */\r\nexport class CurveCurveCloseApproachXY extends RecurseToCurvesGeometryHandler {\r\n  private _geometryB: AnyCurve | undefined;\r\n  private _circularArcB: Arc3d | undefined;\r\n  private _circularRadiusB: number | undefined;\r\n  private setGeometryB(geometryB: AnyCurve | undefined) {\r\n    this._geometryB = geometryB;\r\n    this._circularArcB = undefined;\r\n    this._circularRadiusB = undefined;\r\n    if (geometryB instanceof Arc3d) {\r\n      const r = geometryB.circularRadiusXY();\r\n      if (r !== undefined) {\r\n        this._circularRadiusB = r;\r\n        this._circularArcB = geometryB;\r\n      }\r\n    }\r\n  }\r\n  /**\r\n   * Maximum XY distance (z is ignored). Approach larger than this is not interesting.\r\n   * This is caller defined and can be undefined.\r\n   */\r\n  private _maxDistanceToAccept: number | undefined;\r\n  /** Squared max distance. This is private, and is forced to at least small metric distance squared. */\r\n  private _maxDistanceSquared: number;\r\n  /**\r\n   * Start and end points of line segments that meet closest approach criteria, i.e., they are perpendicular to\r\n   * both curves and their length is smaller than _maxDistanceToAccept.\r\n   */\r\n  private _results: CurveLocationDetailPair[];\r\n\r\n  private static _workPointAA0 = Point3d.create();\r\n  private static _workPointAA1 = Point3d.create();\r\n  private static _workPointBB0 = Point3d.create();\r\n  private static _workPointBB1 = Point3d.create();\r\n  private static _workPointB = Point3d.create();\r\n\r\n  /**\r\n   * Constructor.\r\n   * @param geometryB second curve for intersection. Saved for reference by specific handler methods.\r\n   */\r\n  public constructor(geometryB: AnyCurve | undefined) {\r\n    super();\r\n    this.setGeometryB(geometryB);\r\n    this._maxDistanceSquared = Geometry.smallMetricDistanceSquared;\r\n    this._results = [];\r\n  }\r\n  /** Set the (possibly undefined) max XY distance (z is ignored) to accept. */\r\n  public set maxDistanceToAccept(value: number | undefined) {\r\n    this._maxDistanceToAccept = value;\r\n    if (this._maxDistanceToAccept !== undefined && this._maxDistanceToAccept > 0)\r\n      this._maxDistanceSquared = this._maxDistanceToAccept * this._maxDistanceToAccept;\r\n  }\r\n  /** Access the (possibly undefined) max XY distance (z is ignored) to accept. */\r\n  public get maxDistanceToAccept(): number | undefined {\r\n    return this._maxDistanceToAccept;\r\n  }\r\n  /** Ask if the maxDistanceToAccept value is defined and positive */\r\n  public get isMaxDistanceSet(): boolean {\r\n    return this._maxDistanceToAccept !== undefined && this._maxDistanceToAccept > 0;\r\n  }\r\n  /** Reset the geometry and flags, leaving all other parts unchanged (and preserving accumulated intersections) */\r\n  public resetGeometry(geometryB: AnyCurve) {\r\n    this.setGeometryB(geometryB);\r\n  }\r\n  /** returns true if `fraction` is in [0,1] within tolerance */\r\n  private acceptFraction(fraction: number, fractionTol: number = 1.0e-12) {\r\n    if (fraction < -fractionTol)\r\n      return false;\r\n    if (fraction > 1.0 + fractionTol)\r\n      return false;\r\n    return true;\r\n  }\r\n  /**\r\n   * Return the results structure for the intersection calculation, structured as an array of CurveLocationDetailPair.\r\n   * @param reinitialize if true, a new results structure is created for use by later calls.\r\n   */\r\n  public grabPairedResults(reinitialize: boolean = false): CurveLocationDetailPair[] {\r\n    const result = this._results;\r\n    if (reinitialize)\r\n      this._results = [];\r\n    return result;\r\n  }\r\n  /**\r\n   * If distance between pointA and pointB is less than maxDistance, record CurveLocationDetailPair which is\r\n   * the approach from pointA to pointB.\r\n   */\r\n  private testAndRecordPointPairApproach(\r\n    cpA: CurvePrimitive, fA: number, pointA: Point3d, cpB: CurvePrimitive, fB: number, pointB: Point3d, reversed: boolean,\r\n  ): void {\r\n    const d2 = pointA.distanceSquaredXY(pointB);\r\n    if (d2 < this._maxDistanceSquared) {\r\n      const detailA = CurveLocationDetail.createCurveFractionPoint(cpA, fA, pointA);\r\n      const detailB = CurveLocationDetail.createCurveFractionPoint(cpB, fB, pointB);\r\n      const pair = CurveLocationDetailPair.createCapture(detailA, detailB);\r\n      if (reversed)\r\n        pair.swapDetails();\r\n      this._results.push(pair);\r\n    }\r\n  }\r\n  /**\r\n   * Create a close approach pair if XY distance is within maxDistance.\r\n   * @param localFractionA a fraction on first curve\r\n   * @param cpA the first curve\r\n   * @param fractionA0 start of the first curve\r\n   * @param fractionA1 end of the first curve\r\n   * @param localFractionB a fraction on second curve\r\n   * @param cpB the second curve\r\n   * @param fractionB0 start of the second curve\r\n   * @param fractionB1 end of the second curve\r\n   * @param reversed whether to reverse the details in the pair (e.g., so that detailB refers to geometryB).\r\n   * @param intervalDetails optional CurveLocationDetailPair\r\n   */\r\n  private recordPointWithLocalFractions(\r\n    localFractionA: number,\r\n    cpA: CurvePrimitive,\r\n    fractionA0: number,\r\n    fractionA1: number,\r\n    localFractionB: number,\r\n    cpB: CurvePrimitive,\r\n    fractionB0: number,\r\n    fractionB1: number,\r\n    reversed: boolean,\r\n    intervalDetails?: undefined | CurveLocationDetailPair,\r\n  ): void {\r\n    let globalFractionA, globalFractionB;\r\n    let globalFractionA1, globalFractionB1;\r\n    const isInterval = intervalDetails !== undefined &&\r\n      intervalDetails.detailA.hasFraction1 &&\r\n      intervalDetails.detailB.hasFraction1;\r\n    if (isInterval) {\r\n      globalFractionA = Geometry.interpolate(fractionA0, intervalDetails.detailA.fraction, fractionA1);\r\n      globalFractionB = Geometry.interpolate(fractionB0, intervalDetails.detailB.fraction, fractionB1);\r\n      globalFractionA1 = Geometry.interpolate(fractionA0, intervalDetails.detailA.fraction1!, fractionA1);\r\n      globalFractionB1 = Geometry.interpolate(fractionB0, intervalDetails.detailB.fraction1!, fractionB1);\r\n    } else {\r\n      globalFractionA = globalFractionA1 = Geometry.interpolate(fractionA0, localFractionA, fractionA1);\r\n      globalFractionB = globalFractionB1 = Geometry.interpolate(fractionB0, localFractionB, fractionB1);\r\n    }\r\n    // ignore duplicate of most recent approach\r\n    const numPrevious = this._results.length;\r\n    if (numPrevious > 0 && !isInterval) {\r\n      const oldDetailA = this._results[numPrevious - 1].detailA;\r\n      const oldDetailB = this._results[numPrevious - 1].detailB;\r\n      if (reversed) {\r\n        if (oldDetailB.isSameCurveAndFraction({curve: cpA, fraction: globalFractionA}) &&\r\n          oldDetailA.isSameCurveAndFraction({curve: cpB, fraction: globalFractionB}))\r\n          return;\r\n      } else {\r\n        if (oldDetailA.isSameCurveAndFraction({curve: cpA, fraction: globalFractionA}) &&\r\n          oldDetailB.isSameCurveAndFraction({curve: cpB, fraction: globalFractionB}))\r\n          return;\r\n      }\r\n    }\r\n    const detailA = CurveLocationDetail.createCurveFractionPoint(\r\n      cpA, globalFractionA, cpA.fractionToPoint(globalFractionA),\r\n    );\r\n    const detailB = CurveLocationDetail.createCurveFractionPoint(\r\n      cpB, globalFractionB, cpB.fractionToPoint(globalFractionB),\r\n    );\r\n    if (isInterval) {\r\n      detailA.captureFraction1Point1(globalFractionA1, cpA.fractionToPoint(globalFractionA1));\r\n      detailB.captureFraction1Point1(globalFractionB1, cpB.fractionToPoint(globalFractionB1));\r\n    } else {\r\n      const d2 = detailA.point.distanceSquaredXY(detailB.point);\r\n      if (d2 > this._maxDistanceSquared)\r\n        return;\r\n      detailA.setIntervalRole(CurveIntervalRole.isolated);\r\n      detailB.setIntervalRole(CurveIntervalRole.isolated);\r\n    }\r\n    if (reversed) {\r\n      this._results.push(new CurveLocationDetailPair(detailB, detailA));\r\n    } else {\r\n      this._results.push(new CurveLocationDetailPair(detailA, detailB));\r\n    }\r\n  }\r\n  /**\r\n   * Capture a close approach pair that has point and local fraction but not curve.\r\n   * * Record the pair, each detail modified with global fraction and input curve.\r\n   * * Pair is neither modified nor recorded if it would be a duplicate of the last recorded pair.\r\n   * @param pair details computed with local fractions\r\n   * @param cpA curveA\r\n   * @param fractionA0 global start fraction on curveA\r\n   * @param fractionA1 global end fraction on curveA\r\n   * @param cpB curveB\r\n   * @param fractionB0 global start fraction on curveB\r\n   * @param fractionB1 global end fraction on curveB\r\n   * @param reversed whether to reverse the details in the pair (e.g., so that detailB refers to geometryB).\r\n   */\r\n  private capturePairWithLocalFractions(\r\n    pair: CurveLocationDetailPair,\r\n    cpA: CurvePrimitive,\r\n    fractionA0: number,\r\n    fractionA1: number,\r\n    cpB: CurvePrimitive,\r\n    fractionB0: number,\r\n    fractionB1: number,\r\n    reversed: boolean,\r\n  ) {\r\n    const globalFractionA = Geometry.interpolate(fractionA0, pair.detailA.fraction, fractionA1);\r\n    const globalFractionB = Geometry.interpolate(fractionB0, pair.detailB.fraction, fractionB1);\r\n    // ignore duplicate of most recent pair\r\n    const numPrevious = this._results.length;\r\n    if (numPrevious > 0) {\r\n      const oldDetailA = this._results[numPrevious - 1].detailA;\r\n      const oldDetailB = this._results[numPrevious - 1].detailB;\r\n      if (reversed) {\r\n        if (oldDetailB.isSameCurveAndFraction({curve: cpA, fraction: globalFractionA}) &&\r\n          oldDetailA.isSameCurveAndFraction({curve: cpB, fraction: globalFractionB}))\r\n          return;\r\n      } else {\r\n        if (oldDetailA.isSameCurveAndFraction({curve: cpA, fraction: globalFractionA}) &&\r\n          oldDetailB.isSameCurveAndFraction({curve: cpB, fraction: globalFractionB}))\r\n          return;\r\n      }\r\n    }\r\n    // recompute the points just in case\r\n    CurveLocationDetail.createCurveEvaluatedFraction(cpA, globalFractionA, pair.detailA);\r\n    CurveLocationDetail.createCurveEvaluatedFraction(cpB, globalFractionB, pair.detailB);\r\n    pair.detailA.setIntervalRole(CurveIntervalRole.isolated);\r\n    pair.detailB.setIntervalRole(CurveIntervalRole.isolated);\r\n    if (reversed)\r\n      pair.swapDetails();\r\n    this._results.push(pair);\r\n  }\r\n  /**\r\n   * Emit recordPoint for multiple pairs (on full curve) if within maxDistance.\r\n   * @param cpA first curve primitive (possibly different from curve in detailA, but fraction compatible)\r\n   * @param cpB second curve primitive (possibly different from curve in detailA, but fraction compatible)\r\n   * @param pairs array of pairs\r\n   * @param reversed whether to reverse the details in the pair (e.g., so that detailB refers to geometryB).\r\n   */\r\n  public recordPairs(\r\n    cpA: CurvePrimitive, cpB: CurvePrimitive, pairs: CurveLocationDetailPair[] | undefined, reversed: boolean,\r\n  ): void {\r\n    if (pairs !== undefined) {\r\n      for (const p of pairs) {\r\n        this.recordPointWithLocalFractions(\r\n          p.detailA.fraction, cpA, 0, 1, p.detailB.fraction, cpB, 0, 1, reversed, p,\r\n        );\r\n      }\r\n    }\r\n  }\r\n  /**\r\n   * Record fully assembled (but possibly reversed) detail pair.\r\n   * @param detailA first detail\r\n   * @param detailB second detail\r\n   * @param reversed whether to reverse the details in the pair (e.g., so that detailB refers to geometryB).\r\n   */\r\n  public captureDetailPair(\r\n    detailA: CurveLocationDetail | undefined, detailB: CurveLocationDetail | undefined, reversed: boolean,\r\n  ): void {\r\n    if (detailA && detailB) {\r\n      if (reversed) {\r\n        this._results.push(CurveLocationDetailPair.createCapture(detailB, detailA));\r\n      } else {\r\n        this._results.push(CurveLocationDetailPair.createCapture(detailA, detailB));\r\n      }\r\n    }\r\n  }\r\n  private static updatePointToSegmentDistance(\r\n    fractionA: number,\r\n    pointA: Point3d,\r\n    pointB0: Point3d,\r\n    pointB1: Point3d,\r\n    fractionB: number,\r\n    maxDistanceSquared: number,\r\n    closestApproach: CurveLocationDetailPair,   // modified on return\r\n  ): boolean {\r\n    let updated = false;\r\n    if (fractionB < 0)\r\n      fractionB = 0;\r\n    else if (fractionB > 1)\r\n      fractionB = 1;\r\n    this._workPointB = pointB0.interpolate(fractionB, pointB1, this._workPointB);\r\n    const distanceSquared = this._workPointB.distanceSquaredXY(pointA);\r\n    if (distanceSquared <= Math.min(maxDistanceSquared, closestApproach.detailA.a)) {\r\n      closestApproach.detailA.setFP(fractionA, pointA, undefined, distanceSquared);\r\n      closestApproach.detailB.setFP(fractionB, this._workPointB, undefined, distanceSquared);\r\n      updated = true;\r\n    }\r\n    return updated;\r\n  }\r\n  /**\r\n   * Return fractions of close approach within maxDistance between two line segments (a0,a1) and (b0,b1).\r\n   * * Math details can be found at core/geometry/internaldocs/Curve.md\r\n   * @param a0 start point of line a\r\n   * @param a1 end point of line a\r\n   * @param b0 start point of line b\r\n   * @param b1 end point of line b\r\n   * @param maxDistanceSquared maximum distance squared (assumed to be positive)\r\n   * @returns the fractional (not xy) coordinates in result.x and result.y. result.x is fraction on line a.\r\n   * result.y is fraction on line b.\r\n   */\r\n  private static segmentSegmentBoundedApproach(\r\n    a0: Point3d,\r\n    a1: Point3d,\r\n    b0: Point3d,\r\n    b1: Point3d,\r\n    maxDistanceSquared: number,\r\n  ): CurveLocationDetailPair | undefined {\r\n    const ux = a1.x - a0.x;\r\n    const uy = a1.y - a0.y;\r\n    const vx = b1.x - b0.x;\r\n    const vy = b1.y - b0.y;\r\n    const e00x = b0.x - a0.x;\r\n    const e00y = b0.y - a0.y;\r\n    const e01x = b1.x - a0.x;\r\n    const e01y = b1.y - a0.y;\r\n    const e10x = b0.x - a1.x;\r\n    const e10y = b0.y - a1.y;\r\n    const hab0 = Geometry.crossProductXYXY(ux, uy, e00x, e00y);\r\n    const hab1 = Geometry.crossProductXYXY(ux, uy, e01x, e01y);\r\n    const hba0 = -Geometry.crossProductXYXY(vx, vy, e00x, e00y);\r\n    const hba1 = -Geometry.crossProductXYXY(vx, vy, e10x, e10y);\r\n    if (hab0 * hab1 < 0.0 && hba0 * hba1 < 0.0) { // true intersection, strictly within both segments\r\n      const fractionA = -hba0 / (hba1 - hba0);\r\n      const fractionB = -hab0 / (hab1 - hab0);\r\n      return CurveLocationDetailPair.createCapture(\r\n        CurveLocationDetail.createCurveFractionPoint(undefined, fractionA, a0.interpolate(fractionA, a1)),\r\n        CurveLocationDetail.createCurveFractionPoint(undefined, fractionB, b0.interpolate(fractionB, b1)),\r\n      );\r\n    }\r\n    // there's no intersection, so find the closest approach within maxDistance from an endpoint\r\n    const closestApproach = new CurveLocationDetailPair();\r\n    closestApproach.detailA.a = 2 * maxDistanceSquared; // init to an approach that's too far away\r\n    let reversed = false;\r\n    const uu = Geometry.hypotenuseSquaredXY(ux, uy);\r\n    if (hab0 * hab0 <= maxDistanceSquared * uu) { // test distance of b0 to u\r\n      const fractionA = Geometry.dotProductXYXY(ux, uy, e00x, e00y) / uu;\r\n      if (this.updatePointToSegmentDistance(0, b0, a0, a1, fractionA, maxDistanceSquared, closestApproach))\r\n        reversed = true;\r\n    }\r\n    if (hab1 * hab1 <= maxDistanceSquared * uu) { // test distance of b1 to u\r\n      const fractionA = Geometry.dotProductXYXY(ux, uy, e01x, e01y) / uu;\r\n      if (this.updatePointToSegmentDistance(1, b1, a0, a1, fractionA, maxDistanceSquared, closestApproach))\r\n        reversed = true;\r\n    }\r\n    const vv = Geometry.hypotenuseSquaredXY(vx, vy);\r\n    if (hba0 * hba0 <= maxDistanceSquared * vv) { // test distance of a0 to v\r\n      const fractionB = -Geometry.dotProductXYXY(vx, vy, e00x, e00y) / vv;\r\n      if (this.updatePointToSegmentDistance(0, a0, b0, b1, fractionB, maxDistanceSquared, closestApproach))\r\n        reversed = false;\r\n    }\r\n    if (hba1 * hba1 <= maxDistanceSquared * vv) { // test distance of a1 to v\r\n      const fractionB = -Geometry.dotProductXYXY(vx, vy, e10x, e10y) / vv;\r\n      if (this.updatePointToSegmentDistance(1, a1, b0, b1, fractionB, maxDistanceSquared, closestApproach))\r\n        reversed = false;\r\n    }\r\n    if (closestApproach.detailA.a > maxDistanceSquared)\r\n      return undefined;\r\n    if (reversed)\r\n      closestApproach.swapDetails();\r\n    return closestApproach;\r\n  }\r\n  /**\r\n   * Check different combination of fractions on curveA and curveB. If distance between points at 2 fractions\r\n   * is less than maxDistance, record CurveLocationDetailPair which is the approach between the 2 points.\r\n   * Optionally, record close approaches of one curve's points if they fall between the other curve's points.\r\n   * @param cpA curveA\r\n   * @param fA0 fraction0 on curveA\r\n   * @param fA1 fraction1 on curveA\r\n   * @param testProjectionOnA whether to record projections of the given curveB points onto curveA\r\n   * @param cpB curveB\r\n   * @param fB0 fraction0 on curveB\r\n   * @param fB1 fraction0 on curveB\r\n   * @param testProjectionOnB whether to record projections of the given curveA points onto curveB\r\n   * @param reversed whether to reverse the details in the pair (e.g., so that detailB refers to geometryB).\r\n   */\r\n  private testAndRecordFractionalPairApproach(\r\n    cpA: CurvePrimitive,\r\n    fA0: number,\r\n    fA1: number,\r\n    testProjectionOnA: boolean,\r\n    cpB: CurvePrimitive,\r\n    fB0: number,\r\n    fB1: number,\r\n    testProjectionOnB: boolean,\r\n    reversed: boolean,\r\n  ): void {\r\n    const pointA0 = cpA.fractionToPoint(fA0);\r\n    const pointA1 = cpA.fractionToPoint(fA1);\r\n    const pointB0 = cpB.fractionToPoint(fB0);\r\n    const pointB1 = cpB.fractionToPoint(fB1);\r\n    this.testAndRecordPointPairApproach(cpA, fA0, pointA0, cpB, fB0, pointB0, reversed);\r\n    this.testAndRecordPointPairApproach(cpA, fA1, pointA1, cpB, fB0, pointB0, reversed);\r\n    this.testAndRecordPointPairApproach(cpA, fA0, pointA0, cpB, fB1, pointB1, reversed);\r\n    this.testAndRecordPointPairApproach(cpA, fA1, pointA1, cpB, fB1, pointB1, reversed);\r\n    if (testProjectionOnB) {\r\n      this.testAndRecordProjection(cpA, fA0, pointA0, cpB, fB0, fB1, reversed);\r\n      this.testAndRecordProjection(cpA, fA1, pointA1, cpB, fB0, fB1, reversed);\r\n    }\r\n    if (testProjectionOnA) {\r\n      this.testAndRecordProjection(cpB, fB0, pointB0, cpA, fA0, fA1, !reversed);\r\n      this.testAndRecordProjection(cpB, fB1, pointB1, cpA, fA0, fA1, !reversed);\r\n    }\r\n  }\r\n  /** Find the closest approach between pointA and cpB. Add the approach if it's within fB0 and fB1. */\r\n  private testAndRecordProjection(\r\n    cpA: CurvePrimitive, fA: number, pointA: Point3d, cpB: CurvePrimitive, fB0: number, fB1: number, reversed: boolean,\r\n  ) {\r\n    // NO NO NO -- this is 3D closest point --- need 2d !!\r\n    const detail = cpB.closestPoint(pointA, false);\r\n    if (detail) {\r\n      const fB = Geometry.restrictToInterval(detail.fraction, fB0, fB1);\r\n      if (fB === detail.fraction) { // if fraction is within fB0 and fB1\r\n        this.testAndRecordPointPairApproach(cpA, fA, pointA, cpB, detail.fraction, detail.point, reversed);\r\n      }\r\n    }\r\n  }\r\n  /**\r\n   * Compute intersection of two line segments.\r\n   * Filter by extension rules.\r\n   * Record with fraction mapping.\r\n   * * The fraction mappings allow portions of a linestring to be passed here.\r\n   */\r\n  private computeSegmentSegment3D(\r\n    cpA: CurvePrimitive,\r\n    pointA0: Point3d,\r\n    fractionA0: number,\r\n    pointA1: Point3d,\r\n    fractionA1: number,\r\n    cpB: CurvePrimitive,\r\n    pointB0: Point3d,\r\n    fractionB0: number,\r\n    pointB1: Point3d,\r\n    fractionB1: number,\r\n    reversed: boolean,\r\n  ): void {\r\n    // compute a pair with fractions local to segments\r\n    const approach = CurveCurveCloseApproachXY.segmentSegmentBoundedApproach(\r\n      pointA0, pointA1, pointB0, pointB1, this._maxDistanceSquared,\r\n    );\r\n    // adjust the pair to refer to input curves and global fractions, then record it if new\r\n    if (approach) {\r\n      approach.detailA.setCurve(cpA);\r\n      approach.detailB.setCurve(cpB);\r\n      this.capturePairWithLocalFractions(approach, cpA, fractionA0, fractionA1, cpB, fractionB0, fractionB1, reversed);\r\n    }\r\n  }\r\n  /** Low level dispatch of segment with segment. */\r\n  private dispatchSegmentSegment(\r\n    cpA: CurvePrimitive,\r\n    pointA0: Point3d,\r\n    fractionA0: number,\r\n    pointA1: Point3d,\r\n    fractionA1: number,\r\n    cpB: CurvePrimitive,\r\n    pointB0: Point3d,\r\n    fractionB0: number,\r\n    pointB1: Point3d,\r\n    fractionB1: number,\r\n    reversed: boolean,\r\n  ): void {\r\n    this.computeSegmentSegment3D(\r\n      cpA, pointA0, fractionA0, pointA1, fractionA1,\r\n      cpB, pointB0, fractionB0, pointB1, fractionB1,\r\n      reversed,\r\n    );\r\n  }\r\n  /**\r\n   * Low level dispatch of segment with arc.\r\n   * Find close approaches within maxDistance between a line segments (pointA0, pointA1) and an arc.\r\n   * To consider:\r\n   * 1) intersection between arc and segment.\r\n   * 2) arc endpoints to segment endpoints or arc endpoints projection to the segment.\r\n   * 3) line parallel to arc tangent.\r\n   * @param cpA curve A (line segment or line string)\r\n   * @param pointA0 start point of the segment\r\n   * @param fractionA0 fraction of the start of the segment\r\n   * @param pointA1 end point of the segment\r\n   * @param fractionA1 fraction of the end of the segment\r\n   * @param arc the arc\r\n   * @param reversed whether to reverse the details in the pair (e.g., so that detailB refers to geometryB).\r\n   */\r\n  private dispatchSegmentArc(\r\n    cpA: CurvePrimitive,\r\n    pointA0: Point3d,\r\n    fractionA0: number,\r\n    pointA1: Point3d,\r\n    fractionA1: number,\r\n    arc: Arc3d,\r\n    reversed: boolean,\r\n  ): void {\r\n    // 1) intersection between arc and segment\r\n    // Suppose:\r\n    // Arc: X = C + cU + sV where c = cos(theta) and s = sin(theta)\r\n    // Line: contains points A0 and A1\r\n    // The arc intersects the line at point X if det(A0, A1, X) = 0 with homogeneous xyw points and vectors.\r\n    // With equational X:  det(A0, A1, C) + c*det(A0, A1, U) + s*det(A0, A1, V) = 0.\r\n    // solve for theta.\r\n    // evaluate points.\r\n    // project back to line.\r\n    let intersectionFound = false;\r\n    const data = arc.toTransformedVectors();\r\n    const pointA0Local = pointA0;\r\n    const pointA1Local = pointA1;\r\n    const alpha = Geometry.tripleProductXYW(pointA0Local, 1, pointA1Local, 1, data.center, 1); //  det(A0, A1, C)\r\n    const beta = Geometry.tripleProductXYW(pointA0Local, 1, pointA1Local, 1, data.vector0, 0); // det(A0, A1, U)\r\n    const gamma = Geometry.tripleProductXYW(pointA0Local, 1, pointA1Local, 1, data.vector90, 0); // det(A0, A1, V)\r\n    const cosines = new GrowableFloat64Array(2);\r\n    const sines = new GrowableFloat64Array(2);\r\n    const radians = new GrowableFloat64Array(2);\r\n    const numRoots = AnalyticRoots.appendImplicitLineUnitCircleIntersections( // solve the equation\r\n      alpha, beta, gamma, cosines, sines, radians,\r\n    );\r\n    for (let i = 0; i < numRoots; i++) {\r\n      const arcPoint = data.center.plus2Scaled(\r\n        data.vector0, cosines.atUncheckedIndex(i), data.vector90, sines.atUncheckedIndex(i),\r\n      );\r\n      const arcFraction = data.sweep.radiansToSignedPeriodicFraction(radians.atUncheckedIndex(i));\r\n      const lineFraction = SmallSystem.lineSegment3dXYClosestPointUnbounded(pointA0Local, pointA1Local, arcPoint);\r\n      // only add if the point is within the start and end fractions of both line segment and arc\r\n      if (lineFraction !== undefined && this.acceptFraction(lineFraction) && this.acceptFraction(arcFraction)) {\r\n        this.recordPointWithLocalFractions(\r\n          lineFraction, cpA, fractionA0, fractionA1, arcFraction, arc, 0, 1, reversed,\r\n        );\r\n        intersectionFound = true;\r\n      }\r\n    }\r\n    if (intersectionFound)\r\n      return;\r\n    // 2) endpoints to endpoints or endpoints projection to the other curve\r\n    this.testAndRecordFractionalPairApproach(cpA, fractionA0, fractionA1, true, arc, 0, 1, false, reversed);\r\n    // 3) line parallel to arc tangent.\r\n    // If line does not intersect the arc, then the closest (and/or the furthest) point on arc to the line is a\r\n    // point where the tangent line on arc at that point is parallel to the line.\r\n    const dotUT = data.vector0.crossProductStartEndXY(pointA0, pointA1);\r\n    const dotVT = data.vector90.crossProductStartEndXY(pointA0, pointA1);\r\n    const parallelRadians = Math.atan2(dotVT, dotUT);\r\n    for (const radians1 of [parallelRadians, parallelRadians + Math.PI]) {\r\n      const arcPoint = data.center.plus2Scaled(data.vector0, Math.cos(radians1), data.vector90, Math.sin(radians1));\r\n      const arcFraction = data.sweep.radiansToSignedPeriodicFraction(radians1);\r\n      const lineFraction = SmallSystem.lineSegment3dXYClosestPointUnbounded(pointA0Local, pointA1Local, arcPoint);\r\n      // only add if the point is within the start and end fractions of both line segment and arc\r\n      if (lineFraction !== undefined && this.acceptFraction(lineFraction) && this.acceptFraction(arcFraction)) {\r\n        this.recordPointWithLocalFractions(\r\n          lineFraction, cpA, fractionA0, fractionA1, arcFraction, arc, 0, 1, reversed,\r\n        );\r\n      }\r\n    }\r\n  }\r\n  /** Low level dispatch of circular arc with circular arc. radiusA must be larger than or equal to radiusB. */\r\n  private dispatchCircularCircularOrdered(\r\n    cpA: Arc3d, radiusA: number, cpB: Arc3d, radiusB: number, reversed: boolean,\r\n  ): void {\r\n    const c = cpA.center.distance(cpB.center);\r\n    const e = this._maxDistanceToAccept !== undefined ? this._maxDistanceToAccept : Geometry.smallMetricDistance;\r\n    if (c > radiusA + radiusB + e) // distance between circles is more than max distance\r\n      return;\r\n    // TODO: 1) intersection between arcs\r\n    // 2) endpoints to endpoints\r\n    this.testAndRecordFractionalPairApproach(cpA, 0, 1, false, cpB, 0, 1, false, reversed);\r\n    // 3) line from one arc to another (perpendicular to arc tangents along center-center line)\r\n    if (!Geometry.isSmallMetricDistance(c)) {\r\n      const vectorAB = Vector3d.createStartEnd(cpA.center, cpB.center);\r\n      vectorAB.scaleInPlace(1.0 / c);\r\n      for (const rA of [-radiusA, radiusA]) {\r\n        for (const rB of [-radiusB, radiusB]) {\r\n          const tangentDistance = c - rA + rB;\r\n          if (tangentDistance < e) {\r\n            const detailA = this.resolveDirectionToArcXYFraction(cpA, vectorAB, rA);\r\n            if (detailA) {\r\n              const detailB = this.resolveDirectionToArcXYFraction(cpB, vectorAB, rB);\r\n              if (detailB)\r\n                this.captureDetailPair(detailA, detailB, reversed);\r\n            }\r\n          }\r\n        }\r\n      }\r\n    }\r\n  }\r\n  /** Find the fractional point (if any) on the circular `arc` in the direction of `radialVector`. */\r\n  private resolveDirectionToArcXYFraction(\r\n    arc: Arc3d, radialVector: Vector3d, scale: number,\r\n  ): CurveLocationDetail | undefined {\r\n    // The scale ultimately only affects the direction --- easiest way to use it is two multiplies.\r\n    const c = scale * arc.matrixRef.columnDotXYZ(0, radialVector.x, radialVector.y, 0);\r\n    const s = scale * arc.matrixRef.columnDotXYZ(1, radialVector.x, radialVector.y, 0);\r\n    const radians = Math.atan2(s, c);\r\n    const fraction = arc.sweep.radiansToPositivePeriodicFraction(radians, 0);\r\n    if (fraction < 1.0)\r\n      return CurveLocationDetail.createCurveEvaluatedFraction(arc, fraction);\r\n    return undefined;\r\n  }\r\n  /** Low level dispatch of arc with arc. Only circular arcs are supported. */\r\n  private dispatchArcArc(cpA: Arc3d, cpB: Arc3d, reversed: boolean): void {\r\n    const rangeA = cpA.range();\r\n    const rangeB = cpB.range();\r\n    rangeA.expandInPlace(this._maxDistanceToAccept!);\r\n    if (!rangeB.intersectsRangeXY(rangeA))\r\n      return;\r\n    if (this._circularArcB) {\r\n      const radiusB = this._circularRadiusB!;\r\n      const radiusA = cpA.circularRadiusXY();\r\n      if (radiusA !== undefined) {\r\n        if (radiusA >= radiusB)\r\n          this.dispatchCircularCircularOrdered(cpA, radiusA, cpB, radiusB, reversed);\r\n        else\r\n          this.dispatchCircularCircularOrdered(cpB, radiusB, cpA, radiusA, !reversed);\r\n        return;\r\n      }\r\n    }\r\n  }\r\n  /** Low level dispatch of arc with (beziers of) a bspline curve */\r\n  private dispatchArcBsplineCurve3d(cpA: Arc3d, cpB: BSplineCurve3d, reversed: boolean): void {\r\n    const ls = LineString3d.create();\r\n    cpB.emitStrokes(ls);\r\n    this.computeArcLineString(cpA, ls, reversed);\r\n  }\r\n  /** Low level dispatch of (beziers of) a bspline curve with (beziers of) a bspline curve */\r\n  private dispatchBSplineCurve3dBSplineCurve3d(\r\n    bcurveA: BSplineCurve3dBase, bcurveB: BSplineCurve3dBase, reversed: boolean,\r\n  ): void {\r\n    const lsA = LineString3d.create();\r\n    bcurveA.emitStrokes(lsA);\r\n    const lsB = LineString3d.create();\r\n    bcurveB.emitStrokes(lsB);\r\n    this.computeLineStringLineString(lsA, lsB, reversed);\r\n  }\r\n  /** Low level dispatch of linestring with (beziers of) a bspline curve */\r\n  public dispatchLineStringBSplineCurve(lsA: LineString3d, curveB: BSplineCurve3d, reversed: boolean): any {\r\n    const lsB = LineString3d.create();\r\n    curveB.emitStrokes(lsB);\r\n    this.computeLineStringLineString(lsA, lsB, reversed);\r\n  }\r\n  /** Low level dispatch of segment with (beziers of) a bspline curve */\r\n  public dispatchSegmentBsplineCurve(segA: LineSegment3d, curveB: BSplineCurve3d, reversed: boolean): any {\r\n    const lsB = LineString3d.create();\r\n    curveB.emitStrokes(lsB);\r\n    this.computeSegmentLineString(segA, lsB, reversed);\r\n  }\r\n  /** Detail computation for segment approaching linestring. */\r\n  public computeSegmentLineString(segA: LineSegment3d, lsB: LineString3d, reversed: boolean): void {\r\n    const numB = lsB.numPoints();\r\n    const deltaFracB = Geometry.safeDivideFraction(1, numB - 1, 0);\r\n    const pointA0 = segA.point0Ref;\r\n    const pointA1 = segA.point1Ref;\r\n    const pointB0 = CurveCurveCloseApproachXY._workPointBB0;\r\n    const pointB1 = CurveCurveCloseApproachXY._workPointBB1;\r\n    for (let i = 0; i < numB - 1; ++i) {\r\n      const fB0 = i * deltaFracB; // global linestring fractions\r\n      const fB1 = (i + 1 === numB - 1) ? 1.0 : (i + 1) * deltaFracB;  // make sure we nail the end fraction\r\n      lsB.packedPoints.getPoint3dAtUncheckedPointIndex(i, pointB0);\r\n      lsB.packedPoints.getPoint3dAtUncheckedPointIndex(i + 1, pointB1);\r\n      this.dispatchSegmentSegment(segA, pointA0, 0.0, pointA1, 1.0, lsB, pointB0, fB0, pointB1, fB1, reversed);\r\n    }\r\n  }\r\n  /** Detail computation for arc approaching linestring. */\r\n  public computeArcLineString(arcA: Arc3d, lsB: LineString3d, reversed: boolean): any {\r\n    const rangeA = arcA.range();\r\n    const rangeB = lsB.range();\r\n    rangeA.expandInPlace(this._maxDistanceToAccept!);\r\n    if (!rangeB.intersectsRangeXY(rangeA))\r\n      return;\r\n    const pointB0 = CurveCurveCloseApproachXY._workPointBB0;\r\n    const pointB1 = CurveCurveCloseApproachXY._workPointBB1;\r\n    const numB = lsB.numPoints();\r\n    if (numB > 1) {\r\n      const dfB = 1.0 / (numB - 1);\r\n      let fB0;\r\n      let fB1;\r\n      fB0 = 0.0;\r\n      lsB.pointAt(0, pointB0);\r\n      for (let ib = 1; ib < numB; ib++, pointB0.setFrom(pointB1), fB0 = fB1) {\r\n        lsB.pointAt(ib, pointB1);\r\n        fB1 = ib * dfB;\r\n        this.dispatchSegmentArc(lsB, pointB0, fB0, pointB1, fB1, arcA, !reversed);\r\n      }\r\n    }\r\n    return undefined;\r\n  }\r\n  /** Low level dispatch of curve collection. */\r\n  private dispatchCurveCollection(geomA: AnyCurve, geomAHandler: (geomA: any) => any): void {\r\n    const geomB = this._geometryB;  // save\r\n    if (!geomB || !geomB.children || !(geomB instanceof CurveCollection))\r\n      return;\r\n    for (const child of geomB.children) {\r\n      this.resetGeometry(child);\r\n      geomAHandler(geomA);\r\n    }\r\n    this._geometryB = geomB;  // restore\r\n  }\r\n  /** Low level dispatch to geomA given a CurveChainWithDistanceIndex in geometryB. */\r\n  private dispatchCurveChainWithDistanceIndex(geomA: AnyCurve, geomAHandler: (geomA: any) => any): void {\r\n    if (!this._geometryB || !(this._geometryB instanceof CurveChainWithDistanceIndex))\r\n      return;\r\n    if (geomA instanceof CurveChainWithDistanceIndex) {\r\n      assert(!!\"call handleCurveChainWithDistanceIndex(geomA) instead\");\r\n      return;\r\n    }\r\n    const index0 = this._results.length;\r\n    const geomB = this._geometryB;  // save\r\n    for (const child of geomB.path.children) {\r\n      this.resetGeometry(child);\r\n      geomAHandler(geomA);\r\n    }\r\n    this.resetGeometry(geomB);  // restore\r\n    this._results = CurveChainWithDistanceIndex.convertChildDetailToChainDetail(this._results, index0, undefined, geomB, true);\r\n  }\r\n  /** Double dispatch handler for strongly typed segment. */\r\n  public override handleLineSegment3d(segmentA: LineSegment3d): any {\r\n    if (this._geometryB instanceof LineSegment3d) {\r\n      const segmentB = this._geometryB;\r\n      this.dispatchSegmentSegment(\r\n        segmentA, segmentA.point0Ref, 0.0, segmentA.point1Ref, 1.0,\r\n        segmentB, segmentB.point0Ref, 0.0, segmentB.point1Ref, 1.0,\r\n        false,\r\n      );\r\n    } else if (this._geometryB instanceof LineString3d) {\r\n      this.computeSegmentLineString(segmentA, this._geometryB, false);\r\n    } else if (this._geometryB instanceof Arc3d) {\r\n      this.dispatchSegmentArc(segmentA, segmentA.point0Ref, 0.0, segmentA.point1Ref, 1.0, this._geometryB, false);\r\n    } else if (this._geometryB instanceof BSplineCurve3d) {\r\n      this.dispatchSegmentBsplineCurve(segmentA, this._geometryB, false);\r\n    } else if (this._geometryB instanceof CurveCollection) {\r\n      this.dispatchCurveCollection(segmentA, this.handleLineSegment3d.bind(this));\r\n    } else if (this._geometryB instanceof CurveChainWithDistanceIndex) {\r\n      this.dispatchCurveChainWithDistanceIndex(segmentA, this.handleLineSegment3d.bind(this));\r\n    }\r\n    return undefined;\r\n  }\r\n  /**\r\n   * Set bits for comparison to range xy\r\n   * * bit 0x01 => x smaller than range.low.x\r\n   * * bit 0x02 => x larger than range.high.x\r\n   * * bit 0x04 => y smaller than range.low.y\r\n   * * bit 0x08 => y larger than range.high.y\r\n   * * If we divide XY plane into 9 areas using the range, the function returns 0 for points\r\n   * inside the range. Below is other binary numbers returned by the function for all 9 areas:\r\n   *   1001 | 1000 | 1010\r\n   *   ------------------\r\n   *    1   |  0   |  10\r\n   *   ------------------\r\n   *   101  | 100  | 110\r\n   * @param xy point to test\r\n   * @param range range for comparison\r\n   */\r\n  private classifyBitsPointRangeXY(x: number, y: number, range: Range3d): number {\r\n    let result = 0;\r\n    if (x < range.low.x)\r\n      result = 0x01;\r\n    else if (x > range.high.x)\r\n      result = 0x02;\r\n    // note the OR operation\r\n    if (y < range.low.y)\r\n      result |= 0x04;\r\n    else if (y > range.high.y)\r\n      result |= 0x08;\r\n    return result;\r\n  }\r\n  /** Low level dispatch of line string with line string. */\r\n  private computeLineStringLineString(lsA: LineString3d, lsB: LineString3d, reversed: boolean): void {\r\n    const rangeA = lsA.range();\r\n    const rangeB = lsB.range();\r\n    rangeA.expandInPlace(this._maxDistanceToAccept!);\r\n    if (!rangeB.intersectsRangeXY(rangeA))\r\n      return;\r\n    let bitB0: number;\r\n    let bitB1: number;\r\n    const rangeA1 = Range3d.createNull();\r\n    const pointA0 = CurveCurveCloseApproachXY._workPointAA0;\r\n    const pointA1 = CurveCurveCloseApproachXY._workPointAA1;\r\n    const pointB0 = CurveCurveCloseApproachXY._workPointBB0;\r\n    const pointB1 = CurveCurveCloseApproachXY._workPointBB1;\r\n    const numA = lsA.numPoints();\r\n    const numB = lsB.numPoints();\r\n    if (numA > 1 && numB > 1) {\r\n      const dfA = 1.0 / (numA - 1);\r\n      const dfB = 1.0 / (numB - 1);\r\n      let fA0 = 0.0;\r\n      let fA1, fB0, fB1;\r\n      lsA.pointAt(0, pointA0);\r\n      for (let ia = 1; ia < numA; ia++, pointA0.setFrom(pointA1), fA0 = fA1) {\r\n        fA1 = ia * dfA;\r\n        fB0 = 0.0;\r\n        lsA.pointAt(ia, pointA1);\r\n        // rangeA1 is around line segment [A0,A1] expanded by max distance\r\n        rangeA1.setNull();\r\n        rangeA1.extendPoint(pointA0);\r\n        rangeA1.extendPoint(pointA1);\r\n        rangeA1.expandInPlace(this._maxDistanceToAccept!);\r\n        if (rangeA1.intersectsRangeXY(rangeB)) {\r\n          lsB.pointAt(0, pointB0);\r\n          bitB0 = this.classifyBitsPointRangeXY(pointB0.x, pointB0.y, rangeA1);\r\n          for (let ib = 1; ib < numB; ib++, pointB0.setFrom(pointB1), fB0 = fB1, bitB0 = bitB1) {\r\n            lsB.pointAt(ib, pointB1);\r\n            bitB1 = this.classifyBitsPointRangeXY(pointB1.x, pointB1.y, rangeA1);\r\n            fB1 = ib * dfB;\r\n            // DO NOT study the segment in detail if both bitB bits are on for any of the 4 planes\r\n            // (i.e., no intersection between rangeA1 and the range around line segment [B0,B1])\r\n            if ((bitB0 & bitB1) === 0)\r\n              this.dispatchSegmentSegment(lsA, pointA0, fA0, pointA1, fA1, lsB, pointB0, fB0, pointB1, fB1, reversed);\r\n          }\r\n        }\r\n      }\r\n    }\r\n  }\r\n  /** Double dispatch handler for strongly typed linestring. */\r\n  public override handleLineString3d(lsA: LineString3d): any {\r\n    if (this._geometryB instanceof LineString3d) {\r\n      const lsB = this._geometryB;\r\n      this.computeLineStringLineString(lsA, lsB, false);\r\n    } else if (this._geometryB instanceof LineSegment3d) {\r\n      this.computeSegmentLineString(this._geometryB, lsA, true);\r\n    } else if (this._geometryB instanceof Arc3d) {\r\n      this.computeArcLineString(this._geometryB, lsA, true);\r\n    } else if (this._geometryB instanceof BSplineCurve3d) {\r\n      this.dispatchLineStringBSplineCurve(lsA, this._geometryB, false);\r\n    } else if (this._geometryB instanceof CurveCollection) {\r\n      this.dispatchCurveCollection(lsA, this.handleLineString3d.bind(this));\r\n    } else if (this._geometryB instanceof CurveChainWithDistanceIndex) {\r\n      this.dispatchCurveChainWithDistanceIndex(lsA, this.handleLineString3d.bind(this));\r\n    }\r\n    return undefined;\r\n  }\r\n  /** Double dispatch handler for strongly typed arc. */\r\n  public override handleArc3d(arc0: Arc3d): any {\r\n    if (this._geometryB instanceof LineSegment3d) {\r\n      this.dispatchSegmentArc(\r\n        this._geometryB, this._geometryB.point0Ref, 0.0, this._geometryB.point1Ref, 1.0, arc0, true,\r\n      );\r\n    } else if (this._geometryB instanceof LineString3d) {\r\n      this.computeArcLineString(arc0, this._geometryB, false);\r\n    } else if (this._geometryB instanceof Arc3d) {\r\n      this.dispatchArcArc(arc0, this._geometryB, false);\r\n    } else if (this._geometryB instanceof BSplineCurve3d) {\r\n      this.dispatchArcBsplineCurve3d(arc0, this._geometryB, false);\r\n    } else if (this._geometryB instanceof CurveCollection) {\r\n      this.dispatchCurveCollection(arc0, this.handleArc3d.bind(this));\r\n    } else if (this._geometryB instanceof CurveChainWithDistanceIndex) {\r\n      this.dispatchCurveChainWithDistanceIndex(arc0, this.handleArc3d.bind(this));\r\n    }\r\n    return undefined;\r\n  }\r\n  /** Double dispatch handler for strongly typed bspline curve. */\r\n  public override handleBSplineCurve3d(curve: BSplineCurve3d): any {\r\n    if (this._geometryB instanceof LineSegment3d) {\r\n      this.dispatchSegmentBsplineCurve(this._geometryB, curve, true);\r\n    } else if (this._geometryB instanceof LineString3d) {\r\n      this.dispatchLineStringBSplineCurve(this._geometryB, curve, true);\r\n    } else if (this._geometryB instanceof Arc3d) {\r\n      this.dispatchArcBsplineCurve3d(this._geometryB, curve, true);\r\n    } else if (this._geometryB instanceof BSplineCurve3dBase) {\r\n      this.dispatchBSplineCurve3dBSplineCurve3d(curve, this._geometryB, false);\r\n    } else if (this._geometryB instanceof CurveCollection) {\r\n      this.dispatchCurveCollection(curve, this.handleBSplineCurve3d.bind(this));\r\n    } else if (this._geometryB instanceof CurveChainWithDistanceIndex) {\r\n      this.dispatchCurveChainWithDistanceIndex(curve, this.handleBSplineCurve3d.bind(this));\r\n    }\r\n    return undefined;\r\n  }\r\n  /** Double dispatch handler for strongly typed CurveChainWithDistanceIndex. */\r\n  public override handleCurveChainWithDistanceIndex(chain: CurveChainWithDistanceIndex): any {\r\n    super.handleCurveChainWithDistanceIndex(chain);\r\n    // if _geometryB is also a CurveChainWithDistanceIndex, it will already have been converted by dispatchCurveChainWithDistanceIndex\r\n    this._results = CurveChainWithDistanceIndex.convertChildDetailToChainDetail(this._results, 0, chain, undefined, true);\r\n  }\r\n  /** Double dispatch handler for strongly typed homogeneous bspline curve .. */\r\n  public override handleBSplineCurve3dH(_curve: BSplineCurve3dH): any {\r\n    /*\r\n    //NEEDS WORK -- make \"dispatch\" methods tolerant of both 3d and 3dH.\r\n    // \"easy\" if both present BezierCurve3dH span loaders\r\n    if (this._geometryB instanceof LineSegment3d) {\r\n      this.dispatchSegmentBsplineCurve(\r\n        this._geometryB, this._extendB, this._geometryB.point0Ref, 0.0, this._geometryB.point1Ref, 1.0, this._extendB,\r\n        curve, this._extendA, true);\r\n    } else if (this._geometryB instanceof LineString3d) {\r\n      this.dispatchLineStringBSplineCurve(this._geometryB, this._extendB, curve, this._extendA, true);\r\n    } else if (this._geometryB instanceof Arc3d) {\r\n      this.dispatchArcBsplineCurve3d(this._geometryB, this._extendB, curve, this._extendA, true);\r\n    }\r\n    */\r\n    return undefined;\r\n  }\r\n}\r\n"]}
+{"version":3,"file":"CurveCurveCloseApproachXY.js","sourceRoot":"","sources":["../../../../src/curve/internalContexts/CurveCurveCloseApproachXY.ts"],"names":[],"mappings":";AAAA;;;+FAG+F;;;AAE/F;;GAEG;AAEH,sDAA6C;AAC7C,6DAAgF;AAEhF,6CAA0C;AAC1C,sEAAkF;AAClF,gFAA6E;AAC7E,sEAA2D;AAC3D,kDAAiD;AACjD,kDAG+B;AAC/B,4DAAwE;AACxE,oCAAiC;AACjC,gFAA6E;AAC7E,wDAAqD;AACrD,gEAAyG;AAGzG,oDAAiD;AACjD,kDAA+C;AAE/C,+BAA+B;AAE/B;;;;;;;;;;;;;GAaG;AACH,MAAa,yBAA0B,SAAQ,gDAA8B;IAInE,YAAY,CAAC,SAA+B;QAClD,IAAI,CAAC,UAAU,GAAG,SAAS,CAAC;QAC5B,IAAI,CAAC,aAAa,GAAG,SAAS,CAAC;QAC/B,IAAI,CAAC,gBAAgB,GAAG,SAAS,CAAC;QAClC,IAAI,SAAS,YAAY,aAAK,EAAE;YAC9B,MAAM,CAAC,GAAG,SAAS,CAAC,gBAAgB,EAAE,CAAC;YACvC,IAAI,CAAC,KAAK,SAAS,EAAE;gBACnB,IAAI,CAAC,gBAAgB,GAAG,CAAC,CAAC;gBAC1B,IAAI,CAAC,aAAa,GAAG,SAAS,CAAC;aAChC;SACF;IACH,CAAC;IAoBD;;;OAGG;IACH,YAAmB,SAA+B;QAChD,KAAK,EAAE,CAAC;QACR,IAAI,CAAC,YAAY,CAAC,SAAS,CAAC,CAAC;QAC7B,IAAI,CAAC,mBAAmB,GAAG,mBAAQ,CAAC,0BAA0B,CAAC;QAC/D,IAAI,CAAC,QAAQ,GAAG,EAAE,CAAC;IACrB,CAAC;IACD,6EAA6E;IAC7E,IAAW,mBAAmB,CAAC,KAAyB;QACtD,IAAI,CAAC,oBAAoB,GAAG,KAAK,CAAC;QAClC,IAAI,IAAI,CAAC,oBAAoB,KAAK,SAAS,IAAI,IAAI,CAAC,oBAAoB,GAAG,CAAC;YAC1E,IAAI,CAAC,mBAAmB,GAAG,IAAI,CAAC,oBAAoB,GAAG,IAAI,CAAC,oBAAoB,CAAC;IACrF,CAAC;IACD,gFAAgF;IAChF,IAAW,mBAAmB;QAC5B,OAAO,IAAI,CAAC,oBAAoB,CAAC;IACnC,CAAC;IACD,mEAAmE;IACnE,IAAW,gBAAgB;QACzB,OAAO,IAAI,CAAC,oBAAoB,KAAK,SAAS,IAAI,IAAI,CAAC,oBAAoB,GAAG,CAAC,CAAC;IAClF,CAAC;IACD,iHAAiH;IAC1G,aAAa,CAAC,SAAmB;QACtC,IAAI,CAAC,YAAY,CAAC,SAAS,CAAC,CAAC;IAC/B,CAAC;IACD,8DAA8D;IACtD,cAAc,CAAC,QAAgB,EAAE,cAAsB,OAAO;QACpE,IAAI,QAAQ,GAAG,CAAC,WAAW;YACzB,OAAO,KAAK,CAAC;QACf,IAAI,QAAQ,GAAG,GAAG,GAAG,WAAW;YAC9B,OAAO,KAAK,CAAC;QACf,OAAO,IAAI,CAAC;IACd,CAAC;IACD;;;OAGG;IACI,iBAAiB,CAAC,eAAwB,KAAK;QACpD,MAAM,MAAM,GAAG,IAAI,CAAC,QAAQ,CAAC;QAC7B,IAAI,YAAY;YACd,IAAI,CAAC,QAAQ,GAAG,EAAE,CAAC;QACrB,OAAO,MAAM,CAAC;IAChB,CAAC;IACD;;;OAGG;IACK,8BAA8B,CACpC,GAAmB,EAAE,EAAU,EAAE,MAAe,EAAE,GAAmB,EAAE,EAAU,EAAE,MAAe,EAAE,QAAiB;QAErH,MAAM,EAAE,GAAG,MAAM,CAAC,iBAAiB,CAAC,MAAM,CAAC,CAAC;QAC5C,IAAI,EAAE,GAAG,IAAI,CAAC,mBAAmB,EAAE;YACjC,MAAM,OAAO,GAAG,yCAAmB,CAAC,wBAAwB,CAAC,GAAG,EAAE,EAAE,EAAE,MAAM,CAAC,CAAC;YAC9E,MAAM,OAAO,GAAG,yCAAmB,CAAC,wBAAwB,CAAC,GAAG,EAAE,EAAE,EAAE,MAAM,CAAC,CAAC;YAC9E,MAAM,IAAI,GAAG,6CAAuB,CAAC,aAAa,CAAC,OAAO,EAAE,OAAO,CAAC,CAAC;YACrE,IAAI,QAAQ;gBACV,IAAI,CAAC,WAAW,EAAE,CAAC;YACrB,IAAI,CAAC,QAAQ,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC;SAC1B;IACH,CAAC;IACD;;;;;;;;;;;;OAYG;IACK,6BAA6B,CACnC,cAAsB,EACtB,GAAmB,EACnB,UAAkB,EAClB,UAAkB,EAClB,cAAsB,EACtB,GAAmB,EACnB,UAAkB,EAClB,UAAkB,EAClB,QAAiB,EACjB,eAAqD;QAErD,IAAI,eAAe,EAAE,eAAe,CAAC;QACrC,IAAI,gBAAgB,EAAE,gBAAgB,CAAC;QACvC,MAAM,UAAU,GAAG,eAAe,KAAK,SAAS;YAC9C,eAAe,CAAC,OAAO,CAAC,YAAY;YACpC,eAAe,CAAC,OAAO,CAAC,YAAY,CAAC;QACvC,IAAI,UAAU,EAAE;YACd,eAAe,GAAG,mBAAQ,CAAC,WAAW,CAAC,UAAU,EAAE,eAAe,CAAC,OAAO,CAAC,QAAQ,EAAE,UAAU,CAAC,CAAC;YACjG,eAAe,GAAG,mBAAQ,CAAC,WAAW,CAAC,UAAU,EAAE,eAAe,CAAC,OAAO,CAAC,QAAQ,EAAE,UAAU,CAAC,CAAC;YACjG,gBAAgB,GAAG,mBAAQ,CAAC,WAAW,CAAC,UAAU,EAAE,eAAe,CAAC,OAAO,CAAC,SAAU,EAAE,UAAU,CAAC,CAAC;YACpG,gBAAgB,GAAG,mBAAQ,CAAC,WAAW,CAAC,UAAU,EAAE,eAAe,CAAC,OAAO,CAAC,SAAU,EAAE,UAAU,CAAC,CAAC;SACrG;aAAM;YACL,eAAe,GAAG,gBAAgB,GAAG,mBAAQ,CAAC,WAAW,CAAC,UAAU,EAAE,cAAc,EAAE,UAAU,CAAC,CAAC;YAClG,eAAe,GAAG,gBAAgB,GAAG,mBAAQ,CAAC,WAAW,CAAC,UAAU,EAAE,cAAc,EAAE,UAAU,CAAC,CAAC;SACnG;QACD,2CAA2C;QAC3C,MAAM,WAAW,GAAG,IAAI,CAAC,QAAQ,CAAC,MAAM,CAAC;QACzC,IAAI,WAAW,GAAG,CAAC,IAAI,CAAC,UAAU,EAAE;YAClC,MAAM,UAAU,GAAG,IAAI,CAAC,QAAQ,CAAC,WAAW,GAAG,CAAC,CAAC,CAAC,OAAO,CAAC;YAC1D,MAAM,UAAU,GAAG,IAAI,CAAC,QAAQ,CAAC,WAAW,GAAG,CAAC,CAAC,CAAC,OAAO,CAAC;YAC1D,IAAI,QAAQ,EAAE;gBACZ,IAAI,UAAU,CAAC,sBAAsB,CAAC,EAAE,KAAK,EAAE,GAAG,EAAE,QAAQ,EAAE,eAAe,EAAE,CAAC;oBAC9E,UAAU,CAAC,sBAAsB,CAAC,EAAE,KAAK,EAAE,GAAG,EAAE,QAAQ,EAAE,eAAe,EAAE,CAAC;oBAC5E,OAAO;aACV;iBAAM;gBACL,IAAI,UAAU,CAAC,sBAAsB,CAAC,EAAE,KAAK,EAAE,GAAG,EAAE,QAAQ,EAAE,eAAe,EAAE,CAAC;oBAC9E,UAAU,CAAC,sBAAsB,CAAC,EAAE,KAAK,EAAE,GAAG,EAAE,QAAQ,EAAE,eAAe,EAAE,CAAC;oBAC5E,OAAO;aACV;SACF;QACD,MAAM,OAAO,GAAG,yCAAmB,CAAC,wBAAwB,CAC1D,GAAG,EAAE,eAAe,EAAE,GAAG,CAAC,eAAe,CAAC,eAAe,CAAC,CAC3D,CAAC;QACF,MAAM,OAAO,GAAG,yCAAmB,CAAC,wBAAwB,CAC1D,GAAG,EAAE,eAAe,EAAE,GAAG,CAAC,eAAe,CAAC,eAAe,CAAC,CAC3D,CAAC;QACF,IAAI,UAAU,EAAE;YACd,OAAO,CAAC,sBAAsB,CAAC,gBAAgB,EAAE,GAAG,CAAC,eAAe,CAAC,gBAAgB,CAAC,CAAC,CAAC;YACxF,OAAO,CAAC,sBAAsB,CAAC,gBAAgB,EAAE,GAAG,CAAC,eAAe,CAAC,gBAAgB,CAAC,CAAC,CAAC;SACzF;aAAM;YACL,MAAM,EAAE,GAAG,OAAO,CAAC,KAAK,CAAC,iBAAiB,CAAC,OAAO,CAAC,KAAK,CAAC,CAAC;YAC1D,IAAI,EAAE,GAAG,IAAI,CAAC,mBAAmB;gBAC/B,OAAO;YACT,OAAO,CAAC,eAAe,CAAC,uCAAiB,CAAC,QAAQ,CAAC,CAAC;YACpD,OAAO,CAAC,eAAe,CAAC,uCAAiB,CAAC,QAAQ,CAAC,CAAC;SACrD;QACD,IAAI,QAAQ,EAAE;YACZ,IAAI,CAAC,QAAQ,CAAC,IAAI,CAAC,IAAI,6CAAuB,CAAC,OAAO,EAAE,OAAO,CAAC,CAAC,CAAC;SACnE;aAAM;YACL,IAAI,CAAC,QAAQ,CAAC,IAAI,CAAC,IAAI,6CAAuB,CAAC,OAAO,EAAE,OAAO,CAAC,CAAC,CAAC;SACnE;IACH,CAAC;IACD;;;;;;;;;;;;OAYG;IACK,6BAA6B,CACnC,IAA6B,EAC7B,GAAmB,EACnB,UAAkB,EAClB,UAAkB,EAClB,GAAmB,EACnB,UAAkB,EAClB,UAAkB,EAClB,QAAiB;QAEjB,MAAM,eAAe,GAAG,mBAAQ,CAAC,WAAW,CAAC,UAAU,EAAE,IAAI,CAAC,OAAO,CAAC,QAAQ,EAAE,UAAU,CAAC,CAAC;QAC5F,MAAM,eAAe,GAAG,mBAAQ,CAAC,WAAW,CAAC,UAAU,EAAE,IAAI,CAAC,OAAO,CAAC,QAAQ,EAAE,UAAU,CAAC,CAAC;QAC5F,uCAAuC;QACvC,MAAM,WAAW,GAAG,IAAI,CAAC,QAAQ,CAAC,MAAM,CAAC;QACzC,IAAI,WAAW,GAAG,CAAC,EAAE;YACnB,MAAM,UAAU,GAAG,IAAI,CAAC,QAAQ,CAAC,WAAW,GAAG,CAAC,CAAC,CAAC,OAAO,CAAC;YAC1D,MAAM,UAAU,GAAG,IAAI,CAAC,QAAQ,CAAC,WAAW,GAAG,CAAC,CAAC,CAAC,OAAO,CAAC;YAC1D,IAAI,QAAQ,EAAE;gBACZ,IAAI,UAAU,CAAC,sBAAsB,CAAC,EAAE,KAAK,EAAE,GAAG,EAAE,QAAQ,EAAE,eAAe,EAAE,CAAC;oBAC9E,UAAU,CAAC,sBAAsB,CAAC,EAAE,KAAK,EAAE,GAAG,EAAE,QAAQ,EAAE,eAAe,EAAE,CAAC;oBAC5E,OAAO;aACV;iBAAM;gBACL,IAAI,UAAU,CAAC,sBAAsB,CAAC,EAAE,KAAK,EAAE,GAAG,EAAE,QAAQ,EAAE,eAAe,EAAE,CAAC;oBAC9E,UAAU,CAAC,sBAAsB,CAAC,EAAE,KAAK,EAAE,GAAG,EAAE,QAAQ,EAAE,eAAe,EAAE,CAAC;oBAC5E,OAAO;aACV;SACF;QACD,oCAAoC;QACpC,yCAAmB,CAAC,4BAA4B,CAAC,GAAG,EAAE,eAAe,EAAE,IAAI,CAAC,OAAO,CAAC,CAAC;QACrF,yCAAmB,CAAC,4BAA4B,CAAC,GAAG,EAAE,eAAe,EAAE,IAAI,CAAC,OAAO,CAAC,CAAC;QACrF,IAAI,CAAC,OAAO,CAAC,eAAe,CAAC,uCAAiB,CAAC,QAAQ,CAAC,CAAC;QACzD,IAAI,CAAC,OAAO,CAAC,eAAe,CAAC,uCAAiB,CAAC,QAAQ,CAAC,CAAC;QACzD,IAAI,QAAQ;YACV,IAAI,CAAC,WAAW,EAAE,CAAC;QACrB,IAAI,CAAC,QAAQ,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC;IAC3B,CAAC;IACD;;;;;;OAMG;IACI,WAAW,CAChB,GAAmB,EAAE,GAAmB,EAAE,KAA4C,EAAE,QAAiB;QAEzG,IAAI,KAAK,KAAK,SAAS,EAAE;YACvB,KAAK,MAAM,CAAC,IAAI,KAAK,EAAE;gBACrB,IAAI,CAAC,6BAA6B,CAChC,CAAC,CAAC,OAAO,CAAC,QAAQ,EAAE,GAAG,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,OAAO,CAAC,QAAQ,EAAE,GAAG,EAAE,CAAC,EAAE,CAAC,EAAE,QAAQ,EAAE,CAAC,CAC1E,CAAC;aACH;SACF;IACH,CAAC;IACD;;;;;OAKG;IACI,iBAAiB,CACtB,OAAwC,EAAE,OAAwC,EAAE,QAAiB;QAErG,IAAI,OAAO,IAAI,OAAO,EAAE;YACtB,IAAI,QAAQ,EAAE;gBACZ,IAAI,CAAC,QAAQ,CAAC,IAAI,CAAC,6CAAuB,CAAC,aAAa,CAAC,OAAO,EAAE,OAAO,CAAC,CAAC,CAAC;aAC7E;iBAAM;gBACL,IAAI,CAAC,QAAQ,CAAC,IAAI,CAAC,6CAAuB,CAAC,aAAa,CAAC,OAAO,EAAE,OAAO,CAAC,CAAC,CAAC;aAC7E;SACF;IACH,CAAC;IACO,MAAM,CAAC,4BAA4B,CACzC,SAAiB,EACjB,MAAe,EACf,OAAgB,EAChB,OAAgB,EAChB,SAAiB,EACjB,kBAA0B,EAC1B,eAAwC;QAExC,IAAI,OAAO,GAAG,KAAK,CAAC;QACpB,IAAI,SAAS,GAAG,CAAC;YACf,SAAS,GAAG,CAAC,CAAC;aACX,IAAI,SAAS,GAAG,CAAC;YACpB,SAAS,GAAG,CAAC,CAAC;QAChB,IAAI,CAAC,WAAW,GAAG,OAAO,CAAC,WAAW,CAAC,SAAS,EAAE,OAAO,EAAE,IAAI,CAAC,WAAW,CAAC,CAAC;QAC7E,MAAM,eAAe,GAAG,IAAI,CAAC,WAAW,CAAC,iBAAiB,CAAC,MAAM,CAAC,CAAC;QACnE,IAAI,eAAe,IAAI,IAAI,CAAC,GAAG,CAAC,kBAAkB,EAAE,eAAe,CAAC,OAAO,CAAC,CAAC,CAAC,EAAE;YAC9E,eAAe,CAAC,OAAO,CAAC,KAAK,CAAC,SAAS,EAAE,MAAM,EAAE,SAAS,EAAE,eAAe,CAAC,CAAC;YAC7E,eAAe,CAAC,OAAO,CAAC,KAAK,CAAC,SAAS,EAAE,IAAI,CAAC,WAAW,EAAE,SAAS,EAAE,eAAe,CAAC,CAAC;YACvF,OAAO,GAAG,IAAI,CAAC;SAChB;QACD,OAAO,OAAO,CAAC;IACjB,CAAC;IACD;;;;;;;;;;OAUG;IACK,MAAM,CAAC,6BAA6B,CAC1C,EAAW,EACX,EAAW,EACX,EAAW,EACX,EAAW,EACX,kBAA0B;QAE1B,MAAM,EAAE,GAAG,EAAE,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC;QACvB,MAAM,EAAE,GAAG,EAAE,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC;QACvB,MAAM,EAAE,GAAG,EAAE,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC;QACvB,MAAM,EAAE,GAAG,EAAE,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC;QACvB,MAAM,IAAI,GAAG,EAAE,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC;QACzB,MAAM,IAAI,GAAG,EAAE,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC;QACzB,MAAM,IAAI,GAAG,EAAE,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC;QACzB,MAAM,IAAI,GAAG,EAAE,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC;QACzB,MAAM,IAAI,GAAG,EAAE,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC;QACzB,MAAM,IAAI,GAAG,EAAE,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC;QACzB,MAAM,IAAI,GAAG,mBAAQ,CAAC,gBAAgB,CAAC,EAAE,EAAE,EAAE,EAAE,IAAI,EAAE,IAAI,CAAC,CAAC;QAC3D,MAAM,IAAI,GAAG,mBAAQ,CAAC,gBAAgB,CAAC,EAAE,EAAE,EAAE,EAAE,IAAI,EAAE,IAAI,CAAC,CAAC;QAC3D,MAAM,IAAI,GAAG,CAAC,mBAAQ,CAAC,gBAAgB,CAAC,EAAE,EAAE,EAAE,EAAE,IAAI,EAAE,IAAI,CAAC,CAAC;QAC5D,MAAM,IAAI,GAAG,CAAC,mBAAQ,CAAC,gBAAgB,CAAC,EAAE,EAAE,EAAE,EAAE,IAAI,EAAE,IAAI,CAAC,CAAC;QAC5D,IAAI,IAAI,GAAG,IAAI,GAAG,GAAG,IAAI,IAAI,GAAG,IAAI,GAAG,GAAG,EAAE,EAAE,mDAAmD;YAC/F,MAAM,SAAS,GAAG,CAAC,IAAI,GAAG,CAAC,IAAI,GAAG,IAAI,CAAC,CAAC;YACxC,MAAM,SAAS,GAAG,CAAC,IAAI,GAAG,CAAC,IAAI,GAAG,IAAI,CAAC,CAAC;YACxC,OAAO,6CAAuB,CAAC,aAAa,CAC1C,yCAAmB,CAAC,wBAAwB,CAAC,SAAS,EAAE,SAAS,EAAE,EAAE,CAAC,WAAW,CAAC,SAAS,EAAE,EAAE,CAAC,CAAC,EACjG,yCAAmB,CAAC,wBAAwB,CAAC,SAAS,EAAE,SAAS,EAAE,EAAE,CAAC,WAAW,CAAC,SAAS,EAAE,EAAE,CAAC,CAAC,CAClG,CAAC;SACH;QACD,4FAA4F;QAC5F,MAAM,eAAe,GAAG,IAAI,6CAAuB,EAAE,CAAC;QACtD,eAAe,CAAC,OAAO,CAAC,CAAC,GAAG,CAAC,GAAG,kBAAkB,CAAC,CAAC,0CAA0C;QAC9F,IAAI,QAAQ,GAAG,KAAK,CAAC;QACrB,MAAM,EAAE,GAAG,mBAAQ,CAAC,mBAAmB,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC;QAChD,IAAI,IAAI,GAAG,IAAI,IAAI,kBAAkB,GAAG,EAAE,EAAE,EAAE,2BAA2B;YACvE,MAAM,SAAS,GAAG,mBAAQ,CAAC,cAAc,CAAC,EAAE,EAAE,EAAE,EAAE,IAAI,EAAE,IAAI,CAAC,GAAG,EAAE,CAAC;YACnE,IAAI,IAAI,CAAC,4BAA4B,CAAC,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,SAAS,EAAE,kBAAkB,EAAE,eAAe,CAAC;gBAClG,QAAQ,GAAG,IAAI,CAAC;SACnB;QACD,IAAI,IAAI,GAAG,IAAI,IAAI,kBAAkB,GAAG,EAAE,EAAE,EAAE,2BAA2B;YACvE,MAAM,SAAS,GAAG,mBAAQ,CAAC,cAAc,CAAC,EAAE,EAAE,EAAE,EAAE,IAAI,EAAE,IAAI,CAAC,GAAG,EAAE,CAAC;YACnE,IAAI,IAAI,CAAC,4BAA4B,CAAC,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,SAAS,EAAE,kBAAkB,EAAE,eAAe,CAAC;gBAClG,QAAQ,GAAG,IAAI,CAAC;SACnB;QACD,MAAM,EAAE,GAAG,mBAAQ,CAAC,mBAAmB,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC;QAChD,IAAI,IAAI,GAAG,IAAI,IAAI,kBAAkB,GAAG,EAAE,EAAE,EAAE,2BAA2B;YACvE,MAAM,SAAS,GAAG,CAAC,mBAAQ,CAAC,cAAc,CAAC,EAAE,EAAE,EAAE,EAAE,IAAI,EAAE,IAAI,CAAC,GAAG,EAAE,CAAC;YACpE,IAAI,IAAI,CAAC,4BAA4B,CAAC,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,SAAS,EAAE,kBAAkB,EAAE,eAAe,CAAC;gBAClG,QAAQ,GAAG,KAAK,CAAC;SACpB;QACD,IAAI,IAAI,GAAG,IAAI,IAAI,kBAAkB,GAAG,EAAE,EAAE,EAAE,2BAA2B;YACvE,MAAM,SAAS,GAAG,CAAC,mBAAQ,CAAC,cAAc,CAAC,EAAE,EAAE,EAAE,EAAE,IAAI,EAAE,IAAI,CAAC,GAAG,EAAE,CAAC;YACpE,IAAI,IAAI,CAAC,4BAA4B,CAAC,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,SAAS,EAAE,kBAAkB,EAAE,eAAe,CAAC;gBAClG,QAAQ,GAAG,KAAK,CAAC;SACpB;QACD,IAAI,eAAe,CAAC,OAAO,CAAC,CAAC,GAAG,kBAAkB;YAChD,OAAO,SAAS,CAAC;QACnB,IAAI,QAAQ;YACV,eAAe,CAAC,WAAW,EAAE,CAAC;QAChC,OAAO,eAAe,CAAC;IACzB,CAAC;IACD;;;;;;;;;;;;;OAaG;IACK,mCAAmC,CACzC,GAAmB,EACnB,GAAW,EACX,GAAW,EACX,iBAA0B,EAC1B,GAAmB,EACnB,GAAW,EACX,GAAW,EACX,iBAA0B,EAC1B,QAAiB;QAEjB,MAAM,OAAO,GAAG,GAAG,CAAC,eAAe,CAAC,GAAG,CAAC,CAAC;QACzC,MAAM,OAAO,GAAG,GAAG,CAAC,eAAe,CAAC,GAAG,CAAC,CAAC;QACzC,MAAM,OAAO,GAAG,GAAG,CAAC,eAAe,CAAC,GAAG,CAAC,CAAC;QACzC,MAAM,OAAO,GAAG,GAAG,CAAC,eAAe,CAAC,GAAG,CAAC,CAAC;QACzC,IAAI,CAAC,8BAA8B,CAAC,GAAG,EAAE,GAAG,EAAE,OAAO,EAAE,GAAG,EAAE,GAAG,EAAE,OAAO,EAAE,QAAQ,CAAC,CAAC;QACpF,IAAI,CAAC,8BAA8B,CAAC,GAAG,EAAE,GAAG,EAAE,OAAO,EAAE,GAAG,EAAE,GAAG,EAAE,OAAO,EAAE,QAAQ,CAAC,CAAC;QACpF,IAAI,CAAC,8BAA8B,CAAC,GAAG,EAAE,GAAG,EAAE,OAAO,EAAE,GAAG,EAAE,GAAG,EAAE,OAAO,EAAE,QAAQ,CAAC,CAAC;QACpF,IAAI,CAAC,8BAA8B,CAAC,GAAG,EAAE,GAAG,EAAE,OAAO,EAAE,GAAG,EAAE,GAAG,EAAE,OAAO,EAAE,QAAQ,CAAC,CAAC;QACpF,IAAI,iBAAiB,EAAE;YACrB,IAAI,CAAC,uBAAuB,CAAC,GAAG,EAAE,GAAG,EAAE,OAAO,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,QAAQ,CAAC,CAAC;YACzE,IAAI,CAAC,uBAAuB,CAAC,GAAG,EAAE,GAAG,EAAE,OAAO,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,QAAQ,CAAC,CAAC;SAC1E;QACD,IAAI,iBAAiB,EAAE;YACrB,IAAI,CAAC,uBAAuB,CAAC,GAAG,EAAE,GAAG,EAAE,OAAO,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC,QAAQ,CAAC,CAAC;YAC1E,IAAI,CAAC,uBAAuB,CAAC,GAAG,EAAE,GAAG,EAAE,OAAO,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC,QAAQ,CAAC,CAAC;SAC3E;IACH,CAAC;IACD;;;;OAIG;IACK,oCAAoC,CAAC,MAAsB,EAAE,MAAe;QAClF,IAAI,CAAC,CAAC,MAAM,YAAY,aAAK,CAAC,IAAI,CAAC,CAAC,MAAM,YAAY,6BAAa,CAAC,EAAE;YACpE,IAAA,qBAAM,EAAC,CAAC,CAAC,0EAA0E,CAAC,CAAC;YACrF,OAAO,SAAS,CAAC;SAClB;QACD,MAAM,KAAK,GAAG,CAAC,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,CAAC,CAAC,CAAC,iEAAiE;QACrG,MAAM,eAAe,GAAG,IAAI,uCAA8B,CAAC,MAAM,EAAE,MAAM,CAAC,CAAC;QAC3E,MAAM,cAAc,GAAG,IAAI,0BAAiB,CAAC,eAAe,CAAC,CAAC;QAC9D,IAAI,sBAAsB,GAAG,mBAAQ,CAAC,qBAAqB,CAAC;QAC5D,IAAI,iBAAqC,CAAC;QAC1C,IAAI,SAA8B,CAAC;QACnC,KAAK,MAAM,IAAI,IAAI,KAAK,EAAE;YACxB,cAAc,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC;YAC1B,IAAI,cAAc,CAAC,aAAa,EAAE,EAAE;gBAClC,MAAM,cAAc,GAAG,cAAc,CAAC,IAAI,EAAE,CAAC;gBAC7C,IAAI,IAAI,CAAC,cAAc,CAAC,cAAc,CAAC,EAAE;oBACvC,MAAM,MAAM,GAAG,MAAM,CAAC,eAAe,CAAC,cAAc,CAAC,CAAC;oBACtD,MAAM,mBAAmB,GAAG,MAAM,CAAC,iBAAiB,CAAC,MAAM,CAAC,CAAC;oBAC7D,IAAI,mBAAmB,GAAG,sBAAsB,EAAE;wBAChD,sBAAsB,GAAG,mBAAmB,CAAC;wBAC7C,iBAAiB,GAAG,cAAc,CAAC;wBACnC,SAAS,GAAG,MAAM,CAAC;qBACpB;iBACF;aACF;SACF;QACD,IAAI,iBAAiB,IAAI,SAAS;YAChC,OAAO,yCAAmB,CAAC,wBAAwB,CAAC,MAAM,EAAE,iBAAiB,EAAE,SAAS,CAAC,CAAC;QAC5F,OAAO,SAAS,CAAC;IACnB,CAAC;IACD,qGAAqG;IAC7F,uBAAuB,CAC7B,GAAmB,EAAE,EAAU,EAAE,MAAe,EAAE,GAAmB,EAAE,GAAW,EAAE,GAAW,EAAE,QAAiB;QAElH,MAAM,MAAM,GAAG,IAAI,CAAC,oCAAoC,CAAC,GAAG,EAAE,MAAM,CAAC,CAAC;QACtE,IAAI,MAAM,EAAE;YACV,MAAM,EAAE,GAAG,mBAAQ,CAAC,kBAAkB,CAAC,MAAM,CAAC,QAAQ,EAAE,GAAG,EAAE,GAAG,CAAC,CAAC;YAClE,IAAI,EAAE,KAAK,MAAM,CAAC,QAAQ,EAAE,EAAE,oCAAoC;gBAChE,IAAI,CAAC,8BAA8B,CAAC,GAAG,EAAE,EAAE,EAAE,MAAM,EAAE,GAAG,EAAE,MAAM,CAAC,QAAQ,EAAE,MAAM,CAAC,KAAK,EAAE,QAAQ,CAAC,CAAC;aACpG;SACF;IACH,CAAC;IACD;;;;;OAKG;IACK,uBAAuB,CAC7B,GAAmB,EACnB,OAAgB,EAChB,UAAkB,EAClB,OAAgB,EAChB,UAAkB,EAClB,GAAmB,EACnB,OAAgB,EAChB,UAAkB,EAClB,OAAgB,EAChB,UAAkB,EAClB,QAAiB;QAEjB,kDAAkD;QAClD,MAAM,QAAQ,GAAG,yBAAyB,CAAC,6BAA6B,CACtE,OAAO,EAAE,OAAO,EAAE,OAAO,EAAE,OAAO,EAAE,IAAI,CAAC,mBAAmB,CAC7D,CAAC;QACF,uFAAuF;QACvF,IAAI,QAAQ,EAAE;YACZ,QAAQ,CAAC,OAAO,CAAC,QAAQ,CAAC,GAAG,CAAC,CAAC;YAC/B,QAAQ,CAAC,OAAO,CAAC,QAAQ,CAAC,GAAG,CAAC,CAAC;YAC/B,IAAI,CAAC,6BAA6B,CAAC,QAAQ,EAAE,GAAG,EAAE,UAAU,EAAE,UAAU,EAAE,GAAG,EAAE,UAAU,EAAE,UAAU,EAAE,QAAQ,CAAC,CAAC;SAClH;IACH,CAAC;IACD,kDAAkD;IAC1C,sBAAsB,CAC5B,GAAmB,EACnB,OAAgB,EAChB,UAAkB,EAClB,OAAgB,EAChB,UAAkB,EAClB,GAAmB,EACnB,OAAgB,EAChB,UAAkB,EAClB,OAAgB,EAChB,UAAkB,EAClB,QAAiB;QAEjB,IAAI,CAAC,uBAAuB,CAC1B,GAAG,EAAE,OAAO,EAAE,UAAU,EAAE,OAAO,EAAE,UAAU,EAC7C,GAAG,EAAE,OAAO,EAAE,UAAU,EAAE,OAAO,EAAE,UAAU,EAC7C,QAAQ,CACT,CAAC;IACJ,CAAC;IACD;;;;;;;;;;;;;;OAcG;IACK,kBAAkB,CACxB,GAAmB,EACnB,OAAgB,EAChB,UAAkB,EAClB,OAAgB,EAChB,UAAkB,EAClB,GAAU,EACV,QAAiB;QAEjB,4DAA4D;QAC5D,WAAW;QACX,+DAA+D;QAC/D,kCAAkC;QAClC,wGAAwG;QACxG,gFAAgF;QAChF,mBAAmB;QACnB,mBAAmB;QACnB,wBAAwB;QACxB,IAAI,iBAAiB,GAAG,KAAK,CAAC;QAC9B,MAAM,IAAI,GAAG,GAAG,CAAC,oBAAoB,EAAE,CAAC;QACxC,MAAM,YAAY,GAAG,OAAO,CAAC;QAC7B,MAAM,YAAY,GAAG,OAAO,CAAC;QAC7B,MAAM,KAAK,GAAG,mBAAQ,CAAC,gBAAgB,CAAC,YAAY,EAAE,CAAC,EAAE,YAAY,EAAE,CAAC,EAAE,IAAI,CAAC,MAAM,EAAE,CAAC,CAAC,CAAC,CAAC,kBAAkB;QAC7G,MAAM,IAAI,GAAG,mBAAQ,CAAC,gBAAgB,CAAC,YAAY,EAAE,CAAC,EAAE,YAAY,EAAE,CAAC,EAAE,IAAI,CAAC,OAAO,EAAE,CAAC,CAAC,CAAC,CAAC,iBAAiB;QAC5G,MAAM,KAAK,GAAG,mBAAQ,CAAC,gBAAgB,CAAC,YAAY,EAAE,CAAC,EAAE,YAAY,EAAE,CAAC,EAAE,IAAI,CAAC,QAAQ,EAAE,CAAC,CAAC,CAAC,CAAC,iBAAiB;QAC9G,MAAM,OAAO,GAAG,IAAI,2CAAoB,CAAC,CAAC,CAAC,CAAC;QAC5C,MAAM,KAAK,GAAG,IAAI,2CAAoB,CAAC,CAAC,CAAC,CAAC;QAC1C,MAAM,OAAO,GAAG,IAAI,2CAAoB,CAAC,CAAC,CAAC,CAAC;QAC5C,MAAM,QAAQ,GAAG,2BAAa,CAAC,yCAAyC,CAAE,qBAAqB;QAC7F,KAAK,EAAE,IAAI,EAAE,KAAK,EAAE,OAAO,EAAE,KAAK,EAAE,OAAO,CAC5C,CAAC;QACF,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,QAAQ,EAAE,CAAC,EAAE,EAAE;YACjC,MAAM,QAAQ,GAAG,IAAI,CAAC,MAAM,CAAC,WAAW,CACtC,IAAI,CAAC,OAAO,EAAE,OAAO,CAAC,gBAAgB,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,QAAQ,EAAE,KAAK,CAAC,gBAAgB,CAAC,CAAC,CAAC,CACpF,CAAC;YACF,MAAM,WAAW,GAAG,IAAI,CAAC,KAAK,CAAC,+BAA+B,CAAC,OAAO,CAAC,gBAAgB,CAAC,CAAC,CAAC,CAAC,CAAC;YAC5F,MAAM,YAAY,GAAG,yBAAW,CAAC,oCAAoC,CAAC,YAAY,EAAE,YAAY,EAAE,QAAQ,CAAC,CAAC;YAC5G,2FAA2F;YAC3F,IAAI,YAAY,KAAK,SAAS,IAAI,IAAI,CAAC,cAAc,CAAC,YAAY,CAAC,IAAI,IAAI,CAAC,cAAc,CAAC,WAAW,CAAC,EAAE;gBACvG,IAAI,CAAC,6BAA6B,CAChC,YAAY,EAAE,GAAG,EAAE,UAAU,EAAE,UAAU,EAAE,WAAW,EAAE,GAAG,EAAE,CAAC,EAAE,CAAC,EAAE,QAAQ,CAC5E,CAAC;gBACF,iBAAiB,GAAG,IAAI,CAAC;aAC1B;SACF;QACD,IAAI,iBAAiB;YACnB,OAAO;QACT,wEAAwE;QACxE,IAAI,CAAC,mCAAmC,CAAC,GAAG,EAAE,UAAU,EAAE,UAAU,EAAE,IAAI,EAAE,GAAG,EAAE,CAAC,EAAE,CAAC,EAAE,IAAI,EAAE,QAAQ,CAAC,CAAC;QACvG,uDAAuD;QACvD,2GAA2G;QAC3G,6EAA6E;QAC7E,MAAM,KAAK,GAAG,IAAI,CAAC,OAAO,CAAC,sBAAsB,CAAC,OAAO,EAAE,OAAO,CAAC,CAAC;QACpE,MAAM,KAAK,GAAG,IAAI,CAAC,QAAQ,CAAC,sBAAsB,CAAC,OAAO,EAAE,OAAO,CAAC,CAAC;QACrE,MAAM,eAAe,GAAG,IAAI,CAAC,KAAK,CAAC,KAAK,EAAE,KAAK,CAAC,CAAC;QACjD,KAAK,MAAM,QAAQ,IAAI,CAAC,eAAe,EAAE,eAAe,GAAG,IAAI,CAAC,EAAE,CAAC,EAAE;YACnE,MAAM,QAAQ,GAAG,IAAI,CAAC,MAAM,CAAC,WAAW,CAAC,IAAI,CAAC,OAAO,EAAE,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,EAAE,IAAI,CAAC,QAAQ,EAAE,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,CAAC,CAAC;YAC9G,MAAM,WAAW,GAAG,IAAI,CAAC,KAAK,CAAC,+BAA+B,CAAC,QAAQ,CAAC,CAAC;YACzE,MAAM,YAAY,GAAG,yBAAW,CAAC,oCAAoC,CAAC,YAAY,EAAE,YAAY,EAAE,QAAQ,CAAC,CAAC;YAC5G,2FAA2F;YAC3F,IAAI,YAAY,KAAK,SAAS,IAAI,IAAI,CAAC,cAAc,CAAC,YAAY,CAAC,IAAI,IAAI,CAAC,cAAc,CAAC,WAAW,CAAC,EAAE;gBACvG,IAAI,CAAC,6BAA6B,CAChC,YAAY,EAAE,GAAG,EAAE,UAAU,EAAE,UAAU,EAAE,WAAW,EAAE,GAAG,EAAE,CAAC,EAAE,CAAC,EAAE,QAAQ,CAC5E,CAAC;aACH;SACF;IACH,CAAC;IACD,6DAA6D;IACrD,iBAAiB,CACvB,MAAa,EAAE,MAAa,EAAE,QAAiB,EAAE,eAAuC;QAExF,MAAM,MAAM,GAAG,CAAC,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,CAAC,CAAC,CAAE,8DAA8D;QACpG,MAAM,MAAM,GAAG,CAAC,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,CAAC,CAAC;QACpC,MAAM,cAAc,GAAG,IAAI,wCAA+B,CAAC,eAAe,CAAC,CAAC;QAC5E,KAAK,MAAM,KAAK,IAAI,MAAM,EAAE;YAC1B,KAAK,MAAM,KAAK,IAAI,MAAM,EAAE;gBAC1B,cAAc,CAAC,KAAK,CAAC,KAAK,EAAE,KAAK,CAAC,CAAC;gBACnC,IAAI,cAAc,CAAC,aAAa,EAAE,EAAE;oBAClC,MAAM,cAAc,GAAG,cAAc,CAAC,IAAI,EAAE,CAAC;oBAC7C,MAAM,cAAc,GAAG,cAAc,CAAC,IAAI,EAAE,CAAC;oBAC7C,IAAI,IAAI,CAAC,cAAc,CAAC,cAAc,CAAC,IAAI,IAAI,CAAC,cAAc,CAAC,cAAc,CAAC,EAAE;wBAC9E,IAAI,CAAC,6BAA6B,CAChC,cAAc,EAAE,MAAM,EAAE,CAAC,EAAE,CAAC,EAAE,cAAc,EAAE,MAAM,EAAE,CAAC,EAAE,CAAC,EAAE,QAAQ,CACrE,CAAC;qBACH;iBACF;aACF;SACF;IACH,CAAC;IACD,wDAAwD;IAChD,0BAA0B,CAAC,MAAa,EAAE,MAAa,EAAE,QAAiB;QAChF,MAAM,eAAe,GAAG,IAAI,yCAAgC,CAAC,MAAM,EAAE,MAAM,CAAC,CAAC;QAC7E,IAAI,CAAC,iBAAiB,CAAC,MAAM,EAAE,MAAM,EAAE,QAAQ,EAAE,eAAe,CAAC,CAAC;IACpE,CAAC;IACD,4CAA4C;IACpC,cAAc,CAAC,GAAU,EAAE,GAAU,EAAE,QAAiB;QAC9D,MAAM,MAAM,GAAG,GAAG,CAAC,KAAK,EAAE,CAAC;QAC3B,MAAM,MAAM,GAAG,GAAG,CAAC,KAAK,EAAE,CAAC;QAC3B,MAAM,CAAC,aAAa,CAAC,IAAI,CAAC,oBAAqB,CAAC,CAAC;QACjD,IAAI,CAAC,MAAM,CAAC,iBAAiB,CAAC,MAAM,CAAC;YACnC,OAAO;QACT,uEAAuE;QACvE,IAAI,CAAC,mCAAmC,CAAC,GAAG,EAAE,CAAC,EAAE,CAAC,EAAE,IAAI,EAAE,GAAG,EAAE,CAAC,EAAE,CAAC,EAAE,IAAI,EAAE,QAAQ,CAAC,CAAC;QACrF,gEAAgE;QAChE,IAAI,CAAC,0BAA0B,CAAC,GAAG,EAAE,GAAG,EAAE,QAAQ,CAAC,CAAC;IACtD,CAAC;IACD,kEAAkE;IAC1D,yBAAyB,CAAC,GAAU,EAAE,GAAmB,EAAE,QAAiB;QAClF,MAAM,EAAE,GAAG,2BAAY,CAAC,MAAM,EAAE,CAAC;QACjC,GAAG,CAAC,WAAW,CAAC,EAAE,CAAC,CAAC;QACpB,IAAI,CAAC,oBAAoB,CAAC,GAAG,EAAE,EAAE,EAAE,QAAQ,CAAC,CAAC;IAC/C,CAAC;IACD,2FAA2F;IACnF,oCAAoC,CAC1C,OAA2B,EAAE,OAA2B,EAAE,QAAiB;QAE3E,MAAM,GAAG,GAAG,2BAAY,CAAC,MAAM,EAAE,CAAC;QAClC,OAAO,CAAC,WAAW,CAAC,GAAG,CAAC,CAAC;QACzB,MAAM,GAAG,GAAG,2BAAY,CAAC,MAAM,EAAE,CAAC;QAClC,OAAO,CAAC,WAAW,CAAC,GAAG,CAAC,CAAC;QACzB,IAAI,CAAC,2BAA2B,CAAC,GAAG,EAAE,GAAG,EAAE,QAAQ,CAAC,CAAC;IACvD,CAAC;IACD,yEAAyE;IAClE,8BAA8B,CAAC,GAAiB,EAAE,MAAsB,EAAE,QAAiB;QAChG,MAAM,GAAG,GAAG,2BAAY,CAAC,MAAM,EAAE,CAAC;QAClC,MAAM,CAAC,WAAW,CAAC,GAAG,CAAC,CAAC;QACxB,IAAI,CAAC,2BAA2B,CAAC,GAAG,EAAE,GAAG,EAAE,QAAQ,CAAC,CAAC;IACvD,CAAC;IACD,sEAAsE;IAC/D,2BAA2B,CAAC,IAAmB,EAAE,MAAsB,EAAE,QAAiB;QAC/F,MAAM,GAAG,GAAG,2BAAY,CAAC,MAAM,EAAE,CAAC;QAClC,MAAM,CAAC,WAAW,CAAC,GAAG,CAAC,CAAC;QACxB,IAAI,CAAC,wBAAwB,CAAC,IAAI,EAAE,GAAG,EAAE,QAAQ,CAAC,CAAC;IACrD,CAAC;IACD,6DAA6D;IACtD,wBAAwB,CAAC,IAAmB,EAAE,GAAiB,EAAE,QAAiB;QACvF,MAAM,IAAI,GAAG,GAAG,CAAC,SAAS,EAAE,CAAC;QAC7B,MAAM,UAAU,GAAG,mBAAQ,CAAC,kBAAkB,CAAC,CAAC,EAAE,IAAI,GAAG,CAAC,EAAE,CAAC,CAAC,CAAC;QAC/D,MAAM,OAAO,GAAG,IAAI,CAAC,SAAS,CAAC;QAC/B,MAAM,OAAO,GAAG,IAAI,CAAC,SAAS,CAAC;QAC/B,MAAM,OAAO,GAAG,yBAAyB,CAAC,aAAa,CAAC;QACxD,MAAM,OAAO,GAAG,yBAAyB,CAAC,aAAa,CAAC;QACxD,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,GAAG,CAAC,EAAE,EAAE,CAAC,EAAE;YACjC,MAAM,GAAG,GAAG,CAAC,GAAG,UAAU,CAAC,CAAC,8BAA8B;YAC1D,MAAM,GAAG,GAAG,CAAC,CAAC,GAAG,CAAC,KAAK,IAAI,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,UAAU,CAAC,CAAE,qCAAqC;YACrG,GAAG,CAAC,YAAY,CAAC,+BAA+B,CAAC,CAAC,EAAE,OAAO,CAAC,CAAC;YAC7D,GAAG,CAAC,YAAY,CAAC,+BAA+B,CAAC,CAAC,GAAG,CAAC,EAAE,OAAO,CAAC,CAAC;YACjE,IAAI,CAAC,sBAAsB,CAAC,IAAI,EAAE,OAAO,EAAE,GAAG,EAAE,OAAO,EAAE,GAAG,EAAE,GAAG,EAAE,OAAO,EAAE,GAAG,EAAE,OAAO,EAAE,GAAG,EAAE,QAAQ,CAAC,CAAC;SAC1G;IACH,CAAC;IACD,yDAAyD;IAClD,oBAAoB,CAAC,IAAW,EAAE,GAAiB,EAAE,QAAiB;QAC3E,MAAM,MAAM,GAAG,IAAI,CAAC,KAAK,EAAE,CAAC;QAC5B,MAAM,MAAM,GAAG,GAAG,CAAC,KAAK,EAAE,CAAC;QAC3B,MAAM,CAAC,aAAa,CAAC,IAAI,CAAC,oBAAqB,CAAC,CAAC;QACjD,IAAI,CAAC,MAAM,CAAC,iBAAiB,CAAC,MAAM,CAAC;YACnC,OAAO;QACT,MAAM,OAAO,GAAG,yBAAyB,CAAC,aAAa,CAAC;QACxD,MAAM,OAAO,GAAG,yBAAyB,CAAC,aAAa,CAAC;QACxD,MAAM,IAAI,GAAG,GAAG,CAAC,SAAS,EAAE,CAAC;QAC7B,IAAI,IAAI,GAAG,CAAC,EAAE;YACZ,MAAM,GAAG,GAAG,GAAG,GAAG,CAAC,IAAI,GAAG,CAAC,CAAC,CAAC;YAC7B,IAAI,GAAG,CAAC;YACR,IAAI,GAAG,CAAC;YACR,GAAG,GAAG,GAAG,CAAC;YACV,GAAG,CAAC,OAAO,CAAC,CAAC,EAAE,OAAO,CAAC,CAAC;YACxB,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,IAAI,EAAE,EAAE,EAAE,EAAE,OAAO,CAAC,OAAO,CAAC,OAAO,CAAC,EAAE,GAAG,GAAG,GAAG,EAAE;gBACrE,GAAG,CAAC,OAAO,CAAC,EAAE,EAAE,OAAO,CAAC,CAAC;gBACzB,GAAG,GAAG,EAAE,GAAG,GAAG,CAAC;gBACf,IAAI,CAAC,kBAAkB,CAAC,GAAG,EAAE,OAAO,EAAE,GAAG,EAAE,OAAO,EAAE,GAAG,EAAE,IAAI,EAAE,CAAC,QAAQ,CAAC,CAAC;aAC3E;SACF;QACD,OAAO,SAAS,CAAC;IACnB,CAAC;IACD,8CAA8C;IACtC,uBAAuB,CAAC,KAAe,EAAE,YAAiC;QAChF,MAAM,KAAK,GAAG,IAAI,CAAC,UAAU,CAAC,CAAE,OAAO;QACvC,IAAI,CAAC,KAAK,IAAI,CAAC,KAAK,CAAC,QAAQ,IAAI,CAAC,CAAC,KAAK,YAAY,iCAAe,CAAC;YAClE,OAAO;QACT,KAAK,MAAM,KAAK,IAAI,KAAK,CAAC,QAAQ,EAAE;YAClC,IAAI,CAAC,aAAa,CAAC,KAAK,CAAC,CAAC;YAC1B,YAAY,CAAC,KAAK,CAAC,CAAC;SACrB;QACD,IAAI,CAAC,UAAU,GAAG,KAAK,CAAC,CAAE,UAAU;IACtC,CAAC;IACD,oFAAoF;IAC5E,mCAAmC,CAAC,KAAe,EAAE,YAAiC;QAC5F,IAAI,CAAC,IAAI,CAAC,UAAU,IAAI,CAAC,CAAC,IAAI,CAAC,UAAU,YAAY,yDAA2B,CAAC;YAC/E,OAAO;QACT,IAAI,KAAK,YAAY,yDAA2B,EAAE;YAChD,IAAA,qBAAM,EAAC,CAAC,CAAC,uDAAuD,CAAC,CAAC;YAClE,OAAO;SACR;QACD,MAAM,MAAM,GAAG,IAAI,CAAC,QAAQ,CAAC,MAAM,CAAC;QACpC,MAAM,KAAK,GAAG,IAAI,CAAC,UAAU,CAAC,CAAE,OAAO;QACvC,KAAK,MAAM,KAAK,IAAI,KAAK,CAAC,IAAI,CAAC,QAAQ,EAAE;YACvC,IAAI,CAAC,aAAa,CAAC,KAAK,CAAC,CAAC;YAC1B,YAAY,CAAC,KAAK,CAAC,CAAC;SACrB;QACD,IAAI,CAAC,aAAa,CAAC,KAAK,CAAC,CAAC,CAAE,UAAU;QACtC,IAAI,CAAC,QAAQ,GAAG,yDAA2B,CAAC,+BAA+B,CAAC,IAAI,CAAC,QAAQ,EAAE,MAAM,EAAE,SAAS,EAAE,KAAK,EAAE,IAAI,CAAC,CAAC;IAC7H,CAAC;IACD,0DAA0D;IAC1C,mBAAmB,CAAC,QAAuB;QACzD,IAAI,IAAI,CAAC,UAAU,YAAY,6BAAa,EAAE;YAC5C,MAAM,QAAQ,GAAG,IAAI,CAAC,UAAU,CAAC;YACjC,IAAI,CAAC,sBAAsB,CACzB,QAAQ,EAAE,QAAQ,CAAC,SAAS,EAAE,GAAG,EAAE,QAAQ,CAAC,SAAS,EAAE,GAAG,EAC1D,QAAQ,EAAE,QAAQ,CAAC,SAAS,EAAE,GAAG,EAAE,QAAQ,CAAC,SAAS,EAAE,GAAG,EAC1D,KAAK,CACN,CAAC;SACH;aAAM,IAAI,IAAI,CAAC,UAAU,YAAY,2BAAY,EAAE;YAClD,IAAI,CAAC,wBAAwB,CAAC,QAAQ,EAAE,IAAI,CAAC,UAAU,EAAE,KAAK,CAAC,CAAC;SACjE;aAAM,IAAI,IAAI,CAAC,UAAU,YAAY,aAAK,EAAE;YAC3C,IAAI,CAAC,kBAAkB,CAAC,QAAQ,EAAE,QAAQ,CAAC,SAAS,EAAE,GAAG,EAAE,QAAQ,CAAC,SAAS,EAAE,GAAG,EAAE,IAAI,CAAC,UAAU,EAAE,KAAK,CAAC,CAAC;SAC7G;aAAM,IAAI,IAAI,CAAC,UAAU,YAAY,6BAAc,EAAE;YACpD,IAAI,CAAC,2BAA2B,CAAC,QAAQ,EAAE,IAAI,CAAC,UAAU,EAAE,KAAK,CAAC,CAAC;SACpE;aAAM,IAAI,IAAI,CAAC,UAAU,YAAY,iCAAe,EAAE;YACrD,IAAI,CAAC,uBAAuB,CAAC,QAAQ,EAAE,IAAI,CAAC,mBAAmB,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC;SAC7E;aAAM,IAAI,IAAI,CAAC,UAAU,YAAY,yDAA2B,EAAE;YACjE,IAAI,CAAC,mCAAmC,CAAC,QAAQ,EAAE,IAAI,CAAC,mBAAmB,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC;SACzF;QACD,OAAO,SAAS,CAAC;IACnB,CAAC;IACD;;;;;;;;;;;;;;;OAeG;IACK,wBAAwB,CAAC,CAAS,EAAE,CAAS,EAAE,KAAc;QACnE,IAAI,MAAM,GAAG,CAAC,CAAC;QACf,IAAI,CAAC,GAAG,KAAK,CAAC,GAAG,CAAC,CAAC;YACjB,MAAM,GAAG,IAAI,CAAC;aACX,IAAI,CAAC,GAAG,KAAK,CAAC,IAAI,CAAC,CAAC;YACvB,MAAM,GAAG,IAAI,CAAC;QAChB,wBAAwB;QACxB,IAAI,CAAC,GAAG,KAAK,CAAC,GAAG,CAAC,CAAC;YACjB,MAAM,IAAI,IAAI,CAAC;aACZ,IAAI,CAAC,GAAG,KAAK,CAAC,IAAI,CAAC,CAAC;YACvB,MAAM,IAAI,IAAI,CAAC;QACjB,OAAO,MAAM,CAAC;IAChB,CAAC;IACD,0DAA0D;IAClD,2BAA2B,CAAC,GAAiB,EAAE,GAAiB,EAAE,QAAiB;QACzF,MAAM,MAAM,GAAG,GAAG,CAAC,KAAK,EAAE,CAAC;QAC3B,MAAM,MAAM,GAAG,GAAG,CAAC,KAAK,EAAE,CAAC;QAC3B,MAAM,CAAC,aAAa,CAAC,IAAI,CAAC,oBAAqB,CAAC,CAAC;QACjD,IAAI,CAAC,MAAM,CAAC,iBAAiB,CAAC,MAAM,CAAC;YACnC,OAAO;QACT,IAAI,KAAa,CAAC;QAClB,IAAI,KAAa,CAAC;QAClB,MAAM,OAAO,GAAG,eAAO,CAAC,UAAU,EAAE,CAAC;QACrC,MAAM,OAAO,GAAG,yBAAyB,CAAC,aAAa,CAAC;QACxD,MAAM,OAAO,GAAG,yBAAyB,CAAC,aAAa,CAAC;QACxD,MAAM,OAAO,GAAG,yBAAyB,CAAC,aAAa,CAAC;QACxD,MAAM,OAAO,GAAG,yBAAyB,CAAC,aAAa,CAAC;QACxD,MAAM,IAAI,GAAG,GAAG,CAAC,SAAS,EAAE,CAAC;QAC7B,MAAM,IAAI,GAAG,GAAG,CAAC,SAAS,EAAE,CAAC;QAC7B,IAAI,IAAI,GAAG,CAAC,IAAI,IAAI,GAAG,CAAC,EAAE;YACxB,MAAM,GAAG,GAAG,GAAG,GAAG,CAAC,IAAI,GAAG,CAAC,CAAC,CAAC;YAC7B,MAAM,GAAG,GAAG,GAAG,GAAG,CAAC,IAAI,GAAG,CAAC,CAAC,CAAC;YAC7B,IAAI,GAAG,GAAG,GAAG,CAAC;YACd,IAAI,GAAG,EAAE,GAAG,EAAE,GAAG,CAAC;YAClB,GAAG,CAAC,OAAO,CAAC,CAAC,EAAE,OAAO,CAAC,CAAC;YACxB,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,IAAI,EAAE,EAAE,EAAE,EAAE,OAAO,CAAC,OAAO,CAAC,OAAO,CAAC,EAAE,GAAG,GAAG,GAAG,EAAE;gBACrE,GAAG,GAAG,EAAE,GAAG,GAAG,CAAC;gBACf,GAAG,GAAG,GAAG,CAAC;gBACV,GAAG,CAAC,OAAO,CAAC,EAAE,EAAE,OAAO,CAAC,CAAC;gBACzB,kEAAkE;gBAClE,OAAO,CAAC,OAAO,EAAE,CAAC;gBAClB,OAAO,CAAC,WAAW,CAAC,OAAO,CAAC,CAAC;gBAC7B,OAAO,CAAC,WAAW,CAAC,OAAO,CAAC,CAAC;gBAC7B,OAAO,CAAC,aAAa,CAAC,IAAI,CAAC,oBAAqB,CAAC,CAAC;gBAClD,IAAI,OAAO,CAAC,iBAAiB,CAAC,MAAM,CAAC,EAAE;oBACrC,GAAG,CAAC,OAAO,CAAC,CAAC,EAAE,OAAO,CAAC,CAAC;oBACxB,KAAK,GAAG,IAAI,CAAC,wBAAwB,CAAC,OAAO,CAAC,CAAC,EAAE,OAAO,CAAC,CAAC,EAAE,OAAO,CAAC,CAAC;oBACrE,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,IAAI,EAAE,EAAE,EAAE,EAAE,OAAO,CAAC,OAAO,CAAC,OAAO,CAAC,EAAE,GAAG,GAAG,GAAG,EAAE,KAAK,GAAG,KAAK,EAAE;wBACpF,GAAG,CAAC,OAAO,CAAC,EAAE,EAAE,OAAO,CAAC,CAAC;wBACzB,KAAK,GAAG,IAAI,CAAC,wBAAwB,CAAC,OAAO,CAAC,CAAC,EAAE,OAAO,CAAC,CAAC,EAAE,OAAO,CAAC,CAAC;wBACrE,GAAG,GAAG,EAAE,GAAG,GAAG,CAAC;wBACf,sFAAsF;wBACtF,oFAAoF;wBACpF,IAAI,CAAC,KAAK,GAAG,KAAK,CAAC,KAAK,CAAC;4BACvB,IAAI,CAAC,sBAAsB,CAAC,GAAG,EAAE,OAAO,EAAE,GAAG,EAAE,OAAO,EAAE,GAAG,EAAE,GAAG,EAAE,OAAO,EAAE,GAAG,EAAE,OAAO,EAAE,GAAG,EAAE,QAAQ,CAAC,CAAC;qBAC3G;iBACF;aACF;SACF;IACH,CAAC;IACD,6DAA6D;IAC7C,kBAAkB,CAAC,GAAiB;QAClD,IAAI,IAAI,CAAC,UAAU,YAAY,2BAAY,EAAE;YAC3C,MAAM,GAAG,GAAG,IAAI,CAAC,UAAU,CAAC;YAC5B,IAAI,CAAC,2BAA2B,CAAC,GAAG,EAAE,GAAG,EAAE,KAAK,CAAC,CAAC;SACnD;aAAM,IAAI,IAAI,CAAC,UAAU,YAAY,6BAAa,EAAE;YACnD,IAAI,CAAC,wBAAwB,CAAC,IAAI,CAAC,UAAU,EAAE,GAAG,EAAE,IAAI,CAAC,CAAC;SAC3D;aAAM,IAAI,IAAI,CAAC,UAAU,YAAY,aAAK,EAAE;YAC3C,IAAI,CAAC,oBAAoB,CAAC,IAAI,CAAC,UAAU,EAAE,GAAG,EAAE,IAAI,CAAC,CAAC;SACvD;aAAM,IAAI,IAAI,CAAC,UAAU,YAAY,6BAAc,EAAE;YACpD,IAAI,CAAC,8BAA8B,CAAC,GAAG,EAAE,IAAI,CAAC,UAAU,EAAE,KAAK,CAAC,CAAC;SAClE;aAAM,IAAI,IAAI,CAAC,UAAU,YAAY,iCAAe,EAAE;YACrD,IAAI,CAAC,uBAAuB,CAAC,GAAG,EAAE,IAAI,CAAC,kBAAkB,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC;SACvE;aAAM,IAAI,IAAI,CAAC,UAAU,YAAY,yDAA2B,EAAE;YACjE,IAAI,CAAC,mCAAmC,CAAC,GAAG,EAAE,IAAI,CAAC,kBAAkB,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC;SACnF;QACD,OAAO,SAAS,CAAC;IACnB,CAAC;IACD,sDAAsD;IACtC,WAAW,CAAC,IAAW;QACrC,IAAI,IAAI,CAAC,UAAU,YAAY,6BAAa,EAAE;YAC5C,IAAI,CAAC,kBAAkB,CACrB,IAAI,CAAC,UAAU,EAAE,IAAI,CAAC,UAAU,CAAC,SAAS,EAAE,GAAG,EAAE,IAAI,CAAC,UAAU,CAAC,SAAS,EAAE,GAAG,EAAE,IAAI,EAAE,IAAI,CAC5F,CAAC;SACH;aAAM,IAAI,IAAI,CAAC,UAAU,YAAY,2BAAY,EAAE;YAClD,IAAI,CAAC,oBAAoB,CAAC,IAAI,EAAE,IAAI,CAAC,UAAU,EAAE,KAAK,CAAC,CAAC;SACzD;aAAM,IAAI,IAAI,CAAC,UAAU,YAAY,aAAK,EAAE;YAC3C,IAAI,CAAC,cAAc,CAAC,IAAI,EAAE,IAAI,CAAC,UAAU,EAAE,KAAK,CAAC,CAAC;SACnD;aAAM,IAAI,IAAI,CAAC,UAAU,YAAY,6BAAc,EAAE;YACpD,IAAI,CAAC,yBAAyB,CAAC,IAAI,EAAE,IAAI,CAAC,UAAU,EAAE,KAAK,CAAC,CAAC;SAC9D;aAAM,IAAI,IAAI,CAAC,UAAU,YAAY,iCAAe,EAAE;YACrD,IAAI,CAAC,uBAAuB,CAAC,IAAI,EAAE,IAAI,CAAC,WAAW,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC;SACjE;aAAM,IAAI,IAAI,CAAC,UAAU,YAAY,yDAA2B,EAAE;YACjE,IAAI,CAAC,mCAAmC,CAAC,IAAI,EAAE,IAAI,CAAC,WAAW,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC;SAC7E;QACD,OAAO,SAAS,CAAC;IACnB,CAAC;IACD,gEAAgE;IAChD,oBAAoB,CAAC,KAAqB;QACxD,IAAI,IAAI,CAAC,UAAU,YAAY,6BAAa,EAAE;YAC5C,IAAI,CAAC,2BAA2B,CAAC,IAAI,CAAC,UAAU,EAAE,KAAK,EAAE,IAAI,CAAC,CAAC;SAChE;aAAM,IAAI,IAAI,CAAC,UAAU,YAAY,2BAAY,EAAE;YAClD,IAAI,CAAC,8BAA8B,CAAC,IAAI,CAAC,UAAU,EAAE,KAAK,EAAE,IAAI,CAAC,CAAC;SACnE;aAAM,IAAI,IAAI,CAAC,UAAU,YAAY,aAAK,EAAE;YAC3C,IAAI,CAAC,yBAAyB,CAAC,IAAI,CAAC,UAAU,EAAE,KAAK,EAAE,IAAI,CAAC,CAAC;SAC9D;aAAM,IAAI,IAAI,CAAC,UAAU,YAAY,iCAAkB,EAAE;YACxD,IAAI,CAAC,oCAAoC,CAAC,KAAK,EAAE,IAAI,CAAC,UAAU,EAAE,KAAK,CAAC,CAAC;SAC1E;aAAM,IAAI,IAAI,CAAC,UAAU,YAAY,iCAAe,EAAE;YACrD,IAAI,CAAC,uBAAuB,CAAC,KAAK,EAAE,IAAI,CAAC,oBAAoB,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC;SAC3E;aAAM,IAAI,IAAI,CAAC,UAAU,YAAY,yDAA2B,EAAE;YACjE,IAAI,CAAC,mCAAmC,CAAC,KAAK,EAAE,IAAI,CAAC,oBAAoB,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC;SACvF;QACD,OAAO,SAAS,CAAC;IACnB,CAAC;IACD,8EAA8E;IAC9D,iCAAiC,CAAC,KAAkC;QAClF,KAAK,CAAC,iCAAiC,CAAC,KAAK,CAAC,CAAC;QAC/C,kIAAkI;QAClI,IAAI,CAAC,QAAQ,GAAG,yDAA2B,CAAC,+BAA+B,CAAC,IAAI,CAAC,QAAQ,EAAE,CAAC,EAAE,KAAK,EAAE,SAAS,EAAE,IAAI,CAAC,CAAC;IACxH,CAAC;IACD,8EAA8E;IAC9D,qBAAqB,CAAC,MAAuB;QAC3D;;;;;;;;;;;;UAYE;QACF,OAAO,SAAS,CAAC;IACnB,CAAC;;AAt1Bc,uCAAa,GAAG,yBAAO,CAAC,MAAM,EAAE,CAAC;AACjC,uCAAa,GAAG,yBAAO,CAAC,MAAM,EAAE,CAAC;AACjC,uCAAa,GAAG,yBAAO,CAAC,MAAM,EAAE,CAAC;AACjC,uCAAa,GAAG,yBAAO,CAAC,MAAM,EAAE,CAAC;AACjC,qCAAW,GAAG,yBAAO,CAAC,MAAM,EAAE,CAAC;AAjCnC,8DAAyB","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 Curve\r\n */\r\n\r\nimport { assert } from \"@itwin/core-bentley\";\r\nimport { BSplineCurve3d, BSplineCurve3dBase } from \"../../bspline/BSplineCurve\";\r\nimport { BSplineCurve3dH } from \"../../bspline/BSplineCurve3dH\";\r\nimport { Geometry } from \"../../Geometry\";\r\nimport { RecurseToCurvesGeometryHandler } from \"../../geometry3d/GeometryHandler\";\r\nimport { GrowableFloat64Array } from \"../../geometry3d/GrowableFloat64Array\";\r\nimport { Point3d } from \"../../geometry3d/Point3dVector3d\";\r\nimport { Range3d } from \"../../geometry3d/Range\";\r\nimport {\r\n  CurveCurveCloseApproachXYRRtoRRD, CurvePointCloseApproachXYRtoRD,\r\n  Newton1dUnbounded, Newton2dUnboundedWithDerivative, NewtonEvaluatorRRtoRRD,\r\n} from \"../../numerics/Newton\";\r\nimport { AnalyticRoots, SmallSystem } from \"../../numerics/Polynomials\";\r\nimport { Arc3d } from \"../Arc3d\";\r\nimport { CurveChainWithDistanceIndex } from \"../CurveChainWithDistanceIndex\";\r\nimport { CurveCollection } from \"../CurveCollection\";\r\nimport { CurveIntervalRole, CurveLocationDetail, CurveLocationDetailPair } from \"../CurveLocationDetail\";\r\nimport { CurvePrimitive } from \"../CurvePrimitive\";\r\nimport { AnyCurve } from \"../CurveTypes\";\r\nimport { LineSegment3d } from \"../LineSegment3d\";\r\nimport { LineString3d } from \"../LineString3d\";\r\n\r\n// cspell:word XYRR currentdFdX\r\n\r\n/**\r\n * Handler class for XY close approach between _geometryB and another geometry.\r\n * * Approach means the XY distance (z is ignored) between _geometryB and another geometry.\r\n * * Closest approach is a measure of the proximity of one curve to another. It's the length of the shortest line\r\n * segment perpendicular to both curves; if the curves intersect, the closest approach is zero. In the context of\r\n * this class, z-coordinates are ignored, so the closest approach is as seen in the top view. If you have coplanar\r\n * input curves and want to find closest approach in their plane, rotate them first into a plane parallel to the\r\n * xy-plane, then afterward, rotate the results back as required.\r\n * * Close approach can also be from a curve endpoint perpendicular to another curve or from a curve endpoint to\r\n * another curve endpoint.\r\n * * Instances are initialized and called from CurveCurve.\r\n * * geometryB is saved for later reference.\r\n * @internal\r\n */\r\nexport class CurveCurveCloseApproachXY extends RecurseToCurvesGeometryHandler {\r\n  private _geometryB: AnyCurve | undefined;\r\n  private _circularArcB: Arc3d | undefined;\r\n  private _circularRadiusB: number | undefined;\r\n  private setGeometryB(geometryB: AnyCurve | undefined) {\r\n    this._geometryB = geometryB;\r\n    this._circularArcB = undefined;\r\n    this._circularRadiusB = undefined;\r\n    if (geometryB instanceof Arc3d) {\r\n      const r = geometryB.circularRadiusXY();\r\n      if (r !== undefined) {\r\n        this._circularRadiusB = r;\r\n        this._circularArcB = geometryB;\r\n      }\r\n    }\r\n  }\r\n  /**\r\n   * Maximum XY distance (z is ignored). Approach larger than this is not interesting.\r\n   * This is caller defined and can be undefined.\r\n   */\r\n  private _maxDistanceToAccept: number | undefined;\r\n  /** Squared max distance. This is private, and is forced to at least small metric distance squared. */\r\n  private _maxDistanceSquared: number;\r\n  /**\r\n   * Start and end points of line segments that meet closest approach criteria, i.e., they are perpendicular to\r\n   * both curves and their length is smaller than _maxDistanceToAccept.\r\n   */\r\n  private _results: CurveLocationDetailPair[];\r\n\r\n  private static _workPointAA0 = Point3d.create();\r\n  private static _workPointAA1 = Point3d.create();\r\n  private static _workPointBB0 = Point3d.create();\r\n  private static _workPointBB1 = Point3d.create();\r\n  private static _workPointB = Point3d.create();\r\n\r\n  /**\r\n   * Constructor.\r\n   * @param geometryB second curve for intersection. Saved for reference by specific handler methods.\r\n   */\r\n  public constructor(geometryB: AnyCurve | undefined) {\r\n    super();\r\n    this.setGeometryB(geometryB);\r\n    this._maxDistanceSquared = Geometry.smallMetricDistanceSquared;\r\n    this._results = [];\r\n  }\r\n  /** Set the (possibly undefined) max XY distance (z is ignored) to accept. */\r\n  public set maxDistanceToAccept(value: number | undefined) {\r\n    this._maxDistanceToAccept = value;\r\n    if (this._maxDistanceToAccept !== undefined && this._maxDistanceToAccept > 0)\r\n      this._maxDistanceSquared = this._maxDistanceToAccept * this._maxDistanceToAccept;\r\n  }\r\n  /** Access the (possibly undefined) max XY distance (z is ignored) to accept. */\r\n  public get maxDistanceToAccept(): number | undefined {\r\n    return this._maxDistanceToAccept;\r\n  }\r\n  /** Ask if the maxDistanceToAccept value is defined and positive */\r\n  public get isMaxDistanceSet(): boolean {\r\n    return this._maxDistanceToAccept !== undefined && this._maxDistanceToAccept > 0;\r\n  }\r\n  /** Reset the geometry and flags, leaving all other parts unchanged (and preserving accumulated intersections) */\r\n  public resetGeometry(geometryB: AnyCurve) {\r\n    this.setGeometryB(geometryB);\r\n  }\r\n  /** returns true if `fraction` is in [0,1] within tolerance */\r\n  private acceptFraction(fraction: number, fractionTol: number = 1.0e-12) {\r\n    if (fraction < -fractionTol)\r\n      return false;\r\n    if (fraction > 1.0 + fractionTol)\r\n      return false;\r\n    return true;\r\n  }\r\n  /**\r\n   * Return the results structure for the intersection calculation, structured as an array of CurveLocationDetailPair.\r\n   * @param reinitialize if true, a new results structure is created for use by later calls.\r\n   */\r\n  public grabPairedResults(reinitialize: boolean = false): CurveLocationDetailPair[] {\r\n    const result = this._results;\r\n    if (reinitialize)\r\n      this._results = [];\r\n    return result;\r\n  }\r\n  /**\r\n   * If distance between pointA and pointB is less than maxDistance, record CurveLocationDetailPair which is\r\n   * the approach from pointA to pointB.\r\n   */\r\n  private testAndRecordPointPairApproach(\r\n    cpA: CurvePrimitive, fA: number, pointA: Point3d, cpB: CurvePrimitive, fB: number, pointB: Point3d, reversed: boolean,\r\n  ): void {\r\n    const d2 = pointA.distanceSquaredXY(pointB);\r\n    if (d2 < this._maxDistanceSquared) {\r\n      const detailA = CurveLocationDetail.createCurveFractionPoint(cpA, fA, pointA);\r\n      const detailB = CurveLocationDetail.createCurveFractionPoint(cpB, fB, pointB);\r\n      const pair = CurveLocationDetailPair.createCapture(detailA, detailB);\r\n      if (reversed)\r\n        pair.swapDetails();\r\n      this._results.push(pair);\r\n    }\r\n  }\r\n  /**\r\n   * Create a close approach pair if XY distance is within maxDistance.\r\n   * @param localFractionA a fraction on first curve\r\n   * @param cpA the first curve\r\n   * @param fractionA0 start of the first curve\r\n   * @param fractionA1 end of the first curve\r\n   * @param localFractionB a fraction on second curve\r\n   * @param cpB the second curve\r\n   * @param fractionB0 start of the second curve\r\n   * @param fractionB1 end of the second curve\r\n   * @param reversed whether to reverse the details in the pair (e.g., so that detailB refers to geometryB).\r\n   * @param intervalDetails optional CurveLocationDetailPair\r\n   */\r\n  private recordPointWithLocalFractions(\r\n    localFractionA: number,\r\n    cpA: CurvePrimitive,\r\n    fractionA0: number,\r\n    fractionA1: number,\r\n    localFractionB: number,\r\n    cpB: CurvePrimitive,\r\n    fractionB0: number,\r\n    fractionB1: number,\r\n    reversed: boolean,\r\n    intervalDetails?: undefined | CurveLocationDetailPair,\r\n  ): void {\r\n    let globalFractionA, globalFractionB;\r\n    let globalFractionA1, globalFractionB1;\r\n    const isInterval = intervalDetails !== undefined &&\r\n      intervalDetails.detailA.hasFraction1 &&\r\n      intervalDetails.detailB.hasFraction1;\r\n    if (isInterval) {\r\n      globalFractionA = Geometry.interpolate(fractionA0, intervalDetails.detailA.fraction, fractionA1);\r\n      globalFractionB = Geometry.interpolate(fractionB0, intervalDetails.detailB.fraction, fractionB1);\r\n      globalFractionA1 = Geometry.interpolate(fractionA0, intervalDetails.detailA.fraction1!, fractionA1);\r\n      globalFractionB1 = Geometry.interpolate(fractionB0, intervalDetails.detailB.fraction1!, fractionB1);\r\n    } else {\r\n      globalFractionA = globalFractionA1 = Geometry.interpolate(fractionA0, localFractionA, fractionA1);\r\n      globalFractionB = globalFractionB1 = Geometry.interpolate(fractionB0, localFractionB, fractionB1);\r\n    }\r\n    // ignore duplicate of most recent approach\r\n    const numPrevious = this._results.length;\r\n    if (numPrevious > 0 && !isInterval) {\r\n      const oldDetailA = this._results[numPrevious - 1].detailA;\r\n      const oldDetailB = this._results[numPrevious - 1].detailB;\r\n      if (reversed) {\r\n        if (oldDetailB.isSameCurveAndFraction({ curve: cpA, fraction: globalFractionA }) &&\r\n          oldDetailA.isSameCurveAndFraction({ curve: cpB, fraction: globalFractionB }))\r\n          return;\r\n      } else {\r\n        if (oldDetailA.isSameCurveAndFraction({ curve: cpA, fraction: globalFractionA }) &&\r\n          oldDetailB.isSameCurveAndFraction({ curve: cpB, fraction: globalFractionB }))\r\n          return;\r\n      }\r\n    }\r\n    const detailA = CurveLocationDetail.createCurveFractionPoint(\r\n      cpA, globalFractionA, cpA.fractionToPoint(globalFractionA),\r\n    );\r\n    const detailB = CurveLocationDetail.createCurveFractionPoint(\r\n      cpB, globalFractionB, cpB.fractionToPoint(globalFractionB),\r\n    );\r\n    if (isInterval) {\r\n      detailA.captureFraction1Point1(globalFractionA1, cpA.fractionToPoint(globalFractionA1));\r\n      detailB.captureFraction1Point1(globalFractionB1, cpB.fractionToPoint(globalFractionB1));\r\n    } else {\r\n      const d2 = detailA.point.distanceSquaredXY(detailB.point);\r\n      if (d2 > this._maxDistanceSquared)\r\n        return;\r\n      detailA.setIntervalRole(CurveIntervalRole.isolated);\r\n      detailB.setIntervalRole(CurveIntervalRole.isolated);\r\n    }\r\n    if (reversed) {\r\n      this._results.push(new CurveLocationDetailPair(detailB, detailA));\r\n    } else {\r\n      this._results.push(new CurveLocationDetailPair(detailA, detailB));\r\n    }\r\n  }\r\n  /**\r\n   * Capture a close approach pair that has point and local fraction but not curve.\r\n   * * Record the pair, each detail modified with global fraction and input curve.\r\n   * * Pair is neither modified nor recorded if it would be a duplicate of the last recorded pair.\r\n   * @param pair details computed with local fractions\r\n   * @param cpA curveA\r\n   * @param fractionA0 global start fraction on curveA\r\n   * @param fractionA1 global end fraction on curveA\r\n   * @param cpB curveB\r\n   * @param fractionB0 global start fraction on curveB\r\n   * @param fractionB1 global end fraction on curveB\r\n   * @param reversed whether to reverse the details in the pair (e.g., so that detailB refers to geometryB).\r\n   */\r\n  private capturePairWithLocalFractions(\r\n    pair: CurveLocationDetailPair,\r\n    cpA: CurvePrimitive,\r\n    fractionA0: number,\r\n    fractionA1: number,\r\n    cpB: CurvePrimitive,\r\n    fractionB0: number,\r\n    fractionB1: number,\r\n    reversed: boolean,\r\n  ) {\r\n    const globalFractionA = Geometry.interpolate(fractionA0, pair.detailA.fraction, fractionA1);\r\n    const globalFractionB = Geometry.interpolate(fractionB0, pair.detailB.fraction, fractionB1);\r\n    // ignore duplicate of most recent pair\r\n    const numPrevious = this._results.length;\r\n    if (numPrevious > 0) {\r\n      const oldDetailA = this._results[numPrevious - 1].detailA;\r\n      const oldDetailB = this._results[numPrevious - 1].detailB;\r\n      if (reversed) {\r\n        if (oldDetailB.isSameCurveAndFraction({ curve: cpA, fraction: globalFractionA }) &&\r\n          oldDetailA.isSameCurveAndFraction({ curve: cpB, fraction: globalFractionB }))\r\n          return;\r\n      } else {\r\n        if (oldDetailA.isSameCurveAndFraction({ curve: cpA, fraction: globalFractionA }) &&\r\n          oldDetailB.isSameCurveAndFraction({ curve: cpB, fraction: globalFractionB }))\r\n          return;\r\n      }\r\n    }\r\n    // recompute the points just in case\r\n    CurveLocationDetail.createCurveEvaluatedFraction(cpA, globalFractionA, pair.detailA);\r\n    CurveLocationDetail.createCurveEvaluatedFraction(cpB, globalFractionB, pair.detailB);\r\n    pair.detailA.setIntervalRole(CurveIntervalRole.isolated);\r\n    pair.detailB.setIntervalRole(CurveIntervalRole.isolated);\r\n    if (reversed)\r\n      pair.swapDetails();\r\n    this._results.push(pair);\r\n  }\r\n  /**\r\n   * Emit recordPoint for multiple pairs (on full curve) if within maxDistance.\r\n   * @param cpA first curve primitive (possibly different from curve in detailA, but fraction compatible)\r\n   * @param cpB second curve primitive (possibly different from curve in detailA, but fraction compatible)\r\n   * @param pairs array of pairs\r\n   * @param reversed whether to reverse the details in the pair (e.g., so that detailB refers to geometryB).\r\n   */\r\n  public recordPairs(\r\n    cpA: CurvePrimitive, cpB: CurvePrimitive, pairs: CurveLocationDetailPair[] | undefined, reversed: boolean,\r\n  ): void {\r\n    if (pairs !== undefined) {\r\n      for (const p of pairs) {\r\n        this.recordPointWithLocalFractions(\r\n          p.detailA.fraction, cpA, 0, 1, p.detailB.fraction, cpB, 0, 1, reversed, p,\r\n        );\r\n      }\r\n    }\r\n  }\r\n  /**\r\n   * Record fully assembled (but possibly reversed) detail pair.\r\n   * @param detailA first detail\r\n   * @param detailB second detail\r\n   * @param reversed whether to reverse the details in the pair (e.g., so that detailB refers to geometryB).\r\n   */\r\n  public captureDetailPair(\r\n    detailA: CurveLocationDetail | undefined, detailB: CurveLocationDetail | undefined, reversed: boolean,\r\n  ): void {\r\n    if (detailA && detailB) {\r\n      if (reversed) {\r\n        this._results.push(CurveLocationDetailPair.createCapture(detailB, detailA));\r\n      } else {\r\n        this._results.push(CurveLocationDetailPair.createCapture(detailA, detailB));\r\n      }\r\n    }\r\n  }\r\n  private static updatePointToSegmentDistance(\r\n    fractionA: number,\r\n    pointA: Point3d,\r\n    pointB0: Point3d,\r\n    pointB1: Point3d,\r\n    fractionB: number,\r\n    maxDistanceSquared: number,\r\n    closestApproach: CurveLocationDetailPair,   // modified on return\r\n  ): boolean {\r\n    let updated = false;\r\n    if (fractionB < 0)\r\n      fractionB = 0;\r\n    else if (fractionB > 1)\r\n      fractionB = 1;\r\n    this._workPointB = pointB0.interpolate(fractionB, pointB1, this._workPointB);\r\n    const distanceSquared = this._workPointB.distanceSquaredXY(pointA);\r\n    if (distanceSquared <= Math.min(maxDistanceSquared, closestApproach.detailA.a)) {\r\n      closestApproach.detailA.setFP(fractionA, pointA, undefined, distanceSquared);\r\n      closestApproach.detailB.setFP(fractionB, this._workPointB, undefined, distanceSquared);\r\n      updated = true;\r\n    }\r\n    return updated;\r\n  }\r\n  /**\r\n   * Return fractions of close approach within maxDistance between two line segments (a0,a1) and (b0,b1).\r\n   * * Math details can be found at core/geometry/internaldocs/Curve.md\r\n   * @param a0 start point of line a\r\n   * @param a1 end point of line a\r\n   * @param b0 start point of line b\r\n   * @param b1 end point of line b\r\n   * @param maxDistanceSquared maximum distance squared (assumed to be positive)\r\n   * @returns the fractional (not xy) coordinates in result.x and result.y. result.x is fraction on line a.\r\n   * result.y is fraction on line b.\r\n   */\r\n  private static segmentSegmentBoundedApproach(\r\n    a0: Point3d,\r\n    a1: Point3d,\r\n    b0: Point3d,\r\n    b1: Point3d,\r\n    maxDistanceSquared: number,\r\n  ): CurveLocationDetailPair | undefined {\r\n    const ux = a1.x - a0.x;\r\n    const uy = a1.y - a0.y;\r\n    const vx = b1.x - b0.x;\r\n    const vy = b1.y - b0.y;\r\n    const e00x = b0.x - a0.x;\r\n    const e00y = b0.y - a0.y;\r\n    const e01x = b1.x - a0.x;\r\n    const e01y = b1.y - a0.y;\r\n    const e10x = b0.x - a1.x;\r\n    const e10y = b0.y - a1.y;\r\n    const hab0 = Geometry.crossProductXYXY(ux, uy, e00x, e00y);\r\n    const hab1 = Geometry.crossProductXYXY(ux, uy, e01x, e01y);\r\n    const hba0 = -Geometry.crossProductXYXY(vx, vy, e00x, e00y);\r\n    const hba1 = -Geometry.crossProductXYXY(vx, vy, e10x, e10y);\r\n    if (hab0 * hab1 < 0.0 && hba0 * hba1 < 0.0) { // true intersection, strictly within both segments\r\n      const fractionA = -hba0 / (hba1 - hba0);\r\n      const fractionB = -hab0 / (hab1 - hab0);\r\n      return CurveLocationDetailPair.createCapture(\r\n        CurveLocationDetail.createCurveFractionPoint(undefined, fractionA, a0.interpolate(fractionA, a1)),\r\n        CurveLocationDetail.createCurveFractionPoint(undefined, fractionB, b0.interpolate(fractionB, b1)),\r\n      );\r\n    }\r\n    // there's no intersection, so find the closest approach within maxDistance from an endpoint\r\n    const closestApproach = new CurveLocationDetailPair();\r\n    closestApproach.detailA.a = 2 * maxDistanceSquared; // init to an approach that's too far away\r\n    let reversed = false;\r\n    const uu = Geometry.hypotenuseSquaredXY(ux, uy);\r\n    if (hab0 * hab0 <= maxDistanceSquared * uu) { // test distance of b0 to u\r\n      const fractionA = Geometry.dotProductXYXY(ux, uy, e00x, e00y) / uu;\r\n      if (this.updatePointToSegmentDistance(0, b0, a0, a1, fractionA, maxDistanceSquared, closestApproach))\r\n        reversed = true;\r\n    }\r\n    if (hab1 * hab1 <= maxDistanceSquared * uu) { // test distance of b1 to u\r\n      const fractionA = Geometry.dotProductXYXY(ux, uy, e01x, e01y) / uu;\r\n      if (this.updatePointToSegmentDistance(1, b1, a0, a1, fractionA, maxDistanceSquared, closestApproach))\r\n        reversed = true;\r\n    }\r\n    const vv = Geometry.hypotenuseSquaredXY(vx, vy);\r\n    if (hba0 * hba0 <= maxDistanceSquared * vv) { // test distance of a0 to v\r\n      const fractionB = -Geometry.dotProductXYXY(vx, vy, e00x, e00y) / vv;\r\n      if (this.updatePointToSegmentDistance(0, a0, b0, b1, fractionB, maxDistanceSquared, closestApproach))\r\n        reversed = false;\r\n    }\r\n    if (hba1 * hba1 <= maxDistanceSquared * vv) { // test distance of a1 to v\r\n      const fractionB = -Geometry.dotProductXYXY(vx, vy, e10x, e10y) / vv;\r\n      if (this.updatePointToSegmentDistance(1, a1, b0, b1, fractionB, maxDistanceSquared, closestApproach))\r\n        reversed = false;\r\n    }\r\n    if (closestApproach.detailA.a > maxDistanceSquared)\r\n      return undefined;\r\n    if (reversed)\r\n      closestApproach.swapDetails();\r\n    return closestApproach;\r\n  }\r\n  /**\r\n   * Check different combination of fractions on curveA and curveB. If distance between points at 2 fractions\r\n   * is less than maxDistance, record CurveLocationDetailPair which is the approach between the 2 points.\r\n   * Optionally, record close approaches of one curve's points if they fall between the other curve's points.\r\n   * @param cpA curveA\r\n   * @param fA0 fraction0 on curveA\r\n   * @param fA1 fraction1 on curveA\r\n   * @param testProjectionOnA whether to record projections of the given curveB points onto curveA\r\n   * @param cpB curveB\r\n   * @param fB0 fraction0 on curveB\r\n   * @param fB1 fraction0 on curveB\r\n   * @param testProjectionOnB whether to record projections of the given curveA points onto curveB\r\n   * @param reversed whether to reverse the details in the pair (e.g., so that detailB refers to geometryB).\r\n   */\r\n  private testAndRecordFractionalPairApproach(\r\n    cpA: CurvePrimitive,\r\n    fA0: number,\r\n    fA1: number,\r\n    testProjectionOnA: boolean,\r\n    cpB: CurvePrimitive,\r\n    fB0: number,\r\n    fB1: number,\r\n    testProjectionOnB: boolean,\r\n    reversed: boolean,\r\n  ): void {\r\n    const pointA0 = cpA.fractionToPoint(fA0);\r\n    const pointA1 = cpA.fractionToPoint(fA1);\r\n    const pointB0 = cpB.fractionToPoint(fB0);\r\n    const pointB1 = cpB.fractionToPoint(fB1);\r\n    this.testAndRecordPointPairApproach(cpA, fA0, pointA0, cpB, fB0, pointB0, reversed);\r\n    this.testAndRecordPointPairApproach(cpA, fA1, pointA1, cpB, fB0, pointB0, reversed);\r\n    this.testAndRecordPointPairApproach(cpA, fA0, pointA0, cpB, fB1, pointB1, reversed);\r\n    this.testAndRecordPointPairApproach(cpA, fA1, pointA1, cpB, fB1, pointB1, reversed);\r\n    if (testProjectionOnB) {\r\n      this.testAndRecordProjection(cpA, fA0, pointA0, cpB, fB0, fB1, reversed);\r\n      this.testAndRecordProjection(cpA, fA1, pointA1, cpB, fB0, fB1, reversed);\r\n    }\r\n    if (testProjectionOnA) {\r\n      this.testAndRecordProjection(cpB, fB0, pointB0, cpA, fA0, fA1, !reversed);\r\n      this.testAndRecordProjection(cpB, fB1, pointB1, cpA, fA0, fA1, !reversed);\r\n    }\r\n  }\r\n  /**\r\n   * Return XY closest approach between a curve primitive and a point.\r\n   * Currently, this function only supports Arc3d and LineSegment.\r\n   * Note that this function doesn't handle endpoints.\r\n   */\r\n  private getPointCurveClosestApproachXYNewton(curveP: CurvePrimitive, pointQ: Point3d): CurveLocationDetail | undefined {\r\n    if (!(curveP instanceof Arc3d) && !(curveP instanceof LineSegment3d)) {\r\n      assert(!!\"getPointCurveClosestApproachXYNewton only supports Arc3d and LineSegment\");\r\n      return undefined;\r\n    }\r\n    const seeds = [0.2, 0.4, 0.6, 0.8]; // HEURISTIC: arcs have up to 4 perpendiculars; lines have only 1\r\n    const newtonEvaluator = new CurvePointCloseApproachXYRtoRD(curveP, pointQ);\r\n    const newtonSearcher = new Newton1dUnbounded(newtonEvaluator);\r\n    let minCloseApproachLength = Geometry.largeCoordinateResult;\r\n    let minCurvePFraction: number | undefined;\r\n    let minPointP: Point3d | undefined;\r\n    for (const seed of seeds) {\r\n      newtonSearcher.setX(seed);\r\n      if (newtonSearcher.runIterations()) {\r\n        const curvePFraction = newtonSearcher.getX();\r\n        if (this.acceptFraction(curvePFraction)) {\r\n          const pointP = curveP.fractionToPoint(curvePFraction);\r\n          const closeApproachLength = pointP.distanceSquaredXY(pointQ);\r\n          if (closeApproachLength < minCloseApproachLength) {\r\n            minCloseApproachLength = closeApproachLength;\r\n            minCurvePFraction = curvePFraction;\r\n            minPointP = pointP;\r\n          }\r\n        }\r\n      }\r\n    }\r\n    if (minCurvePFraction && minPointP)\r\n      return CurveLocationDetail.createCurveFractionPoint(curveP, minCurvePFraction, minPointP);\r\n    return undefined;\r\n  }\r\n  /** Find the closest approach between pointA and cpB. Add the approach if it's within fB0 and fB1. */\r\n  private testAndRecordProjection(\r\n    cpA: CurvePrimitive, fA: number, pointA: Point3d, cpB: CurvePrimitive, fB0: number, fB1: number, reversed: boolean,\r\n  ) {\r\n    const detail = this.getPointCurveClosestApproachXYNewton(cpB, pointA);\r\n    if (detail) {\r\n      const fB = Geometry.restrictToInterval(detail.fraction, fB0, fB1);\r\n      if (fB === detail.fraction) { // if fraction is within fB0 and fB1\r\n        this.testAndRecordPointPairApproach(cpA, fA, pointA, cpB, detail.fraction, detail.point, reversed);\r\n      }\r\n    }\r\n  }\r\n  /**\r\n   * Compute intersection of two line segments.\r\n   * Filter by extension rules.\r\n   * Record with fraction mapping.\r\n   * * The fraction mappings allow portions of a linestring to be passed here.\r\n   */\r\n  private computeSegmentSegment3D(\r\n    cpA: CurvePrimitive,\r\n    pointA0: Point3d,\r\n    fractionA0: number,\r\n    pointA1: Point3d,\r\n    fractionA1: number,\r\n    cpB: CurvePrimitive,\r\n    pointB0: Point3d,\r\n    fractionB0: number,\r\n    pointB1: Point3d,\r\n    fractionB1: number,\r\n    reversed: boolean,\r\n  ): void {\r\n    // compute a pair with fractions local to segments\r\n    const approach = CurveCurveCloseApproachXY.segmentSegmentBoundedApproach(\r\n      pointA0, pointA1, pointB0, pointB1, this._maxDistanceSquared,\r\n    );\r\n    // adjust the pair to refer to input curves and global fractions, then record it if new\r\n    if (approach) {\r\n      approach.detailA.setCurve(cpA);\r\n      approach.detailB.setCurve(cpB);\r\n      this.capturePairWithLocalFractions(approach, cpA, fractionA0, fractionA1, cpB, fractionB0, fractionB1, reversed);\r\n    }\r\n  }\r\n  /** Low level dispatch of segment with segment. */\r\n  private dispatchSegmentSegment(\r\n    cpA: CurvePrimitive,\r\n    pointA0: Point3d,\r\n    fractionA0: number,\r\n    pointA1: Point3d,\r\n    fractionA1: number,\r\n    cpB: CurvePrimitive,\r\n    pointB0: Point3d,\r\n    fractionB0: number,\r\n    pointB1: Point3d,\r\n    fractionB1: number,\r\n    reversed: boolean,\r\n  ): void {\r\n    this.computeSegmentSegment3D(\r\n      cpA, pointA0, fractionA0, pointA1, fractionA1,\r\n      cpB, pointB0, fractionB0, pointB1, fractionB1,\r\n      reversed,\r\n    );\r\n  }\r\n  /**\r\n   * Low level dispatch of line segment with arc.\r\n   * Find close approaches within maxDistance between a line segments (pointA0, pointA1) and an arc.\r\n   * To consider:\r\n   * 1) intersection between arc and segment.\r\n   * 2) endpoints to endpoints or endpoints projection to the other curve.\r\n   * 3) arc tangent parallel to line segment (or line string).\r\n   * @param cpA curve A (line segment or line string; if it is a line string, then the fractions must specify a segment)\r\n   * @param pointA0 start point of the segment\r\n   * @param fractionA0 fraction of the start of the segment\r\n   * @param pointA1 end point of the segment\r\n   * @param fractionA1 fraction of the end of the segment\r\n   * @param arc the arc\r\n   * @param reversed whether to reverse the details in the pair (e.g., so that detailB refers to geometryB).\r\n   */\r\n  private dispatchSegmentArc(\r\n    cpA: CurvePrimitive,\r\n    pointA0: Point3d,\r\n    fractionA0: number,\r\n    pointA1: Point3d,\r\n    fractionA1: number,\r\n    arc: Arc3d,\r\n    reversed: boolean,\r\n  ): void {\r\n    // 1) intersection between arc and line segment (or string).\r\n    // Suppose:\r\n    // Arc: X = C + cU + sV where c = cos(theta) and s = sin(theta)\r\n    // Line: contains points A0 and A1\r\n    // The arc intersects the line at point X if det(A0, A1, X) = 0 with homogeneous xyw points and vectors.\r\n    // With equational X:  det(A0, A1, C) + c*det(A0, A1, U) + s*det(A0, A1, V) = 0.\r\n    // solve for theta.\r\n    // evaluate points.\r\n    // project back to line.\r\n    let intersectionFound = false;\r\n    const data = arc.toTransformedVectors();\r\n    const pointA0Local = pointA0;\r\n    const pointA1Local = pointA1;\r\n    const alpha = Geometry.tripleProductXYW(pointA0Local, 1, pointA1Local, 1, data.center, 1); //  det(A0, A1, C)\r\n    const beta = Geometry.tripleProductXYW(pointA0Local, 1, pointA1Local, 1, data.vector0, 0); // det(A0, A1, U)\r\n    const gamma = Geometry.tripleProductXYW(pointA0Local, 1, pointA1Local, 1, data.vector90, 0); // det(A0, A1, V)\r\n    const cosines = new GrowableFloat64Array(2);\r\n    const sines = new GrowableFloat64Array(2);\r\n    const radians = new GrowableFloat64Array(2);\r\n    const numRoots = AnalyticRoots.appendImplicitLineUnitCircleIntersections( // solve the equation\r\n      alpha, beta, gamma, cosines, sines, radians,\r\n    );\r\n    for (let i = 0; i < numRoots; i++) {\r\n      const arcPoint = data.center.plus2Scaled(\r\n        data.vector0, cosines.atUncheckedIndex(i), data.vector90, sines.atUncheckedIndex(i),\r\n      );\r\n      const arcFraction = data.sweep.radiansToSignedPeriodicFraction(radians.atUncheckedIndex(i));\r\n      const lineFraction = SmallSystem.lineSegment3dXYClosestPointUnbounded(pointA0Local, pointA1Local, arcPoint);\r\n      // only add if the point is within the start and end fractions of both line segment and arc\r\n      if (lineFraction !== undefined && this.acceptFraction(lineFraction) && this.acceptFraction(arcFraction)) {\r\n        this.recordPointWithLocalFractions(\r\n          lineFraction, cpA, fractionA0, fractionA1, arcFraction, arc, 0, 1, reversed,\r\n        );\r\n        intersectionFound = true;\r\n      }\r\n    }\r\n    if (intersectionFound)\r\n      return;\r\n    // 2) endpoints to endpoints or endpoints projection to the other curve.\r\n    this.testAndRecordFractionalPairApproach(cpA, fractionA0, fractionA1, true, arc, 0, 1, true, reversed);\r\n    // 3) arc tangent parallel to line segment (or string).\r\n    // If line does not intersect the arc, then the closest (and/or the furthest) point on arc to the line is a\r\n    // point where the tangent line on arc at that point is parallel to the line.\r\n    const dotUT = data.vector0.crossProductStartEndXY(pointA0, pointA1);\r\n    const dotVT = data.vector90.crossProductStartEndXY(pointA0, pointA1);\r\n    const parallelRadians = Math.atan2(dotVT, dotUT);\r\n    for (const radians1 of [parallelRadians, parallelRadians + Math.PI]) {\r\n      const arcPoint = data.center.plus2Scaled(data.vector0, Math.cos(radians1), data.vector90, Math.sin(radians1));\r\n      const arcFraction = data.sweep.radiansToSignedPeriodicFraction(radians1);\r\n      const lineFraction = SmallSystem.lineSegment3dXYClosestPointUnbounded(pointA0Local, pointA1Local, arcPoint);\r\n      // only add if the point is within the start and end fractions of both line segment and arc\r\n      if (lineFraction !== undefined && this.acceptFraction(lineFraction) && this.acceptFraction(arcFraction)) {\r\n        this.recordPointWithLocalFractions(\r\n          lineFraction, cpA, fractionA0, fractionA1, arcFraction, arc, 0, 1, reversed,\r\n        );\r\n      }\r\n    }\r\n  }\r\n  /** Solve Newton for 2 arcs and the given newtonEvaluator. */\r\n  private solveArcArcNewton(\r\n    curveP: Arc3d, curveQ: Arc3d, reversed: boolean, newtonEvaluator: NewtonEvaluatorRRtoRRD,\r\n  ): void {\r\n    const seedsU = [0.2, 0.4, 0.6, 0.8];  // HEURISTIC: 2 arcs have up to 4 perpendiculars/intersections\r\n    const seedsV = [0.2, 0.4, 0.6, 0.8];\r\n    const newtonSearcher = new Newton2dUnboundedWithDerivative(newtonEvaluator);\r\n    for (const seedU of seedsU) {\r\n      for (const seedV of seedsV) {\r\n        newtonSearcher.setUV(seedU, seedV);\r\n        if (newtonSearcher.runIterations()) {\r\n          const curvePFraction = newtonSearcher.getU();\r\n          const curveQFraction = newtonSearcher.getV();\r\n          if (this.acceptFraction(curvePFraction) && this.acceptFraction(curveQFraction)) {\r\n            this.recordPointWithLocalFractions(\r\n              curvePFraction, curveP, 0, 1, curveQFraction, curveQ, 0, 1, reversed,\r\n            );\r\n          }\r\n        }\r\n      }\r\n    }\r\n  }\r\n  /** Find and store perpendicular line between 2 arcs. */\r\n  private findPerpLineXYArcArcNewton(curveP: Arc3d, curveQ: Arc3d, reversed: boolean): void {\r\n    const newtonEvaluator = new CurveCurveCloseApproachXYRRtoRRD(curveP, curveQ);\r\n    this.solveArcArcNewton(curveP, curveQ, reversed, newtonEvaluator);\r\n  }\r\n  /** Low level dispatch of arc with Arc3d. */\r\n  private dispatchArcArc(cpA: Arc3d, cpB: Arc3d, reversed: boolean): void {\r\n    const rangeA = cpA.range();\r\n    const rangeB = cpB.range();\r\n    rangeA.expandInPlace(this._maxDistanceToAccept!);\r\n    if (!rangeB.intersectsRangeXY(rangeA))\r\n      return;\r\n    // 1) endpoints to endpoints or endpoints projection to the other curve\r\n    this.testAndRecordFractionalPairApproach(cpA, 0, 1, true, cpB, 0, 1, true, reversed);\r\n    // 2) perpendicular line between 2 arcs (includes intersections)\r\n    this.findPerpLineXYArcArcNewton(cpA, cpB, reversed);\r\n  }\r\n  /** Low level dispatch of arc with (beziers of) a bspline curve */\r\n  private dispatchArcBsplineCurve3d(cpA: Arc3d, cpB: BSplineCurve3d, reversed: boolean): void {\r\n    const ls = LineString3d.create();\r\n    cpB.emitStrokes(ls);\r\n    this.computeArcLineString(cpA, ls, reversed);\r\n  }\r\n  /** Low level dispatch of (beziers of) a bspline curve with (beziers of) a bspline curve */\r\n  private dispatchBSplineCurve3dBSplineCurve3d(\r\n    bcurveA: BSplineCurve3dBase, bcurveB: BSplineCurve3dBase, reversed: boolean,\r\n  ): void {\r\n    const lsA = LineString3d.create();\r\n    bcurveA.emitStrokes(lsA);\r\n    const lsB = LineString3d.create();\r\n    bcurveB.emitStrokes(lsB);\r\n    this.computeLineStringLineString(lsA, lsB, reversed);\r\n  }\r\n  /** Low level dispatch of linestring with (beziers of) a bspline curve */\r\n  public dispatchLineStringBSplineCurve(lsA: LineString3d, curveB: BSplineCurve3d, reversed: boolean): any {\r\n    const lsB = LineString3d.create();\r\n    curveB.emitStrokes(lsB);\r\n    this.computeLineStringLineString(lsA, lsB, reversed);\r\n  }\r\n  /** Low level dispatch of segment with (beziers of) a bspline curve */\r\n  public dispatchSegmentBsplineCurve(segA: LineSegment3d, curveB: BSplineCurve3d, reversed: boolean): any {\r\n    const lsB = LineString3d.create();\r\n    curveB.emitStrokes(lsB);\r\n    this.computeSegmentLineString(segA, lsB, reversed);\r\n  }\r\n  /** Detail computation for segment approaching linestring. */\r\n  public computeSegmentLineString(segA: LineSegment3d, lsB: LineString3d, reversed: boolean): void {\r\n    const numB = lsB.numPoints();\r\n    const deltaFracB = Geometry.safeDivideFraction(1, numB - 1, 0);\r\n    const pointA0 = segA.point0Ref;\r\n    const pointA1 = segA.point1Ref;\r\n    const pointB0 = CurveCurveCloseApproachXY._workPointBB0;\r\n    const pointB1 = CurveCurveCloseApproachXY._workPointBB1;\r\n    for (let i = 0; i < numB - 1; ++i) {\r\n      const fB0 = i * deltaFracB; // global linestring fractions\r\n      const fB1 = (i + 1 === numB - 1) ? 1.0 : (i + 1) * deltaFracB;  // make sure we nail the end fraction\r\n      lsB.packedPoints.getPoint3dAtUncheckedPointIndex(i, pointB0);\r\n      lsB.packedPoints.getPoint3dAtUncheckedPointIndex(i + 1, pointB1);\r\n      this.dispatchSegmentSegment(segA, pointA0, 0.0, pointA1, 1.0, lsB, pointB0, fB0, pointB1, fB1, reversed);\r\n    }\r\n  }\r\n  /** Detail computation for arc approaching linestring. */\r\n  public computeArcLineString(arcA: Arc3d, lsB: LineString3d, reversed: boolean): any {\r\n    const rangeA = arcA.range();\r\n    const rangeB = lsB.range();\r\n    rangeA.expandInPlace(this._maxDistanceToAccept!);\r\n    if (!rangeB.intersectsRangeXY(rangeA))\r\n      return;\r\n    const pointB0 = CurveCurveCloseApproachXY._workPointBB0;\r\n    const pointB1 = CurveCurveCloseApproachXY._workPointBB1;\r\n    const numB = lsB.numPoints();\r\n    if (numB > 1) {\r\n      const dfB = 1.0 / (numB - 1);\r\n      let fB0;\r\n      let fB1;\r\n      fB0 = 0.0;\r\n      lsB.pointAt(0, pointB0);\r\n      for (let ib = 1; ib < numB; ib++, pointB0.setFrom(pointB1), fB0 = fB1) {\r\n        lsB.pointAt(ib, pointB1);\r\n        fB1 = ib * dfB;\r\n        this.dispatchSegmentArc(lsB, pointB0, fB0, pointB1, fB1, arcA, !reversed);\r\n      }\r\n    }\r\n    return undefined;\r\n  }\r\n  /** Low level dispatch of curve collection. */\r\n  private dispatchCurveCollection(geomA: AnyCurve, geomAHandler: (geomA: any) => any): void {\r\n    const geomB = this._geometryB;  // save\r\n    if (!geomB || !geomB.children || !(geomB instanceof CurveCollection))\r\n      return;\r\n    for (const child of geomB.children) {\r\n      this.resetGeometry(child);\r\n      geomAHandler(geomA);\r\n    }\r\n    this._geometryB = geomB;  // restore\r\n  }\r\n  /** Low level dispatch to geomA given a CurveChainWithDistanceIndex in geometryB. */\r\n  private dispatchCurveChainWithDistanceIndex(geomA: AnyCurve, geomAHandler: (geomA: any) => any): void {\r\n    if (!this._geometryB || !(this._geometryB instanceof CurveChainWithDistanceIndex))\r\n      return;\r\n    if (geomA instanceof CurveChainWithDistanceIndex) {\r\n      assert(!!\"call handleCurveChainWithDistanceIndex(geomA) instead\");\r\n      return;\r\n    }\r\n    const index0 = this._results.length;\r\n    const geomB = this._geometryB;  // save\r\n    for (const child of geomB.path.children) {\r\n      this.resetGeometry(child);\r\n      geomAHandler(geomA);\r\n    }\r\n    this.resetGeometry(geomB);  // restore\r\n    this._results = CurveChainWithDistanceIndex.convertChildDetailToChainDetail(this._results, index0, undefined, geomB, true);\r\n  }\r\n  /** Double dispatch handler for strongly typed segment. */\r\n  public override handleLineSegment3d(segmentA: LineSegment3d): any {\r\n    if (this._geometryB instanceof LineSegment3d) {\r\n      const segmentB = this._geometryB;\r\n      this.dispatchSegmentSegment(\r\n        segmentA, segmentA.point0Ref, 0.0, segmentA.point1Ref, 1.0,\r\n        segmentB, segmentB.point0Ref, 0.0, segmentB.point1Ref, 1.0,\r\n        false,\r\n      );\r\n    } else if (this._geometryB instanceof LineString3d) {\r\n      this.computeSegmentLineString(segmentA, this._geometryB, false);\r\n    } else if (this._geometryB instanceof Arc3d) {\r\n      this.dispatchSegmentArc(segmentA, segmentA.point0Ref, 0.0, segmentA.point1Ref, 1.0, this._geometryB, false);\r\n    } else if (this._geometryB instanceof BSplineCurve3d) {\r\n      this.dispatchSegmentBsplineCurve(segmentA, this._geometryB, false);\r\n    } else if (this._geometryB instanceof CurveCollection) {\r\n      this.dispatchCurveCollection(segmentA, this.handleLineSegment3d.bind(this));\r\n    } else if (this._geometryB instanceof CurveChainWithDistanceIndex) {\r\n      this.dispatchCurveChainWithDistanceIndex(segmentA, this.handleLineSegment3d.bind(this));\r\n    }\r\n    return undefined;\r\n  }\r\n  /**\r\n   * Set bits for comparison to range xy\r\n   * * bit 0x01 => x smaller than range.low.x\r\n   * * bit 0x02 => x larger than range.high.x\r\n   * * bit 0x04 => y smaller than range.low.y\r\n   * * bit 0x08 => y larger than range.high.y\r\n   * * If we divide XY plane into 9 areas using the range, the function returns 0 for points\r\n   * inside the range. Below is other binary numbers returned by the function for all 9 areas:\r\n   *   1001 | 1000 | 1010\r\n   *   ------------------\r\n   *    1   |  0   |  10\r\n   *   ------------------\r\n   *   101  | 100  | 110\r\n   * @param xy point to test\r\n   * @param range range for comparison\r\n   */\r\n  private classifyBitsPointRangeXY(x: number, y: number, range: Range3d): number {\r\n    let result = 0;\r\n    if (x < range.low.x)\r\n      result = 0x01;\r\n    else if (x > range.high.x)\r\n      result = 0x02;\r\n    // note the OR operation\r\n    if (y < range.low.y)\r\n      result |= 0x04;\r\n    else if (y > range.high.y)\r\n      result |= 0x08;\r\n    return result;\r\n  }\r\n  /** Low level dispatch of line string with line string. */\r\n  private computeLineStringLineString(lsA: LineString3d, lsB: LineString3d, reversed: boolean): void {\r\n    const rangeA = lsA.range();\r\n    const rangeB = lsB.range();\r\n    rangeA.expandInPlace(this._maxDistanceToAccept!);\r\n    if (!rangeB.intersectsRangeXY(rangeA))\r\n      return;\r\n    let bitB0: number;\r\n    let bitB1: number;\r\n    const rangeA1 = Range3d.createNull();\r\n    const pointA0 = CurveCurveCloseApproachXY._workPointAA0;\r\n    const pointA1 = CurveCurveCloseApproachXY._workPointAA1;\r\n    const pointB0 = CurveCurveCloseApproachXY._workPointBB0;\r\n    const pointB1 = CurveCurveCloseApproachXY._workPointBB1;\r\n    const numA = lsA.numPoints();\r\n    const numB = lsB.numPoints();\r\n    if (numA > 1 && numB > 1) {\r\n      const dfA = 1.0 / (numA - 1);\r\n      const dfB = 1.0 / (numB - 1);\r\n      let fA0 = 0.0;\r\n      let fA1, fB0, fB1;\r\n      lsA.pointAt(0, pointA0);\r\n      for (let ia = 1; ia < numA; ia++, pointA0.setFrom(pointA1), fA0 = fA1) {\r\n        fA1 = ia * dfA;\r\n        fB0 = 0.0;\r\n        lsA.pointAt(ia, pointA1);\r\n        // rangeA1 is around line segment [A0,A1] expanded by max distance\r\n        rangeA1.setNull();\r\n        rangeA1.extendPoint(pointA0);\r\n        rangeA1.extendPoint(pointA1);\r\n        rangeA1.expandInPlace(this._maxDistanceToAccept!);\r\n        if (rangeA1.intersectsRangeXY(rangeB)) {\r\n          lsB.pointAt(0, pointB0);\r\n          bitB0 = this.classifyBitsPointRangeXY(pointB0.x, pointB0.y, rangeA1);\r\n          for (let ib = 1; ib < numB; ib++, pointB0.setFrom(pointB1), fB0 = fB1, bitB0 = bitB1) {\r\n            lsB.pointAt(ib, pointB1);\r\n            bitB1 = this.classifyBitsPointRangeXY(pointB1.x, pointB1.y, rangeA1);\r\n            fB1 = ib * dfB;\r\n            // DO NOT study the segment in detail if both bitB bits are on for any of the 4 planes\r\n            // (i.e., no intersection between rangeA1 and the range around line segment [B0,B1])\r\n            if ((bitB0 & bitB1) === 0)\r\n              this.dispatchSegmentSegment(lsA, pointA0, fA0, pointA1, fA1, lsB, pointB0, fB0, pointB1, fB1, reversed);\r\n          }\r\n        }\r\n      }\r\n    }\r\n  }\r\n  /** Double dispatch handler for strongly typed linestring. */\r\n  public override handleLineString3d(lsA: LineString3d): any {\r\n    if (this._geometryB instanceof LineString3d) {\r\n      const lsB = this._geometryB;\r\n      this.computeLineStringLineString(lsA, lsB, false);\r\n    } else if (this._geometryB instanceof LineSegment3d) {\r\n      this.computeSegmentLineString(this._geometryB, lsA, true);\r\n    } else if (this._geometryB instanceof Arc3d) {\r\n      this.computeArcLineString(this._geometryB, lsA, true);\r\n    } else if (this._geometryB instanceof BSplineCurve3d) {\r\n      this.dispatchLineStringBSplineCurve(lsA, this._geometryB, false);\r\n    } else if (this._geometryB instanceof CurveCollection) {\r\n      this.dispatchCurveCollection(lsA, this.handleLineString3d.bind(this));\r\n    } else if (this._geometryB instanceof CurveChainWithDistanceIndex) {\r\n      this.dispatchCurveChainWithDistanceIndex(lsA, this.handleLineString3d.bind(this));\r\n    }\r\n    return undefined;\r\n  }\r\n  /** Double dispatch handler for strongly typed arc. */\r\n  public override handleArc3d(arc0: Arc3d): any {\r\n    if (this._geometryB instanceof LineSegment3d) {\r\n      this.dispatchSegmentArc(\r\n        this._geometryB, this._geometryB.point0Ref, 0.0, this._geometryB.point1Ref, 1.0, arc0, true,\r\n      );\r\n    } else if (this._geometryB instanceof LineString3d) {\r\n      this.computeArcLineString(arc0, this._geometryB, false);\r\n    } else if (this._geometryB instanceof Arc3d) {\r\n      this.dispatchArcArc(arc0, this._geometryB, false);\r\n    } else if (this._geometryB instanceof BSplineCurve3d) {\r\n      this.dispatchArcBsplineCurve3d(arc0, this._geometryB, false);\r\n    } else if (this._geometryB instanceof CurveCollection) {\r\n      this.dispatchCurveCollection(arc0, this.handleArc3d.bind(this));\r\n    } else if (this._geometryB instanceof CurveChainWithDistanceIndex) {\r\n      this.dispatchCurveChainWithDistanceIndex(arc0, this.handleArc3d.bind(this));\r\n    }\r\n    return undefined;\r\n  }\r\n  /** Double dispatch handler for strongly typed bspline curve. */\r\n  public override handleBSplineCurve3d(curve: BSplineCurve3d): any {\r\n    if (this._geometryB instanceof LineSegment3d) {\r\n      this.dispatchSegmentBsplineCurve(this._geometryB, curve, true);\r\n    } else if (this._geometryB instanceof LineString3d) {\r\n      this.dispatchLineStringBSplineCurve(this._geometryB, curve, true);\r\n    } else if (this._geometryB instanceof Arc3d) {\r\n      this.dispatchArcBsplineCurve3d(this._geometryB, curve, true);\r\n    } else if (this._geometryB instanceof BSplineCurve3dBase) {\r\n      this.dispatchBSplineCurve3dBSplineCurve3d(curve, this._geometryB, false);\r\n    } else if (this._geometryB instanceof CurveCollection) {\r\n      this.dispatchCurveCollection(curve, this.handleBSplineCurve3d.bind(this));\r\n    } else if (this._geometryB instanceof CurveChainWithDistanceIndex) {\r\n      this.dispatchCurveChainWithDistanceIndex(curve, this.handleBSplineCurve3d.bind(this));\r\n    }\r\n    return undefined;\r\n  }\r\n  /** Double dispatch handler for strongly typed CurveChainWithDistanceIndex. */\r\n  public override handleCurveChainWithDistanceIndex(chain: CurveChainWithDistanceIndex): any {\r\n    super.handleCurveChainWithDistanceIndex(chain);\r\n    // if _geometryB is also a CurveChainWithDistanceIndex, it will already have been converted by dispatchCurveChainWithDistanceIndex\r\n    this._results = CurveChainWithDistanceIndex.convertChildDetailToChainDetail(this._results, 0, chain, undefined, true);\r\n  }\r\n  /** Double dispatch handler for strongly typed homogeneous bspline curve .. */\r\n  public override handleBSplineCurve3dH(_curve: BSplineCurve3dH): any {\r\n    /*\r\n    //NEEDS WORK -- make \"dispatch\" methods tolerant of both 3d and 3dH.\r\n    // \"easy\" if both present BezierCurve3dH span loaders\r\n    if (this._geometryB instanceof LineSegment3d) {\r\n      this.dispatchSegmentBsplineCurve(\r\n        this._geometryB, this._extendB, this._geometryB.point0Ref, 0.0, this._geometryB.point1Ref, 1.0, this._extendB,\r\n        curve, this._extendA, true);\r\n    } else if (this._geometryB instanceof LineString3d) {\r\n      this.dispatchLineStringBSplineCurve(this._geometryB, this._extendB, curve, this._extendA, true);\r\n    } else if (this._geometryB instanceof Arc3d) {\r\n      this.dispatchArcBsplineCurve3d(this._geometryB, this._extendB, curve, this._extendA, true);\r\n    }\r\n    */\r\n    return undefined;\r\n  }\r\n}\r\n"]}

package/lib/cjs/curve/internalContexts/CurveCurveIntersectXY.d.ts

@@ -1,9 +1,7 @@
-import { BezierCurveBase } from "../../bspline/BezierCurveBase";
 import { BSplineCurve3d } from "../../bspline/BSplineCurve";
 import { BSplineCurve3dH } from "../../bspline/BSplineCurve3dH";
 import { RecurseToCurvesGeometryHandler } from "../../geometry3d/GeometryHandler";
 import { Matrix4d } from "../../geometry4d/Matrix4d";
-import { NewtonEvaluatorRRtoRRD } from "../../numerics/Newton";
 import { Arc3d } from "../Arc3d";
 import { CurveChainWithDistanceIndex } from "../CurveChainWithDistanceIndex";
 import { CurveLocationDetailPair } from "../CurveLocationDetail";
@@ -11,19 +9,6 @@ import { CurvePrimitive } from "../CurvePrimitive";
 import { AnyCurve } from "../CurveTypes";
 import { LineSegment3d } from "../LineSegment3d";
 import { LineString3d } from "../LineString3d";
-/**
- * Private class for refining bezier-bezier intersections.
- * * The inputs are assumed pre-transformed so that the target condition is to match x and y coordinates.
- * @internal
- */
-export declare class BezierBezierIntersectionXYRRToRRD extends NewtonEvaluatorRRtoRRD {
-    private _curveA;
-    private _curveB;
-    private _rayA;
-    private _rayB;
-    constructor(curveA: BezierCurveBase, curveB: BezierCurveBase);
-    evaluate(fractionA: number, fractionB: number): boolean;
-}
 /**
  * Handler class for XY intersections between _geometryB and another geometry.
  * * Instances are initialized and called from CurveCurve.
@@ -105,7 +90,20 @@ export declare class CurveCurveIntersectXY extends RecurseToCurvesGeometryHandle
     private computeSegmentSegment3DH;
     private dispatchSegmentSegment;
     private dispatchSegmentArc;
+    /**
+     * Compute the intersection of two xy-arcs.
+     * * Each matrix has [U V C] in (x,y,w) form from homogeneous projection (local to world).
+     * * Arcs are ordered so that matrixA is better conditioned.
+     */
     private dispatchArcArcThisOrder;
+    /**
+     * We have 2 xy-arcs.
+     * 1- We pick the arc that is closest to circular (larger condition number is closer to circular).
+     * 2- Transform it to local coords, where it becomes the unit xy-circle.
+     * 3- Use the same map to transform the other arc.
+     * 4- Find the intersection of arc and unit circle.
+     * 5- Convert intersection angles to fractions and record intersections.
+     */
     private dispatchArcArc;
     private dispatchArcBsplineCurve3d;
     /** Apply the transformation to bezier curves. Optionally construct ranges. */
@@ -132,7 +130,7 @@ export declare class CurveCurveIntersectXY extends RecurseToCurvesGeometryHandle
     private static setTransformedWorkPoints;
     /** Low level dispatch of curve collection. */
     private dispatchCurveCollection;
-    /** Low level dispatch to geomA given a CurveChainWithDistanceIndex in geometryB. */
+    /** Low level dispatch of CurveChainWithDistanceIndex. */
     private dispatchCurveChainWithDistanceIndex;
     /** Double dispatch handler for strongly typed segment. */
     handleLineSegment3d(segmentA: LineSegment3d): any;