@itwin/core-geometry 5.3.0-dev.2 → 5.3.0-dev.23

package/lib/esm/curve/spiral/IntegratedSpiral3d.js

@@ -21,35 +21,40 @@ import { Ray3d } from "../../geometry3d/Ray3d";
 import { Plane3dByOriginAndVectors } from "../../geometry3d/Plane3dByOriginAndVectors";
 /**
  * An IntegratedSpiral3d is a curve defined by integrating its curvature.
- * * The first integral of curvature (with respect to distance along the curve) is the bearing angle (in radians)
+ * * The first integral of curvature (with respect to distance along the curve) is the bearing angle (in radians).
  * * Integrating (cos(theta), sin(theta)) gives displacement from the start point, and thus the actual curve position.
- * * The curvature functions of interest are all symmetric snap functions in the NormalizedTransition class.
- * * `TransitionConditionalProperties` implements the computations of the interrelationship of radii, bearing, and length.
+ * * The curvature functions of interest are all symmetric snap functions in the [[NormalizedTransition]] class.
+ * * [[TransitionConditionalProperties]] implements the computations of the interrelationship of radii, bearing, and length.
  * @public
  */
 export class IntegratedSpiral3d extends TransitionSpiral3d {
-    /** String name for schema properties */
+    /** String name for schema properties. */
     curvePrimitiveType = "transitionSpiral";
-    /** start and end radii as a Segment1d */
+    /** Start and end radii as a Segment1d. */
     radius01;
-    /** start and end bearings as an AngleSweep */
+    /** Start and end bearings as an AngleSweep. */
     bearing01;
-    /** stroked approximation of entire spiral. */
+    /** Stroked approximation of entire spiral. */
     _globalStrokes;
-    /** stroked approximation of active spiral.
-     * * Same count as global -- possibly overly fine, but it gives some consistency between same clothoid constructed as partial versus complete.
-     * * If no trimming, this points to the same place as the _globalStrokes !!!  Don't double transform!!!
+    /**
+     * Stroked approximation of active spiral.
+     * * Same count as global -- possibly overly fine, but it gives some consistency between same clothoid constructed
+     * as partial versus complete.
+     * * If no trimming, this points to the same place as the _globalStrokes.  DO NOT double transform.
      */
     _activeStrokes;
-    /** Return the internal stroked form of the (possibly partial) spiral   */
-    get activeStrokes() { return this._activeStrokes !== undefined ? this._activeStrokes : this._globalStrokes; }
+    /** Return the internal stroked form of the (possibly partial) spiral. */
+    get activeStrokes() {
+        return this._activeStrokes !== undefined ? this._activeStrokes : this._globalStrokes;
+    }
+    /** Evaluator for transition. */
     _evaluator;
-    /** Total curve arc length (computed) */
+    /** Total curve arc length (computed). */
     _arcLength01;
-    /** Curvatures (inverse radii) at start and end */
+    /** Curvatures (inverse radii) at start and end. */
     _curvature01;
-    /** evaluator for transition */
-    // constructor demands all bearing, radius, and length data -- caller determines usual dependency of "any 4 determine the 5th"
+    // constructor demands all bearing, radius, and length data.
+    // caller determines usual dependency of "any 4 determine the 5th".
     constructor(spiralType, evaluator, radius01, bearing01, activeFractionInterval, localToWorld, arcLength, properties) {
         super(spiralType, localToWorld, activeFractionInterval, properties);
         this._evaluator = evaluator;
@@ -57,29 +62,33 @@ export class IntegratedSpiral3d extends TransitionSpiral3d {
         this.bearing01 = bearing01;
         this._arcLength01 = arcLength;
         this._globalStrokes = LineString3d.create();
-        // initialize for compiler -- but this will be recomputed in refreshComputeProperties ...
+        // initialize for compiler; but this will be recomputed in refreshComputeProperties
         this._curvature01 = Segment1d.create(0, 1);
         this.refreshComputedProperties();
     }
-    /** default spiral type name. (clothoid) */
+    /** Default spiral type name: clothoid. */
     static defaultSpiralType = "clothoid";
-    /** use the integrated function to return an angle at fractional position. */
+    /** Use the integrated function to return an angle at fractional position. */
     globalFractionToBearingRadians(fraction) {
         const areaFraction = this._evaluator.fractionToArea(fraction);
         const dx = this._arcLength01;
-        return this.bearing01.startRadians + areaFraction * dx * this._curvature01.signedDelta() + fraction * this._curvature01.x0 * dx;
+        return this.bearing01.startRadians
+            + areaFraction * dx * this._curvature01.signedDelta()
+            + fraction * this._curvature01.x0 * dx;
     }
-    /** use the integrated function to return an angle at fractional position. */
+    /** Use the integrated function to return an angle at fractional position. */
     globalFractionToCurvature(fraction) {
         const f = this._evaluator.fractionToCurvatureFraction(fraction);
         return this._curvature01.fractionToPoint(f);
     }
-    /** Return the bearing at given fraction of the active interval .... */
+    /** Return the bearing at given fraction of the active interval. */
     fractionToBearingRadians(activeFraction) {
         const fraction = this.activeFractionInterval.fractionToPoint(activeFraction);
-        return this.bearing01.startRadians + fraction * this._arcLength01 * (this._curvature01.x0 + 0.5 * fraction * (this._curvature01.x1 - this._curvature01.x0));
+        return this.bearing01.startRadians
+            + fraction * this._arcLength01 * (this._curvature01.x0 + 0.5 * fraction * (this._curvature01.x1 - this._curvature01.x0));
     }
-    /** Return the curvature at given fraction of the active interval ...
+    /**
+     * Return the curvature at given fraction of the active interval.
      * * The `undefined` result is to match the abstract class -- it cannot actually occur.
      */
     fractionToCurvature(activeFraction) {
@@ -96,14 +105,15 @@ export class IntegratedSpiral3d extends TransitionSpiral3d {
         IntegratedSpiral3d._gaussWeight = new Float64Array(5);
         IntegratedSpiral3d._gaussMapper = (xA, xB, xMapped, wMapped) => Quadrature.setupGauss5(xA, xB, xMapped, wMapped);
     }
-    /** Evaluate and sum the gauss quadrature formulas to integrate cos(theta), sin(theta) fractional subset of a reference length.
-     * (recall that theta is a nonlinear function of the fraction.)
+    /**
+     * Evaluate and sum the gauss quadrature formulas to integrate cos(theta), sin(theta) fractional subset of a reference
+     * length (recall that theta is a nonlinear function of the fraction).
      * * This is a single interval of gaussian integration.
-     * * The fraction is on the full spiral (not in the mapped active interval)
+     * * The fraction is on the full spiral (not in the mapped active interval).
      * @param xyz advancing integrated point.
-     * @param fractionA fraction at start of interval
+     * @param fractionA fraction at start of interval.
      * @param fractionB fraction at end of interval.
-     * @param unitArcLength length of curve for 0 to 1 fractional
+     * @param unitArcLength length of curve for 0 to 1 fractional.
      */
     fullSpiralIncrementalIntegral(xyz, fractionA, fractionB, applyMatrix) {
         const gaussFraction = IntegratedSpiral3d._gaussFraction;
@@ -124,7 +134,7 @@ export class IntegratedSpiral3d extends TransitionSpiral3d {
         else
             xyz.addXYZInPlace(dx, dy, 0.0);
     }
-    /** Recompute strokes */
+    /** Recompute strokes. */
     refreshComputedProperties() {
         this._curvature01 = Segment1d.create(TransitionSpiral3d.radiusToCurvature(this.radius01.x0), TransitionSpiral3d.radiusToCurvature(this.radius01.x1));
         this._globalStrokes.clear();
@@ -143,8 +153,9 @@ export class IntegratedSpiral3d extends TransitionSpiral3d {
             if (this._activeStrokes === undefined)
                 this._activeStrokes = LineString3d.create();
             this._activeStrokes.clear();
-            // finer strokes in the active interval ... same fraction step, but mapped
-            // This assumes factionToPoint acts normally within refreshComputedProperties -- that depends on the global strokes we just computed, but not on the active strokes
+            // The active interval has finer strokes: it's the same fraction step but mapped to a sub-interval.
+            // Below assumes fractionToPoint depends upon the global strokes we just computed, and not
+            // on the active strokes we are computing.
             for (let i = 0; i <= numInterval; i++) {
                 const localFraction = i * fractionStep;
                 this._activeStrokes.addPoint(this.fractionToPoint(localFraction));
@@ -153,10 +164,11 @@ export class IntegratedSpiral3d extends TransitionSpiral3d {
     }
     /**
      * Create a transition spiral with radius and bearing conditions.
-     * @param radius01 radius (inverse curvature) at start and end. (radius of zero means straight line)
-     * @param bearing01 bearing angles at start and end.  bearings are measured from the x axis, positive clockwise towards y axis
+     * @param radius01 radius (inverse curvature) at start and end (radius of zero means straight line).
+     * @param bearing01 bearing angles at start and end.  bearings are measured from the x axis, positive clockwise
+     * towards y axis.
      * @param activeFractionInterval fractional limits of the active portion of the spiral.
-     * @param localToWorld placement frame.  Fractional coordinate 0 is at the origin.
+     * @param localToWorld placement frame. Fractional coordinate 0 is at the origin.
      */
     static createRadiusRadiusBearingBearing(radius01, bearing01, activeFractionInterval, localToWorld, typeName) {
         const arcLength = TransitionSpiral3d.radiusRadiusSweepRadiansToArcLength(radius01.x0, radius01.x1, bearing01.sweepRadians);
@@ -169,14 +181,15 @@ export class IntegratedSpiral3d extends TransitionSpiral3d {
     }
     /**
      * Create a transition spiral.
-     * * Inputs must provide exactly 4 of the 5 values `[radius0,radius1,bearing0,bearing1,length`.
+     * * Inputs must provide exactly 4 of the 5 values `[radius0,radius1,bearing0,bearing1,length]`.
      * @param spiralType one of "clothoid", "bloss", "biquadratic", "cosine", "sine".  If undefined, "clothoid" is used.
-     * @param radius0 radius (or 0 for tangent to line) at start
-     * @param radius1 radius (or 0 for tangent to line) at end
+     * @param radius0 radius (or 0 for tangent to line) at start.
+     * @param radius1 radius (or 0 for tangent to line) at end.
      * @param bearing0 bearing, measured CCW from x axis at start.
      * @param bearing1 bearing, measured CCW from x axis at end.
-     * @param fractionInterval optional fractional interval for an "active" portion of the curve.   if omitted, the full [0,1] is used.
-     * @param localToWorld placement transform
+     * @param fractionInterval optional fractional interval for an "active" portion of the curve. If omitted, the full
+     * [0,1] is used.
+     * @param localToWorld placement transform.
      */
     static createFrom4OutOf5(spiralType, radius0, radius1, bearing0, bearing1, arcLength, fractionInterval, localToWorld) {
         if (spiralType === undefined)
@@ -203,11 +216,11 @@ export class IntegratedSpiral3d extends TransitionSpiral3d {
         this._arcLength01 = other._arcLength01;
         return this;
     }
-    /** Deep clone of this spiral */
+    /** Deep clone of this spiral. */
     clone() {
         return new IntegratedSpiral3d(this._spiralType, this._evaluator, this.radius01.clone(), this.bearing01.clone(), this.activeFractionInterval.clone(), this.localToWorld.clone(), this._arcLength01, this._designProperties?.clone());
     }
-    /** apply `transform` to this spiral's local to world transform. */
+    /** Apply `transform` to this spiral's local to world transform. */
     tryTransformInPlace(transformA) {
         const rigidData = this.applyRigidPartOfTransform(transformA);
         if (rigidData !== undefined) {
@@ -224,32 +237,50 @@ export class IntegratedSpiral3d extends TransitionSpiral3d {
     startPoint(result) {
         return this.activeStrokes.startPoint(result);
     }
-    /** return the spiral end point. */
+    /** Return the spiral end point. */
     endPoint(result) {
         return this.activeStrokes.endPoint(result);
     }
-    /** test if the local to world transform places the spiral xy plane into `plane` */
+    /** Test if the local to world transform places the spiral xy plane into `plane`. */
     isInPlane(plane) {
         return plane.isPointInPlane(this.localToWorld.origin)
             && Geometry.isSameCoordinate(0.0, this.localToWorld.matrix.dotColumnX(plane.getNormalRef()))
             && Geometry.isSameCoordinate(0.0, this.localToWorld.matrix.dotColumnY(plane.getNormalRef()));
     }
-    /** Return length of the spiral.  Because TransitionSpiral is parameterized directly in terms of distance along, this is a simple return value. */
-    quickLength() { return this.curveLength(); }
-    /** Return length of the spiral.  Because TransitionSpiral is parameterized directly in terms of distance along, this is a simple return value. */
-    curveLength() { return this._arcLength01 * (this._activeFractionInterval.absoluteDelta()); }
-    /** Return (unsigned) length of the spiral between fractions.  Because TransitionSpiral is parameterized directly in terms of distance along, this is a simple return value. */
+    /**
+     * Return length of the spiral.
+     * * Because TransitionSpiral is parameterized directly in terms of distance along, this is a simple return value.
+     */
+    quickLength() {
+        return this.curveLength();
+    }
+    /**
+     * Return length of the spiral.
+     * * Because TransitionSpiral is parameterized directly in terms of distance along, this is a simple return value.
+     */
+    curveLength() {
+        return this._arcLength01 * (this._activeFractionInterval.absoluteDelta());
+    }
+    /**
+     * Return (unsigned) length of the spiral between fractions.
+     * * Because TransitionSpiral is parameterized directly in terms of distance along, this is a simple return value.
+     */
     curveLengthBetweenFractions(fraction0, fraction1) {
         return this._arcLength01 * (this._activeFractionInterval.absoluteDelta() * Math.abs(fraction1 - fraction0));
     }
-    /** Test if `other` is an instance of `TransitionSpiral3d` */
-    isSameGeometryClass(other) { return other instanceof TransitionSpiral3d; }
-    /** Add strokes from this spiral to `dest`.
+    /** Test if `other` is an instance of `TransitionSpiral3d`. */
+    isSameGeometryClass(other) {
+        return other instanceof TransitionSpiral3d;
+    }
+    /**
+     * Add strokes from this spiral to `dest`.
      * * Linestrings will usually stroke as just their points.
      * * If maxEdgeLength is given, this will sub-stroke within the linestring -- not what we want.
      */
-    emitStrokes(dest, options) { this.activeStrokes.emitStrokes(dest, options); }
-    /** emit stroke fragments to `dest` handler. */
+    emitStrokes(dest, options) {
+        this.activeStrokes.emitStrokes(dest, options);
+    }
+    /** Emit stroke fragments to `dest` handler. */
     emitStrokableParts(dest, options) {
         const n = this.computeStrokeCountForOptions(options);
         dest.startParentCurvePrimitive(this);
@@ -264,8 +295,8 @@ export class IntegratedSpiral3d extends TransitionSpiral3d {
         dest.endParentCurvePrimitive(this);
     }
     /**
-     * return the stroke count required for given options.
-     * @param options StrokeOptions that determine count
+     * Return the stroke count required for given options.
+     * @param options StrokeOptions that determine count.
      */
     computeStrokeCountForOptions(options) {
         let numStroke;
@@ -280,8 +311,9 @@ export class IntegratedSpiral3d extends TransitionSpiral3d {
         }
         return numStroke;
     }
-    /** Reverse the active interval and active strokes.
-     * * Primary defining data remains unchanged !!!
+    /**
+     * Reverse the active interval and active strokes.
+     * * Primary defining data remains unchanged.
      */
     reverseInPlace() {
         this.activeFractionInterval.reverseInPlace();
@@ -319,10 +351,10 @@ export class IntegratedSpiral3d extends TransitionSpiral3d {
         }
         else {
             const clampedGlobalFraction = Geometry.clampToStartEnd(targetGlobalFraction, 0, 1);
-            const index0 = Math.trunc(clampedGlobalFraction * numStrokes); // This indexes the point to the left of the query.
+            const index0 = Math.trunc(clampedGlobalFraction * numStrokes); // this indexes the point to the left of the query
             const globalFraction0 = index0 / numStrokes;
             result = this._globalStrokes.packedPoints.getPoint3dAtUncheckedPointIndex(index0, result);
-            // GeometryCoreTestIO.consoleLog(" fractionToPoint ", activeFraction, this.activeFractionInterval, "( global integration " + globalFraction0 + " to " + globalFraction + ")", index0);
+            // GeometryCoreTestIO.consoleLog("fractionToPoint ", activeFraction, this.activeFractionInterval, "( global integration " + globalFraction0 + " to " + globalFraction + ")", index0);
             this.fullSpiralIncrementalIntegral(result, globalFraction0, targetGlobalFraction, true);
         }
         return result;
@@ -349,10 +381,11 @@ export class IntegratedSpiral3d extends TransitionSpiral3d {
         result.matrix.applyGivensColumnOp(0, 1, c, s);
         return result;
     }
-    /** Return a plane with
-     *
-     * * origin at fractional position along the curve
-     * * vectorU is the first derivative, i.e. tangent vector with length equal to the rate of change with respect to the fraction.
+    /**
+     * Return a plane with
+     * * origin at fractional position along the curve.
+     * * vectorU is the first derivative, i.e. tangent vector with length equal to the rate of change with respect to
+     * the fraction.
      * * vectorV is the second derivative, i.e.derivative of vectorU.
      */
     fractionToPointAnd2Derivatives(activeFraction, result) {
@@ -369,11 +402,11 @@ export class IntegratedSpiral3d extends TransitionSpiral3d {
         vectorY.scaleInPlace(this.globalFractionToCurvature(globalFraction));
         return Plane3dByOriginAndVectors.createCapture(origin, vectorX, vectorY, result);
     }
-    /** Second step of double dispatch:  call `handler.handleTransitionSpiral(this)` */
+    /** Second step of double dispatch:  call `handler.handleTransitionSpiral(this)`. */
     dispatchToGeometryHandler(handler) {
         return handler.handleTransitionSpiral(this);
     }
-    /** compare various coordinate quantities */
+    /** Compare various coordinate quantities. */
     isAlmostEqual(other) {
         if (other instanceof IntegratedSpiral3d) {
             return this.radius01.isAlmostEqual(other.radius01)

package/lib/esm/curve/spiral/IntegratedSpiral3d.js.map

@@ -1 +1 @@
-{"version":3,"file":"IntegratedSpiral3d.js","sourceRoot":"","sources":["../../../../src/curve/spiral/IntegratedSpiral3d.ts"],"names":[],"mappings":"AAAA;;;+FAG+F;AAC/F;;GAEG;AAEH,OAAO,EAAE,kBAAkB,EAAE,MAAM,sBAAsB,CAAC;AAC1D,OAAO,EAAE,SAAS,EAAE,MAAM,4BAA4B,CAAC;AACvD,OAAO,EAAE,UAAU,EAAE,MAAM,6BAA6B,CAAC;AACzD,OAAO,EAAE,SAAS,EAAE,MAAM,4BAA4B,CAAC;AACvD,OAAO,EAAE,YAAY,EAAE,MAAM,iBAAiB,CAAC;AAC/C,OAAO,EAAE,oBAAoB,EAAE,MAAM,wBAAwB,CAAC;AAC9D,OAAO,EAAE,+BAA+B,EAAE,MAAM,mCAAmC,CAAC;AACpF,OAAO,EAAE,UAAU,EAAE,MAAM,2BAA2B,CAAC;AACvD,OAAO,EAAE,OAAO,EAAE,MAAM,kCAAkC,CAAC;AAC3D,OAAO,EAAE,QAAQ,EAAE,MAAM,2BAA2B,CAAC;AAGrD,OAAO,EAAE,SAAS,EAAE,QAAQ,EAAE,MAAM,gBAAgB,CAAC;AACrD,OAAO,EAAE,aAAa,EAAE,MAAM,kBAAkB,CAAC;AAEjD,OAAO,EAAE,KAAK,EAAE,MAAM,wBAAwB,CAAC;AAC/C,OAAO,EAAE,yBAAyB,EAAE,MAAM,4CAA4C,CAAC;AAEvF;;;;;;;GAOG;AACH,MAAM,OAAO,kBAAmB,SAAQ,kBAAkB;IACxD,wCAAwC;IACxB,kBAAkB,GAAG,kBAAkB,CAAC;IAExD,yCAAyC;IAClC,QAAQ,CAAY;IAC3B,8CAA8C;IACvC,SAAS,CAAa;IAC7B,8CAA8C;IACtC,cAAc,CAAe;IACrC;;;OAGG;IACK,cAAc,CAAgB;IACtC,0EAA0E;IAC1E,IAAW,aAAa,KAAmB,OAAO,IAAI,CAAC,cAAc,KAAK,SAAS,CAAC,CAAC,CAAC,IAAI,CAAC,cAAc,CAAC,CAAC,CAAC,IAAI,CAAC,cAAc,CAAC,CAAC,CAAC;IAC1H,UAAU,CAAuB;IACzC,wCAAwC;IAChC,YAAY,CAAS;IAC7B,kDAAkD;IAC1C,YAAY,CAAY;IAChC,+BAA+B;IAC/B,8HAA8H;IAC9H,YAAoB,UAA8B,EAChD,SAA+B,EAC/B,QAAmB,EACnB,SAAqB,EACrB,sBAAiC,EACjC,YAAuB,EACvB,SAAiB,EACjB,UAAuD;QACvD,KAAK,CAAC,UAAU,EAAE,YAAY,EAAE,sBAAsB,EAAE,UAAU,CAAC,CAAC;QACpE,IAAI,CAAC,UAAU,GAAG,SAAS,CAAC;QAC5B,IAAI,CAAC,QAAQ,GAAG,QAAQ,CAAC;QACzB,IAAI,CAAC,SAAS,GAAG,SAAS,CAAC;QAC3B,IAAI,CAAC,YAAY,GAAG,SAAS,CAAC;QAC9B,IAAI,CAAC,cAAc,GAAG,YAAY,CAAC,MAAM,EAAE,CAAC;QAC5C,yFAAyF;QACzF,IAAI,CAAC,YAAY,GAAG,SAAS,CAAC,MAAM,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;QAC3C,IAAI,CAAC,yBAAyB,EAAE,CAAC;IACnC,CAAC;IACD,2CAA2C;IACpC,MAAM,CAAU,iBAAiB,GAAG,UAAU,CAAC;IAEtD,6EAA6E;IAEtE,8BAA8B,CAAC,QAAgB;QACpD,MAAM,YAAY,GAAG,IAAI,CAAC,UAAU,CAAC,cAAc,CAAC,QAAQ,CAAC,CAAC;QAC9D,MAAM,EAAE,GAAG,IAAI,CAAC,YAAY,CAAC;QAC7B,OAAO,IAAI,CAAC,SAAS,CAAC,YAAY,GAAG,YAAY,GAAG,EAAE,GAAG,IAAI,CAAC,YAAY,CAAC,WAAW,EAAE,GAAG,QAAQ,GAAG,IAAI,CAAC,YAAY,CAAC,EAAE,GAAG,EAAE,CAAC;IAClI,CAAC;IACD,6EAA6E;IACtE,yBAAyB,CAAC,QAAgB;QAC/C,MAAM,CAAC,GAAG,IAAI,CAAC,UAAU,CAAC,2BAA2B,CAAC,QAAQ,CAAC,CAAC;QAChE,OAAO,IAAI,CAAC,YAAY,CAAC,eAAe,CAAC,CAAC,CAAC,CAAC;IAC9C,CAAC;IAED,uEAAuE;IAEhE,wBAAwB,CAAC,cAAsB;QACpD,MAAM,QAAQ,GAAG,IAAI,CAAC,sBAAsB,CAAC,eAAe,CAAC,cAAc,CAAC,CAAC;QAC7E,OAAO,IAAI,CAAC,SAAS,CAAC,YAAY,GAAG,QAAQ,GAAG,IAAI,CAAC,YAAY,GAAG,CAAC,IAAI,CAAC,YAAY,CAAC,EAAE,GAAG,GAAG,GAAG,QAAQ,GAAG,CAAC,IAAI,CAAC,YAAY,CAAC,EAAE,GAAG,IAAI,CAAC,YAAY,CAAC,EAAE,CAAC,CAAC,CAAC;IAC9J,CAAC;IACD;;OAEG;IACa,mBAAmB,CAAC,cAAsB;QACxD,uBAAuB;QACvB,OAAO,IAAI,CAAC,YAAY,CAAC,eAAe,CAAC,IAAI,CAAC,sBAAsB,CAAC,eAAe,CAAC,cAAc,CAAC,CAAC,CAAC;IACxG,CAAC;IAED,2EAA2E;IAEnE,MAAM,CAAC,cAAc,CAAe;IACpC,MAAM,CAAC,YAAY,CAAe;IAClC,MAAM,CAAC,YAAY,CAAiF;IAC5G,0CAA0C;IACnC,MAAM,CAAC,aAAa;QACzB,kBAAkB,CAAC,cAAc,GAAG,IAAI,YAAY,CAAC,CAAC,CAAC,CAAC;QACxD,kBAAkB,CAAC,YAAY,GAAG,IAAI,YAAY,CAAC,CAAC,CAAC,CAAC;QACtD,kBAAkB,CAAC,YAAY,GAAG,CAAC,EAAE,EAAE,EAAE,EAAE,OAAO,EAAE,OAAO,EAAE,EAAE,CAAC,UAAU,CAAC,WAAW,CAAC,EAAE,EAAE,EAAE,EAAE,OAAO,EAAE,OAAO,CAAC,CAAC;IACnH,CAAC;IACD;;;;;;;;OAQG;IACK,6BAA6B,CAAC,GAAY,EAAE,SAAiB,EAAE,SAAiB,EAAE,WAAoB;QAC5G,MAAM,aAAa,GAAG,kBAAkB,CAAC,cAAc,CAAC;QACxD,MAAM,WAAW,GAAG,kBAAkB,CAAC,YAAY,CAAC;QACpD,MAAM,OAAO,GAAG,kBAAkB,CAAC,YAAY,CAAC,SAAS,EAAE,SAAS,EAAE,aAAa,EAAE,WAAW,CAAC,CAAC;QAClG,MAAM,MAAM,GAAG,IAAI,CAAC,YAAY,CAAC;QACjC,IAAI,CAAC,GAAG,CAAC,CAAC;QACV,IAAI,EAAE,GAAG,GAAG,CAAC;QACb,IAAI,EAAE,GAAG,GAAG,CAAC;QACb,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,OAAO,EAAE,CAAC,EAAE,EAAE,CAAC;YACjC,MAAM,OAAO,GAAG,IAAI,CAAC,8BAA8B,CAAC,aAAa,CAAC,CAAC,CAAC,CAAC,CAAC;YACtE,CAAC,GAAG,WAAW,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC;YAC5B,EAAE,IAAI,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,OAAO,CAAC,CAAC;YAC5B,EAAE,IAAI,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,OAAO,CAAC,CAAC;QAC9B,CAAC;QACD,IAAI,WAAW;YACb,QAAQ,CAAC,qBAAqB,CAAC,GAAG,EAAE,IAAI,CAAC,YAAY,CAAC,MAAM,EAAE,EAAE,CAAC,EAAE,EAAE,EAAE,CAAC,EAAE,EAAE,EAAE,CAAC,EAAE,GAAG,EAAE,EAAE,GAAG,CAAC,CAAC;;YAE7F,GAAG,CAAC,aAAa,CAAC,EAAE,EAAE,EAAE,EAAE,GAAG,CAAC,CAAC;IAEnC,CAAC;IACD,wBAAwB;IACR,yBAAyB;QACvC,IAAI,CAAC,YAAY,GAAG,SAAS,CAAC,MAAM,CAClC,kBAAkB,CAAC,iBAAiB,CAAC,IAAI,CAAC,QAAQ,CAAC,EAAE,CAAC,EACtD,kBAAkB,CAAC,iBAAiB,CAAC,IAAI,CAAC,QAAQ,CAAC,EAAE,CAAC,CAAC,CAAC;QAC1D,IAAI,CAAC,cAAc,CAAC,KAAK,EAAE,CAAC;QAC5B,MAAM,YAAY,GAAG,OAAO,CAAC,MAAM,EAAE,CAAC;QACtC,IAAI,CAAC,cAAc,CAAC,iBAAiB,CAAC,YAAY,CAAC,CAAC;QACpD,MAAM,WAAW,GAAG,EAAE,CAAC;QAEvB,MAAM,YAAY,GAAG,GAAG,GAAG,WAAW,CAAC;QACvC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,WAAW,EAAE,CAAC,EAAE,EAAE,CAAC;YACtC,MAAM,SAAS,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,YAAY,CAAC;YACzC,MAAM,SAAS,GAAG,CAAC,GAAG,YAAY,CAAC;YACnC,IAAI,CAAC,6BAA6B,CAAC,YAAY,EAAE,SAAS,EAAE,SAAS,EAAE,KAAK,CAAC,CAAC;YAC9E,IAAI,CAAC,cAAc,CAAC,iBAAiB,CAAC,YAAY,CAAC,CAAC;QACtD,CAAC;QACD,IAAI,CAAC,cAAc,CAAC,mBAAmB,CAAC,IAAI,CAAC,YAAY,CAAC,CAAC;QAC3D,IAAI,CAAC,IAAI,CAAC,sBAAsB,CAAC,SAAS,EAAE,CAAC;YAC3C,IAAI,IAAI,CAAC,cAAc,KAAK,SAAS;gBACnC,IAAI,CAAC,cAAc,GAAG,YAAY,CAAC,MAAM,EAAE,CAAC;YAC9C,IAAI,CAAC,cAAc,CAAC,KAAK,EAAE,CAAC;YAC5B,0EAA0E;YAC1E,mKAAmK;YACnK,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,WAAW,EAAE,CAAC,EAAE,EAAE,CAAC;gBACtC,MAAM,aAAa,GAAG,CAAC,GAAG,YAAY,CAAC;gBACvC,IAAI,CAAC,cAAc,CAAC,QAAQ,CAAC,IAAI,CAAC,eAAe,CAAC,aAAa,CAAC,CAAC,CAAC;YACpE,CAAC;QACH,CAAC;IACH,CAAC;IACD;;;;;;OAMG;IACI,MAAM,CAAC,gCAAgC,CAAC,QAAmB,EAAE,SAAqB,EAAE,sBAAiC,EAAE,YAAuB,EAAE,QAAiB;QACtK,MAAM,SAAS,GAAG,kBAAkB,CAAC,mCAAmC,CAAC,QAAQ,CAAC,EAAE,EAAE,QAAQ,CAAC,EAAE,EAAE,SAAS,CAAC,YAAY,CAAC,CAAC;QAC3H,IAAI,QAAQ,KAAK,SAAS;YACxB,QAAQ,GAAG,UAAU,CAAC;QACxB,MAAM,SAAS,GAAG,oBAAoB,CAAC,aAAa,CAAC,QAAQ,CAAC,CAAC;QAC/D,IAAI,CAAC,SAAS;YACZ,OAAO,SAAS,CAAC;QACnB,OAAO,IAAI,kBAAkB,CAAC,QAAQ,EAAE,SAAS,EAC/C,QAAQ,CAAC,KAAK,EAAE,EAChB,SAAS,CAAC,KAAK,EAAE,EAAE,sBAAsB,CAAC,KAAK,EAAE,EAAE,YAAY,CAAC,KAAK,EAAE,EAAE,SAAS,EAClF,IAAI,+BAA+B,CAAC,QAAQ,CAAC,EAAE,EAAE,QAAQ,CAAC,EAAE,EAC1D,SAAS,CAAC,UAAU,CAAC,KAAK,EAAE,EAAE,SAAS,CAAC,QAAQ,CAAC,KAAK,EAAE,EACxD,SAAS,CAAC,CAAC,CAAC;IAClB,CAAC;IACD;;;;;;;;;;OAUG;IACI,MAAM,CAAC,iBAAiB,CAC7B,UAA8B,EAC9B,OAA2B,EAC3B,OAA2B,EAC3B,QAA2B,EAC3B,QAA2B,EAC3B,SAA6B,EAC7B,gBAAuC,EACvC,YAAuB;QACvB,IAAI,UAAU,KAAK,SAAS;YAC1B,UAAU,GAAG,UAAU,CAAC;QAC1B,MAAM,SAAS,GAAG,oBAAoB,CAAC,aAAa,CAAC,UAAU,CAAC,CAAC;QACjE,IAAI,CAAC,SAAS;YACZ,OAAO,SAAS,CAAC;QACnB,MAAM,IAAI,GAAG,IAAI,+BAA+B,CAAC,OAAO,EAAE,OAAO,EAAE,QAAQ,EAAE,QAAQ,EAAE,SAAS,CAAC,CAAC;QAClG,MAAM,KAAK,GAAG,IAAI,CAAC,KAAK,EAAE,CAAC;QAC3B,IAAI,CAAC,IAAI,CAAC,0BAA0B,EAAE;YACpC,OAAO,SAAS,CAAC;QACnB,IAAI,gBAAgB,KAAK,SAAS;YAChC,gBAAgB,GAAG,SAAS,CAAC,MAAM,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;QAC5C,OAAO,IAAI,kBAAkB,CAC3B,UAAU,EACV,SAAS,EACT,SAAS,CAAC,MAAM,CAAC,IAAI,CAAC,OAAO,EAAE,IAAI,CAAC,OAAO,CAAC,EAC5C,UAAU,CAAC,cAAc,CAAC,IAAI,CAAC,QAAS,EAAE,IAAI,CAAC,QAAS,CAAC,EACzD,gBAAgB,CAAC,CAAC,CAAC,gBAAgB,CAAC,KAAK,EAAE,CAAC,CAAC,CAAC,SAAS,CAAC,MAAM,CAAC,CAAC,EAAE,CAAC,CAAC,EACpE,YAAY,EAAE,IAAI,CAAC,WAAY,EAAE,KAAK,CAAC,CAAC;IAC5C,CAAC;IACD,kDAAkD;IAC3C,OAAO,CAAC,KAAyB;QACtC,IAAI,CAAC,YAAY,CAAC,OAAO,CAAC,KAAK,CAAC,YAAY,CAAC,CAAC;QAC9C,IAAI,CAAC,QAAQ,CAAC,OAAO,CAAC,KAAK,CAAC,QAAQ,CAAC,CAAC;QACtC,IAAI,CAAC,YAAY,CAAC,OAAO,CAAC,KAAK,CAAC,YAAY,CAAC,CAAC;QAC9C,IAAI,CAAC,SAAS,CAAC,OAAO,CAAC,KAAK,CAAC,SAAS,CAAC,CAAC;QACxC,IAAI,CAAC,YAAY,CAAC,OAAO,CAAC,KAAK,CAAC,YAAY,CAAC,CAAC;QAC9C,IAAI,CAAC,sBAAsB,CAAC,OAAO,CAAC,KAAK,CAAC,sBAAsB,CAAC,CAAC;QAClE,IAAI,CAAC,YAAY,GAAG,KAAK,CAAC,YAAY,CAAC;QACvC,OAAO,IAAI,CAAC;IACd,CAAC;IACD,gCAAgC;IAChB,KAAK;QACnB,OAAO,IAAI,kBAAkB,CAAC,IAAI,CAAC,WAAW,EAAE,IAAI,CAAC,UAAU,EAC7D,IAAI,CAAC,QAAQ,CAAC,KAAK,EAAE,EAAE,IAAI,CAAC,SAAS,CAAC,KAAK,EAAE,EAC7C,IAAI,CAAC,sBAAsB,CAAC,KAAK,EAAE,EAAE,IAAI,CAAC,YAAY,CAAC,KAAK,EAAE,EAAE,IAAI,CAAC,YAAY,EACjF,IAAI,CAAC,iBAAiB,EAAE,KAAK,EAAE,CAAC,CAAC;IACrC,CAAC;IAED,mEAAmE;IAC5D,mBAAmB,CAAC,UAAqB;QAE9C,MAAM,SAAS,GAAG,IAAI,CAAC,yBAAyB,CAAC,UAAU,CAAC,CAAC;QAC7D,IAAI,SAAS,KAAK,SAAS,EAAE,CAAC;YAC5B,IAAI,CAAC,YAAY,CAAC,EAAE,IAAI,SAAS,CAAC,KAAK,CAAC;YACxC,IAAI,CAAC,YAAY,CAAC,EAAE,IAAI,SAAS,CAAC,KAAK,CAAC;YACxC,IAAI,CAAC,QAAQ,CAAC,EAAE,IAAI,SAAS,CAAC,KAAK,CAAC;YACpC,IAAI,CAAC,QAAQ,CAAC,EAAE,IAAI,SAAS,CAAC,KAAK,CAAC;YACpC,IAAI,CAAC,YAAY,IAAI,SAAS,CAAC,KAAK,CAAC;QACvC,CAAC;QACD,IAAI,CAAC,yBAAyB,EAAE,CAAC;QACjC,OAAO,IAAI,CAAC;IACd,CAAC;IAED,qCAAqC;IACrB,UAAU,CAAC,MAAgB;QACzC,OAAO,IAAI,CAAC,aAAa,CAAC,UAAU,CAAC,MAAM,CAAC,CAAC;IAC/C,CAAC;IACD,mCAAmC;IACnB,QAAQ,CAAC,MAAgB;QACvC,OAAO,IAAI,CAAC,aAAa,CAAC,QAAQ,CAAC,MAAM,CAAC,CAAC;IAC7C,CAAC;IACD,mFAAmF;IAC5E,SAAS,CAAC,KAAmC;QAClD,OAAO,KAAK,CAAC,cAAc,CAAC,IAAI,CAAC,YAAY,CAAC,MAAiB,CAAC;eAC3D,QAAQ,CAAC,gBAAgB,CAAC,GAAG,EAAE,IAAI,CAAC,YAAY,CAAC,MAAM,CAAC,UAAU,CAAC,KAAK,CAAC,YAAY,EAAE,CAAC,CAAC;eACzF,QAAQ,CAAC,gBAAgB,CAAC,GAAG,EAAE,IAAI,CAAC,YAAY,CAAC,MAAM,CAAC,UAAU,CAAC,KAAK,CAAC,YAAY,EAAE,CAAC,CAAC,CAAC;IACjG,CAAC;IACD,kJAAkJ;IAC3I,WAAW,KAAK,OAAO,IAAI,CAAC,WAAW,EAAE,CAAC,CAAC,CAAC;IACnD,kJAAkJ;IAClI,WAAW,KAAK,OAAO,IAAI,CAAC,YAAY,GAAG,CAAC,IAAI,CAAC,uBAAuB,CAAC,aAAa,EAAE,CAAC,CAAC,CAAC,CAAC;IAC5G,+KAA+K;IAC/J,2BAA2B,CAAC,SAAiB,EAAE,SAAiB;QAC9E,OAAO,IAAI,CAAC,YAAY,GAAG,CAAC,IAAI,CAAC,uBAAuB,CAAC,aAAa,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,SAAS,GAAG,SAAS,CAAC,CAAC,CAAC;IAC9G,CAAC;IACD,6DAA6D;IACtD,mBAAmB,CAAC,KAAU,IAAa,OAAO,KAAK,YAAY,kBAAkB,CAAC,CAAC,CAAC;IAC/F;;;OAGG;IACI,WAAW,CAAC,IAAkB,EAAE,OAAuB,IAAU,IAAI,CAAC,aAAa,CAAC,WAAW,CAAC,IAAI,EAAE,OAAO,CAAC,CAAC,CAAC,CAAC;IACxH,+CAA+C;IACxC,kBAAkB,CAAC,IAAoB,EAAE,OAAuB;QACrE,MAAM,CAAC,GAAG,IAAI,CAAC,4BAA4B,CAAC,OAAO,CAAC,CAAC;QACrD,IAAI,CAAC,yBAAyB,CAAC,IAAI,CAAC,CAAC;QACrC,MAAM,aAAa,GAAG,IAAI,CAAC,aAAa,CAAC;QACzC,MAAM,aAAa,GAAG,IAAI,CAAC,6BAA6B,KAAK,SAAS,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,IAAI,CAAC,6BAA6B,EAAE,CAAC;QACtH,IAAI,CAAC,aAAa,IAAI,CAAC,IAAI,aAAa,CAAC,SAAS,EAAE,EAAE,CAAC;YACrD,IAAI,CAAC,aAAa,CAAC,kBAAkB,CAAC,IAAI,EAAE,OAAO,CAAC,CAAC;QACvD,CAAC;aAAM,CAAC;YACN,IAAI,CAAC,qCAAqC,CAAC,IAAI,EAAE,CAAC,EAAE,GAAG,EAAE,GAAG,CAAC,CAAC;QAChE,CAAC;QACD,IAAI,CAAC,uBAAuB,CAAC,IAAI,CAAC,CAAC;IACrC,CAAC;IAED;;;OAGG;IAEI,4BAA4B,CAAC,OAAuB;QACzD,IAAI,SAAS,CAAC;QACd,IAAI,OAAO,EAAE,CAAC;YACZ,MAAM,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,QAAQ,CAAC,EAAE,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,QAAQ,CAAC,EAAE,CAAC,CAAC,CAAC;YAC9E,SAAS,GAAG,OAAO,CAAC,oBAAoB,CAAC,IAAI,EAAE,IAAI,CAAC,SAAS,CAAC,YAAY,CAAC,CAAC;YAC5E,SAAS,GAAG,OAAO,CAAC,kBAAkB,CAAC,SAAS,EAAE,IAAI,CAAC,WAAW,EAAE,CAAC,CAAC;YACtE,SAAS,GAAG,OAAO,CAAC,2BAA2B,CAAC,SAAS,CAAC,CAAC;QAC7D,CAAC;aAAM,CAAC;YACN,SAAS,GAAG,aAAa,CAAC,aAAa,CAAC,SAAS,EAAE,CAAC,EAAE,IAAI,CAAC,SAAS,CAAC,YAAY,CAAC,CAAC;QACrF,CAAC;QACD,OAAO,SAAS,CAAC;IACnB,CAAC;IAED;;OAEG;IAEI,cAAc;QACnB,IAAI,CAAC,sBAAsB,CAAC,cAAc,EAAE,CAAC;QAC7C,IAAI,IAAI,CAAC,cAAc,KAAK,SAAS;YACnC,IAAI,CAAC,cAAc,GAAG,IAAI,CAAC,cAAc,CAAC,KAAK,EAAE,CAAC;QACpD,IAAI,CAAC,cAAc,CAAC,cAAc,EAAE,CAAC;IACvC,CAAC;IACD,qDAAqD;IAC9C,eAAe,CAAC,cAAsB,EAAE,MAAgB;QAC7D,MAAM,oBAAoB,GAAG,IAAI,CAAC,sBAAsB,CAAC,eAAe,CAAC,cAAc,CAAC,CAAC;QACzF,MAAM,UAAU,GAAG,IAAI,CAAC,cAAc,CAAC,YAAY,CAAC,MAAM,GAAG,CAAC,CAAC;QAC/D,IAAI,cAAc,GAAG,GAAG,EAAE,CAAC;YACzB,MAAM,GAAG,IAAI,CAAC,cAAc,CAAC,YAAY,CAAC,IAAI,CAAC,MAAM,CAAE,CAAC;YACxD,MAAM,eAAe,GAAG,GAAG,GAAG,UAAU,CAAC;YACzC,IAAI,qBAAqB,GAAG,GAAG,CAAC;YAChC,IAAI,kBAAkB,GAAG,qBAAqB,GAAG,eAAe,CAAC;YACjE,OAAO,kBAAkB,GAAG,oBAAoB,EAAE,CAAC;gBACjD,IAAI,CAAC,6BAA6B,CAAC,MAAM,EAAE,qBAAqB,EAAE,kBAAkB,EAAE,IAAI,CAAC,CAAC;gBAC5F,qBAAqB,GAAG,kBAAkB,CAAC;gBAC3C,kBAAkB,IAAI,eAAe,CAAC;YACxC,CAAC;YACD,IAAI,CAAC,6BAA6B,CAAC,MAAM,EAAE,qBAAqB,EAAE,oBAAoB,EAAE,IAAI,CAAC,CAAC;QAChG,CAAC;aAAM,IAAI,cAAc,GAAG,GAAG,EAAE,CAAC;YAChC,MAAM,GAAG,IAAI,CAAC,cAAc,CAAC,YAAY,CAAC,KAAK,CAAC,MAAM,CAAE,CAAC;YACzD,MAAM,eAAe,GAAG,GAAG,GAAG,UAAU,CAAC;YACzC,IAAI,qBAAqB,GAAG,GAAG,CAAC;YAChC,IAAI,kBAAkB,GAAG,qBAAqB,GAAG,eAAe,CAAC;YACjE,OAAO,kBAAkB,GAAG,oBAAoB,EAAE,CAAC;gBACjD,IAAI,CAAC,6BAA6B,CAAC,MAAM,EAAE,qBAAqB,EAAE,kBAAkB,EAAE,IAAI,CAAC,CAAC;gBAC5F,qBAAqB,GAAG,kBAAkB,CAAC;gBAC3C,kBAAkB,IAAI,eAAe,CAAC;YACxC,CAAC;YACD,IAAI,CAAC,6BAA6B,CAAC,MAAM,EAAE,qBAAqB,EAAE,oBAAoB,EAAE,IAAI,CAAC,CAAC;QAChG,CAAC;aAAM,CAAC;YACN,MAAM,qBAAqB,GAAG,QAAQ,CAAC,eAAe,CAAC,oBAAoB,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;YACnF,MAAM,MAAM,GAAG,IAAI,CAAC,KAAK,CAAC,qBAAqB,GAAG,UAAU,CAAC,CAAC,CAAC,mDAAmD;YAClH,MAAM,eAAe,GAAG,MAAM,GAAG,UAAU,CAAC;YAC5C,MAAM,GAAG,IAAI,CAAC,cAAc,CAAC,YAAY,CAAC,+BAA+B,CAAC,MAAM,EAAE,MAAM,CAAC,CAAC;YAC1F,sLAAsL;YACtL,IAAI,CAAC,6BAA6B,CAAC,MAAM,EAAE,eAAe,EAAE,oBAAoB,EAAE,IAAI,CAAC,CAAC;QAE1F,CAAC;QACD,OAAO,MAAM,CAAC;IAChB,CAAC;IACD,oEAAoE;IAC7D,4BAA4B,CAAC,cAAsB,EAAE,MAAc;QACxE,MAAM,cAAc,GAAG,IAAI,CAAC,sBAAsB,CAAC,eAAe,CAAC,cAAc,CAAC,CAAC;QACnF,MAAM,GAAG,MAAM,CAAC,CAAC,CAAC,MAAM,CAAC,CAAC,CAAC,KAAK,CAAC,UAAU,EAAE,CAAC;QAC9C,IAAI,CAAC,eAAe,CAAC,cAAc,EAAE,MAAM,CAAC,MAAM,CAAC,CAAC;QACpD,MAAM,OAAO,GAAG,IAAI,CAAC,8BAA8B,CAAC,cAAc,CAAC,CAAC;QACpE,MAAM,CAAC,GAAG,IAAI,CAAC,YAAY,GAAG,IAAI,CAAC,sBAAsB,CAAC,WAAW,EAAE,CAAC;QACxE,IAAI,CAAC,YAAY,CAAC,MAAM,CAAC,UAAU,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,OAAO,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,OAAO,CAAC,EAAE,MAAM,CAAC,SAAS,CAAC,CAAC;QACpG,OAAO,MAAM,CAAC;IAChB,CAAC;IAED,sDAAsD;IACtC,qBAAqB,CAAC,cAAsB,EAAE,MAAkB;QAC9E,MAAM,cAAc,GAAG,IAAI,CAAC,sBAAsB,CAAC,eAAe,CAAC,cAAc,CAAC,CAAC;QACnF,MAAM,GAAG,MAAM,CAAC,CAAC,CAAC,MAAM,CAAC,CAAC,CAAC,SAAS,CAAC,cAAc,EAAE,CAAC;QACtD,MAAM,CAAC,MAAM,CAAC,OAAO,CAAC,IAAI,CAAC,eAAe,CAAC,cAAc,CAAC,CAAC,CAAC;QAC5D,QAAQ,CAAC,uBAAuB,CAAC,IAAI,CAAC,YAAY,CAAC,MAAM,EAAE,SAAS,CAAC,GAAG,EAAE,MAAM,CAAC,MAAM,CAAC,CAAC;QAEzF,MAAM,OAAO,GAAG,IAAI,CAAC,8BAA8B,CAAC,cAAc,CAAC,CAAC;QACpE,MAAM,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,OAAO,CAAC,CAAC;QAC5B,MAAM,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,OAAO,CAAC,CAAC;QAC5B,MAAM,CAAC,MAAM,CAAC,mBAAmB,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;QAC9C,OAAO,MAAM,CAAC;IAChB,CAAC;IACD;;;;;OAKG;IACI,8BAA8B,CAAC,cAAsB,EAAE,MAAkC;QAC9F,MAAM,cAAc,GAAG,IAAI,CAAC,sBAAsB,CAAC,eAAe,CAAC,cAAc,CAAC,CAAC;QACnF,MAAM,MAAM,GAAG,IAAI,CAAC,eAAe,CAAC,cAAc,CAAC,CAAC;QACpD,MAAM,OAAO,GAAG,IAAI,CAAC,8BAA8B,CAAC,cAAc,CAAC,CAAC;QACpE,MAAM,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,OAAO,CAAC,CAAC;QAC5B,MAAM,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,OAAO,CAAC,CAAC;QAC5B,MAAM,KAAK,GAAG,IAAI,CAAC,sBAAsB,CAAC,WAAW,EAAE,CAAC;QACxD,MAAM,CAAC,GAAG,KAAK,CAAC;QAChB,MAAM,CAAC,GAAG,CAAC,GAAG,KAAK,CAAC;QACpB,MAAM,OAAO,GAAG,IAAI,CAAC,YAAY,CAAC,MAAM,CAAC,UAAU,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC;QAClE,MAAM,OAAO,GAAG,IAAI,CAAC,YAAY,CAAC,MAAM,CAAC,UAAU,CAAC,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC;QACnE,OAAO,CAAC,YAAY,CAAC,IAAI,CAAC,yBAAyB,CAAC,cAAc,CAAC,CAAC,CAAC;QACrE,OAAO,yBAAyB,CAAC,aAAa,CAAC,MAAM,EAAE,OAAO,EAAE,OAAO,EAAE,MAAM,CAAC,CAAC;IACnF,CAAC;IACD,mFAAmF;IAC5E,yBAAyB,CAAC,OAAwB;QACvD,OAAO,OAAO,CAAC,sBAAsB,CAAC,IAAI,CAAC,CAAC;IAC9C,CAAC;IAED,4CAA4C;IAC5B,aAAa,CAAC,KAAqB;QACjD,IAAI,KAAK,YAAY,kBAAkB,EAAE,CAAC;YACxC,OAAO,IAAI,CAAC,QAAQ,CAAC,aAAa,CAAC,KAAK,CAAC,QAAQ,CAAC;mBAC7C,IAAI,CAAC,SAAS,CAAC,6BAA6B,CAAC,KAAK,CAAC,SAAS,CAAC;mBAC7D,IAAI,CAAC,YAAY,CAAC,aAAa,CAAC,KAAK,CAAC,YAAY,CAAC;mBACnD,QAAQ,CAAC,gBAAgB,CAAC,IAAI,CAAC,YAAY,EAAE,KAAK,CAAC,YAAY,CAAC;mBAChE,IAAI,CAAC,sBAAsB,CAAC,aAAa,CAAC,KAAK,CAAC,sBAAsB,CAAC;mBACvE,IAAI,CAAC,YAAY,CAAC,aAAa,CAAC,KAAK,CAAC,YAAY,CAAC,CAAC;QAC3D,CAAC;QACD,OAAO,KAAK,CAAC;IACf,CAAC;;AAEH,sDAAsD;AACtD,kBAAkB,CAAC,aAAa,EAAE,CAAC","sourcesContent":["/*---------------------------------------------------------------------------------------------\n* Copyright (c) Bentley Systems, Incorporated. All rights reserved.\n* See LICENSE.md in the project root for license terms and full copyright notice.\n*--------------------------------------------------------------------------------------------*/\n/** @packageDocumentation\n * @module Curve\n */\n\nimport { TransitionSpiral3d } from \"./TransitionSpiral3d\";\nimport { Segment1d } from \"../../geometry3d/Segment1d\";\nimport { AngleSweep } from \"../../geometry3d/AngleSweep\";\nimport { Transform } from \"../../geometry3d/Transform\";\nimport { LineString3d } from \"../LineString3d\";\nimport { NormalizedTransition } from \"./NormalizedTransition\";\nimport { TransitionConditionalProperties } from \"./TransitionConditionalProperties\";\nimport { Quadrature } from \"../../numerics/Quadrature\";\nimport { Point3d } from \"../../geometry3d/Point3dVector3d\";\nimport { Matrix3d } from \"../../geometry3d/Matrix3d\";\nimport { Angle } from \"../../geometry3d/Angle\";\nimport { Plane3dByOriginAndUnitNormal } from \"../../geometry3d/Plane3dByOriginAndUnitNormal\";\nimport { AxisOrder, Geometry } from \"../../Geometry\";\nimport { StrokeOptions } from \"../StrokeOptions\";\nimport { GeometryHandler, IStrokeHandler } from \"../../geometry3d/GeometryHandler\";\nimport { Ray3d } from \"../../geometry3d/Ray3d\";\nimport { Plane3dByOriginAndVectors } from \"../../geometry3d/Plane3dByOriginAndVectors\";\nimport { GeometryQuery } from \"../GeometryQuery\";\n/**\n * An IntegratedSpiral3d is a curve defined by integrating its curvature.\n * * The first integral of curvature (with respect to distance along the curve) is the bearing angle (in radians)\n * * Integrating (cos(theta), sin(theta)) gives displacement from the start point, and thus the actual curve position.\n * * The curvature functions of interest are all symmetric snap functions in the NormalizedTransition class.\n * * `TransitionConditionalProperties` implements the computations of the interrelationship of radii, bearing, and length.\n * @public\n */\nexport class IntegratedSpiral3d extends TransitionSpiral3d {\n  /** String name for schema properties */\n  public readonly curvePrimitiveType = \"transitionSpiral\";\n\n  /** start and end radii as a Segment1d */\n  public radius01: Segment1d;\n  /** start and end bearings as an AngleSweep */\n  public bearing01: AngleSweep;\n  /** stroked approximation of entire spiral. */\n  private _globalStrokes: LineString3d;\n  /** stroked approximation of active spiral.\n   * * Same count as global -- possibly overly fine, but it gives some consistency between same clothoid constructed as partial versus complete.\n   * * If no trimming, this points to the same place as the _globalStrokes !!!  Don't double transform!!!\n   */\n  private _activeStrokes?: LineString3d;\n  /** Return the internal stroked form of the (possibly partial) spiral   */\n  public get activeStrokes(): LineString3d { return this._activeStrokes !== undefined ? this._activeStrokes : this._globalStrokes; }\n  private _evaluator: NormalizedTransition;\n  /** Total curve arc length (computed) */\n  private _arcLength01: number;\n  /** Curvatures (inverse radii) at start and end */\n  private _curvature01: Segment1d;\n  /** evaluator for transition */\n  // constructor demands all bearing, radius, and length data -- caller determines usual dependency of \"any 4 determine the 5th\"\n  private constructor(spiralType: string | undefined,\n    evaluator: NormalizedTransition,\n    radius01: Segment1d,\n    bearing01: AngleSweep,\n    activeFractionInterval: Segment1d,\n    localToWorld: Transform,\n    arcLength: number,\n    properties: TransitionConditionalProperties | undefined) {\n    super(spiralType, localToWorld, activeFractionInterval, properties);\n    this._evaluator = evaluator;\n    this.radius01 = radius01;\n    this.bearing01 = bearing01;\n    this._arcLength01 = arcLength;\n    this._globalStrokes = LineString3d.create();\n    // initialize for compiler -- but this will be recomputed in refreshComputeProperties ...\n    this._curvature01 = Segment1d.create(0, 1);\n    this.refreshComputedProperties();\n  }\n  /** default spiral type name. (clothoid) */\n  public static readonly defaultSpiralType = \"clothoid\";\n\n  /** use the integrated function to return an angle at fractional position. */\n\n  public globalFractionToBearingRadians(fraction: number): number {\n    const areaFraction = this._evaluator.fractionToArea(fraction);\n    const dx = this._arcLength01;\n    return this.bearing01.startRadians + areaFraction * dx * this._curvature01.signedDelta() + fraction * this._curvature01.x0 * dx;\n  }\n  /** use the integrated function to return an angle at fractional position. */\n  public globalFractionToCurvature(fraction: number): number {\n    const f = this._evaluator.fractionToCurvatureFraction(fraction);\n    return this._curvature01.fractionToPoint(f);\n  }\n\n  /** Return the bearing at given fraction of the active interval .... */\n\n  public fractionToBearingRadians(activeFraction: number): number {\n    const fraction = this.activeFractionInterval.fractionToPoint(activeFraction);\n    return this.bearing01.startRadians + fraction * this._arcLength01 * (this._curvature01.x0 + 0.5 * fraction * (this._curvature01.x1 - this._curvature01.x0));\n  }\n  /** Return the curvature at given fraction of the active interval ...\n   * * The `undefined` result is to match the abstract class -- it cannot actually occur.\n   */\n  public override fractionToCurvature(activeFraction: number): number | undefined {\n    // BUG? active interval\n    return this._curvature01.fractionToPoint(this.activeFractionInterval.fractionToPoint(activeFraction));\n  }\n\n  // These static variables are reused on calls to integrateFromStartFraction\n\n  private static _gaussFraction: Float64Array;\n  private static _gaussWeight: Float64Array;\n  private static _gaussMapper: (xA: number, xB: number, arrayX: Float64Array, arrayW: Float64Array) => number;\n  /** Initialize class level work arrays. */\n  public static initWorkSpace() {\n    IntegratedSpiral3d._gaussFraction = new Float64Array(5);\n    IntegratedSpiral3d._gaussWeight = new Float64Array(5);\n    IntegratedSpiral3d._gaussMapper = (xA, xB, xMapped, wMapped) => Quadrature.setupGauss5(xA, xB, xMapped, wMapped);\n  }\n  /** Evaluate and sum the gauss quadrature formulas to integrate cos(theta), sin(theta) fractional subset of a reference length.\n   * (recall that theta is a nonlinear function of the fraction.)\n   * * This is a single interval of gaussian integration.\n   * * The fraction is on the full spiral (not in the mapped active interval)\n   * @param xyz advancing integrated point.\n   * @param fractionA fraction at start of interval\n   * @param fractionB fraction at end of interval.\n   * @param unitArcLength length of curve for 0 to 1 fractional\n   */\n  private fullSpiralIncrementalIntegral(xyz: Point3d, fractionA: number, fractionB: number, applyMatrix: boolean) {\n    const gaussFraction = IntegratedSpiral3d._gaussFraction;\n    const gaussWeight = IntegratedSpiral3d._gaussWeight;\n    const numEval = IntegratedSpiral3d._gaussMapper(fractionA, fractionB, gaussFraction, gaussWeight);\n    const deltaL = this._arcLength01;\n    let w = 0;\n    let dx = 0.0;\n    let dy = 0.0;\n    for (let k = 0; k < numEval; k++) {\n      const radians = this.globalFractionToBearingRadians(gaussFraction[k]);\n      w = gaussWeight[k] * deltaL;\n      dx += w * Math.cos(radians);\n      dy += w * Math.sin(radians);\n    }\n    if (applyMatrix)\n      Matrix3d.xyzPlusMatrixTimesXYZ(xyz, this.localToWorld.matrix, { x: dx, y: dy, z: 0.0 }, xyz);\n    else\n      xyz.addXYZInPlace(dx, dy, 0.0);\n\n  }\n  /** Recompute strokes */\n  public override refreshComputedProperties() {\n    this._curvature01 = Segment1d.create(\n      TransitionSpiral3d.radiusToCurvature(this.radius01.x0),\n      TransitionSpiral3d.radiusToCurvature(this.radius01.x1));\n    this._globalStrokes.clear();\n    const currentPoint = Point3d.create();\n    this._globalStrokes.appendStrokePoint(currentPoint);\n    const numInterval = 16;\n\n    const fractionStep = 1.0 / numInterval;\n    for (let i = 1; i <= numInterval; i++) {\n      const fraction0 = (i - 1) * fractionStep;\n      const fraction1 = i * fractionStep;\n      this.fullSpiralIncrementalIntegral(currentPoint, fraction0, fraction1, false);\n      this._globalStrokes.appendStrokePoint(currentPoint);\n    }\n    this._globalStrokes.tryTransformInPlace(this.localToWorld);\n    if (!this.activeFractionInterval.isExact01) {\n      if (this._activeStrokes === undefined)\n        this._activeStrokes = LineString3d.create();\n      this._activeStrokes.clear();\n      // finer strokes in the active interval ... same fraction step, but mapped\n      // This assumes factionToPoint acts normally within refreshComputedProperties -- that depends on the global strokes we just computed, but not on the active strokes\n      for (let i = 0; i <= numInterval; i++) {\n        const localFraction = i * fractionStep;\n        this._activeStrokes.addPoint(this.fractionToPoint(localFraction));\n      }\n    }\n  }\n  /**\n   * Create a transition spiral with radius and bearing conditions.\n   * @param radius01 radius (inverse curvature) at start and end. (radius of zero means straight line)\n   * @param bearing01 bearing angles at start and end.  bearings are measured from the x axis, positive clockwise towards y axis\n   * @param activeFractionInterval fractional limits of the active portion of the spiral.\n   * @param localToWorld placement frame.  Fractional coordinate 0 is at the origin.\n   */\n  public static createRadiusRadiusBearingBearing(radius01: Segment1d, bearing01: AngleSweep, activeFractionInterval: Segment1d, localToWorld: Transform, typeName?: string) {\n    const arcLength = TransitionSpiral3d.radiusRadiusSweepRadiansToArcLength(radius01.x0, radius01.x1, bearing01.sweepRadians);\n    if (typeName === undefined)\n      typeName = \"clothoid\";\n    const evaluator = NormalizedTransition.findEvaluator(typeName);\n    if (!evaluator)\n      return undefined;\n    return new IntegratedSpiral3d(typeName, evaluator,\n      radius01.clone(),\n      bearing01.clone(), activeFractionInterval.clone(), localToWorld.clone(), arcLength,\n      new TransitionConditionalProperties(radius01.x0, radius01.x1,\n        bearing01.startAngle.clone(), bearing01.endAngle.clone(),\n        undefined));\n  }\n  /**\n   * Create a transition spiral.\n   * * Inputs must provide exactly 4 of the 5 values `[radius0,radius1,bearing0,bearing1,length`.\n   * @param spiralType one of \"clothoid\", \"bloss\", \"biquadratic\", \"cosine\", \"sine\".  If undefined, \"clothoid\" is used.\n   * @param radius0 radius (or 0 for tangent to line) at start\n   * @param radius1 radius (or 0 for tangent to line) at end\n   * @param bearing0 bearing, measured CCW from x axis at start.\n   * @param bearing1 bearing, measured CCW from x axis at end.\n   * @param fractionInterval optional fractional interval for an \"active\" portion of the curve.   if omitted, the full [0,1] is used.\n   * @param localToWorld placement transform\n   */\n  public static createFrom4OutOf5(\n    spiralType: string | undefined,\n    radius0: number | undefined,\n    radius1: number | undefined,\n    bearing0: Angle | undefined,\n    bearing1: Angle | undefined,\n    arcLength: number | undefined,\n    fractionInterval: undefined | Segment1d,\n    localToWorld: Transform): IntegratedSpiral3d | undefined {\n    if (spiralType === undefined)\n      spiralType = \"clothoid\";\n    const evaluator = NormalizedTransition.findEvaluator(spiralType);\n    if (!evaluator)\n      return undefined;\n    const data = new TransitionConditionalProperties(radius0, radius1, bearing0, bearing1, arcLength);\n    const data1 = data.clone();\n    if (!data.tryResolveAnySingleUnknown())\n      return undefined;\n    if (fractionInterval === undefined)\n      fractionInterval = Segment1d.create(0, 1);\n    return new IntegratedSpiral3d(\n      spiralType,\n      evaluator,\n      Segment1d.create(data.radius0, data.radius1),\n      AngleSweep.createStartEnd(data.bearing0!, data.bearing1!),\n      fractionInterval ? fractionInterval.clone() : Segment1d.create(0, 1),\n      localToWorld, data.curveLength!, data1);\n  }\n  /** Copy all defining data from another spiral. */\n  public setFrom(other: IntegratedSpiral3d): IntegratedSpiral3d {\n    this.localToWorld.setFrom(other.localToWorld);\n    this.radius01.setFrom(other.radius01);\n    this._curvature01.setFrom(other._curvature01);\n    this.bearing01.setFrom(other.bearing01);\n    this.localToWorld.setFrom(other.localToWorld);\n    this.activeFractionInterval.setFrom(other.activeFractionInterval);\n    this._arcLength01 = other._arcLength01;\n    return this;\n  }\n  /** Deep clone of this spiral */\n  public override clone(): IntegratedSpiral3d {\n    return new IntegratedSpiral3d(this._spiralType, this._evaluator,\n      this.radius01.clone(), this.bearing01.clone(),\n      this.activeFractionInterval.clone(), this.localToWorld.clone(), this._arcLength01,\n      this._designProperties?.clone());\n  }\n\n  /** apply `transform` to this spiral's local to world transform. */\n  public tryTransformInPlace(transformA: Transform): boolean {\n\n    const rigidData = this.applyRigidPartOfTransform(transformA);\n    if (rigidData !== undefined) {\n      this._curvature01.x0 /= rigidData.scale;\n      this._curvature01.x1 /= rigidData.scale;\n      this.radius01.x0 *= rigidData.scale;\n      this.radius01.x1 *= rigidData.scale;\n      this._arcLength01 *= rigidData.scale;\n    }\n    this.refreshComputedProperties();\n    return true;\n  }\n\n  /** Return the spiral start point. */\n  public override startPoint(result?: Point3d): Point3d {\n    return this.activeStrokes.startPoint(result);\n  }\n  /** return the spiral end point. */\n  public override endPoint(result?: Point3d): Point3d {\n    return this.activeStrokes.endPoint(result);\n  }\n  /** test if the local to world transform places the spiral xy plane into `plane` */\n  public isInPlane(plane: Plane3dByOriginAndUnitNormal): boolean {\n    return plane.isPointInPlane(this.localToWorld.origin as Point3d)\n      && Geometry.isSameCoordinate(0.0, this.localToWorld.matrix.dotColumnX(plane.getNormalRef()))\n      && Geometry.isSameCoordinate(0.0, this.localToWorld.matrix.dotColumnY(plane.getNormalRef()));\n  }\n  /** Return length of the spiral.  Because TransitionSpiral is parameterized directly in terms of distance along, this is a simple return value. */\n  public quickLength() { return this.curveLength(); }\n  /** Return length of the spiral.  Because TransitionSpiral is parameterized directly in terms of distance along, this is a simple return value. */\n  public override curveLength() { return this._arcLength01 * (this._activeFractionInterval.absoluteDelta()); }\n  /** Return (unsigned) length of the spiral between fractions.  Because TransitionSpiral is parameterized directly in terms of distance along, this is a simple return value. */\n  public override curveLengthBetweenFractions(fraction0: number, fraction1: number) {\n    return this._arcLength01 * (this._activeFractionInterval.absoluteDelta() * Math.abs(fraction1 - fraction0));\n  }\n  /** Test if `other` is an instance of `TransitionSpiral3d` */\n  public isSameGeometryClass(other: any): boolean { return other instanceof TransitionSpiral3d; }\n  /** Add strokes from this spiral to `dest`.\n   * * Linestrings will usually stroke as just their points.\n   * * If maxEdgeLength is given, this will sub-stroke within the linestring -- not what we want.\n   */\n  public emitStrokes(dest: LineString3d, options?: StrokeOptions): void { this.activeStrokes.emitStrokes(dest, options); }\n  /** emit stroke fragments to `dest` handler. */\n  public emitStrokableParts(dest: IStrokeHandler, options?: StrokeOptions): void {\n    const n = this.computeStrokeCountForOptions(options);\n    dest.startParentCurvePrimitive(this);\n    const activeStrokes = this.activeStrokes;\n    const preferPrimary = dest.needPrimaryGeometryForStrokes === undefined ? false : dest.needPrimaryGeometryForStrokes();\n    if (!preferPrimary && n <= activeStrokes.numPoints()) {\n      this.activeStrokes.emitStrokableParts(dest, options);\n    } else {\n      dest.announceIntervalForUniformStepStrokes(this, n, 0.0, 1.0);\n    }\n    dest.endParentCurvePrimitive(this);\n  }\n\n  /**\n   * return the stroke count required for given options.\n   * @param options StrokeOptions that determine count\n   */\n\n  public computeStrokeCountForOptions(options?: StrokeOptions): number {\n    let numStroke;\n    if (options) {\n      const rMin = Math.min(Math.abs(this.radius01.x0), Math.abs(this.radius01.x1));\n      numStroke = options.applyTolerancesToArc(rMin, this.bearing01.sweepRadians);\n      numStroke = options.applyMaxEdgeLength(numStroke, this.curveLength());\n      numStroke = options.applyMinStrokesPerPrimitive(numStroke);\n    } else {\n      numStroke = StrokeOptions.applyAngleTol(undefined, 4, this.bearing01.sweepRadians);\n    }\n    return numStroke;\n  }\n\n  /** Reverse the active interval and active strokes.\n   * * Primary defining data remains unchanged !!!\n   */\n\n  public reverseInPlace(): void {\n    this.activeFractionInterval.reverseInPlace();\n    if (this._activeStrokes === undefined)\n      this._activeStrokes = this._globalStrokes.clone();\n    this._activeStrokes.reverseInPlace();\n  }\n  /** Evaluate curve point with respect to fraction. */\n  public fractionToPoint(activeFraction: number, result?: Point3d): Point3d {\n    const targetGlobalFraction = this.activeFractionInterval.fractionToPoint(activeFraction);\n    const numStrokes = this._globalStrokes.packedPoints.length - 1;\n    if (activeFraction > 1.0) {\n      result = this._globalStrokes.packedPoints.back(result)!;\n      const integrationStep = 1.0 / numStrokes;\n      let currentGlobalFraction = 1.0;\n      let nextGlobalFraction = currentGlobalFraction + integrationStep;\n      while (nextGlobalFraction < targetGlobalFraction) {\n        this.fullSpiralIncrementalIntegral(result, currentGlobalFraction, nextGlobalFraction, true);\n        currentGlobalFraction = nextGlobalFraction;\n        nextGlobalFraction += integrationStep;\n      }\n      this.fullSpiralIncrementalIntegral(result, currentGlobalFraction, targetGlobalFraction, true);\n    } else if (activeFraction < 0.0) {\n      result = this._globalStrokes.packedPoints.front(result)!;\n      const integrationStep = 1.0 / numStrokes;\n      let currentGlobalFraction = 0.0;\n      let nextGlobalFraction = currentGlobalFraction - integrationStep;\n      while (nextGlobalFraction > targetGlobalFraction) {\n        this.fullSpiralIncrementalIntegral(result, currentGlobalFraction, nextGlobalFraction, true);\n        currentGlobalFraction = nextGlobalFraction;\n        nextGlobalFraction -= integrationStep;\n      }\n      this.fullSpiralIncrementalIntegral(result, currentGlobalFraction, targetGlobalFraction, true);\n    } else {\n      const clampedGlobalFraction = Geometry.clampToStartEnd(targetGlobalFraction, 0, 1);\n      const index0 = Math.trunc(clampedGlobalFraction * numStrokes); // This indexes the point to the left of the query.\n      const globalFraction0 = index0 / numStrokes;\n      result = this._globalStrokes.packedPoints.getPoint3dAtUncheckedPointIndex(index0, result);\n      // GeometryCoreTestIO.consoleLog(\" fractionToPoint \", activeFraction, this.activeFractionInterval, \"( global integration \" + globalFraction0 + \" to \" + globalFraction + \")\", index0);\n      this.fullSpiralIncrementalIntegral(result, globalFraction0, targetGlobalFraction, true);\n\n    }\n    return result;\n  }\n  /** Evaluate curve point and derivative with respect to fraction. */\n  public fractionToPointAndDerivative(activeFraction: number, result?: Ray3d): Ray3d {\n    const globalFraction = this.activeFractionInterval.fractionToPoint(activeFraction);\n    result = result ? result : Ray3d.createZero();\n    this.fractionToPoint(activeFraction, result.origin);\n    const radians = this.globalFractionToBearingRadians(globalFraction);\n    const a = this._arcLength01 * this.activeFractionInterval.signedDelta();\n    this.localToWorld.matrix.multiplyXY(a * Math.cos(radians), a * Math.sin(radians), result.direction);\n    return result;\n  }\n\n  /** Return the frenet frame at fractional position. */\n  public override fractionToFrenetFrame(activeFraction: number, result?: Transform): Transform {\n    const globalFraction = this.activeFractionInterval.fractionToPoint(activeFraction);\n    result = result ? result : Transform.createIdentity();\n    result.origin.setFrom(this.fractionToPoint(activeFraction));\n    Matrix3d.createRigidFromMatrix3d(this.localToWorld.matrix, AxisOrder.XYZ, result.matrix);\n\n    const radians = this.globalFractionToBearingRadians(globalFraction);\n    const c = Math.cos(radians);\n    const s = Math.sin(radians);\n    result.matrix.applyGivensColumnOp(0, 1, c, s);\n    return result;\n  }\n  /** Return a plane with\n   *\n   * * origin at fractional position along the curve\n   * * vectorU is the first derivative, i.e. tangent vector with length equal to the rate of change with respect to the fraction.\n   * * vectorV is the second derivative, i.e.derivative of vectorU.\n   */\n  public fractionToPointAnd2Derivatives(activeFraction: number, result?: Plane3dByOriginAndVectors): Plane3dByOriginAndVectors | undefined {\n    const globalFraction = this.activeFractionInterval.fractionToPoint(activeFraction);\n    const origin = this.fractionToPoint(activeFraction);\n    const radians = this.globalFractionToBearingRadians(globalFraction);\n    const c = Math.cos(radians);\n    const s = Math.sin(radians);\n    const delta = this.activeFractionInterval.signedDelta();\n    const a = delta;\n    const b = a * delta;\n    const vectorX = this.localToWorld.matrix.multiplyXY(a * c, a * s);\n    const vectorY = this.localToWorld.matrix.multiplyXY(-b * s, b * c);\n    vectorY.scaleInPlace(this.globalFractionToCurvature(globalFraction));\n    return Plane3dByOriginAndVectors.createCapture(origin, vectorX, vectorY, result);\n  }\n  /** Second step of double dispatch:  call `handler.handleTransitionSpiral(this)` */\n  public dispatchToGeometryHandler(handler: GeometryHandler): any {\n    return handler.handleTransitionSpiral(this);\n  }\n\n  /** compare various coordinate quantities */\n  public override isAlmostEqual(other?: GeometryQuery): boolean {\n    if (other instanceof IntegratedSpiral3d) {\n      return this.radius01.isAlmostEqual(other.radius01)\n        && this.bearing01.isAlmostEqualAllowPeriodShift(other.bearing01)\n        && this.localToWorld.isAlmostEqual(other.localToWorld)\n        && Geometry.isSameCoordinate(this._arcLength01, other._arcLength01)\n        && this.activeFractionInterval.isAlmostEqual(other.activeFractionInterval)\n        && this._curvature01.isAlmostEqual(other._curvature01);\n    }\n    return false;\n  }\n}\n// at load time, initialize gauss quadrature workspace\nIntegratedSpiral3d.initWorkSpace();\n"]}
+{"version":3,"file":"IntegratedSpiral3d.js","sourceRoot":"","sources":["../../../../src/curve/spiral/IntegratedSpiral3d.ts"],"names":[],"mappings":"AAAA;;;+FAG+F;AAC/F;;GAEG;AAEH,OAAO,EAAE,kBAAkB,EAAE,MAAM,sBAAsB,CAAC;AAC1D,OAAO,EAAE,SAAS,EAAE,MAAM,4BAA4B,CAAC;AACvD,OAAO,EAAE,UAAU,EAAE,MAAM,6BAA6B,CAAC;AACzD,OAAO,EAAE,SAAS,EAAE,MAAM,4BAA4B,CAAC;AACvD,OAAO,EAAE,YAAY,EAAE,MAAM,iBAAiB,CAAC;AAC/C,OAAO,EAAE,oBAAoB,EAAE,MAAM,wBAAwB,CAAC;AAC9D,OAAO,EAAE,+BAA+B,EAAE,MAAM,mCAAmC,CAAC;AACpF,OAAO,EAAE,UAAU,EAAE,MAAM,2BAA2B,CAAC;AACvD,OAAO,EAAE,OAAO,EAAE,MAAM,kCAAkC,CAAC;AAC3D,OAAO,EAAE,QAAQ,EAAE,MAAM,2BAA2B,CAAC;AAGrD,OAAO,EAAE,SAAS,EAAE,QAAQ,EAAE,MAAM,gBAAgB,CAAC;AACrD,OAAO,EAAE,aAAa,EAAE,MAAM,kBAAkB,CAAC;AAEjD,OAAO,EAAE,KAAK,EAAE,MAAM,wBAAwB,CAAC;AAC/C,OAAO,EAAE,yBAAyB,EAAE,MAAM,4CAA4C,CAAC;AAGvF;;;;;;;GAOG;AACH,MAAM,OAAO,kBAAmB,SAAQ,kBAAkB;IACxD,yCAAyC;IACzB,kBAAkB,GAAG,kBAAkB,CAAC;IACxD,0CAA0C;IACnC,QAAQ,CAAY;IAC3B,+CAA+C;IACxC,SAAS,CAAa;IAC7B,8CAA8C;IACtC,cAAc,CAAe;IACrC;;;;;OAKG;IACK,cAAc,CAAgB;IACtC,yEAAyE;IACzE,IAAW,aAAa;QACtB,OAAO,IAAI,CAAC,cAAc,KAAK,SAAS,CAAC,CAAC,CAAC,IAAI,CAAC,cAAc,CAAC,CAAC,CAAC,IAAI,CAAC,cAAc,CAAC;IACvF,CAAC;IACD,gCAAgC;IACxB,UAAU,CAAuB;IACzC,yCAAyC;IACjC,YAAY,CAAS;IAC7B,mDAAmD;IAC3C,YAAY,CAAY;IAChC,4DAA4D;IAC5D,mEAAmE;IACnE,YACE,UAA8B,EAC9B,SAA+B,EAC/B,QAAmB,EACnB,SAAqB,EACrB,sBAAiC,EACjC,YAAuB,EACvB,SAAiB,EACjB,UAAuD;QAEvD,KAAK,CAAC,UAAU,EAAE,YAAY,EAAE,sBAAsB,EAAE,UAAU,CAAC,CAAC;QACpE,IAAI,CAAC,UAAU,GAAG,SAAS,CAAC;QAC5B,IAAI,CAAC,QAAQ,GAAG,QAAQ,CAAC;QACzB,IAAI,CAAC,SAAS,GAAG,SAAS,CAAC;QAC3B,IAAI,CAAC,YAAY,GAAG,SAAS,CAAC;QAC9B,IAAI,CAAC,cAAc,GAAG,YAAY,CAAC,MAAM,EAAE,CAAC;QAC5C,mFAAmF;QACnF,IAAI,CAAC,YAAY,GAAG,SAAS,CAAC,MAAM,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;QAC3C,IAAI,CAAC,yBAAyB,EAAE,CAAC;IACnC,CAAC;IACD,0CAA0C;IACnC,MAAM,CAAU,iBAAiB,GAAG,UAAU,CAAC;IACtD,6EAA6E;IACtE,8BAA8B,CAAC,QAAgB;QACpD,MAAM,YAAY,GAAG,IAAI,CAAC,UAAU,CAAC,cAAc,CAAC,QAAQ,CAAC,CAAC;QAC9D,MAAM,EAAE,GAAG,IAAI,CAAC,YAAY,CAAC;QAC7B,OAAO,IAAI,CAAC,SAAS,CAAC,YAAY;cAC9B,YAAY,GAAG,EAAE,GAAG,IAAI,CAAC,YAAY,CAAC,WAAW,EAAE;cACnD,QAAQ,GAAG,IAAI,CAAC,YAAY,CAAC,EAAE,GAAG,EAAE,CAAC;IAC3C,CAAC;IACD,6EAA6E;IACtE,yBAAyB,CAAC,QAAgB;QAC/C,MAAM,CAAC,GAAG,IAAI,CAAC,UAAU,CAAC,2BAA2B,CAAC,QAAQ,CAAC,CAAC;QAChE,OAAO,IAAI,CAAC,YAAY,CAAC,eAAe,CAAC,CAAC,CAAC,CAAC;IAC9C,CAAC;IACD,mEAAmE;IAC5D,wBAAwB,CAAC,cAAsB;QACpD,MAAM,QAAQ,GAAG,IAAI,CAAC,sBAAsB,CAAC,eAAe,CAAC,cAAc,CAAC,CAAC;QAC7E,OAAO,IAAI,CAAC,SAAS,CAAC,YAAY;cAC9B,QAAQ,GAAG,IAAI,CAAC,YAAY,GAAG,CAAC,IAAI,CAAC,YAAY,CAAC,EAAE,GAAG,GAAG,GAAG,QAAQ,GAAG,CAAC,IAAI,CAAC,YAAY,CAAC,EAAE,GAAG,IAAI,CAAC,YAAY,CAAC,EAAE,CAAC,CAAC,CAAC;IAC7H,CAAC;IACD;;;OAGG;IACa,mBAAmB,CAAC,cAAsB;QACxD,uBAAuB;QACvB,OAAO,IAAI,CAAC,YAAY,CAAC,eAAe,CAAC,IAAI,CAAC,sBAAsB,CAAC,eAAe,CAAC,cAAc,CAAC,CAAC,CAAC;IACxG,CAAC;IACD,2EAA2E;IACnE,MAAM,CAAC,cAAc,CAAe;IACpC,MAAM,CAAC,YAAY,CAAe;IAClC,MAAM,CAAC,YAAY,CAAiF;IAC5G,0CAA0C;IACnC,MAAM,CAAC,aAAa;QACzB,kBAAkB,CAAC,cAAc,GAAG,IAAI,YAAY,CAAC,CAAC,CAAC,CAAC;QACxD,kBAAkB,CAAC,YAAY,GAAG,IAAI,YAAY,CAAC,CAAC,CAAC,CAAC;QACtD,kBAAkB,CAAC,YAAY,GAAG,CAAC,EAAE,EAAE,EAAE,EAAE,OAAO,EAAE,OAAO,EAAE,EAAE,CAAC,UAAU,CAAC,WAAW,CAAC,EAAE,EAAE,EAAE,EAAE,OAAO,EAAE,OAAO,CAAC,CAAC;IACnH,CAAC;IACD;;;;;;;;;OASG;IACK,6BAA6B,CAAC,GAAY,EAAE,SAAiB,EAAE,SAAiB,EAAE,WAAoB;QAC5G,MAAM,aAAa,GAAG,kBAAkB,CAAC,cAAc,CAAC;QACxD,MAAM,WAAW,GAAG,kBAAkB,CAAC,YAAY,CAAC;QACpD,MAAM,OAAO,GAAG,kBAAkB,CAAC,YAAY,CAAC,SAAS,EAAE,SAAS,EAAE,aAAa,EAAE,WAAW,CAAC,CAAC;QAClG,MAAM,MAAM,GAAG,IAAI,CAAC,YAAY,CAAC;QACjC,IAAI,CAAC,GAAG,CAAC,CAAC;QACV,IAAI,EAAE,GAAG,GAAG,CAAC;QACb,IAAI,EAAE,GAAG,GAAG,CAAC;QACb,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,OAAO,EAAE,CAAC,EAAE,EAAE,CAAC;YACjC,MAAM,OAAO,GAAG,IAAI,CAAC,8BAA8B,CAAC,aAAa,CAAC,CAAC,CAAC,CAAC,CAAC;YACtE,CAAC,GAAG,WAAW,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC;YAC5B,EAAE,IAAI,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,OAAO,CAAC,CAAC;YAC5B,EAAE,IAAI,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,OAAO,CAAC,CAAC;QAC9B,CAAC;QACD,IAAI,WAAW;YACb,QAAQ,CAAC,qBAAqB,CAAC,GAAG,EAAE,IAAI,CAAC,YAAY,CAAC,MAAM,EAAE,EAAE,CAAC,EAAE,EAAE,EAAE,CAAC,EAAE,EAAE,EAAE,CAAC,EAAE,GAAG,EAAE,EAAE,GAAG,CAAC,CAAC;;YAE7F,GAAG,CAAC,aAAa,CAAC,EAAE,EAAE,EAAE,EAAE,GAAG,CAAC,CAAC;IACnC,CAAC;IACD,yBAAyB;IACT,yBAAyB;QACvC,IAAI,CAAC,YAAY,GAAG,SAAS,CAAC,MAAM,CAClC,kBAAkB,CAAC,iBAAiB,CAAC,IAAI,CAAC,QAAQ,CAAC,EAAE,CAAC,EACtD,kBAAkB,CAAC,iBAAiB,CAAC,IAAI,CAAC,QAAQ,CAAC,EAAE,CAAC,CACvD,CAAC;QACF,IAAI,CAAC,cAAc,CAAC,KAAK,EAAE,CAAC;QAC5B,MAAM,YAAY,GAAG,OAAO,CAAC,MAAM,EAAE,CAAC;QACtC,IAAI,CAAC,cAAc,CAAC,iBAAiB,CAAC,YAAY,CAAC,CAAC;QACpD,MAAM,WAAW,GAAG,EAAE,CAAC;QACvB,MAAM,YAAY,GAAG,GAAG,GAAG,WAAW,CAAC;QACvC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,WAAW,EAAE,CAAC,EAAE,EAAE,CAAC;YACtC,MAAM,SAAS,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,YAAY,CAAC;YACzC,MAAM,SAAS,GAAG,CAAC,GAAG,YAAY,CAAC;YACnC,IAAI,CAAC,6BAA6B,CAAC,YAAY,EAAE,SAAS,EAAE,SAAS,EAAE,KAAK,CAAC,CAAC;YAC9E,IAAI,CAAC,cAAc,CAAC,iBAAiB,CAAC,YAAY,CAAC,CAAC;QACtD,CAAC;QACD,IAAI,CAAC,cAAc,CAAC,mBAAmB,CAAC,IAAI,CAAC,YAAY,CAAC,CAAC;QAC3D,IAAI,CAAC,IAAI,CAAC,sBAAsB,CAAC,SAAS,EAAE,CAAC;YAC3C,IAAI,IAAI,CAAC,cAAc,KAAK,SAAS;gBACnC,IAAI,CAAC,cAAc,GAAG,YAAY,CAAC,MAAM,EAAE,CAAC;YAC9C,IAAI,CAAC,cAAc,CAAC,KAAK,EAAE,CAAC;YAC5B,mGAAmG;YACnG,0FAA0F;YAC1F,0CAA0C;YAC1C,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,WAAW,EAAE,CAAC,EAAE,EAAE,CAAC;gBACtC,MAAM,aAAa,GAAG,CAAC,GAAG,YAAY,CAAC;gBACvC,IAAI,CAAC,cAAc,CAAC,QAAQ,CAAC,IAAI,CAAC,eAAe,CAAC,aAAa,CAAC,CAAC,CAAC;YACpE,CAAC;QACH,CAAC;IACH,CAAC;IACD;;;;;;;OAOG;IACI,MAAM,CAAC,gCAAgC,CAC5C,QAAmB,EACnB,SAAqB,EACrB,sBAAiC,EACjC,YAAuB,EACvB,QAAiB;QAEjB,MAAM,SAAS,GAAG,kBAAkB,CAAC,mCAAmC,CACtE,QAAQ,CAAC,EAAE,EAAE,QAAQ,CAAC,EAAE,EAAE,SAAS,CAAC,YAAY,CACjD,CAAC;QACF,IAAI,QAAQ,KAAK,SAAS;YACxB,QAAQ,GAAG,UAAU,CAAC;QACxB,MAAM,SAAS,GAAG,oBAAoB,CAAC,aAAa,CAAC,QAAQ,CAAC,CAAC;QAC/D,IAAI,CAAC,SAAS;YACZ,OAAO,SAAS,CAAC;QACnB,OAAO,IAAI,kBAAkB,CAC3B,QAAQ,EACR,SAAS,EACT,QAAQ,CAAC,KAAK,EAAE,EAChB,SAAS,CAAC,KAAK,EAAE,EACjB,sBAAsB,CAAC,KAAK,EAAE,EAC9B,YAAY,CAAC,KAAK,EAAE,EACpB,SAAS,EACT,IAAI,+BAA+B,CACjC,QAAQ,CAAC,EAAE,EAAE,QAAQ,CAAC,EAAE,EAAE,SAAS,CAAC,UAAU,CAAC,KAAK,EAAE,EAAE,SAAS,CAAC,QAAQ,CAAC,KAAK,EAAE,EAAE,SAAS,CAC9F,CACF,CAAC;IACJ,CAAC;IACD;;;;;;;;;;;OAWG;IACI,MAAM,CAAC,iBAAiB,CAC7B,UAA8B,EAC9B,OAA2B,EAC3B,OAA2B,EAC3B,QAA2B,EAC3B,QAA2B,EAC3B,SAA6B,EAC7B,gBAAuC,EACvC,YAAuB;QAEvB,IAAI,UAAU,KAAK,SAAS;YAC1B,UAAU,GAAG,UAAU,CAAC;QAC1B,MAAM,SAAS,GAAG,oBAAoB,CAAC,aAAa,CAAC,UAAU,CAAC,CAAC;QACjE,IAAI,CAAC,SAAS;YACZ,OAAO,SAAS,CAAC;QACnB,MAAM,IAAI,GAAG,IAAI,+BAA+B,CAAC,OAAO,EAAE,OAAO,EAAE,QAAQ,EAAE,QAAQ,EAAE,SAAS,CAAC,CAAC;QAClG,MAAM,KAAK,GAAG,IAAI,CAAC,KAAK,EAAE,CAAC;QAC3B,IAAI,CAAC,IAAI,CAAC,0BAA0B,EAAE;YACpC,OAAO,SAAS,CAAC;QACnB,IAAI,gBAAgB,KAAK,SAAS;YAChC,gBAAgB,GAAG,SAAS,CAAC,MAAM,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;QAC5C,OAAO,IAAI,kBAAkB,CAC3B,UAAU,EACV,SAAS,EACT,SAAS,CAAC,MAAM,CAAC,IAAI,CAAC,OAAO,EAAE,IAAI,CAAC,OAAO,CAAC,EAC5C,UAAU,CAAC,cAAc,CAAC,IAAI,CAAC,QAAS,EAAE,IAAI,CAAC,QAAS,CAAC,EACzD,gBAAgB,CAAC,CAAC,CAAC,gBAAgB,CAAC,KAAK,EAAE,CAAC,CAAC,CAAC,SAAS,CAAC,MAAM,CAAC,CAAC,EAAE,CAAC,CAAC,EACpE,YAAY,EACZ,IAAI,CAAC,WAAY,EACjB,KAAK,CACN,CAAC;IACJ,CAAC;IACD,kDAAkD;IAC3C,OAAO,CAAC,KAAyB;QACtC,IAAI,CAAC,YAAY,CAAC,OAAO,CAAC,KAAK,CAAC,YAAY,CAAC,CAAC;QAC9C,IAAI,CAAC,QAAQ,CAAC,OAAO,CAAC,KAAK,CAAC,QAAQ,CAAC,CAAC;QACtC,IAAI,CAAC,YAAY,CAAC,OAAO,CAAC,KAAK,CAAC,YAAY,CAAC,CAAC;QAC9C,IAAI,CAAC,SAAS,CAAC,OAAO,CAAC,KAAK,CAAC,SAAS,CAAC,CAAC;QACxC,IAAI,CAAC,YAAY,CAAC,OAAO,CAAC,KAAK,CAAC,YAAY,CAAC,CAAC;QAC9C,IAAI,CAAC,sBAAsB,CAAC,OAAO,CAAC,KAAK,CAAC,sBAAsB,CAAC,CAAC;QAClE,IAAI,CAAC,YAAY,GAAG,KAAK,CAAC,YAAY,CAAC;QACvC,OAAO,IAAI,CAAC;IACd,CAAC;IACD,iCAAiC;IACjB,KAAK;QACnB,OAAO,IAAI,kBAAkB,CAC3B,IAAI,CAAC,WAAW,EAChB,IAAI,CAAC,UAAU,EACf,IAAI,CAAC,QAAQ,CAAC,KAAK,EAAE,EACrB,IAAI,CAAC,SAAS,CAAC,KAAK,EAAE,EACtB,IAAI,CAAC,sBAAsB,CAAC,KAAK,EAAE,EACnC,IAAI,CAAC,YAAY,CAAC,KAAK,EAAE,EACzB,IAAI,CAAC,YAAY,EACjB,IAAI,CAAC,iBAAiB,EAAE,KAAK,EAAE,CAChC,CAAC;IACJ,CAAC;IACD,mEAAmE;IAC5D,mBAAmB,CAAC,UAAqB;QAC9C,MAAM,SAAS,GAAG,IAAI,CAAC,yBAAyB,CAAC,UAAU,CAAC,CAAC;QAC7D,IAAI,SAAS,KAAK,SAAS,EAAE,CAAC;YAC5B,IAAI,CAAC,YAAY,CAAC,EAAE,IAAI,SAAS,CAAC,KAAK,CAAC;YACxC,IAAI,CAAC,YAAY,CAAC,EAAE,IAAI,SAAS,CAAC,KAAK,CAAC;YACxC,IAAI,CAAC,QAAQ,CAAC,EAAE,IAAI,SAAS,CAAC,KAAK,CAAC;YACpC,IAAI,CAAC,QAAQ,CAAC,EAAE,IAAI,SAAS,CAAC,KAAK,CAAC;YACpC,IAAI,CAAC,YAAY,IAAI,SAAS,CAAC,KAAK,CAAC;QACvC,CAAC;QACD,IAAI,CAAC,yBAAyB,EAAE,CAAC;QACjC,OAAO,IAAI,CAAC;IACd,CAAC;IACD,qCAAqC;IACrB,UAAU,CAAC,MAAgB;QACzC,OAAO,IAAI,CAAC,aAAa,CAAC,UAAU,CAAC,MAAM,CAAC,CAAC;IAC/C,CAAC;IACD,mCAAmC;IACnB,QAAQ,CAAC,MAAgB;QACvC,OAAO,IAAI,CAAC,aAAa,CAAC,QAAQ,CAAC,MAAM,CAAC,CAAC;IAC7C,CAAC;IACD,oFAAoF;IAC7E,SAAS,CAAC,KAAmC;QAClD,OAAO,KAAK,CAAC,cAAc,CAAC,IAAI,CAAC,YAAY,CAAC,MAAiB,CAAC;eAC3D,QAAQ,CAAC,gBAAgB,CAAC,GAAG,EAAE,IAAI,CAAC,YAAY,CAAC,MAAM,CAAC,UAAU,CAAC,KAAK,CAAC,YAAY,EAAE,CAAC,CAAC;eACzF,QAAQ,CAAC,gBAAgB,CAAC,GAAG,EAAE,IAAI,CAAC,YAAY,CAAC,MAAM,CAAC,UAAU,CAAC,KAAK,CAAC,YAAY,EAAE,CAAC,CAAC,CAAC;IACjG,CAAC;IACD;;;OAGG;IACI,WAAW;QAChB,OAAO,IAAI,CAAC,WAAW,EAAE,CAAC;IAC5B,CAAC;IACD;;;OAGG;IACa,WAAW;QACzB,OAAO,IAAI,CAAC,YAAY,GAAG,CAAC,IAAI,CAAC,uBAAuB,CAAC,aAAa,EAAE,CAAC,CAAC;IAC5E,CAAC;IACD;;;OAGG;IACa,2BAA2B,CAAC,SAAiB,EAAE,SAAiB;QAC9E,OAAO,IAAI,CAAC,YAAY,GAAG,CAAC,IAAI,CAAC,uBAAuB,CAAC,aAAa,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,SAAS,GAAG,SAAS,CAAC,CAAC,CAAC;IAC9G,CAAC;IACD,8DAA8D;IACvD,mBAAmB,CAAC,KAAU;QACnC,OAAO,KAAK,YAAY,kBAAkB,CAAC;IAC7C,CAAC;IACD;;;;OAIG;IACI,WAAW,CAAC,IAAkB,EAAE,OAAuB;QAC5D,IAAI,CAAC,aAAa,CAAC,WAAW,CAAC,IAAI,EAAE,OAAO,CAAC,CAAC;IAChD,CAAC;IACD,+CAA+C;IACxC,kBAAkB,CAAC,IAAoB,EAAE,OAAuB;QACrE,MAAM,CAAC,GAAG,IAAI,CAAC,4BAA4B,CAAC,OAAO,CAAC,CAAC;QACrD,IAAI,CAAC,yBAAyB,CAAC,IAAI,CAAC,CAAC;QACrC,MAAM,aAAa,GAAG,IAAI,CAAC,aAAa,CAAC;QACzC,MAAM,aAAa,GAAG,IAAI,CAAC,6BAA6B,KAAK,SAAS,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,IAAI,CAAC,6BAA6B,EAAE,CAAC;QACtH,IAAI,CAAC,aAAa,IAAI,CAAC,IAAI,aAAa,CAAC,SAAS,EAAE,EAAE,CAAC;YACrD,IAAI,CAAC,aAAa,CAAC,kBAAkB,CAAC,IAAI,EAAE,OAAO,CAAC,CAAC;QACvD,CAAC;aAAM,CAAC;YACN,IAAI,CAAC,qCAAqC,CAAC,IAAI,EAAE,CAAC,EAAE,GAAG,EAAE,GAAG,CAAC,CAAC;QAChE,CAAC;QACD,IAAI,CAAC,uBAAuB,CAAC,IAAI,CAAC,CAAC;IACrC,CAAC;IACD;;;OAGG;IACI,4BAA4B,CAAC,OAAuB;QACzD,IAAI,SAAS,CAAC;QACd,IAAI,OAAO,EAAE,CAAC;YACZ,MAAM,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,QAAQ,CAAC,EAAE,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,QAAQ,CAAC,EAAE,CAAC,CAAC,CAAC;YAC9E,SAAS,GAAG,OAAO,CAAC,oBAAoB,CAAC,IAAI,EAAE,IAAI,CAAC,SAAS,CAAC,YAAY,CAAC,CAAC;YAC5E,SAAS,GAAG,OAAO,CAAC,kBAAkB,CAAC,SAAS,EAAE,IAAI,CAAC,WAAW,EAAE,CAAC,CAAC;YACtE,SAAS,GAAG,OAAO,CAAC,2BAA2B,CAAC,SAAS,CAAC,CAAC;QAC7D,CAAC;aAAM,CAAC;YACN,SAAS,GAAG,aAAa,CAAC,aAAa,CAAC,SAAS,EAAE,CAAC,EAAE,IAAI,CAAC,SAAS,CAAC,YAAY,CAAC,CAAC;QACrF,CAAC;QACD,OAAO,SAAS,CAAC;IACnB,CAAC;IACD;;;OAGG;IACI,cAAc;QACnB,IAAI,CAAC,sBAAsB,CAAC,cAAc,EAAE,CAAC;QAC7C,IAAI,IAAI,CAAC,cAAc,KAAK,SAAS;YACnC,IAAI,CAAC,cAAc,GAAG,IAAI,CAAC,cAAc,CAAC,KAAK,EAAE,CAAC;QACpD,IAAI,CAAC,cAAc,CAAC,cAAc,EAAE,CAAC;IACvC,CAAC;IACD,qDAAqD;IAC9C,eAAe,CAAC,cAAsB,EAAE,MAAgB;QAC7D,MAAM,oBAAoB,GAAG,IAAI,CAAC,sBAAsB,CAAC,eAAe,CAAC,cAAc,CAAC,CAAC;QACzF,MAAM,UAAU,GAAG,IAAI,CAAC,cAAc,CAAC,YAAY,CAAC,MAAM,GAAG,CAAC,CAAC;QAC/D,IAAI,cAAc,GAAG,GAAG,EAAE,CAAC;YACzB,MAAM,GAAG,IAAI,CAAC,cAAc,CAAC,YAAY,CAAC,IAAI,CAAC,MAAM,CAAE,CAAC;YACxD,MAAM,eAAe,GAAG,GAAG,GAAG,UAAU,CAAC;YACzC,IAAI,qBAAqB,GAAG,GAAG,CAAC;YAChC,IAAI,kBAAkB,GAAG,qBAAqB,GAAG,eAAe,CAAC;YACjE,OAAO,kBAAkB,GAAG,oBAAoB,EAAE,CAAC;gBACjD,IAAI,CAAC,6BAA6B,CAAC,MAAM,EAAE,qBAAqB,EAAE,kBAAkB,EAAE,IAAI,CAAC,CAAC;gBAC5F,qBAAqB,GAAG,kBAAkB,CAAC;gBAC3C,kBAAkB,IAAI,eAAe,CAAC;YACxC,CAAC;YACD,IAAI,CAAC,6BAA6B,CAAC,MAAM,EAAE,qBAAqB,EAAE,oBAAoB,EAAE,IAAI,CAAC,CAAC;QAChG,CAAC;aAAM,IAAI,cAAc,GAAG,GAAG,EAAE,CAAC;YAChC,MAAM,GAAG,IAAI,CAAC,cAAc,CAAC,YAAY,CAAC,KAAK,CAAC,MAAM,CAAE,CAAC;YACzD,MAAM,eAAe,GAAG,GAAG,GAAG,UAAU,CAAC;YACzC,IAAI,qBAAqB,GAAG,GAAG,CAAC;YAChC,IAAI,kBAAkB,GAAG,qBAAqB,GAAG,eAAe,CAAC;YACjE,OAAO,kBAAkB,GAAG,oBAAoB,EAAE,CAAC;gBACjD,IAAI,CAAC,6BAA6B,CAAC,MAAM,EAAE,qBAAqB,EAAE,kBAAkB,EAAE,IAAI,CAAC,CAAC;gBAC5F,qBAAqB,GAAG,kBAAkB,CAAC;gBAC3C,kBAAkB,IAAI,eAAe,CAAC;YACxC,CAAC;YACD,IAAI,CAAC,6BAA6B,CAAC,MAAM,EAAE,qBAAqB,EAAE,oBAAoB,EAAE,IAAI,CAAC,CAAC;QAChG,CAAC;aAAM,CAAC;YACN,MAAM,qBAAqB,GAAG,QAAQ,CAAC,eAAe,CAAC,oBAAoB,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;YACnF,MAAM,MAAM,GAAG,IAAI,CAAC,KAAK,CAAC,qBAAqB,GAAG,UAAU,CAAC,CAAC,CAAC,kDAAkD;YACjH,MAAM,eAAe,GAAG,MAAM,GAAG,UAAU,CAAC;YAC5C,MAAM,GAAG,IAAI,CAAC,cAAc,CAAC,YAAY,CAAC,+BAA+B,CAAC,MAAM,EAAE,MAAM,CAAC,CAAC;YAC1F,qLAAqL;YACrL,IAAI,CAAC,6BAA6B,CAAC,MAAM,EAAE,eAAe,EAAE,oBAAoB,EAAE,IAAI,CAAC,CAAC;QAE1F,CAAC;QACD,OAAO,MAAM,CAAC;IAChB,CAAC;IACD,oEAAoE;IAC7D,4BAA4B,CAAC,cAAsB,EAAE,MAAc;QACxE,MAAM,cAAc,GAAG,IAAI,CAAC,sBAAsB,CAAC,eAAe,CAAC,cAAc,CAAC,CAAC;QACnF,MAAM,GAAG,MAAM,CAAC,CAAC,CAAC,MAAM,CAAC,CAAC,CAAC,KAAK,CAAC,UAAU,EAAE,CAAC;QAC9C,IAAI,CAAC,eAAe,CAAC,cAAc,EAAE,MAAM,CAAC,MAAM,CAAC,CAAC;QACpD,MAAM,OAAO,GAAG,IAAI,CAAC,8BAA8B,CAAC,cAAc,CAAC,CAAC;QACpE,MAAM,CAAC,GAAG,IAAI,CAAC,YAAY,GAAG,IAAI,CAAC,sBAAsB,CAAC,WAAW,EAAE,CAAC;QACxE,IAAI,CAAC,YAAY,CAAC,MAAM,CAAC,UAAU,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,OAAO,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,OAAO,CAAC,EAAE,MAAM,CAAC,SAAS,CAAC,CAAC;QACpG,OAAO,MAAM,CAAC;IAChB,CAAC;IACD,sDAAsD;IACtC,qBAAqB,CAAC,cAAsB,EAAE,MAAkB;QAC9E,MAAM,cAAc,GAAG,IAAI,CAAC,sBAAsB,CAAC,eAAe,CAAC,cAAc,CAAC,CAAC;QACnF,MAAM,GAAG,MAAM,CAAC,CAAC,CAAC,MAAM,CAAC,CAAC,CAAC,SAAS,CAAC,cAAc,EAAE,CAAC;QACtD,MAAM,CAAC,MAAM,CAAC,OAAO,CAAC,IAAI,CAAC,eAAe,CAAC,cAAc,CAAC,CAAC,CAAC;QAC5D,QAAQ,CAAC,uBAAuB,CAAC,IAAI,CAAC,YAAY,CAAC,MAAM,EAAE,SAAS,CAAC,GAAG,EAAE,MAAM,CAAC,MAAM,CAAC,CAAC;QACzF,MAAM,OAAO,GAAG,IAAI,CAAC,8BAA8B,CAAC,cAAc,CAAC,CAAC;QACpE,MAAM,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,OAAO,CAAC,CAAC;QAC5B,MAAM,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,OAAO,CAAC,CAAC;QAC5B,MAAM,CAAC,MAAM,CAAC,mBAAmB,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;QAC9C,OAAO,MAAM,CAAC;IAChB,CAAC;IACD;;;;;;OAMG;IACI,8BAA8B,CACnC,cAAsB,EAAE,MAAkC;QAE1D,MAAM,cAAc,GAAG,IAAI,CAAC,sBAAsB,CAAC,eAAe,CAAC,cAAc,CAAC,CAAC;QACnF,MAAM,MAAM,GAAG,IAAI,CAAC,eAAe,CAAC,cAAc,CAAC,CAAC;QACpD,MAAM,OAAO,GAAG,IAAI,CAAC,8BAA8B,CAAC,cAAc,CAAC,CAAC;QACpE,MAAM,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,OAAO,CAAC,CAAC;QAC5B,MAAM,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,OAAO,CAAC,CAAC;QAC5B,MAAM,KAAK,GAAG,IAAI,CAAC,sBAAsB,CAAC,WAAW,EAAE,CAAC;QACxD,MAAM,CAAC,GAAG,KAAK,CAAC;QAChB,MAAM,CAAC,GAAG,CAAC,GAAG,KAAK,CAAC;QACpB,MAAM,OAAO,GAAG,IAAI,CAAC,YAAY,CAAC,MAAM,CAAC,UAAU,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC;QAClE,MAAM,OAAO,GAAG,IAAI,CAAC,YAAY,CAAC,MAAM,CAAC,UAAU,CAAC,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC;QACnE,OAAO,CAAC,YAAY,CAAC,IAAI,CAAC,yBAAyB,CAAC,cAAc,CAAC,CAAC,CAAC;QACrE,OAAO,yBAAyB,CAAC,aAAa,CAAC,MAAM,EAAE,OAAO,EAAE,OAAO,EAAE,MAAM,CAAC,CAAC;IACnF,CAAC;IACD,oFAAoF;IAC7E,yBAAyB,CAAC,OAAwB;QACvD,OAAO,OAAO,CAAC,sBAAsB,CAAC,IAAI,CAAC,CAAC;IAC9C,CAAC;IACD,6CAA6C;IAC7B,aAAa,CAAC,KAAqB;QACjD,IAAI,KAAK,YAAY,kBAAkB,EAAE,CAAC;YACxC,OAAO,IAAI,CAAC,QAAQ,CAAC,aAAa,CAAC,KAAK,CAAC,QAAQ,CAAC;mBAC7C,IAAI,CAAC,SAAS,CAAC,6BAA6B,CAAC,KAAK,CAAC,SAAS,CAAC;mBAC7D,IAAI,CAAC,YAAY,CAAC,aAAa,CAAC,KAAK,CAAC,YAAY,CAAC;mBACnD,QAAQ,CAAC,gBAAgB,CAAC,IAAI,CAAC,YAAY,EAAE,KAAK,CAAC,YAAY,CAAC;mBAChE,IAAI,CAAC,sBAAsB,CAAC,aAAa,CAAC,KAAK,CAAC,sBAAsB,CAAC;mBACvE,IAAI,CAAC,YAAY,CAAC,aAAa,CAAC,KAAK,CAAC,YAAY,CAAC,CAAC;QAC3D,CAAC;QACD,OAAO,KAAK,CAAC;IACf,CAAC;;AAEH,sDAAsD;AACtD,kBAAkB,CAAC,aAAa,EAAE,CAAC","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/** @packageDocumentation\r\n * @module Curve\r\n */\r\n\r\nimport { TransitionSpiral3d } from \"./TransitionSpiral3d\";\r\nimport { Segment1d } from \"../../geometry3d/Segment1d\";\r\nimport { AngleSweep } from \"../../geometry3d/AngleSweep\";\r\nimport { Transform } from \"../../geometry3d/Transform\";\r\nimport { LineString3d } from \"../LineString3d\";\r\nimport { NormalizedTransition } from \"./NormalizedTransition\";\r\nimport { TransitionConditionalProperties } from \"./TransitionConditionalProperties\";\r\nimport { Quadrature } from \"../../numerics/Quadrature\";\r\nimport { Point3d } from \"../../geometry3d/Point3dVector3d\";\r\nimport { Matrix3d } from \"../../geometry3d/Matrix3d\";\r\nimport { Angle } from \"../../geometry3d/Angle\";\r\nimport { Plane3dByOriginAndUnitNormal } from \"../../geometry3d/Plane3dByOriginAndUnitNormal\";\r\nimport { AxisOrder, Geometry } from \"../../Geometry\";\r\nimport { StrokeOptions } from \"../StrokeOptions\";\r\nimport { GeometryHandler, IStrokeHandler } from \"../../geometry3d/GeometryHandler\";\r\nimport { Ray3d } from \"../../geometry3d/Ray3d\";\r\nimport { Plane3dByOriginAndVectors } from \"../../geometry3d/Plane3dByOriginAndVectors\";\r\nimport { GeometryQuery } from \"../GeometryQuery\";\r\n\r\n/**\r\n * An IntegratedSpiral3d is a curve defined by integrating its curvature.\r\n * * The first integral of curvature (with respect to distance along the curve) is the bearing angle (in radians).\r\n * * Integrating (cos(theta), sin(theta)) gives displacement from the start point, and thus the actual curve position.\r\n * * The curvature functions of interest are all symmetric snap functions in the [[NormalizedTransition]] class.\r\n * * [[TransitionConditionalProperties]] implements the computations of the interrelationship of radii, bearing, and length.\r\n * @public\r\n */\r\nexport class IntegratedSpiral3d extends TransitionSpiral3d {\r\n  /** String name for schema properties. */\r\n  public readonly curvePrimitiveType = \"transitionSpiral\";\r\n  /** Start and end radii as a Segment1d. */\r\n  public radius01: Segment1d;\r\n  /** Start and end bearings as an AngleSweep. */\r\n  public bearing01: AngleSweep;\r\n  /** Stroked approximation of entire spiral. */\r\n  private _globalStrokes: LineString3d;\r\n  /**\r\n   * Stroked approximation of active spiral.\r\n   * * Same count as global -- possibly overly fine, but it gives some consistency between same clothoid constructed\r\n   * as partial versus complete.\r\n   * * If no trimming, this points to the same place as the _globalStrokes.  DO NOT double transform.\r\n   */\r\n  private _activeStrokes?: LineString3d;\r\n  /** Return the internal stroked form of the (possibly partial) spiral. */\r\n  public get activeStrokes(): LineString3d {\r\n    return this._activeStrokes !== undefined ? this._activeStrokes : this._globalStrokes;\r\n  }\r\n  /** Evaluator for transition. */\r\n  private _evaluator: NormalizedTransition;\r\n  /** Total curve arc length (computed). */\r\n  private _arcLength01: number;\r\n  /** Curvatures (inverse radii) at start and end. */\r\n  private _curvature01: Segment1d;\r\n  // constructor demands all bearing, radius, and length data.\r\n  // caller determines usual dependency of \"any 4 determine the 5th\".\r\n  private constructor(\r\n    spiralType: string | undefined,\r\n    evaluator: NormalizedTransition,\r\n    radius01: Segment1d,\r\n    bearing01: AngleSweep,\r\n    activeFractionInterval: Segment1d,\r\n    localToWorld: Transform,\r\n    arcLength: number,\r\n    properties: TransitionConditionalProperties | undefined,\r\n  ) {\r\n    super(spiralType, localToWorld, activeFractionInterval, properties);\r\n    this._evaluator = evaluator;\r\n    this.radius01 = radius01;\r\n    this.bearing01 = bearing01;\r\n    this._arcLength01 = arcLength;\r\n    this._globalStrokes = LineString3d.create();\r\n    // initialize for compiler; but this will be recomputed in refreshComputeProperties\r\n    this._curvature01 = Segment1d.create(0, 1);\r\n    this.refreshComputedProperties();\r\n  }\r\n  /** Default spiral type name: clothoid. */\r\n  public static readonly defaultSpiralType = \"clothoid\";\r\n  /** Use the integrated function to return an angle at fractional position. */\r\n  public globalFractionToBearingRadians(fraction: number): number {\r\n    const areaFraction = this._evaluator.fractionToArea(fraction);\r\n    const dx = this._arcLength01;\r\n    return this.bearing01.startRadians\r\n      + areaFraction * dx * this._curvature01.signedDelta()\r\n      + fraction * this._curvature01.x0 * dx;\r\n  }\r\n  /** Use the integrated function to return an angle at fractional position. */\r\n  public globalFractionToCurvature(fraction: number): number {\r\n    const f = this._evaluator.fractionToCurvatureFraction(fraction);\r\n    return this._curvature01.fractionToPoint(f);\r\n  }\r\n  /** Return the bearing at given fraction of the active interval. */\r\n  public fractionToBearingRadians(activeFraction: number): number {\r\n    const fraction = this.activeFractionInterval.fractionToPoint(activeFraction);\r\n    return this.bearing01.startRadians\r\n      + fraction * this._arcLength01 * (this._curvature01.x0 + 0.5 * fraction * (this._curvature01.x1 - this._curvature01.x0));\r\n  }\r\n  /**\r\n   * Return the curvature at given fraction of the active interval.\r\n   * * The `undefined` result is to match the abstract class -- it cannot actually occur.\r\n   */\r\n  public override fractionToCurvature(activeFraction: number): number | undefined {\r\n    // BUG? active interval\r\n    return this._curvature01.fractionToPoint(this.activeFractionInterval.fractionToPoint(activeFraction));\r\n  }\r\n  // These static variables are reused on calls to integrateFromStartFraction\r\n  private static _gaussFraction: Float64Array;\r\n  private static _gaussWeight: Float64Array;\r\n  private static _gaussMapper: (xA: number, xB: number, arrayX: Float64Array, arrayW: Float64Array) => number;\r\n  /** Initialize class level work arrays. */\r\n  public static initWorkSpace() {\r\n    IntegratedSpiral3d._gaussFraction = new Float64Array(5);\r\n    IntegratedSpiral3d._gaussWeight = new Float64Array(5);\r\n    IntegratedSpiral3d._gaussMapper = (xA, xB, xMapped, wMapped) => Quadrature.setupGauss5(xA, xB, xMapped, wMapped);\r\n  }\r\n  /**\r\n   * Evaluate and sum the gauss quadrature formulas to integrate cos(theta), sin(theta) fractional subset of a reference\r\n   * length (recall that theta is a nonlinear function of the fraction).\r\n   * * This is a single interval of gaussian integration.\r\n   * * The fraction is on the full spiral (not in the mapped active interval).\r\n   * @param xyz advancing integrated point.\r\n   * @param fractionA fraction at start of interval.\r\n   * @param fractionB fraction at end of interval.\r\n   * @param unitArcLength length of curve for 0 to 1 fractional.\r\n   */\r\n  private fullSpiralIncrementalIntegral(xyz: Point3d, fractionA: number, fractionB: number, applyMatrix: boolean) {\r\n    const gaussFraction = IntegratedSpiral3d._gaussFraction;\r\n    const gaussWeight = IntegratedSpiral3d._gaussWeight;\r\n    const numEval = IntegratedSpiral3d._gaussMapper(fractionA, fractionB, gaussFraction, gaussWeight);\r\n    const deltaL = this._arcLength01;\r\n    let w = 0;\r\n    let dx = 0.0;\r\n    let dy = 0.0;\r\n    for (let k = 0; k < numEval; k++) {\r\n      const radians = this.globalFractionToBearingRadians(gaussFraction[k]);\r\n      w = gaussWeight[k] * deltaL;\r\n      dx += w * Math.cos(radians);\r\n      dy += w * Math.sin(radians);\r\n    }\r\n    if (applyMatrix)\r\n      Matrix3d.xyzPlusMatrixTimesXYZ(xyz, this.localToWorld.matrix, { x: dx, y: dy, z: 0.0 }, xyz);\r\n    else\r\n      xyz.addXYZInPlace(dx, dy, 0.0);\r\n  }\r\n  /** Recompute strokes. */\r\n  public override refreshComputedProperties() {\r\n    this._curvature01 = Segment1d.create(\r\n      TransitionSpiral3d.radiusToCurvature(this.radius01.x0),\r\n      TransitionSpiral3d.radiusToCurvature(this.radius01.x1),\r\n    );\r\n    this._globalStrokes.clear();\r\n    const currentPoint = Point3d.create();\r\n    this._globalStrokes.appendStrokePoint(currentPoint);\r\n    const numInterval = 16;\r\n    const fractionStep = 1.0 / numInterval;\r\n    for (let i = 1; i <= numInterval; i++) {\r\n      const fraction0 = (i - 1) * fractionStep;\r\n      const fraction1 = i * fractionStep;\r\n      this.fullSpiralIncrementalIntegral(currentPoint, fraction0, fraction1, false);\r\n      this._globalStrokes.appendStrokePoint(currentPoint);\r\n    }\r\n    this._globalStrokes.tryTransformInPlace(this.localToWorld);\r\n    if (!this.activeFractionInterval.isExact01) {\r\n      if (this._activeStrokes === undefined)\r\n        this._activeStrokes = LineString3d.create();\r\n      this._activeStrokes.clear();\r\n      // The active interval has finer strokes: it's the same fraction step but mapped to a sub-interval.\r\n      // Below assumes fractionToPoint depends upon the global strokes we just computed, and not\r\n      // on the active strokes we are computing.\r\n      for (let i = 0; i <= numInterval; i++) {\r\n        const localFraction = i * fractionStep;\r\n        this._activeStrokes.addPoint(this.fractionToPoint(localFraction));\r\n      }\r\n    }\r\n  }\r\n  /**\r\n   * Create a transition spiral with radius and bearing conditions.\r\n   * @param radius01 radius (inverse curvature) at start and end (radius of zero means straight line).\r\n   * @param bearing01 bearing angles at start and end.  bearings are measured from the x axis, positive clockwise\r\n   * towards y axis.\r\n   * @param activeFractionInterval fractional limits of the active portion of the spiral.\r\n   * @param localToWorld placement frame. Fractional coordinate 0 is at the origin.\r\n   */\r\n  public static createRadiusRadiusBearingBearing(\r\n    radius01: Segment1d,\r\n    bearing01: AngleSweep,\r\n    activeFractionInterval: Segment1d,\r\n    localToWorld: Transform,\r\n    typeName?: string,\r\n  ): IntegratedSpiral3d | undefined {\r\n    const arcLength = TransitionSpiral3d.radiusRadiusSweepRadiansToArcLength(\r\n      radius01.x0, radius01.x1, bearing01.sweepRadians,\r\n    );\r\n    if (typeName === undefined)\r\n      typeName = \"clothoid\";\r\n    const evaluator = NormalizedTransition.findEvaluator(typeName);\r\n    if (!evaluator)\r\n      return undefined;\r\n    return new IntegratedSpiral3d(\r\n      typeName,\r\n      evaluator,\r\n      radius01.clone(),\r\n      bearing01.clone(),\r\n      activeFractionInterval.clone(),\r\n      localToWorld.clone(),\r\n      arcLength,\r\n      new TransitionConditionalProperties(\r\n        radius01.x0, radius01.x1, bearing01.startAngle.clone(), bearing01.endAngle.clone(), undefined,\r\n      ),\r\n    );\r\n  }\r\n  /**\r\n   * Create a transition spiral.\r\n   * * Inputs must provide exactly 4 of the 5 values `[radius0,radius1,bearing0,bearing1,length]`.\r\n   * @param spiralType one of \"clothoid\", \"bloss\", \"biquadratic\", \"cosine\", \"sine\".  If undefined, \"clothoid\" is used.\r\n   * @param radius0 radius (or 0 for tangent to line) at start.\r\n   * @param radius1 radius (or 0 for tangent to line) at end.\r\n   * @param bearing0 bearing, measured CCW from x axis at start.\r\n   * @param bearing1 bearing, measured CCW from x axis at end.\r\n   * @param fractionInterval optional fractional interval for an \"active\" portion of the curve. If omitted, the full\r\n   * [0,1] is used.\r\n   * @param localToWorld placement transform.\r\n   */\r\n  public static createFrom4OutOf5(\r\n    spiralType: string | undefined,\r\n    radius0: number | undefined,\r\n    radius1: number | undefined,\r\n    bearing0: Angle | undefined,\r\n    bearing1: Angle | undefined,\r\n    arcLength: number | undefined,\r\n    fractionInterval: undefined | Segment1d,\r\n    localToWorld: Transform,\r\n  ): IntegratedSpiral3d | undefined {\r\n    if (spiralType === undefined)\r\n      spiralType = \"clothoid\";\r\n    const evaluator = NormalizedTransition.findEvaluator(spiralType);\r\n    if (!evaluator)\r\n      return undefined;\r\n    const data = new TransitionConditionalProperties(radius0, radius1, bearing0, bearing1, arcLength);\r\n    const data1 = data.clone();\r\n    if (!data.tryResolveAnySingleUnknown())\r\n      return undefined;\r\n    if (fractionInterval === undefined)\r\n      fractionInterval = Segment1d.create(0, 1);\r\n    return new IntegratedSpiral3d(\r\n      spiralType,\r\n      evaluator,\r\n      Segment1d.create(data.radius0, data.radius1),\r\n      AngleSweep.createStartEnd(data.bearing0!, data.bearing1!),\r\n      fractionInterval ? fractionInterval.clone() : Segment1d.create(0, 1),\r\n      localToWorld,\r\n      data.curveLength!,\r\n      data1,\r\n    );\r\n  }\r\n  /** Copy all defining data from another spiral. */\r\n  public setFrom(other: IntegratedSpiral3d): IntegratedSpiral3d {\r\n    this.localToWorld.setFrom(other.localToWorld);\r\n    this.radius01.setFrom(other.radius01);\r\n    this._curvature01.setFrom(other._curvature01);\r\n    this.bearing01.setFrom(other.bearing01);\r\n    this.localToWorld.setFrom(other.localToWorld);\r\n    this.activeFractionInterval.setFrom(other.activeFractionInterval);\r\n    this._arcLength01 = other._arcLength01;\r\n    return this;\r\n  }\r\n  /** Deep clone of this spiral. */\r\n  public override clone(): IntegratedSpiral3d {\r\n    return new IntegratedSpiral3d(\r\n      this._spiralType,\r\n      this._evaluator,\r\n      this.radius01.clone(),\r\n      this.bearing01.clone(),\r\n      this.activeFractionInterval.clone(),\r\n      this.localToWorld.clone(),\r\n      this._arcLength01,\r\n      this._designProperties?.clone(),\r\n    );\r\n  }\r\n  /** Apply `transform` to this spiral's local to world transform. */\r\n  public tryTransformInPlace(transformA: Transform): boolean {\r\n    const rigidData = this.applyRigidPartOfTransform(transformA);\r\n    if (rigidData !== undefined) {\r\n      this._curvature01.x0 /= rigidData.scale;\r\n      this._curvature01.x1 /= rigidData.scale;\r\n      this.radius01.x0 *= rigidData.scale;\r\n      this.radius01.x1 *= rigidData.scale;\r\n      this._arcLength01 *= rigidData.scale;\r\n    }\r\n    this.refreshComputedProperties();\r\n    return true;\r\n  }\r\n  /** Return the spiral start point. */\r\n  public override startPoint(result?: Point3d): Point3d {\r\n    return this.activeStrokes.startPoint(result);\r\n  }\r\n  /** Return the spiral end point. */\r\n  public override endPoint(result?: Point3d): Point3d {\r\n    return this.activeStrokes.endPoint(result);\r\n  }\r\n  /** Test if the local to world transform places the spiral xy plane into `plane`. */\r\n  public isInPlane(plane: Plane3dByOriginAndUnitNormal): boolean {\r\n    return plane.isPointInPlane(this.localToWorld.origin as Point3d)\r\n      && Geometry.isSameCoordinate(0.0, this.localToWorld.matrix.dotColumnX(plane.getNormalRef()))\r\n      && Geometry.isSameCoordinate(0.0, this.localToWorld.matrix.dotColumnY(plane.getNormalRef()));\r\n  }\r\n  /**\r\n   * Return length of the spiral.\r\n   * * Because TransitionSpiral is parameterized directly in terms of distance along, this is a simple return value.\r\n   */\r\n  public quickLength() {\r\n    return this.curveLength();\r\n  }\r\n  /**\r\n   * Return length of the spiral.\r\n   * * Because TransitionSpiral is parameterized directly in terms of distance along, this is a simple return value.\r\n   */\r\n  public override curveLength() {\r\n    return this._arcLength01 * (this._activeFractionInterval.absoluteDelta());\r\n  }\r\n  /**\r\n   * Return (unsigned) length of the spiral between fractions.\r\n   * * Because TransitionSpiral is parameterized directly in terms of distance along, this is a simple return value.\r\n   */\r\n  public override curveLengthBetweenFractions(fraction0: number, fraction1: number) {\r\n    return this._arcLength01 * (this._activeFractionInterval.absoluteDelta() * Math.abs(fraction1 - fraction0));\r\n  }\r\n  /** Test if `other` is an instance of `TransitionSpiral3d`. */\r\n  public isSameGeometryClass(other: any): boolean {\r\n    return other instanceof TransitionSpiral3d;\r\n  }\r\n  /**\r\n   * Add strokes from this spiral to `dest`.\r\n   * * Linestrings will usually stroke as just their points.\r\n   * * If maxEdgeLength is given, this will sub-stroke within the linestring -- not what we want.\r\n   */\r\n  public emitStrokes(dest: LineString3d, options?: StrokeOptions): void {\r\n    this.activeStrokes.emitStrokes(dest, options);\r\n  }\r\n  /** Emit stroke fragments to `dest` handler. */\r\n  public emitStrokableParts(dest: IStrokeHandler, options?: StrokeOptions): void {\r\n    const n = this.computeStrokeCountForOptions(options);\r\n    dest.startParentCurvePrimitive(this);\r\n    const activeStrokes = this.activeStrokes;\r\n    const preferPrimary = dest.needPrimaryGeometryForStrokes === undefined ? false : dest.needPrimaryGeometryForStrokes();\r\n    if (!preferPrimary && n <= activeStrokes.numPoints()) {\r\n      this.activeStrokes.emitStrokableParts(dest, options);\r\n    } else {\r\n      dest.announceIntervalForUniformStepStrokes(this, n, 0.0, 1.0);\r\n    }\r\n    dest.endParentCurvePrimitive(this);\r\n  }\r\n  /**\r\n   * Return the stroke count required for given options.\r\n   * @param options StrokeOptions that determine count.\r\n   */\r\n  public computeStrokeCountForOptions(options?: StrokeOptions): number {\r\n    let numStroke;\r\n    if (options) {\r\n      const rMin = Math.min(Math.abs(this.radius01.x0), Math.abs(this.radius01.x1));\r\n      numStroke = options.applyTolerancesToArc(rMin, this.bearing01.sweepRadians);\r\n      numStroke = options.applyMaxEdgeLength(numStroke, this.curveLength());\r\n      numStroke = options.applyMinStrokesPerPrimitive(numStroke);\r\n    } else {\r\n      numStroke = StrokeOptions.applyAngleTol(undefined, 4, this.bearing01.sweepRadians);\r\n    }\r\n    return numStroke;\r\n  }\r\n  /**\r\n   * Reverse the active interval and active strokes.\r\n   * * Primary defining data remains unchanged.\r\n   */\r\n  public reverseInPlace(): void {\r\n    this.activeFractionInterval.reverseInPlace();\r\n    if (this._activeStrokes === undefined)\r\n      this._activeStrokes = this._globalStrokes.clone();\r\n    this._activeStrokes.reverseInPlace();\r\n  }\r\n  /** Evaluate curve point with respect to fraction. */\r\n  public fractionToPoint(activeFraction: number, result?: Point3d): Point3d {\r\n    const targetGlobalFraction = this.activeFractionInterval.fractionToPoint(activeFraction);\r\n    const numStrokes = this._globalStrokes.packedPoints.length - 1;\r\n    if (activeFraction > 1.0) {\r\n      result = this._globalStrokes.packedPoints.back(result)!;\r\n      const integrationStep = 1.0 / numStrokes;\r\n      let currentGlobalFraction = 1.0;\r\n      let nextGlobalFraction = currentGlobalFraction + integrationStep;\r\n      while (nextGlobalFraction < targetGlobalFraction) {\r\n        this.fullSpiralIncrementalIntegral(result, currentGlobalFraction, nextGlobalFraction, true);\r\n        currentGlobalFraction = nextGlobalFraction;\r\n        nextGlobalFraction += integrationStep;\r\n      }\r\n      this.fullSpiralIncrementalIntegral(result, currentGlobalFraction, targetGlobalFraction, true);\r\n    } else if (activeFraction < 0.0) {\r\n      result = this._globalStrokes.packedPoints.front(result)!;\r\n      const integrationStep = 1.0 / numStrokes;\r\n      let currentGlobalFraction = 0.0;\r\n      let nextGlobalFraction = currentGlobalFraction - integrationStep;\r\n      while (nextGlobalFraction > targetGlobalFraction) {\r\n        this.fullSpiralIncrementalIntegral(result, currentGlobalFraction, nextGlobalFraction, true);\r\n        currentGlobalFraction = nextGlobalFraction;\r\n        nextGlobalFraction -= integrationStep;\r\n      }\r\n      this.fullSpiralIncrementalIntegral(result, currentGlobalFraction, targetGlobalFraction, true);\r\n    } else {\r\n      const clampedGlobalFraction = Geometry.clampToStartEnd(targetGlobalFraction, 0, 1);\r\n      const index0 = Math.trunc(clampedGlobalFraction * numStrokes); // this indexes the point to the left of the query\r\n      const globalFraction0 = index0 / numStrokes;\r\n      result = this._globalStrokes.packedPoints.getPoint3dAtUncheckedPointIndex(index0, result);\r\n      // GeometryCoreTestIO.consoleLog(\"fractionToPoint \", activeFraction, this.activeFractionInterval, \"( global integration \" + globalFraction0 + \" to \" + globalFraction + \")\", index0);\r\n      this.fullSpiralIncrementalIntegral(result, globalFraction0, targetGlobalFraction, true);\r\n\r\n    }\r\n    return result;\r\n  }\r\n  /** Evaluate curve point and derivative with respect to fraction. */\r\n  public fractionToPointAndDerivative(activeFraction: number, result?: Ray3d): Ray3d {\r\n    const globalFraction = this.activeFractionInterval.fractionToPoint(activeFraction);\r\n    result = result ? result : Ray3d.createZero();\r\n    this.fractionToPoint(activeFraction, result.origin);\r\n    const radians = this.globalFractionToBearingRadians(globalFraction);\r\n    const a = this._arcLength01 * this.activeFractionInterval.signedDelta();\r\n    this.localToWorld.matrix.multiplyXY(a * Math.cos(radians), a * Math.sin(radians), result.direction);\r\n    return result;\r\n  }\r\n  /** Return the frenet frame at fractional position. */\r\n  public override fractionToFrenetFrame(activeFraction: number, result?: Transform): Transform {\r\n    const globalFraction = this.activeFractionInterval.fractionToPoint(activeFraction);\r\n    result = result ? result : Transform.createIdentity();\r\n    result.origin.setFrom(this.fractionToPoint(activeFraction));\r\n    Matrix3d.createRigidFromMatrix3d(this.localToWorld.matrix, AxisOrder.XYZ, result.matrix);\r\n    const radians = this.globalFractionToBearingRadians(globalFraction);\r\n    const c = Math.cos(radians);\r\n    const s = Math.sin(radians);\r\n    result.matrix.applyGivensColumnOp(0, 1, c, s);\r\n    return result;\r\n  }\r\n  /**\r\n   * Return a plane with\r\n   * * origin at fractional position along the curve.\r\n   * * vectorU is the first derivative, i.e. tangent vector with length equal to the rate of change with respect to\r\n   * the fraction.\r\n   * * vectorV is the second derivative, i.e.derivative of vectorU.\r\n   */\r\n  public fractionToPointAnd2Derivatives(\r\n    activeFraction: number, result?: Plane3dByOriginAndVectors,\r\n  ): Plane3dByOriginAndVectors | undefined {\r\n    const globalFraction = this.activeFractionInterval.fractionToPoint(activeFraction);\r\n    const origin = this.fractionToPoint(activeFraction);\r\n    const radians = this.globalFractionToBearingRadians(globalFraction);\r\n    const c = Math.cos(radians);\r\n    const s = Math.sin(radians);\r\n    const delta = this.activeFractionInterval.signedDelta();\r\n    const a = delta;\r\n    const b = a * delta;\r\n    const vectorX = this.localToWorld.matrix.multiplyXY(a * c, a * s);\r\n    const vectorY = this.localToWorld.matrix.multiplyXY(-b * s, b * c);\r\n    vectorY.scaleInPlace(this.globalFractionToCurvature(globalFraction));\r\n    return Plane3dByOriginAndVectors.createCapture(origin, vectorX, vectorY, result);\r\n  }\r\n  /** Second step of double dispatch:  call `handler.handleTransitionSpiral(this)`. */\r\n  public dispatchToGeometryHandler(handler: GeometryHandler): any {\r\n    return handler.handleTransitionSpiral(this);\r\n  }\r\n  /** Compare various coordinate quantities. */\r\n  public override isAlmostEqual(other?: GeometryQuery): boolean {\r\n    if (other instanceof IntegratedSpiral3d) {\r\n      return this.radius01.isAlmostEqual(other.radius01)\r\n        && this.bearing01.isAlmostEqualAllowPeriodShift(other.bearing01)\r\n        && this.localToWorld.isAlmostEqual(other.localToWorld)\r\n        && Geometry.isSameCoordinate(this._arcLength01, other._arcLength01)\r\n        && this.activeFractionInterval.isAlmostEqual(other.activeFractionInterval)\r\n        && this._curvature01.isAlmostEqual(other._curvature01);\r\n    }\r\n    return false;\r\n  }\r\n}\r\n// at load time, initialize gauss quadrature workspace\r\nIntegratedSpiral3d.initWorkSpace();\r\n"]}

package/lib/esm/curve/spiral/MXCubicAlongArcSpiralEvaluator.d.ts

@@ -1,17 +1,12 @@
 import { CubicEvaluator } from "./CubicEvaluator";
 /**
  * MX Cubic along arc.
- * This is y= m*x^3 with
- * * x any point on the x axis
- * * `fraction` along the spiral goes to `x = fraction * Lx`
- * * m is (1/6RL)
+ * This is y = m*x^3 with
+ * * x any point on the x axis.
+ * * `fraction` along the spiral goes to `x = fraction * Lx`.
+ * * m is (1/6RL).
  * * construction length L is nominal along the curve.
  * * x length Lx is along the axis, determined by two terms of the clothoid x series.
- * *
- * @param localToWorld
- * @param nominalL1
- * @param nominalR1
- * @param activeInterval
  * @internal
  */
 export declare class MXCubicAlongArcEvaluator extends CubicEvaluator {
@@ -23,14 +18,15 @@ export declare class MXCubicAlongArcEvaluator extends CubicEvaluator {
     static computeCubicM(length1: number, radius1: number): number | undefined;
     static create(length1: number, radius1: number): MXCubicAlongArcEvaluator | undefined;
     scaleInPlace(scaleFactor: number): void;
-    /** return a deep copy of the evaluator */
+    /** Return a deep copy of the evaluator. */
     clone(): MXCubicAlongArcEvaluator;
-    /** Member by member matchup ... */
+    /** Member by member matchup. */
     isAlmostEqual(other: any): boolean;
     /**
      * Return a (fast but mediocre) approximation of spiral x position as function of approximate distance along the curve.
      * * This x-to-distance relation is not as precise as the CurvePrimitive method moveSignedDistanceFromFraction.
-     * * It is supported here for users interested in replicating the Czech distance mapping rather than the more accurate CurvePrimitive measurements.
+     * * It is supported here for users interested in replicating the Czech distance mapping rather than the more accurate
+     * CurvePrimitive measurements.
      * @param x distance along the x axis.
      */
     static approximateDistanceAlongToX(nominalLength1: number, nominalRadius1: number, nominalDistanceAlong: number): number;

package/lib/esm/curve/spiral/MXCubicAlongArcSpiralEvaluator.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"MXCubicAlongArcSpiralEvaluator.d.ts","sourceRoot":"","sources":["../../../../src/curve/spiral/MXCubicAlongArcSpiralEvaluator.ts"],"names":[],"mappings":"AASA,OAAO,EAAE,cAAc,EAAE,MAAM,kBAAkB,CAAC;AAClD;;;;;;;;;;;;;;GAcG;AACH,qBAAa,wBAAyB,SAAQ,cAAc;IACnD,cAAc,EAAE,MAAM,CAAC;IACvB,cAAc,EAAE,MAAM,CAAC;IAE9B,uEAAuE;IACvE,OAAO;IAMP,kCAAkC;WACpB,aAAa,CAAC,OAAO,EAAE,MAAM,EAAE,OAAO,EAAE,MAAM,GAAG,MAAM,GAAG,SAAS;WAKnE,MAAM,CAAC,OAAO,EAAE,MAAM,EAAE,OAAO,EAAE,MAAM,GAAG,wBAAwB,GAAG,SAAS;IAQ5E,YAAY,CAAC,WAAW,EAAE,MAAM;IAKhD,0CAA0C;IACnC,KAAK,IAAI,wBAAwB;IACxC,mCAAmC;IAC5B,aAAa,CAAC,KAAK,EAAE,GAAG,GAAG,OAAO;IAOzC;;;;;OAKG;WACW,2BAA2B,CAAC,cAAc,EAAE,MAAM,EAAE,cAAc,EAAE,MAAM,EAAE,oBAAoB,EAAE,MAAM,GAAG,MAAM;CAMhI"}
+{"version":3,"file":"MXCubicAlongArcSpiralEvaluator.d.ts","sourceRoot":"","sources":["../../../../src/curve/spiral/MXCubicAlongArcSpiralEvaluator.ts"],"names":[],"mappings":"AASA,OAAO,EAAE,cAAc,EAAE,MAAM,kBAAkB,CAAC;AAElD;;;;;;;;;GASG;AACH,qBAAa,wBAAyB,SAAQ,cAAc;IACnD,cAAc,EAAE,MAAM,CAAC;IACvB,cAAc,EAAE,MAAM,CAAC;IAC9B,uEAAuE;IACvE,OAAO;IAKP,kCAAkC;WACpB,aAAa,CAAC,OAAO,EAAE,MAAM,EAAE,OAAO,EAAE,MAAM,GAAG,MAAM,GAAG,SAAS;WAInE,MAAM,CAAC,OAAO,EAAE,MAAM,EAAE,OAAO,EAAE,MAAM,GAAG,wBAAwB,GAAG,SAAS;IAO5E,YAAY,CAAC,WAAW,EAAE,MAAM;IAKhD,2CAA2C;IACpC,KAAK,IAAI,wBAAwB;IAGxC,gCAAgC;IACzB,aAAa,CAAC,KAAK,EAAE,GAAG,GAAG,OAAO;IAOzC;;;;;;OAMG;WACW,2BAA2B,CACvC,cAAc,EAAE,MAAM,EAAE,cAAc,EAAE,MAAM,EAAE,oBAAoB,EAAE,MAAM,GAC3E,MAAM;CAMV"}

package/lib/esm/curve/spiral/MXCubicAlongArcSpiralEvaluator.js

@@ -9,17 +9,12 @@ import { Geometry } from "../../Geometry";
 import { CubicEvaluator } from "./CubicEvaluator";
 /**
  * MX Cubic along arc.
- * This is y= m*x^3 with
- * * x any point on the x axis
- * * `fraction` along the spiral goes to `x = fraction * Lx`
- * * m is (1/6RL)
+ * This is y = m*x^3 with
+ * * x any point on the x axis.
+ * * `fraction` along the spiral goes to `x = fraction * Lx`.
+ * * m is (1/6RL).
  * * construction length L is nominal along the curve.
  * * x length Lx is along the axis, determined by two terms of the clothoid x series.
- * *
- * @param localToWorld
- * @param nominalL1
- * @param nominalR1
- * @param activeInterval
  * @internal
  */
 export class MXCubicAlongArcEvaluator extends CubicEvaluator {
@@ -48,9 +43,11 @@ export class MXCubicAlongArcEvaluator extends CubicEvaluator {
         this.nominalRadius1 *= scaleFactor;
         super.scaleInPlace(scaleFactor);
     }
-    /** return a deep copy of the evaluator */
-    clone() { return new MXCubicAlongArcEvaluator(this.nominalLength1, this.nominalRadius1, this.axisLength, this.cubicM); }
-    /** Member by member matchup ... */
+    /** Return a deep copy of the evaluator. */
+    clone() {
+        return new MXCubicAlongArcEvaluator(this.nominalLength1, this.nominalRadius1, this.axisLength, this.cubicM);
+    }
+    /** Member by member matchup. */
     isAlmostEqual(other) {
         if (other instanceof MXCubicAlongArcEvaluator) {
             return Geometry.isSameCoordinate(this.nominalLength1, other.nominalLength1)
@@ -61,7 +58,8 @@ export class MXCubicAlongArcEvaluator extends CubicEvaluator {
     /**
      * Return a (fast but mediocre) approximation of spiral x position as function of approximate distance along the curve.
      * * This x-to-distance relation is not as precise as the CurvePrimitive method moveSignedDistanceFromFraction.
-     * * It is supported here for users interested in replicating the Czech distance mapping rather than the more accurate CurvePrimitive measurements.
+     * * It is supported here for users interested in replicating the Czech distance mapping rather than the more accurate
+     * CurvePrimitive measurements.
      * @param x distance along the x axis.
      */
     static approximateDistanceAlongToX(nominalLength1, nominalRadius1, nominalDistanceAlong) {