@itwin/core-geometry 5.5.0-dev.4 → 5.5.0-dev.5

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

@@ -5,20 +5,21 @@
 /** @packageDocumentation
  * @module Curve
  */
-import { TransitionSpiral3d } from "./TransitionSpiral3d";
-import { Segment1d } from "../../geometry3d/Segment1d";
+import { assert } from "@itwin/core-bentley";
+import { AxisOrder, Geometry } from "../../Geometry";
 import { AngleSweep } from "../../geometry3d/AngleSweep";
+import { Matrix3d } from "../../geometry3d/Matrix3d";
+import { Plane3dByOriginAndVectors } from "../../geometry3d/Plane3dByOriginAndVectors";
+import { Point3d } from "../../geometry3d/Point3dVector3d";
+import { Ray3d } from "../../geometry3d/Ray3d";
+import { Segment1d } from "../../geometry3d/Segment1d";
 import { Transform } from "../../geometry3d/Transform";
+import { Quadrature } from "../../numerics/Quadrature";
 import { LineString3d } from "../LineString3d";
+import { StrokeOptions } from "../StrokeOptions";
 import { NormalizedTransition } from "./NormalizedTransition";
 import { TransitionConditionalProperties } from "./TransitionConditionalProperties";
-import { Quadrature } from "../../numerics/Quadrature";
-import { Point3d } from "../../geometry3d/Point3dVector3d";
-import { Matrix3d } from "../../geometry3d/Matrix3d";
-import { AxisOrder, Geometry } from "../../Geometry";
-import { StrokeOptions } from "../StrokeOptions";
-import { Ray3d } from "../../geometry3d/Ray3d";
-import { Plane3dByOriginAndVectors } from "../../geometry3d/Plane3dByOriginAndVectors";
+import { TransitionSpiral3d } from "./TransitionSpiral3d";
 /**
  * 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).
@@ -90,12 +91,8 @@ export class IntegratedSpiral3d extends TransitionSpiral3d {
         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.
-     * * The `undefined` result is to match the abstract class -- it cannot actually occur.
-     */
+    /** Return the curvature at given fraction of the active interval. */
     fractionToCurvature(activeFraction) {
-        // BUG? active interval
         return this._curvature01.fractionToPoint(this.activeFractionInterval.fractionToPoint(activeFraction));
     }
     // These static variables are reused on calls to integrateFromStartFraction
@@ -329,8 +326,11 @@ export class IntegratedSpiral3d extends TransitionSpiral3d {
     fractionToPoint(activeFraction, result) {
         const targetGlobalFraction = this.activeFractionInterval.fractionToPoint(activeFraction);
         const numStrokes = this._globalStrokes.packedPoints.length - 1;
+        assert(numStrokes > 0, "expect constructor to populate nonzero global strokes");
+        if (numStrokes < 1)
+            return Point3d.createZero(result); // trust but verify
         if (activeFraction > 1.0) {
-            result = this._globalStrokes.packedPoints.back(result);
+            result = this._globalStrokes.packedPoints.getPoint3dAtUncheckedPointIndex(numStrokes, result);
             const integrationStep = 1.0 / numStrokes;
             let currentGlobalFraction = 1.0;
             let nextGlobalFraction = currentGlobalFraction + integrationStep;
@@ -342,7 +342,7 @@ export class IntegratedSpiral3d extends TransitionSpiral3d {
             this.fullSpiralIncrementalIntegral(result, currentGlobalFraction, targetGlobalFraction, true);
         }
         else if (activeFraction < 0.0) {
-            result = this._globalStrokes.packedPoints.front(result);
+            result = this._globalStrokes.packedPoints.getPoint3dAtUncheckedPointIndex(0, result);
             const integrationStep = 1.0 / numStrokes;
             let currentGlobalFraction = 0.0;
             let nextGlobalFraction = currentGlobalFraction - integrationStep;
@@ -355,10 +355,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);
+            assert(0 <= index0 && index0 <= numStrokes, "expect valid index0, at left of containing stroke interval");
             const globalFraction0 = index0 / numStrokes;
             result = this._globalStrokes.packedPoints.getPoint3dAtUncheckedPointIndex(index0, result);
-            // GeometryCoreTestIO.consoleLog("fractionToPoint ", activeFraction, this.activeFractionInterval, "( global integration " + globalFraction0 + " to " + globalFraction + ")", index0);
             this.fullSpiralIncrementalIntegral(result, globalFraction0, targetGlobalFraction, true);
         }
         return result;

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;AAGvF;;;;;;;GAOG;AACH,2CAA2C;AAC3C,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,uHAAuH;QACvH,2CAA2C;QAC3C,MAAM,YAAY,GAAG,IAAI,CAAC,UAAU,CAAC,cAAc,CAAC,QAAQ,CAAC,CAAC;QAC9D,MAAM,SAAS,GAAG,IAAI,CAAC,YAAY,CAAC;QACpC,OAAO,IAAI,CAAC,SAAS,CAAC,YAAY;cAC9B,QAAQ,GAAG,SAAS,GAAG,IAAI,CAAC,YAAY,CAAC,EAAE;cAC3C,YAAY,GAAG,SAAS,GAAG,IAAI,CAAC,YAAY,CAAC,WAAW,EAAE,CAAC;IACjE,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,2CAA2C;QAC3C,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\n// see internaldocs/Spiral.md for more info\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    // calculate area under curvature curve from 0 to fraction and add it to start angle to get angle at the given fraction\r\n    // see internaldocs/Spiral.md for more info\r\n    const areaFraction = this._evaluator.fractionToArea(fraction);\r\n    const arcLength = this._arcLength01;\r\n    return this.bearing01.startRadians\r\n      + fraction * arcLength * this._curvature01.x0\r\n      + areaFraction * arcLength * this._curvature01.signedDelta();\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 applyMatrix if true, apply the localToWorld matrix to the computed delta before adding to xyz.\r\n   */\r\n  private fullSpiralIncrementalIntegral(xyz: Point3d, fractionA: number, fractionB: number, applyMatrix: boolean): void {\r\n    // see internaldocs/Spiral.md for more info\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"]}
+{"version":3,"file":"IntegratedSpiral3d.js","sourceRoot":"","sources":["../../../../src/curve/spiral/IntegratedSpiral3d.ts"],"names":[],"mappings":"AAAA;;;+FAG+F;AAC/F;;GAEG;AAEH,OAAO,EAAE,MAAM,EAAE,MAAM,qBAAqB,CAAC;AAC7C,OAAO,EAAE,SAAS,EAAE,QAAQ,EAAE,MAAM,gBAAgB,CAAC;AAErD,OAAO,EAAE,UAAU,EAAE,MAAM,6BAA6B,CAAC;AAEzD,OAAO,EAAE,QAAQ,EAAE,MAAM,2BAA2B,CAAC;AAErD,OAAO,EAAE,yBAAyB,EAAE,MAAM,4CAA4C,CAAC;AACvF,OAAO,EAAE,OAAO,EAAE,MAAM,kCAAkC,CAAC;AAC3D,OAAO,EAAE,KAAK,EAAE,MAAM,wBAAwB,CAAC;AAC/C,OAAO,EAAE,SAAS,EAAE,MAAM,4BAA4B,CAAC;AACvD,OAAO,EAAE,SAAS,EAAE,MAAM,4BAA4B,CAAC;AACvD,OAAO,EAAE,UAAU,EAAE,MAAM,2BAA2B,CAAC;AAEvD,OAAO,EAAE,YAAY,EAAE,MAAM,iBAAiB,CAAC;AAC/C,OAAO,EAAE,aAAa,EAAE,MAAM,kBAAkB,CAAC;AACjD,OAAO,EAAE,oBAAoB,EAAE,MAAM,wBAAwB,CAAC;AAC9D,OAAO,EAAE,+BAA+B,EAAE,MAAM,mCAAmC,CAAC;AACpF,OAAO,EAAE,kBAAkB,EAAE,MAAM,sBAAsB,CAAC;AAE1D;;;;;;;GAOG;AACH,2CAA2C;AAC3C,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,uHAAuH;QACvH,2CAA2C;QAC3C,MAAM,YAAY,GAAG,IAAI,CAAC,UAAU,CAAC,cAAc,CAAC,QAAQ,CAAC,CAAC;QAC9D,MAAM,SAAS,GAAG,IAAI,CAAC,YAAY,CAAC;QACpC,OAAO,IAAI,CAAC,SAAS,CAAC,YAAY;cAC9B,QAAQ,GAAG,SAAS,GAAG,IAAI,CAAC,YAAY,CAAC,EAAE;cAC3C,YAAY,GAAG,SAAS,GAAG,IAAI,CAAC,YAAY,CAAC,WAAW,EAAE,CAAC;IACjE,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,qEAAqE;IACrD,mBAAmB,CAAC,cAAsB;QACxD,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,2CAA2C;QAC3C,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,QAAQ,EAAE,IAAI,CAAC,QAAQ,CAAC,EACvD,gBAAgB,CAAC,CAAC,CAAC,gBAAgB,CAAC,KAAK,EAAE,CAAC,CAAC,CAAC,SAAS,CAAC,MAAM,CAAC,CAAC,EAAE,CAAC,CAAC,EACpE,YAAY,EAAE,IAAI,CAAC,WAAW,EAAE,KAAK,CAAC,CAAC;IAC3C,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,MAAM,CAAC,UAAU,GAAG,CAAC,EAAE,uDAAuD,CAAC,CAAC;QAChF,IAAI,UAAU,GAAG,CAAC;YAChB,OAAO,OAAO,CAAC,UAAU,CAAC,MAAM,CAAC,CAAC,CAAC,mBAAmB;QACxD,IAAI,cAAc,GAAG,GAAG,EAAE,CAAC;YACzB,MAAM,GAAG,IAAI,CAAC,cAAc,CAAC,YAAY,CAAC,+BAA+B,CAAC,UAAU,EAAE,MAAM,CAAC,CAAC;YAC9F,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,+BAA+B,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC;YACrF,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;YAC9D,MAAM,CAAC,CAAC,IAAI,MAAM,IAAI,MAAM,IAAI,UAAU,EAAE,4DAA4D,CAAC,CAAC;YAC1G,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,IAAI,CAAC,6BAA6B,CAAC,MAAM,EAAE,eAAe,EAAE,oBAAoB,EAAE,IAAI,CAAC,CAAC;QAC1F,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 { assert } from \"@itwin/core-bentley\";\r\nimport { AxisOrder, Geometry } from \"../../Geometry\";\r\nimport { Angle } from \"../../geometry3d/Angle\";\r\nimport { AngleSweep } from \"../../geometry3d/AngleSweep\";\r\nimport { GeometryHandler, IStrokeHandler } from \"../../geometry3d/GeometryHandler\";\r\nimport { Matrix3d } from \"../../geometry3d/Matrix3d\";\r\nimport { Plane3dByOriginAndUnitNormal } from \"../../geometry3d/Plane3dByOriginAndUnitNormal\";\r\nimport { Plane3dByOriginAndVectors } from \"../../geometry3d/Plane3dByOriginAndVectors\";\r\nimport { Point3d } from \"../../geometry3d/Point3dVector3d\";\r\nimport { Ray3d } from \"../../geometry3d/Ray3d\";\r\nimport { Segment1d } from \"../../geometry3d/Segment1d\";\r\nimport { Transform } from \"../../geometry3d/Transform\";\r\nimport { Quadrature } from \"../../numerics/Quadrature\";\r\nimport { GeometryQuery } from \"../GeometryQuery\";\r\nimport { LineString3d } from \"../LineString3d\";\r\nimport { StrokeOptions } from \"../StrokeOptions\";\r\nimport { NormalizedTransition } from \"./NormalizedTransition\";\r\nimport { TransitionConditionalProperties } from \"./TransitionConditionalProperties\";\r\nimport { TransitionSpiral3d } from \"./TransitionSpiral3d\";\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\n// see internaldocs/Spiral.md for more info\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    // calculate area under curvature curve from 0 to fraction and add it to start angle to get angle at the given fraction\r\n    // see internaldocs/Spiral.md for more info\r\n    const areaFraction = this._evaluator.fractionToArea(fraction);\r\n    const arcLength = this._arcLength01;\r\n    return this.bearing01.startRadians\r\n      + fraction * arcLength * this._curvature01.x0\r\n      + areaFraction * arcLength * this._curvature01.signedDelta();\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  /** Return the curvature at given fraction of the active interval. */\r\n  public override fractionToCurvature(activeFraction: number): number {\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 applyMatrix if true, apply the localToWorld matrix to the computed delta before adding to xyz.\r\n   */\r\n  private fullSpiralIncrementalIntegral(xyz: Point3d, fractionA: number, fractionB: number, applyMatrix: boolean): void {\r\n    // see internaldocs/Spiral.md for more info\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, data.curveLength, data1);\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    assert(numStrokes > 0, \"expect constructor to populate nonzero global strokes\");\r\n    if (numStrokes < 1)\r\n      return Point3d.createZero(result); // trust but verify\r\n    if (activeFraction > 1.0) {\r\n      result = this._globalStrokes.packedPoints.getPoint3dAtUncheckedPointIndex(numStrokes, 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.getPoint3dAtUncheckedPointIndex(0, 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);\r\n      assert(0 <= index0 && index0 <= numStrokes, \"expect valid index0, at left of containing stroke interval\");\r\n      const globalFraction0 = index0 / numStrokes;\r\n      result = this._globalStrokes.packedPoints.getPoint3dAtUncheckedPointIndex(index0, result);\r\n      this.fullSpiralIncrementalIntegral(result, globalFraction0, targetGlobalFraction, true);\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/TransitionConditionalProperties.d.ts

@@ -54,7 +54,13 @@ export declare class TransitionConditionalProperties {
      * Examine which properties are defined and compute the (single) undefined.
      * @returns true if the input state had precisely one undefined member.
      */
-    tryResolveAnySingleUnknown(): boolean;
+    tryResolveAnySingleUnknown(): this is {
+        radius0: number;
+        radius1: number;
+        bearing0: Angle;
+        bearing1: Angle;
+        curveLength: number;
+    };
     private almostEqualCoordinate;
     private almostEqualBearing;
     /** Test if `this` and `other` have matching numeric and undefined members. */

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

@@ -1 +1 @@
-{"version":3,"file":"TransitionConditionalProperties.d.ts","sourceRoot":"","sources":["../../../../src/curve/spiral/TransitionConditionalProperties.ts"],"names":[],"mappings":"AASA,OAAO,EAAE,KAAK,EAAE,MAAM,wBAAwB,CAAC;AAG/C;;;;;;;;;;;;GAYG;AAEH,qBAAa,+BAA+B;IAC1C;;;OAGG;IACI,OAAO,EAAE,MAAM,GAAG,SAAS,CAAC;IACnC;;;OAGG;IACI,OAAO,EAAE,MAAM,GAAG,SAAS,CAAC;IACnC;;;OAGG;IACI,QAAQ,EAAE,KAAK,GAAG,SAAS,CAAC;IACnC;;;OAGG;IACI,QAAQ,EAAE,KAAK,GAAG,SAAS,CAAC;IACnC,oBAAoB;IACb,WAAW,EAAE,MAAM,GAAG,SAAS,CAAC;IACvC;;;;;;;OAOG;gBAED,OAAO,EAAE,MAAM,GAAG,SAAS,EAC3B,OAAO,EAAE,MAAM,GAAG,SAAS,EAC3B,QAAQ,EAAE,KAAK,GAAG,SAAS,EAC3B,QAAQ,EAAE,KAAK,GAAG,SAAS,EAC3B,SAAS,EAAE,MAAM,GAAG,SAAS;IAQ/B,kEAAkE;IAC3D,oBAAoB;IAO3B,qEAAqE;IAC9D,KAAK,IAAI,+BAA+B;IAS/C,0EAA0E;IACnE,qBAAqB,IAAI,OAAO;IAUvC;;;OAGG;IACI,0BAA0B,IAAI,OAAO;IAsC5C,OAAO,CAAC,qBAAqB;IAO7B,OAAO,CAAC,kBAAkB;IAO1B,8EAA8E;IACvE,aAAa,CAAC,KAAK,CAAC,EAAE,+BAA+B,GAAG,OAAO;IAetE,iFAAiF;IAC1E,gBAAgB,CAAC,CAAC,EAAE,MAAM,GAAG,IAAI;IAUxC,uEAAuE;WACzD,cAAc,CAC1B,CAAC,EAAE,+BAA+B,GAAG,SAAS,EAAE,CAAC,EAAE,+BAA+B,GAAG,SAAS,GAC7F,OAAO;CAKX"}
+{"version":3,"file":"TransitionConditionalProperties.d.ts","sourceRoot":"","sources":["../../../../src/curve/spiral/TransitionConditionalProperties.ts"],"names":[],"mappings":"AASA,OAAO,EAAE,KAAK,EAAE,MAAM,wBAAwB,CAAC;AAG/C;;;;;;;;;;;;GAYG;AAEH,qBAAa,+BAA+B;IAC1C;;;OAGG;IACI,OAAO,EAAE,MAAM,GAAG,SAAS,CAAC;IACnC;;;OAGG;IACI,OAAO,EAAE,MAAM,GAAG,SAAS,CAAC;IACnC;;;OAGG;IACI,QAAQ,EAAE,KAAK,GAAG,SAAS,CAAC;IACnC;;;OAGG;IACI,QAAQ,EAAE,KAAK,GAAG,SAAS,CAAC;IACnC,oBAAoB;IACb,WAAW,EAAE,MAAM,GAAG,SAAS,CAAC;IACvC;;;;;;;OAOG;gBAED,OAAO,EAAE,MAAM,GAAG,SAAS,EAC3B,OAAO,EAAE,MAAM,GAAG,SAAS,EAC3B,QAAQ,EAAE,KAAK,GAAG,SAAS,EAC3B,QAAQ,EAAE,KAAK,GAAG,SAAS,EAC3B,SAAS,EAAE,MAAM,GAAG,SAAS;IAQ/B,kEAAkE;IAC3D,oBAAoB;IAO3B,qEAAqE;IAC9D,KAAK,IAAI,+BAA+B;IAS/C,0EAA0E;IACnE,qBAAqB,IAAI,OAAO;IAUvC;;;OAGG;IACI,0BAA0B,IAAI,IAAI,IAAI;QAAC,OAAO,EAAE,MAAM,CAAC;QAAC,OAAO,EAAE,MAAM,CAAC;QAAC,QAAQ,EAAE,KAAK,CAAC;QAAC,QAAQ,EAAE,KAAK,CAAC;QAAC,WAAW,EAAE,MAAM,CAAA;KAAC;IAsCtI,OAAO,CAAC,qBAAqB;IAO7B,OAAO,CAAC,kBAAkB;IAO1B,8EAA8E;IACvE,aAAa,CAAC,KAAK,CAAC,EAAE,+BAA+B,GAAG,OAAO;IAetE,iFAAiF;IAC1E,gBAAgB,CAAC,CAAC,EAAE,MAAM,GAAG,IAAI;IAUxC,uEAAuE;WACzD,cAAc,CAC1B,CAAC,EAAE,+BAA+B,GAAG,SAAS,EAAE,CAAC,EAAE,+BAA+B,GAAG,SAAS,GAC7F,OAAO;CAKX"}

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

@@ -1 +1 @@
-{"version":3,"file":"TransitionConditionalProperties.js","sourceRoot":"","sources":["../../../../src/curve/spiral/TransitionConditionalProperties.ts"],"names":[],"mappings":"AAAA;;;+FAG+F;AAE/F;;GAEG;AACH,OAAO,EAAE,QAAQ,EAAE,MAAM,gBAAgB,CAAC;AAC1C,OAAO,EAAE,KAAK,EAAE,MAAM,wBAAwB,CAAC;AAC/C,OAAO,EAAE,kBAAkB,EAAE,MAAM,sBAAsB,CAAC;AAE1D;;;;;;;;;;;;GAYG;AACH,2CAA2C;AAC3C,MAAM,OAAO,+BAA+B;IAC1C;;;OAGG;IACI,OAAO,CAAqB;IACnC;;;OAGG;IACI,OAAO,CAAqB;IACnC;;;OAGG;IACI,QAAQ,CAAoB;IACnC;;;OAGG;IACI,QAAQ,CAAoB;IACnC,oBAAoB;IACb,WAAW,CAAqB;IACvC;;;;;;;OAOG;IACH,YACE,OAA2B,EAC3B,OAA2B,EAC3B,QAA2B,EAC3B,QAA2B,EAC3B,SAA6B;QAE7B,IAAI,CAAC,OAAO,GAAG,OAAO,CAAC;QACvB,IAAI,CAAC,OAAO,GAAG,OAAO,CAAC;QACvB,IAAI,CAAC,QAAQ,GAAG,QAAQ,CAAC;QACzB,IAAI,CAAC,QAAQ,GAAG,QAAQ,CAAC;QACzB,IAAI,CAAC,WAAW,GAAG,SAAS,CAAC;IAC/B,CAAC;IACD,kEAAkE;IAC3D,oBAAoB;QACzB,OAAO,QAAQ,CAAC,SAAS,CAAC,IAAI,CAAC,OAAO,CAAC;cACnC,QAAQ,CAAC,SAAS,CAAC,IAAI,CAAC,OAAO,CAAC;cAChC,QAAQ,CAAC,SAAS,CAAC,IAAI,CAAC,QAAQ,CAAC;cACjC,QAAQ,CAAC,SAAS,CAAC,IAAI,CAAC,QAAQ,CAAC;cACjC,QAAQ,CAAC,SAAS,CAAC,IAAI,CAAC,WAAW,CAAC,CAAC;IAC3C,CAAC;IACD,qEAAqE;IAC9D,KAAK;QACV,OAAO,IAAI,+BAA+B,CACxC,IAAI,CAAC,OAAO,EACZ,IAAI,CAAC,OAAO,EACZ,IAAI,CAAC,QAAQ,KAAK,SAAS,CAAC,CAAC,CAAC,SAAS,CAAC,CAAC,CAAC,IAAI,CAAC,QAAQ,CAAC,KAAK,EAAE,EAC/D,IAAI,CAAC,QAAQ,KAAK,SAAS,CAAC,CAAC,CAAC,SAAS,CAAC,CAAC,CAAC,IAAI,CAAC,QAAQ,CAAC,KAAK,EAAE,EAC/D,IAAI,CAAC,WAAW,CACjB,CAAC;IACJ,CAAC;IACD,0EAA0E;IACnE,qBAAqB;QAC1B,IAAI,IAAI,CAAC,WAAW,KAAK,SAAS,IAAI,IAAI,CAAC,QAAQ,KAAK,SAAS,IAAI,IAAI,CAAC,QAAQ,KAAK,SAAS;eAC3F,IAAI,CAAC,OAAO,KAAK,SAAS,IAAI,IAAI,CAAC,OAAO,KAAK,SAAS,EAAE,CAAC;YAC9D,MAAM,SAAS,GAAG,kBAAkB,CAAC,mCAAmC,CACtE,IAAI,CAAC,OAAO,EAAE,IAAI,CAAC,OAAO,EAAE,IAAI,CAAC,QAAQ,CAAC,OAAO,GAAG,IAAI,CAAC,QAAQ,CAAC,OAAO,CAC1E,CAAC;YACF,OAAO,QAAQ,CAAC,gBAAgB,CAAC,IAAI,CAAC,WAAW,EAAE,SAAS,CAAC,CAAC;QAChE,CAAC;QACD,OAAO,KAAK,CAAC;IACf,CAAC;IACD;;;OAGG;IACI,0BAA0B;QAC/B,IAAI,IAAI,CAAC,qBAAqB,EAAE;YAC9B,OAAO,IAAI,CAAC;QACd,IAAI,IAAI,CAAC,QAAQ,IAAI,IAAI,CAAC,QAAQ,EAAE,CAAC;YACnC,MAAM,YAAY,GAAG,IAAI,CAAC,QAAQ,CAAC,OAAO,GAAG,IAAI,CAAC,QAAQ,CAAC,OAAO,CAAC;YACnE,IAAI,IAAI,CAAC,WAAW,KAAK,SAAS,IAAI,IAAI,CAAC,OAAO,KAAK,SAAS,IAAI,IAAI,CAAC,OAAO,KAAK,SAAS,EAAE,CAAC;gBAC/F,IAAI,CAAC,WAAW,GAAG,kBAAkB,CAAC,mCAAmC,CAAC,IAAI,CAAC,OAAO,EAAE,IAAI,CAAC,OAAO,EAAE,YAAY,CAAC,CAAC;gBACpH,OAAO,IAAI,CAAC;YACd,CAAC;YACD,IAAI,IAAI,CAAC,WAAW,KAAK,SAAS,IAAI,IAAI,CAAC,OAAO,KAAK,SAAS,IAAI,IAAI,CAAC,OAAO,KAAK,SAAS,EAAE,CAAC;gBAC/F,IAAI,CAAC,OAAO,GAAG,kBAAkB,CAAC,kCAAkC,CAAC,IAAI,CAAC,OAAO,EAAE,IAAI,CAAC,WAAW,EAAE,YAAY,CAAC,CAAC;gBACnH,OAAO,IAAI,CAAC;YACd,CAAC;YACD,IAAI,IAAI,CAAC,WAAW,KAAK,SAAS,IAAI,IAAI,CAAC,OAAO,KAAK,SAAS,IAAI,IAAI,CAAC,OAAO,KAAK,SAAS,EAAE,CAAC;gBAC/F,IAAI,CAAC,OAAO,GAAG,kBAAkB,CAAC,kCAAkC,CAAC,IAAI,CAAC,OAAO,EAAE,IAAI,CAAC,WAAW,EAAE,YAAY,CAAC,CAAC;gBACnH,OAAO,IAAI,CAAC;YACd,CAAC;YACD,OAAO,KAAK,CAAC;QACf,CAAC;QACD,oCAAoC;QACpC,IAAI,IAAI,CAAC,WAAW,KAAK,SAAS,IAAI,IAAI,CAAC,OAAO,KAAK,SAAS,IAAI,IAAI,CAAC,OAAO,KAAK,SAAS;YAC5F,OAAO,KAAK,CAAC;QACf,IAAI,IAAI,CAAC,QAAQ,EAAE,CAAC,CAAC,yBAAyB;YAC5C,IAAI,CAAC,QAAQ,GAAG,KAAK,CAAC,aAAa,CACjC,IAAI,CAAC,QAAQ,CAAC,OAAO;gBACrB,kBAAkB,CAAC,gCAAgC,CAAC,IAAI,CAAC,OAAO,EAAE,IAAI,CAAC,OAAO,EAAE,IAAI,CAAC,WAAW,CAAC,CAClG,CAAC;YACF,OAAO,IAAI,CAAC;QACd,CAAC;QACD,IAAI,IAAI,CAAC,QAAQ,EAAE,CAAC,CAAC,yBAAyB;YAC5C,IAAI,CAAC,QAAQ,GAAG,KAAK,CAAC,aAAa,CACjC,IAAI,CAAC,QAAQ,CAAC,OAAO;gBACrB,kBAAkB,CAAC,gCAAgC,CAAC,IAAI,CAAC,OAAO,EAAE,IAAI,CAAC,OAAO,EAAE,IAAI,CAAC,WAAW,CAAC,CAClG,CAAC;YACF,OAAO,IAAI,CAAC;QACd,CAAC;QACD,OAAO,KAAK,CAAC;IACf,CAAC;IACO,qBAAqB,CAAC,CAAqB,EAAE,CAAqB;QACxE,IAAI,CAAC,KAAK,SAAS,IAAI,CAAC,KAAK,SAAS;YACpC,OAAO,IAAI,CAAC;QACd,IAAI,CAAC,KAAK,SAAS,IAAI,CAAC,KAAK,SAAS;YACpC,OAAO,QAAQ,CAAC,gBAAgB,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;QACzC,OAAO,KAAK,CAAC;IACf,CAAC;IACO,kBAAkB,CAAC,CAAoB,EAAE,CAAoB;QACnE,IAAI,CAAC,KAAK,SAAS,IAAI,CAAC,KAAK,SAAS;YACpC,OAAO,IAAI,CAAC;QACd,IAAI,CAAC,KAAK,SAAS,IAAI,CAAC,KAAK,SAAS;YACpC,OAAO,CAAC,CAAC,0BAA0B,CAAC,CAAC,CAAC,CAAC;QACzC,OAAO,KAAK,CAAC;IACf,CAAC;IACD,8EAA8E;IACvE,aAAa,CAAC,KAAuC;QAC1D,IAAI,CAAC,KAAK;YACR,OAAO,KAAK,CAAC;QACf,IAAI,CAAC,IAAI,CAAC,qBAAqB,CAAC,IAAI,CAAC,OAAO,EAAE,KAAK,CAAC,OAAO,CAAC;YAC1D,OAAO,KAAK,CAAC;QACf,IAAI,CAAC,IAAI,CAAC,qBAAqB,CAAC,IAAI,CAAC,OAAO,EAAE,KAAK,CAAC,OAAO,CAAC;YAC1D,OAAO,KAAK,CAAC;QACf,IAAI,CAAC,IAAI,CAAC,kBAAkB,CAAC,IAAI,CAAC,QAAQ,EAAE,KAAK,CAAC,QAAQ,CAAC;YACzD,OAAO,KAAK,CAAC;QACf,IAAI,CAAC,IAAI,CAAC,kBAAkB,CAAC,IAAI,CAAC,QAAQ,EAAE,KAAK,CAAC,QAAQ,CAAC;YACzD,OAAO,KAAK,CAAC;QACf,IAAI,CAAC,IAAI,CAAC,qBAAqB,CAAC,IAAI,CAAC,WAAW,EAAE,KAAK,CAAC,WAAW,CAAC;YAClE,OAAO,KAAK,CAAC;QACf,OAAO,IAAI,CAAC;IACd,CAAC;IACD,iFAAiF;IAC1E,gBAAgB,CAAC,CAAS;QAC/B,IAAI,CAAC,KAAK,CAAC;YACT,OAAO;QACT,IAAI,IAAI,CAAC,OAAO,KAAK,SAAS;YAC5B,IAAI,CAAC,OAAO,IAAI,CAAC,CAAC;QACpB,IAAI,IAAI,CAAC,OAAO,KAAK,SAAS;YAC5B,IAAI,CAAC,OAAO,IAAI,CAAC,CAAC;QACpB,IAAI,IAAI,CAAC,WAAW,KAAK,SAAS;YAChC,IAAI,CAAC,WAAW,IAAI,CAAC,CAAC;IAC1B,CAAC;IACD,uEAAuE;IAChE,MAAM,CAAC,cAAc,CAC1B,CAA8C,EAAE,CAA8C;QAE9F,IAAI,CAAC,KAAK,SAAS;YACjB,OAAO,CAAC,KAAK,SAAS,CAAC;QACzB,OAAO,CAAC,CAAC,aAAa,CAAC,CAAC,CAAC,CAAC;IAC5B,CAAC;CACF","sourcesContent":["/*---------------------------------------------------------------------------------------------\r\n* Copyright (c) Bentley Systems, Incorporated. All rights reserved.\r\n* See LICENSE.md in the project root for license terms and full copyright notice.\r\n*--------------------------------------------------------------------------------------------*/\r\n\r\n/** @packageDocumentation\r\n * @module Curve\r\n */\r\nimport { Geometry } from \"../../Geometry\";\r\nimport { Angle } from \"../../geometry3d/Angle\";\r\nimport { TransitionSpiral3d } from \"./TransitionSpiral3d\";\r\n\r\n/**\r\n * A true transition spiral is a curve defined by its curvature, with the curvature function symmetric about midpoint.\r\n * * The symmetry condition creates a relationship among the following 4 quantities:\r\n * ** curvature0 = curvature (i.e. 1/radius) at start.\r\n * ** curvature1 = curvature (i.e. 1/radius) at end.\r\n * ** sweepRadians = signed turning angle from start to end.\r\n * ** arcLength = length of curve.\r\n * * The relationship is the equation:\r\n * ** `sweepRadians = arcLength * average Curvature = arcLength * 0.5 * (curvature0 + curvature1)`\r\n * * That is, regardless of any curvature properties other than symmetry, specifying any 3 of the quantities fully\r\n * determines the remaining one.\r\n * @public\r\n */\r\n// see internaldocs/Spiral.md for more info\r\nexport class TransitionConditionalProperties {\r\n  /**\r\n   * Radius at start (radius of the osculating circle at the spiral segment's start).\r\n   * radius0 === 0 means straight line at start.\r\n   */\r\n  public radius0: number | undefined;\r\n  /**\r\n   * Radius at end (radius of the osculating circle at the spiral segment's end).\r\n   * radius1 === 0 means straight line at end.\r\n   */\r\n  public radius1: number | undefined;\r\n  /**\r\n   * Bearing at start (the signed angle from global x-axis to the start tangent of the spiral).\r\n   * * If the direction from x-axis to start tangent is CCW, the start bearing is positive.\r\n   */\r\n  public bearing0: Angle | undefined;\r\n  /**\r\n   * Bearing at end (the signed angle from global x-axis to the end tangent of the spiral).\r\n   * * If the direction from x-axis to end tangent is CCW, the end bearing is positive.\r\n   */\r\n  public bearing1: Angle | undefined;\r\n  /** Curve length. */\r\n  public curveLength: number | undefined;\r\n  /**\r\n   * Capture numeric or undefined values.\r\n   * @param radius0 start radius or undefined.\r\n   * @param radius1 end radius or undefined.\r\n   * @param bearing0 start bearing or undefined.\r\n   * @param bearing1 end bearing or undefined.\r\n   * @param arcLength arc length or undefined.\r\n   */\r\n  public constructor(\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  ) {\r\n    this.radius0 = radius0;\r\n    this.radius1 = radius1;\r\n    this.bearing0 = bearing0;\r\n    this.bearing1 = bearing1;\r\n    this.curveLength = arcLength;\r\n  }\r\n  /** Return the number of defined values among the 5 properties. */\r\n  public numDefinedProperties() {\r\n    return Geometry.defined01(this.radius0)\r\n      + Geometry.defined01(this.radius1)\r\n      + Geometry.defined01(this.bearing0)\r\n      + Geometry.defined01(this.bearing1)\r\n      + Geometry.defined01(this.curveLength);\r\n  }\r\n  /** Clone with all 5 properties (i.e., preserve undefined states). */\r\n  public clone(): TransitionConditionalProperties {\r\n    return new TransitionConditionalProperties(\r\n      this.radius0,\r\n      this.radius1,\r\n      this.bearing0 === undefined ? undefined : this.bearing0.clone(),\r\n      this.bearing1 === undefined ? undefined : this.bearing1.clone(),\r\n      this.curveLength,\r\n    );\r\n  }\r\n  /** Return true if all 5 properties are defined and agree equationally. */\r\n  public getIsValidCompleteSet(): boolean {\r\n    if (this.curveLength !== undefined && this.bearing0 !== undefined && this.bearing1 !== undefined\r\n      && this.radius0 !== undefined && this.radius1 !== undefined) {\r\n      const arcLength = TransitionSpiral3d.radiusRadiusSweepRadiansToArcLength(\r\n        this.radius0, this.radius1, this.bearing1.radians - this.bearing0.radians,\r\n      );\r\n      return Geometry.isSameCoordinate(this.curveLength, arcLength);\r\n    }\r\n    return false;\r\n  }\r\n  /**\r\n   * Examine which properties are defined and compute the (single) undefined.\r\n   * @returns true if the input state had precisely one undefined member.\r\n   */\r\n  public tryResolveAnySingleUnknown(): boolean {\r\n    if (this.getIsValidCompleteSet())\r\n      return true;\r\n    if (this.bearing0 && this.bearing1) {\r\n      const sweepRadians = this.bearing1.radians - this.bearing0.radians;\r\n      if (this.curveLength === undefined && this.radius0 !== undefined && this.radius1 !== undefined) {\r\n        this.curveLength = TransitionSpiral3d.radiusRadiusSweepRadiansToArcLength(this.radius0, this.radius1, sweepRadians);\r\n        return true;\r\n      }\r\n      if (this.curveLength !== undefined && this.radius0 === undefined && this.radius1 !== undefined) {\r\n        this.radius0 = TransitionSpiral3d.radius1LengthSweepRadiansToRadius0(this.radius1, this.curveLength, sweepRadians);\r\n        return true;\r\n      }\r\n      if (this.curveLength !== undefined && this.radius0 !== undefined && this.radius1 === undefined) {\r\n        this.radius1 = TransitionSpiral3d.radius0LengthSweepRadiansToRadius1(this.radius0, this.curveLength, sweepRadians);\r\n        return true;\r\n      }\r\n      return false;\r\n    }\r\n    // at least one bearing is undefined\r\n    if (this.curveLength === undefined || this.radius0 === undefined || this.radius1 === undefined)\r\n      return false;\r\n    if (this.bearing0) { // bearing 1 is undefined\r\n      this.bearing1 = Angle.createRadians(\r\n        this.bearing0.radians +\r\n        TransitionSpiral3d.radiusRadiusLengthToSweepRadians(this.radius0, this.radius1, this.curveLength),\r\n      );\r\n      return true;\r\n    }\r\n    if (this.bearing1) { // bearing 0 is undefined\r\n      this.bearing0 = Angle.createRadians(\r\n        this.bearing1.radians -\r\n        TransitionSpiral3d.radiusRadiusLengthToSweepRadians(this.radius0, this.radius1, this.curveLength),\r\n      );\r\n      return true;\r\n    }\r\n    return false;\r\n  }\r\n  private almostEqualCoordinate(a: number | undefined, b: number | undefined): boolean {\r\n    if (a === undefined && b === undefined)\r\n      return true;\r\n    if (a !== undefined && b !== undefined)\r\n      return Geometry.isSameCoordinate(a, b);\r\n    return false;\r\n  }\r\n  private almostEqualBearing(a: Angle | undefined, b: Angle | undefined): boolean {\r\n    if (a === undefined && b === undefined)\r\n      return true;\r\n    if (a !== undefined && b !== undefined)\r\n      return a.isAlmostEqualNoPeriodShift(b);\r\n    return false;\r\n  }\r\n  /** Test if `this` and `other` have matching numeric and undefined members. */\r\n  public isAlmostEqual(other?: TransitionConditionalProperties): boolean {\r\n    if (!other)\r\n      return false;\r\n    if (!this.almostEqualCoordinate(this.radius0, other.radius0))\r\n      return false;\r\n    if (!this.almostEqualCoordinate(this.radius1, other.radius1))\r\n      return false;\r\n    if (!this.almostEqualBearing(this.bearing0, other.bearing0))\r\n      return false;\r\n    if (!this.almostEqualBearing(this.bearing1, other.bearing1))\r\n      return false;\r\n    if (!this.almostEqualCoordinate(this.curveLength, other.curveLength))\r\n      return false;\r\n    return true;\r\n  }\r\n  /** Apply a NONZERO scale factor to all distances. If `a` is zero, do nothing. */\r\n  public applyScaleFactor(a: number): void {\r\n    if (a === 0)\r\n      return;\r\n    if (this.radius0 !== undefined)\r\n      this.radius0 *= a;\r\n    if (this.radius1 !== undefined)\r\n      this.radius1 *= a;\r\n    if (this.curveLength !== undefined)\r\n      this.curveLength *= a;\r\n  }\r\n  /** Test if `a` and `b` have matching numeric and undefined members. */\r\n  public static areAlmostEqual(\r\n    a: TransitionConditionalProperties | undefined, b: TransitionConditionalProperties | undefined,\r\n  ): boolean {\r\n    if (a === undefined)\r\n      return b === undefined;\r\n    return a.isAlmostEqual(b);\r\n  }\r\n}\r\n"]}
+{"version":3,"file":"TransitionConditionalProperties.js","sourceRoot":"","sources":["../../../../src/curve/spiral/TransitionConditionalProperties.ts"],"names":[],"mappings":"AAAA;;;+FAG+F;AAE/F;;GAEG;AACH,OAAO,EAAE,QAAQ,EAAE,MAAM,gBAAgB,CAAC;AAC1C,OAAO,EAAE,KAAK,EAAE,MAAM,wBAAwB,CAAC;AAC/C,OAAO,EAAE,kBAAkB,EAAE,MAAM,sBAAsB,CAAC;AAE1D;;;;;;;;;;;;GAYG;AACH,2CAA2C;AAC3C,MAAM,OAAO,+BAA+B;IAC1C;;;OAGG;IACI,OAAO,CAAqB;IACnC;;;OAGG;IACI,OAAO,CAAqB;IACnC;;;OAGG;IACI,QAAQ,CAAoB;IACnC;;;OAGG;IACI,QAAQ,CAAoB;IACnC,oBAAoB;IACb,WAAW,CAAqB;IACvC;;;;;;;OAOG;IACH,YACE,OAA2B,EAC3B,OAA2B,EAC3B,QAA2B,EAC3B,QAA2B,EAC3B,SAA6B;QAE7B,IAAI,CAAC,OAAO,GAAG,OAAO,CAAC;QACvB,IAAI,CAAC,OAAO,GAAG,OAAO,CAAC;QACvB,IAAI,CAAC,QAAQ,GAAG,QAAQ,CAAC;QACzB,IAAI,CAAC,QAAQ,GAAG,QAAQ,CAAC;QACzB,IAAI,CAAC,WAAW,GAAG,SAAS,CAAC;IAC/B,CAAC;IACD,kEAAkE;IAC3D,oBAAoB;QACzB,OAAO,QAAQ,CAAC,SAAS,CAAC,IAAI,CAAC,OAAO,CAAC;cACnC,QAAQ,CAAC,SAAS,CAAC,IAAI,CAAC,OAAO,CAAC;cAChC,QAAQ,CAAC,SAAS,CAAC,IAAI,CAAC,QAAQ,CAAC;cACjC,QAAQ,CAAC,SAAS,CAAC,IAAI,CAAC,QAAQ,CAAC;cACjC,QAAQ,CAAC,SAAS,CAAC,IAAI,CAAC,WAAW,CAAC,CAAC;IAC3C,CAAC;IACD,qEAAqE;IAC9D,KAAK;QACV,OAAO,IAAI,+BAA+B,CACxC,IAAI,CAAC,OAAO,EACZ,IAAI,CAAC,OAAO,EACZ,IAAI,CAAC,QAAQ,KAAK,SAAS,CAAC,CAAC,CAAC,SAAS,CAAC,CAAC,CAAC,IAAI,CAAC,QAAQ,CAAC,KAAK,EAAE,EAC/D,IAAI,CAAC,QAAQ,KAAK,SAAS,CAAC,CAAC,CAAC,SAAS,CAAC,CAAC,CAAC,IAAI,CAAC,QAAQ,CAAC,KAAK,EAAE,EAC/D,IAAI,CAAC,WAAW,CACjB,CAAC;IACJ,CAAC;IACD,0EAA0E;IACnE,qBAAqB;QAC1B,IAAI,IAAI,CAAC,WAAW,KAAK,SAAS,IAAI,IAAI,CAAC,QAAQ,KAAK,SAAS,IAAI,IAAI,CAAC,QAAQ,KAAK,SAAS;eAC3F,IAAI,CAAC,OAAO,KAAK,SAAS,IAAI,IAAI,CAAC,OAAO,KAAK,SAAS,EAAE,CAAC;YAC9D,MAAM,SAAS,GAAG,kBAAkB,CAAC,mCAAmC,CACtE,IAAI,CAAC,OAAO,EAAE,IAAI,CAAC,OAAO,EAAE,IAAI,CAAC,QAAQ,CAAC,OAAO,GAAG,IAAI,CAAC,QAAQ,CAAC,OAAO,CAC1E,CAAC;YACF,OAAO,QAAQ,CAAC,gBAAgB,CAAC,IAAI,CAAC,WAAW,EAAE,SAAS,CAAC,CAAC;QAChE,CAAC;QACD,OAAO,KAAK,CAAC;IACf,CAAC;IACD;;;OAGG;IACI,0BAA0B;QAC/B,IAAI,IAAI,CAAC,qBAAqB,EAAE;YAC9B,OAAO,IAAI,CAAC;QACd,IAAI,IAAI,CAAC,QAAQ,IAAI,IAAI,CAAC,QAAQ,EAAE,CAAC;YACnC,MAAM,YAAY,GAAG,IAAI,CAAC,QAAQ,CAAC,OAAO,GAAG,IAAI,CAAC,QAAQ,CAAC,OAAO,CAAC;YACnE,IAAI,IAAI,CAAC,WAAW,KAAK,SAAS,IAAI,IAAI,CAAC,OAAO,KAAK,SAAS,IAAI,IAAI,CAAC,OAAO,KAAK,SAAS,EAAE,CAAC;gBAC/F,IAAI,CAAC,WAAW,GAAG,kBAAkB,CAAC,mCAAmC,CAAC,IAAI,CAAC,OAAO,EAAE,IAAI,CAAC,OAAO,EAAE,YAAY,CAAC,CAAC;gBACpH,OAAO,IAAI,CAAC;YACd,CAAC;YACD,IAAI,IAAI,CAAC,WAAW,KAAK,SAAS,IAAI,IAAI,CAAC,OAAO,KAAK,SAAS,IAAI,IAAI,CAAC,OAAO,KAAK,SAAS,EAAE,CAAC;gBAC/F,IAAI,CAAC,OAAO,GAAG,kBAAkB,CAAC,kCAAkC,CAAC,IAAI,CAAC,OAAO,EAAE,IAAI,CAAC,WAAW,EAAE,YAAY,CAAC,CAAC;gBACnH,OAAO,IAAI,CAAC;YACd,CAAC;YACD,IAAI,IAAI,CAAC,WAAW,KAAK,SAAS,IAAI,IAAI,CAAC,OAAO,KAAK,SAAS,IAAI,IAAI,CAAC,OAAO,KAAK,SAAS,EAAE,CAAC;gBAC/F,IAAI,CAAC,OAAO,GAAG,kBAAkB,CAAC,kCAAkC,CAAC,IAAI,CAAC,OAAO,EAAE,IAAI,CAAC,WAAW,EAAE,YAAY,CAAC,CAAC;gBACnH,OAAO,IAAI,CAAC;YACd,CAAC;YACD,OAAO,KAAK,CAAC;QACf,CAAC;QACD,oCAAoC;QACpC,IAAI,IAAI,CAAC,WAAW,KAAK,SAAS,IAAI,IAAI,CAAC,OAAO,KAAK,SAAS,IAAI,IAAI,CAAC,OAAO,KAAK,SAAS;YAC5F,OAAO,KAAK,CAAC;QACf,IAAI,IAAI,CAAC,QAAQ,EAAE,CAAC,CAAC,yBAAyB;YAC5C,IAAI,CAAC,QAAQ,GAAG,KAAK,CAAC,aAAa,CACjC,IAAI,CAAC,QAAQ,CAAC,OAAO;gBACrB,kBAAkB,CAAC,gCAAgC,CAAC,IAAI,CAAC,OAAO,EAAE,IAAI,CAAC,OAAO,EAAE,IAAI,CAAC,WAAW,CAAC,CAClG,CAAC;YACF,OAAO,IAAI,CAAC;QACd,CAAC;QACD,IAAI,IAAI,CAAC,QAAQ,EAAE,CAAC,CAAC,yBAAyB;YAC5C,IAAI,CAAC,QAAQ,GAAG,KAAK,CAAC,aAAa,CACjC,IAAI,CAAC,QAAQ,CAAC,OAAO;gBACrB,kBAAkB,CAAC,gCAAgC,CAAC,IAAI,CAAC,OAAO,EAAE,IAAI,CAAC,OAAO,EAAE,IAAI,CAAC,WAAW,CAAC,CAClG,CAAC;YACF,OAAO,IAAI,CAAC;QACd,CAAC;QACD,OAAO,KAAK,CAAC;IACf,CAAC;IACO,qBAAqB,CAAC,CAAqB,EAAE,CAAqB;QACxE,IAAI,CAAC,KAAK,SAAS,IAAI,CAAC,KAAK,SAAS;YACpC,OAAO,IAAI,CAAC;QACd,IAAI,CAAC,KAAK,SAAS,IAAI,CAAC,KAAK,SAAS;YACpC,OAAO,QAAQ,CAAC,gBAAgB,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;QACzC,OAAO,KAAK,CAAC;IACf,CAAC;IACO,kBAAkB,CAAC,CAAoB,EAAE,CAAoB;QACnE,IAAI,CAAC,KAAK,SAAS,IAAI,CAAC,KAAK,SAAS;YACpC,OAAO,IAAI,CAAC;QACd,IAAI,CAAC,KAAK,SAAS,IAAI,CAAC,KAAK,SAAS;YACpC,OAAO,CAAC,CAAC,0BAA0B,CAAC,CAAC,CAAC,CAAC;QACzC,OAAO,KAAK,CAAC;IACf,CAAC;IACD,8EAA8E;IACvE,aAAa,CAAC,KAAuC;QAC1D,IAAI,CAAC,KAAK;YACR,OAAO,KAAK,CAAC;QACf,IAAI,CAAC,IAAI,CAAC,qBAAqB,CAAC,IAAI,CAAC,OAAO,EAAE,KAAK,CAAC,OAAO,CAAC;YAC1D,OAAO,KAAK,CAAC;QACf,IAAI,CAAC,IAAI,CAAC,qBAAqB,CAAC,IAAI,CAAC,OAAO,EAAE,KAAK,CAAC,OAAO,CAAC;YAC1D,OAAO,KAAK,CAAC;QACf,IAAI,CAAC,IAAI,CAAC,kBAAkB,CAAC,IAAI,CAAC,QAAQ,EAAE,KAAK,CAAC,QAAQ,CAAC;YACzD,OAAO,KAAK,CAAC;QACf,IAAI,CAAC,IAAI,CAAC,kBAAkB,CAAC,IAAI,CAAC,QAAQ,EAAE,KAAK,CAAC,QAAQ,CAAC;YACzD,OAAO,KAAK,CAAC;QACf,IAAI,CAAC,IAAI,CAAC,qBAAqB,CAAC,IAAI,CAAC,WAAW,EAAE,KAAK,CAAC,WAAW,CAAC;YAClE,OAAO,KAAK,CAAC;QACf,OAAO,IAAI,CAAC;IACd,CAAC;IACD,iFAAiF;IAC1E,gBAAgB,CAAC,CAAS;QAC/B,IAAI,CAAC,KAAK,CAAC;YACT,OAAO;QACT,IAAI,IAAI,CAAC,OAAO,KAAK,SAAS;YAC5B,IAAI,CAAC,OAAO,IAAI,CAAC,CAAC;QACpB,IAAI,IAAI,CAAC,OAAO,KAAK,SAAS;YAC5B,IAAI,CAAC,OAAO,IAAI,CAAC,CAAC;QACpB,IAAI,IAAI,CAAC,WAAW,KAAK,SAAS;YAChC,IAAI,CAAC,WAAW,IAAI,CAAC,CAAC;IAC1B,CAAC;IACD,uEAAuE;IAChE,MAAM,CAAC,cAAc,CAC1B,CAA8C,EAAE,CAA8C;QAE9F,IAAI,CAAC,KAAK,SAAS;YACjB,OAAO,CAAC,KAAK,SAAS,CAAC;QACzB,OAAO,CAAC,CAAC,aAAa,CAAC,CAAC,CAAC,CAAC;IAC5B,CAAC;CACF","sourcesContent":["/*---------------------------------------------------------------------------------------------\r\n* Copyright (c) Bentley Systems, Incorporated. All rights reserved.\r\n* See LICENSE.md in the project root for license terms and full copyright notice.\r\n*--------------------------------------------------------------------------------------------*/\r\n\r\n/** @packageDocumentation\r\n * @module Curve\r\n */\r\nimport { Geometry } from \"../../Geometry\";\r\nimport { Angle } from \"../../geometry3d/Angle\";\r\nimport { TransitionSpiral3d } from \"./TransitionSpiral3d\";\r\n\r\n/**\r\n * A true transition spiral is a curve defined by its curvature, with the curvature function symmetric about midpoint.\r\n * * The symmetry condition creates a relationship among the following 4 quantities:\r\n * ** curvature0 = curvature (i.e. 1/radius) at start.\r\n * ** curvature1 = curvature (i.e. 1/radius) at end.\r\n * ** sweepRadians = signed turning angle from start to end.\r\n * ** arcLength = length of curve.\r\n * * The relationship is the equation:\r\n * ** `sweepRadians = arcLength * average Curvature = arcLength * 0.5 * (curvature0 + curvature1)`\r\n * * That is, regardless of any curvature properties other than symmetry, specifying any 3 of the quantities fully\r\n * determines the remaining one.\r\n * @public\r\n */\r\n// see internaldocs/Spiral.md for more info\r\nexport class TransitionConditionalProperties {\r\n  /**\r\n   * Radius at start (radius of the osculating circle at the spiral segment's start).\r\n   * radius0 === 0 means straight line at start.\r\n   */\r\n  public radius0: number | undefined;\r\n  /**\r\n   * Radius at end (radius of the osculating circle at the spiral segment's end).\r\n   * radius1 === 0 means straight line at end.\r\n   */\r\n  public radius1: number | undefined;\r\n  /**\r\n   * Bearing at start (the signed angle from global x-axis to the start tangent of the spiral).\r\n   * * If the direction from x-axis to start tangent is CCW, the start bearing is positive.\r\n   */\r\n  public bearing0: Angle | undefined;\r\n  /**\r\n   * Bearing at end (the signed angle from global x-axis to the end tangent of the spiral).\r\n   * * If the direction from x-axis to end tangent is CCW, the end bearing is positive.\r\n   */\r\n  public bearing1: Angle | undefined;\r\n  /** Curve length. */\r\n  public curveLength: number | undefined;\r\n  /**\r\n   * Capture numeric or undefined values.\r\n   * @param radius0 start radius or undefined.\r\n   * @param radius1 end radius or undefined.\r\n   * @param bearing0 start bearing or undefined.\r\n   * @param bearing1 end bearing or undefined.\r\n   * @param arcLength arc length or undefined.\r\n   */\r\n  public constructor(\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  ) {\r\n    this.radius0 = radius0;\r\n    this.radius1 = radius1;\r\n    this.bearing0 = bearing0;\r\n    this.bearing1 = bearing1;\r\n    this.curveLength = arcLength;\r\n  }\r\n  /** Return the number of defined values among the 5 properties. */\r\n  public numDefinedProperties() {\r\n    return Geometry.defined01(this.radius0)\r\n      + Geometry.defined01(this.radius1)\r\n      + Geometry.defined01(this.bearing0)\r\n      + Geometry.defined01(this.bearing1)\r\n      + Geometry.defined01(this.curveLength);\r\n  }\r\n  /** Clone with all 5 properties (i.e., preserve undefined states). */\r\n  public clone(): TransitionConditionalProperties {\r\n    return new TransitionConditionalProperties(\r\n      this.radius0,\r\n      this.radius1,\r\n      this.bearing0 === undefined ? undefined : this.bearing0.clone(),\r\n      this.bearing1 === undefined ? undefined : this.bearing1.clone(),\r\n      this.curveLength,\r\n    );\r\n  }\r\n  /** Return true if all 5 properties are defined and agree equationally. */\r\n  public getIsValidCompleteSet(): boolean {\r\n    if (this.curveLength !== undefined && this.bearing0 !== undefined && this.bearing1 !== undefined\r\n      && this.radius0 !== undefined && this.radius1 !== undefined) {\r\n      const arcLength = TransitionSpiral3d.radiusRadiusSweepRadiansToArcLength(\r\n        this.radius0, this.radius1, this.bearing1.radians - this.bearing0.radians,\r\n      );\r\n      return Geometry.isSameCoordinate(this.curveLength, arcLength);\r\n    }\r\n    return false;\r\n  }\r\n  /**\r\n   * Examine which properties are defined and compute the (single) undefined.\r\n   * @returns true if the input state had precisely one undefined member.\r\n   */\r\n  public tryResolveAnySingleUnknown(): this is {radius0: number, radius1: number, bearing0: Angle, bearing1: Angle, curveLength: number} {\r\n    if (this.getIsValidCompleteSet())\r\n      return true;\r\n    if (this.bearing0 && this.bearing1) {\r\n      const sweepRadians = this.bearing1.radians - this.bearing0.radians;\r\n      if (this.curveLength === undefined && this.radius0 !== undefined && this.radius1 !== undefined) {\r\n        this.curveLength = TransitionSpiral3d.radiusRadiusSweepRadiansToArcLength(this.radius0, this.radius1, sweepRadians);\r\n        return true;\r\n      }\r\n      if (this.curveLength !== undefined && this.radius0 === undefined && this.radius1 !== undefined) {\r\n        this.radius0 = TransitionSpiral3d.radius1LengthSweepRadiansToRadius0(this.radius1, this.curveLength, sweepRadians);\r\n        return true;\r\n      }\r\n      if (this.curveLength !== undefined && this.radius0 !== undefined && this.radius1 === undefined) {\r\n        this.radius1 = TransitionSpiral3d.radius0LengthSweepRadiansToRadius1(this.radius0, this.curveLength, sweepRadians);\r\n        return true;\r\n      }\r\n      return false;\r\n    }\r\n    // at least one bearing is undefined\r\n    if (this.curveLength === undefined || this.radius0 === undefined || this.radius1 === undefined)\r\n      return false;\r\n    if (this.bearing0) { // bearing 1 is undefined\r\n      this.bearing1 = Angle.createRadians(\r\n        this.bearing0.radians +\r\n        TransitionSpiral3d.radiusRadiusLengthToSweepRadians(this.radius0, this.radius1, this.curveLength),\r\n      );\r\n      return true;\r\n    }\r\n    if (this.bearing1) { // bearing 0 is undefined\r\n      this.bearing0 = Angle.createRadians(\r\n        this.bearing1.radians -\r\n        TransitionSpiral3d.radiusRadiusLengthToSweepRadians(this.radius0, this.radius1, this.curveLength),\r\n      );\r\n      return true;\r\n    }\r\n    return false;\r\n  }\r\n  private almostEqualCoordinate(a: number | undefined, b: number | undefined): boolean {\r\n    if (a === undefined && b === undefined)\r\n      return true;\r\n    if (a !== undefined && b !== undefined)\r\n      return Geometry.isSameCoordinate(a, b);\r\n    return false;\r\n  }\r\n  private almostEqualBearing(a: Angle | undefined, b: Angle | undefined): boolean {\r\n    if (a === undefined && b === undefined)\r\n      return true;\r\n    if (a !== undefined && b !== undefined)\r\n      return a.isAlmostEqualNoPeriodShift(b);\r\n    return false;\r\n  }\r\n  /** Test if `this` and `other` have matching numeric and undefined members. */\r\n  public isAlmostEqual(other?: TransitionConditionalProperties): boolean {\r\n    if (!other)\r\n      return false;\r\n    if (!this.almostEqualCoordinate(this.radius0, other.radius0))\r\n      return false;\r\n    if (!this.almostEqualCoordinate(this.radius1, other.radius1))\r\n      return false;\r\n    if (!this.almostEqualBearing(this.bearing0, other.bearing0))\r\n      return false;\r\n    if (!this.almostEqualBearing(this.bearing1, other.bearing1))\r\n      return false;\r\n    if (!this.almostEqualCoordinate(this.curveLength, other.curveLength))\r\n      return false;\r\n    return true;\r\n  }\r\n  /** Apply a NONZERO scale factor to all distances. If `a` is zero, do nothing. */\r\n  public applyScaleFactor(a: number): void {\r\n    if (a === 0)\r\n      return;\r\n    if (this.radius0 !== undefined)\r\n      this.radius0 *= a;\r\n    if (this.radius1 !== undefined)\r\n      this.radius1 *= a;\r\n    if (this.curveLength !== undefined)\r\n      this.curveLength *= a;\r\n  }\r\n  /** Test if `a` and `b` have matching numeric and undefined members. */\r\n  public static areAlmostEqual(\r\n    a: TransitionConditionalProperties | undefined, b: TransitionConditionalProperties | undefined,\r\n  ): boolean {\r\n    if (a === undefined)\r\n      return b === undefined;\r\n    return a.isAlmostEqual(b);\r\n  }\r\n}\r\n"]}

package/lib/esm/geometry3d/Ellipsoid.d.ts

@@ -114,7 +114,7 @@ export declare class Ellipsoid implements Clipper {
      * * For a given pair of points on an ellipsoid, construct an arc (possibly elliptical) which
      *   * passes through both points
      *   * is completely within the ellipsoid surface
-     *   * has its centerEvaluate a point on the ellipsoid at angles give in radians.
+     *   * has its center
      * * If the ellipsoid is a sphere, this is the shortest great-circle arc between the two points.
      * * If the ellipsoid is not a sphere, this is close to but not precisely the shortest path.
      * @param thetaARadians longitude, in radians, for pointA
@@ -122,11 +122,10 @@ export declare class Ellipsoid implements Clipper {
      * @param thetaBRadians longitude, in radians, for pointB
      * @param phiBRadians latitude, in radians, for pointB
      * @param result optional preallocated result
+     * @param returns great arc, or undefined if input points are identical or diametrically opposite.
      */
     radiansPairToGreatArc(thetaARadians: number, phiARadians: number, thetaBRadians: number, phiBRadians: number, result?: Arc3d): Arc3d | undefined;
-    /**
-     * See radiansPairToGreatArc, which does this computation with positions from `angleA` and `angleB` directly as radians
-     */
+    /** See [[radiansPairToGreatArc]], which does this computation with positions from `angleA` and `angleB` directly as radians. */
     anglePairToGreatArc(angleA: LongitudeLatitudeNumber, angleB: LongitudeLatitudeNumber, result?: Arc3d): Arc3d | undefined;
     /**
      * Construct an arc for the section cut of a plane with the ellipsoid.
@@ -167,11 +166,24 @@ export declare class Ellipsoid implements Clipper {
      */
     constantLatitudeArc(longitudeSweep: AngleSweep, latitude: Angle, result?: Arc3d): Arc3d | undefined;
     /**
-     * * create a section arc with and end at positions A and B, and in plane with the normal at a fractional
-     *    interpolation between.
+     * Create a section arc with start and end at positions A and B, and in plane with the normal at a fractional
+     * interpolation between.
+     * @param angleA start point of arc (given as angles on this ellipsoid)
+     * @param intermediateNormalFraction fraction at which to interpolate normals at A and B to define arc plane
+     * @param angleB end point of arc (given as angles on this ellipsoid)
+     * @param result optional pre-allocated object to populate and return
+     * @returns arc in the plane defined by the normal at the intermediate point, or `undefined` if calculation fails.
+     */
+    sectionArcInPlaneOfInterpolatedNormal(angleA: LongitudeLatitudeNumber, intermediateNormalFraction: number, angleB: LongitudeLatitudeNumber, result?: Arc3d): Arc3d | undefined;
+    /**
+     * Create a section arc with start and end at positions A and B, and in plane with the normal at a fractional
+     * interpolation between.
      * @param angleA start point of arc (given as angles on this ellipsoid)
-     * @param intermediateNormalFraction
+     * @param intermediateNormalFraction fraction at which to interpolate normals at A and B to define arc plane
      * @param angleB end point of arc (given as angles on this ellipsoid)
+     * @returns arc in the plane defined by the normal at the intermediate point. If calculation fails, return an
+     * arc with zero matrix.
+     * @deprecated Prefer [[sectionArcInPlaneOfInterpolatedNormal]], which has expanded return type.
      */
     sectionArcWithIntermediateNormal(angleA: LongitudeLatitudeNumber, intermediateNormalFraction: number, angleB: LongitudeLatitudeNumber): Arc3d;
     /**

package/lib/esm/geometry3d/Ellipsoid.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"Ellipsoid.d.ts","sourceRoot":"","sources":["../../../src/geometry3d/Ellipsoid.ts"],"names":[],"mappings":"AAKA;;GAEG;AAEH,OAAO,EAAE,6BAA6B,EAA2B,MAAM,kCAAkC,CAAC;AAC1G,OAAO,EAAE,OAAO,EAAE,MAAM,uBAAuB,CAAC;AAChD,OAAO,EAAE,KAAK,EAAE,MAAM,gBAAgB,CAAC;AAEvC,OAAO,EAAE,oBAAoB,EAAE,kCAAkC,EAAE,MAAM,yBAAyB,CAAC;AAEnG,OAAO,EAAE,OAAO,EAAE,MAAM,uBAAuB,CAAC;AAKhD,OAAO,EAAE,KAAK,EAAE,MAAM,SAAS,CAAC;AAChC,OAAO,EAAE,UAAU,EAAE,MAAM,cAAc,CAAC;AAC1C,OAAO,EAAE,SAAS,EAAE,MAAM,mBAAmB,CAAC;AAC9C,OAAO,EAAE,uBAAuB,EAAE,MAAM,6BAA6B,CAAC;AACtE,OAAO,EAAE,QAAQ,EAAE,MAAM,YAAY,CAAC;AACtC,OAAO,EAAE,4BAA4B,EAAE,MAAM,gCAAgC,CAAC;AAC9E,OAAO,EAAE,yBAAyB,EAAE,MAAM,6BAA6B,CAAC;AACxE,OAAO,EAAE,OAAO,EAAE,QAAQ,EAAE,MAAM,mBAAmB,CAAC;AACtD,OAAO,EAAW,OAAO,EAAE,MAAM,SAAS,CAAC;AAC3C,OAAO,EAAE,KAAK,EAAE,MAAM,SAAS,CAAC;AAChC,OAAO,EAAE,SAAS,EAAE,MAAM,aAAa,CAAC;AACxC,OAAO,EAAE,MAAM,EAAE,MAAM,YAAY,CAAC;AA4GpC;;;;;;;;;GASG;AACH,qBAAa,SAAU,YAAW,OAAO;IACvC,OAAO,CAAC,UAAU,CAAY;IAC9B,OAAO,CAAC,gBAAgB,CAAW;IACnC,OAAO,CAAC,gBAAgB,CAAW;IACnC,OAAO,CAAC,WAAW,CAAU;IAC7B,OAAO,CAAC,WAAW,CAAU;IAC7B,OAAO;IAOP;;OAEG;WACW,MAAM,CAAC,iBAAiB,CAAC,EAAE,SAAS,GAAG,QAAQ,GAAG,SAAS;IAQzE;;;;;;;OAOG;WACW,uBAAuB,CAAC,MAAM,EAAE,OAAO,EAAE,IAAI,EAAE,QAAQ,GAAG,SAAS,EAAE,OAAO,EAAE,MAAM,EAAE,OAAO,EAAE,MAAM,EAAE,OAAO,EAAE,MAAM,GAAG,SAAS;IAQhJ;;;;;;OAMG;IACH,IAAW,YAAY,IAAI,SAAS,CAA4B;IAChE;;;;;;OAMG;IACI,YAAY,CAAC,UAAU,EAAE,MAAM,EAAE,MAAM,CAAC,EAAE,OAAO,GAAG,OAAO,GAAG,SAAS;IAG9E;;;;;OAKG;IACI,YAAY,CAAC,UAAU,EAAE,MAAM,EAAE,MAAM,CAAC,EAAE,OAAO,GAAG,OAAO;IAIlE,2CAA2C;IACpC,KAAK,IAAI,SAAS;IAGzB,oCAAoC;IAC7B,aAAa,CAAC,KAAK,EAAE,SAAS,GAAG,OAAO;IAG/C,wCAAwC;IACjC,mBAAmB,CAAC,SAAS,EAAE,SAAS,GAAG,OAAO;IAIzD;;;OAGG;IACI,gBAAgB,CAAC,SAAS,EAAE,SAAS,GAAG,SAAS,GAAG,SAAS;IAKpE;;;OAGG;IACI,qBAAqB,CAAC,UAAU,EAAE,OAAO,GAAG,uBAAuB,GAAG,SAAS;IAKtF;;OAEG;IACI,aAAa,CAAC,QAAQ,EAAE,OAAO,GAAG,KAAK,GAAG,SAAS;IAoB1D;;;;;;;;OAQG;IACI,YAAY,CAAC,GAAG,EAAE,KAAK,EAAE,YAAY,EAAE,MAAM,EAAE,GAAG,SAAS,EAAE,GAAG,EAAE,OAAO,EAAE,GAAG,SAAS,EAAE,eAAe,EAAE,uBAAuB,EAAE,GAAG,SAAS,GAAG,MAAM;IA+B/J,4DAA4D;IACrD,yBAAyB,CAAC,aAAa,EAAE,MAAM,EAAE,aAAa,EAAE,MAAM,EAAE,WAAW,EAAE,MAAM,EAAE,WAAW,EAAE,MAAM,EAAE,MAAM,CAAC,EAAE,OAAO,GAAG,OAAO;IAqCnJ;;;;;OAKG;IACI,cAAc,CAAC,YAAY,EAAE,MAAM,EAAE,UAAU,EAAE,MAAM,EAAE,MAAM,CAAC,EAAE,OAAO,GAAG,OAAO;IAQ1F;;;;;;;;;;;;OAYG;IACI,qBAAqB,CAC1B,aAAa,EAAE,MAAM,EAAE,WAAW,EAAE,MAAM,EAC1C,aAAa,EAAE,MAAM,EAAE,WAAW,EAAE,MAAM,EAC1C,MAAM,CAAC,EAAE,KAAK,GAAG,KAAK,GAAG,SAAS;IAapC;;OAEG;IACI,mBAAmB,CAAC,MAAM,EAAE,uBAAuB,EAAE,MAAM,EAAE,uBAAuB,EAAE,MAAM,CAAC,EAAE,KAAK,GAAG,KAAK,GAAG,SAAS;IAI/H;;;OAGG;IACI,kBAAkB,CAAC,KAAK,EAAE,4BAA4B,GAAG,KAAK,GAAG,SAAS;IAsBjF;;;;;OAKG;IACI,uCAAuC,CAAC,YAAY,EAAE,uBAAuB,EAAE,YAAY,EAAE,uBAAuB,EAAE,aAAa,EAAE,QAAQ,EAClJ,MAAM,CAAC,EAAE,KAAK,GAAG,KAAK,GAAG,SAAS;IA6BpC;;;;;;;;;;OAUG;IACI,gCAAgC,CACrC,aAAa,EAAE,MAAM,EAAE,WAAW,EAAE,MAAM,EAC1C,aAAa,EAAE,MAAM,EAAE,WAAW,EAAE,MAAM,EAC1C,MAAM,CAAC,EAAE,SAAS,GAAG,SAAS,GAAG,SAAS;IAc5C;;;;;OAKG;IACI,oBAAoB,CAAC,SAAS,EAAE,KAAK,EAAE,aAAa,EAAE,UAAU,EAAE,MAAM,CAAC,EAAE,KAAK,GAAG,KAAK,GAAG,SAAS;IAU3G;;;;;OAKG;IACI,mBAAmB,CAAC,cAAc,EAAE,UAAU,EAAE,QAAQ,EAAE,KAAK,EAAE,MAAM,CAAC,EAAE,KAAK,GAAG,KAAK,GAAG,SAAS;IAY1G;;;;;;OAMG;IACI,gCAAgC,CACrC,MAAM,EAAE,uBAAuB,EAC/B,0BAA0B,EAAE,MAAM,EAClC,MAAM,EAAE,uBAAuB,GAAG,KAAK;IAQzC;;;;;;;;;;OAUG;IACI,4BAA4B,CAAC,YAAY,EAAE,MAAM,EAAE,UAAU,EAAE,MAAM,EAAE,iBAAiB,UAAO,EAAE,MAAM,CAAC,EAAE,yBAAyB,GAAG,yBAAyB;IAkBtK;;;;;;;;;;;;OAYG;IACI,6BAA6B,CAAC,YAAY,EAAE,MAAM,EAAE,UAAU,EAAE,MAAM,EAC3E,KAAK,EAAE,OAAO,EACd,OAAO,EAAE,QAAQ,EACjB,KAAK,EAAE,QAAQ,EACf,YAAY,EAAE,QAAQ,EACtB,QAAQ,EAAE,QAAQ,EAClB,UAAU,EAAE,QAAQ;IAmBtB;;;;;;;OAOG;IACI,oBAAoB,CAAC,YAAY,EAAE,MAAM,EAAE,UAAU,EAAE,MAAM,EAAE,MAAM,CAAC,EAAE,SAAS,GAAG,SAAS,GAAG,SAAS;IAIhH;;;;;;OAMG;IACI,sBAAsB,CAAC,YAAY,EAAE,MAAM,EAAE,UAAU,EAAE,MAAM,EAAE,MAAM,CAAC,EAAE,KAAK,GAAG,KAAK,GAAG,SAAS;IAK1G;;OAEG;IACI,qBAAqB,CAAC,MAAM,EAAE,QAAQ,EAAE,MAAM,CAAC,EAAE,uBAAuB,GAAG,uBAAuB;IAWzG;;;;OAIG;IACI,yCAAyC,CAAC,cAAc,EAAE,SAAS,GAAG,SAAS,EAAE,WAAW,EAAE,uBAAuB,EAAE,MAAM,CAAC,EAAE,uBAAuB,GAAG,uBAAuB,GAAG,SAAS;IAMpM;;;OAGG;WACW,sBAAsB,CAAC,SAAS,EAAE,SAAS,GAAG,SAAS,EAAE,YAAY,EAAE,MAAM,EAAE,UAAU,EAAE,MAAM,EAAE,MAAM,CAAC,EAAE,KAAK,GAAG,KAAK,GAAG,SAAS;IAWnJ,uDAAuD;IAChD,iBAAiB,CAAC,KAAK,EAAE,OAAO,GAAG,OAAO;IAMjD,+GAA+G;IACxG,+BAA+B,CAAC,EAAE,EAAE,MAAM,EAAE,EAAE,EAAE,MAAM,EAAE,MAAM,EAAE,OAAO,EAAE,MAAM,EAAE,OAAO,EAAE,QAAQ,CAAC,EAAE,oBAAoB,GAAG,OAAO;IAqC1I,0GAA0G;IACnG,2BAA2B,CAAC,GAAG,EAAE,KAAK,EAAE,QAAQ,CAAC,EAAE,kCAAkC,GAAG,OAAO;CA6CvG;AACD;;;;;;;GAOG;AACH,qBAAa,cAAe,YAAW,SAAS;IACvC,SAAS,EAAE,SAAS,CAAC;IACrB,cAAc,EAAE,UAAU,CAAC;IAC3B,aAAa,EAAE,UAAU,CAAC;IACjC;;;;;OAKG;IACH,OAAO;IAKP;;;;;OAKG;WACW,aAAa,CAAC,SAAS,EAAE,SAAS,EAAE,cAAc,EAAE,UAAU,EAAE,aAAa,EAAE,UAAU,GAAG,cAAc;IAGxH,uFAAuF;IAChF,iBAAiB,CAAC,iBAAiB,EAAE,MAAM,EAAE,gBAAgB,EAAE,MAAM,EAAE,MAAM,CAAC,EAAE,OAAO,GAAG,OAAO;IAGxG;;OAEG;IACI,4BAA4B,CAAC,iBAAiB,EAAE,MAAM,EAAE,gBAAgB,EAAE,MAAM,EAAE,MAAM,CAAC,EAAE,yBAAyB,GAAG,yBAAyB;IAUvJ,qFAAqF;IAC9E,KAAK,CAAC,MAAM,CAAC,EAAE,OAAO,GAAG,OAAO;IAGvC;;OAEG;IACI,YAAY,CAAC,GAAG,EAAE,KAAK,EAAE,eAAe,EAAE,OAAO,EAAE,qCAAqC,GAAE,OAAe,GAAG,6BAA6B,EAAE;IA4BlJ;;;;;OAKG;IACI,cAAc,CAAC,QAAQ,EAAE,uBAAuB,EAAE,sBAAsB,GAAE,OAAc,GAAG,OAAO;IAKzG;;;;;;;;OAQG;IACI,qBAAqB,CAAC,QAAQ,EAAE,uBAAuB,EAAE,MAAM,CAAC,EAAE,KAAK,GAAG,KAAK,GAAG,SAAS;IAOlG;;;;;;OAMG;IACI,kBAAkB,CAAC,iBAAiB,EAAE,MAAM,EAAE,WAAW,EAAE,MAAM,EAAE,CAAC,GAAE,MAAU,EAAE,MAAM,CAAC,EAAE,uBAAuB,GAAG,uBAAuB;IAGnJ,8DAA8D;IACvD,qBAAqB,CAAC,UAAU,EAAE,OAAO,GAAG,uBAAuB,GAAG,SAAS;CAGvF;AA4DD;;;GAGG;AACH,qBAAa,iBAAiB;IAC5B,8BAA8B;IACvB,YAAY,EAAE,MAAM,CAAC;IAC5B,+BAA+B;IACxB,UAAU,EAAE,MAAM,CAAC;IACnB,KAAK,EAAE,OAAO,CAAC;IACf,MAAM,EAAE,QAAQ,CAAC;IACjB,IAAI,EAAE,QAAQ,CAAC;IACf,OAAO,EAAE,QAAQ,CAAC;IAClB,KAAK,EAAE,QAAQ,CAAC;IAChB,UAAU,EAAE,QAAQ,CAAC;IACrB,OAAO,EAAE,QAAQ,CAAC;;IAYzB,mDAAmD;IAC5C,kCAAkC,CAAC,SAAS,EAAE,SAAS;IAI9D,OAAO,CAAC,MAAM,CAAC,SAAS,CAAC,CAAW;IACpC,OAAO,CAAC,MAAM,CAAC,SAAS,CAAC,CAAW;IACpC,OAAO,CAAC,MAAM,CAAC,YAAY,CAAC,CAAW;IACvC;;;;;;;;OAQG;WACW,sBAAsB,CAAC,MAAM,EAAE,iBAAiB,EAAE,MAAM,EAAE,iBAAiB,EAAE,MAAM,EAAE,iBAAiB,EAAE,MAAM,EAAE,YAAY;IAkB1I;;OAEG;IACI,QAAQ,IAAI,uBAAuB;CAG3C;AACD;;;;GAIG;AACH,qBAAa,kBAAkB;IAC7B,OAAO,CAAC,WAAW,CAAQ;IAC3B,OAAO,CAAC,WAAW,CAAsB;IACzC,OAAO,CAAC,kBAAkB,CAAqB;IAC/C,OAAO;IAIP;;;;;;;OAOG;WACW,kBAAkB,CAAC,iBAAiB,EAAE,SAAS,EAC3D,WAAW,EAAE,uBAAuB,EAAE,SAAS,EAAE,uBAAuB,EAAE,OAAO,EAAE,MAAM,GAAG,KAAK,GAAG,iBAAiB,EAAE,GAAG,SAAS;IAsDrI,OAAO,CAAC,uCAAuC;IAc/C,OAAO,CAAC,WAAW;IAUnB;;;;OAIG;IACH,OAAO,CAAC,SAAS;IAqBjB,OAAO,CAAC,KAAK;IAGb;;;;;;;;;;OAUG;WACW,kCAAkC,CAAC,SAAS,EAAE,SAAS,EACnE,MAAM,EAAE,uBAAuB,EAC/B,MAAM,EAAE,uBAAuB,EAAE,SAAS,EAAE,MAAM,EAClD,4BAA4B,EAAE,MAAM,EACpC,4BAA4B,EAAE,MAAM,GAAG;QAAE,YAAY,EAAE,KAAK,CAAC;QAAC,oCAAoC,EAAE,MAAM,CAAA;KAAE,GAAG,SAAS;CA6B3H"}
+{"version":3,"file":"Ellipsoid.d.ts","sourceRoot":"","sources":["../../../src/geometry3d/Ellipsoid.ts"],"names":[],"mappings":"AAKA;;GAEG;AAEH,OAAO,EAAE,6BAA6B,EAA2B,MAAM,kCAAkC,CAAC;AAC1G,OAAO,EAAE,OAAO,EAAE,MAAM,uBAAuB,CAAC;AAChD,OAAO,EAAE,KAAK,EAAE,MAAM,gBAAgB,CAAC;AAEvC,OAAO,EAAE,oBAAoB,EAAE,kCAAkC,EAAE,MAAM,yBAAyB,CAAC;AAEnG,OAAO,EAAE,OAAO,EAAE,MAAM,uBAAuB,CAAC;AAKhD,OAAO,EAAE,KAAK,EAAE,MAAM,SAAS,CAAC;AAChC,OAAO,EAAE,UAAU,EAAE,MAAM,cAAc,CAAC;AAC1C,OAAO,EAAE,SAAS,EAAE,MAAM,mBAAmB,CAAC;AAC9C,OAAO,EAAE,uBAAuB,EAAE,MAAM,6BAA6B,CAAC;AACtE,OAAO,EAAE,QAAQ,EAAE,MAAM,YAAY,CAAC;AACtC,OAAO,EAAE,4BAA4B,EAAE,MAAM,gCAAgC,CAAC;AAC9E,OAAO,EAAE,yBAAyB,EAAE,MAAM,6BAA6B,CAAC;AACxE,OAAO,EAAE,OAAO,EAAE,QAAQ,EAAE,MAAM,mBAAmB,CAAC;AACtD,OAAO,EAAW,OAAO,EAAE,MAAM,SAAS,CAAC;AAC3C,OAAO,EAAE,KAAK,EAAE,MAAM,SAAS,CAAC;AAChC,OAAO,EAAE,SAAS,EAAE,MAAM,aAAa,CAAC;AACxC,OAAO,EAAE,MAAM,EAAE,MAAM,YAAY,CAAC;AA4GpC;;;;;;;;;GASG;AACH,qBAAa,SAAU,YAAW,OAAO;IACvC,OAAO,CAAC,UAAU,CAAY;IAC9B,OAAO,CAAC,gBAAgB,CAAW;IACnC,OAAO,CAAC,gBAAgB,CAAW;IACnC,OAAO,CAAC,WAAW,CAAU;IAC7B,OAAO,CAAC,WAAW,CAAU;IAC7B,OAAO;IAOP;;OAEG;WACW,MAAM,CAAC,iBAAiB,CAAC,EAAE,SAAS,GAAG,QAAQ,GAAG,SAAS;IAQzE;;;;;;;OAOG;WACW,uBAAuB,CAAC,MAAM,EAAE,OAAO,EAAE,IAAI,EAAE,QAAQ,GAAG,SAAS,EAAE,OAAO,EAAE,MAAM,EAAE,OAAO,EAAE,MAAM,EAAE,OAAO,EAAE,MAAM,GAAG,SAAS;IAQhJ;;;;;;OAMG;IACH,IAAW,YAAY,IAAI,SAAS,CAA4B;IAChE;;;;;;OAMG;IACI,YAAY,CAAC,UAAU,EAAE,MAAM,EAAE,MAAM,CAAC,EAAE,OAAO,GAAG,OAAO,GAAG,SAAS;IAG9E;;;;;OAKG;IACI,YAAY,CAAC,UAAU,EAAE,MAAM,EAAE,MAAM,CAAC,EAAE,OAAO,GAAG,OAAO;IAIlE,2CAA2C;IACpC,KAAK,IAAI,SAAS;IAGzB,oCAAoC;IAC7B,aAAa,CAAC,KAAK,EAAE,SAAS,GAAG,OAAO;IAG/C,wCAAwC;IACjC,mBAAmB,CAAC,SAAS,EAAE,SAAS,GAAG,OAAO;IAIzD;;;OAGG;IACI,gBAAgB,CAAC,SAAS,EAAE,SAAS,GAAG,SAAS,GAAG,SAAS;IAKpE;;;OAGG;IACI,qBAAqB,CAAC,UAAU,EAAE,OAAO,GAAG,uBAAuB,GAAG,SAAS;IAKtF;;OAEG;IACI,aAAa,CAAC,QAAQ,EAAE,OAAO,GAAG,KAAK,GAAG,SAAS;IAoB1D;;;;;;;;OAQG;IACI,YAAY,CAAC,GAAG,EAAE,KAAK,EAAE,YAAY,EAAE,MAAM,EAAE,GAAG,SAAS,EAAE,GAAG,EAAE,OAAO,EAAE,GAAG,SAAS,EAAE,eAAe,EAAE,uBAAuB,EAAE,GAAG,SAAS,GAAG,MAAM;IA+B/J,4DAA4D;IACrD,yBAAyB,CAAC,aAAa,EAAE,MAAM,EAAE,aAAa,EAAE,MAAM,EAAE,WAAW,EAAE,MAAM,EAAE,WAAW,EAAE,MAAM,EAAE,MAAM,CAAC,EAAE,OAAO,GAAG,OAAO;IAqCnJ;;;;;OAKG;IACI,cAAc,CAAC,YAAY,EAAE,MAAM,EAAE,UAAU,EAAE,MAAM,EAAE,MAAM,CAAC,EAAE,OAAO,GAAG,OAAO;IAQ1F;;;;;;;;;;;;;OAaG;IACI,qBAAqB,CAC1B,aAAa,EAAE,MAAM,EAAE,WAAW,EAAE,MAAM,EAC1C,aAAa,EAAE,MAAM,EAAE,WAAW,EAAE,MAAM,EAC1C,MAAM,CAAC,EAAE,KAAK,GAAG,KAAK,GAAG,SAAS;IAYpC,gIAAgI;IACzH,mBAAmB,CAAC,MAAM,EAAE,uBAAuB,EAAE,MAAM,EAAE,uBAAuB,EAAE,MAAM,CAAC,EAAE,KAAK,GAAG,KAAK,GAAG,SAAS;IAI/H;;;OAGG;IACI,kBAAkB,CAAC,KAAK,EAAE,4BAA4B,GAAG,KAAK,GAAG,SAAS;IAsBjF;;;;;OAKG;IACI,uCAAuC,CAAC,YAAY,EAAE,uBAAuB,EAAE,YAAY,EAAE,uBAAuB,EAAE,aAAa,EAAE,QAAQ,EAClJ,MAAM,CAAC,EAAE,KAAK,GAAG,KAAK,GAAG,SAAS;IA6BpC;;;;;;;;;;OAUG;IACI,gCAAgC,CACrC,aAAa,EAAE,MAAM,EAAE,WAAW,EAAE,MAAM,EAC1C,aAAa,EAAE,MAAM,EAAE,WAAW,EAAE,MAAM,EAC1C,MAAM,CAAC,EAAE,SAAS,GAAG,SAAS,GAAG,SAAS;IAc5C;;;;;OAKG;IACI,oBAAoB,CAAC,SAAS,EAAE,KAAK,EAAE,aAAa,EAAE,UAAU,EAAE,MAAM,CAAC,EAAE,KAAK,GAAG,KAAK,GAAG,SAAS;IAU3G;;;;;OAKG;IACI,mBAAmB,CAAC,cAAc,EAAE,UAAU,EAAE,QAAQ,EAAE,KAAK,EAAE,MAAM,CAAC,EAAE,KAAK,GAAG,KAAK,GAAG,SAAS;IAY1G;;;;;;;;OAQG;IACI,qCAAqC,CAAC,MAAM,EAAE,uBAAuB,EAAE,0BAA0B,EAAE,MAAM,EAAE,MAAM,EAAE,uBAAuB,EAAE,MAAM,CAAC,EAAE,KAAK,GAAG,KAAK,GAAG,SAAS;IAQrL;;;;;;;;;OASG;IACI,gCAAgC,CAAC,MAAM,EAAE,uBAAuB,EAAE,0BAA0B,EAAE,MAAM,EAAE,MAAM,EAAE,uBAAuB,GAAG,KAAK;IAKpJ;;;;;;;;;;OAUG;IACI,4BAA4B,CAAC,YAAY,EAAE,MAAM,EAAE,UAAU,EAAE,MAAM,EAAE,iBAAiB,UAAO,EAAE,MAAM,CAAC,EAAE,yBAAyB,GAAG,yBAAyB;IAkBtK;;;;;;;;;;;;OAYG;IACI,6BAA6B,CAAC,YAAY,EAAE,MAAM,EAAE,UAAU,EAAE,MAAM,EAC3E,KAAK,EAAE,OAAO,EACd,OAAO,EAAE,QAAQ,EACjB,KAAK,EAAE,QAAQ,EACf,YAAY,EAAE,QAAQ,EACtB,QAAQ,EAAE,QAAQ,EAClB,UAAU,EAAE,QAAQ;IAmBtB;;;;;;;OAOG;IACI,oBAAoB,CAAC,YAAY,EAAE,MAAM,EAAE,UAAU,EAAE,MAAM,EAAE,MAAM,CAAC,EAAE,SAAS,GAAG,SAAS,GAAG,SAAS;IAIhH;;;;;;OAMG;IACI,sBAAsB,CAAC,YAAY,EAAE,MAAM,EAAE,UAAU,EAAE,MAAM,EAAE,MAAM,CAAC,EAAE,KAAK,GAAG,KAAK,GAAG,SAAS;IAK1G;;OAEG;IACI,qBAAqB,CAAC,MAAM,EAAE,QAAQ,EAAE,MAAM,CAAC,EAAE,uBAAuB,GAAG,uBAAuB;IAWzG;;;;OAIG;IACI,yCAAyC,CAAC,cAAc,EAAE,SAAS,GAAG,SAAS,EAAE,WAAW,EAAE,uBAAuB,EAAE,MAAM,CAAC,EAAE,uBAAuB,GAAG,uBAAuB,GAAG,SAAS;IAMpM;;;OAGG;WACW,sBAAsB,CAAC,SAAS,EAAE,SAAS,GAAG,SAAS,EAAE,YAAY,EAAE,MAAM,EAAE,UAAU,EAAE,MAAM,EAAE,MAAM,CAAC,EAAE,KAAK,GAAG,KAAK,GAAG,SAAS;IAWnJ,uDAAuD;IAChD,iBAAiB,CAAC,KAAK,EAAE,OAAO,GAAG,OAAO;IAMjD,+GAA+G;IACxG,+BAA+B,CAAC,EAAE,EAAE,MAAM,EAAE,EAAE,EAAE,MAAM,EAAE,MAAM,EAAE,OAAO,EAAE,MAAM,EAAE,OAAO,EAAE,QAAQ,CAAC,EAAE,oBAAoB,GAAG,OAAO;IAqC1I,0GAA0G;IACnG,2BAA2B,CAAC,GAAG,EAAE,KAAK,EAAE,QAAQ,CAAC,EAAE,kCAAkC,GAAG,OAAO;CA6CvG;AACD;;;;;;;GAOG;AACH,qBAAa,cAAe,YAAW,SAAS;IACvC,SAAS,EAAE,SAAS,CAAC;IACrB,cAAc,EAAE,UAAU,CAAC;IAC3B,aAAa,EAAE,UAAU,CAAC;IACjC;;;;;OAKG;IACH,OAAO;IAKP;;;;;OAKG;WACW,aAAa,CAAC,SAAS,EAAE,SAAS,EAAE,cAAc,EAAE,UAAU,EAAE,aAAa,EAAE,UAAU,GAAG,cAAc;IAGxH,uFAAuF;IAChF,iBAAiB,CAAC,iBAAiB,EAAE,MAAM,EAAE,gBAAgB,EAAE,MAAM,EAAE,MAAM,CAAC,EAAE,OAAO,GAAG,OAAO;IAGxG;;OAEG;IACI,4BAA4B,CAAC,iBAAiB,EAAE,MAAM,EAAE,gBAAgB,EAAE,MAAM,EAAE,MAAM,CAAC,EAAE,yBAAyB,GAAG,yBAAyB;IAUvJ,qFAAqF;IAC9E,KAAK,CAAC,MAAM,CAAC,EAAE,OAAO,GAAG,OAAO;IAGvC;;OAEG;IACI,YAAY,CAAC,GAAG,EAAE,KAAK,EAAE,eAAe,EAAE,OAAO,EAAE,qCAAqC,GAAE,OAAe,GAAG,6BAA6B,EAAE;IA4BlJ;;;;;OAKG;IACI,cAAc,CAAC,QAAQ,EAAE,uBAAuB,EAAE,sBAAsB,GAAE,OAAc,GAAG,OAAO;IAKzG;;;;;;;;OAQG;IACI,qBAAqB,CAAC,QAAQ,EAAE,uBAAuB,EAAE,MAAM,CAAC,EAAE,KAAK,GAAG,KAAK,GAAG,SAAS;IAOlG;;;;;;OAMG;IACI,kBAAkB,CAAC,iBAAiB,EAAE,MAAM,EAAE,WAAW,EAAE,MAAM,EAAE,CAAC,GAAE,MAAU,EAAE,MAAM,CAAC,EAAE,uBAAuB,GAAG,uBAAuB;IAGnJ,8DAA8D;IACvD,qBAAqB,CAAC,UAAU,EAAE,OAAO,GAAG,uBAAuB,GAAG,SAAS;CAGvF;AA4DD;;;GAGG;AACH,qBAAa,iBAAiB;IAC5B,8BAA8B;IACvB,YAAY,EAAE,MAAM,CAAC;IAC5B,+BAA+B;IACxB,UAAU,EAAE,MAAM,CAAC;IACnB,KAAK,EAAE,OAAO,CAAC;IACf,MAAM,EAAE,QAAQ,CAAC;IACjB,IAAI,EAAE,QAAQ,CAAC;IACf,OAAO,EAAE,QAAQ,CAAC;IAClB,KAAK,EAAE,QAAQ,CAAC;IAChB,UAAU,EAAE,QAAQ,CAAC;IACrB,OAAO,EAAE,QAAQ,CAAC;;IAYzB,mDAAmD;IAC5C,kCAAkC,CAAC,SAAS,EAAE,SAAS;IAI9D,OAAO,CAAC,MAAM,CAAC,SAAS,CAAC,CAAW;IACpC,OAAO,CAAC,MAAM,CAAC,SAAS,CAAC,CAAW;IACpC,OAAO,CAAC,MAAM,CAAC,YAAY,CAAC,CAAW;IACvC;;;;;;;;OAQG;WACW,sBAAsB,CAAC,MAAM,EAAE,iBAAiB,EAAE,MAAM,EAAE,iBAAiB,EAAE,MAAM,EAAE,iBAAiB,EAAE,MAAM,EAAE,YAAY;IAkB1I;;OAEG;IACI,QAAQ,IAAI,uBAAuB;CAG3C;AACD;;;;GAIG;AACH,qBAAa,kBAAkB;IAC7B,OAAO,CAAC,WAAW,CAAQ;IAC3B,OAAO,CAAC,WAAW,CAAsB;IACzC,OAAO,CAAC,kBAAkB,CAAqB;IAC/C,OAAO;IAIP;;;;;;;OAOG;WACW,kBAAkB,CAAC,iBAAiB,EAAE,SAAS,EAC3D,WAAW,EAAE,uBAAuB,EAAE,SAAS,EAAE,uBAAuB,EAAE,OAAO,EAAE,MAAM,GAAG,KAAK,GAAG,iBAAiB,EAAE,GAAG,SAAS;IAsDrI,OAAO,CAAC,uCAAuC;IAc/C,OAAO,CAAC,WAAW;IAUnB;;;;OAIG;IACH,OAAO,CAAC,SAAS;IAqBjB,OAAO,CAAC,KAAK;IAGb;;;;;;;;;;OAUG;WACW,kCAAkC,CAAC,SAAS,EAAE,SAAS,EACnE,MAAM,EAAE,uBAAuB,EAC/B,MAAM,EAAE,uBAAuB,EAAE,SAAS,EAAE,MAAM,EAClD,4BAA4B,EAAE,MAAM,EACpC,4BAA4B,EAAE,MAAM,GAAG;QAAE,YAAY,EAAE,KAAK,CAAC;QAAC,oCAAoC,EAAE,MAAM,CAAA;KAAE,GAAG,SAAS;CA6B3H"}

package/lib/esm/geometry3d/Ellipsoid.js

@@ -348,7 +348,7 @@ export class Ellipsoid {
      * * For a given pair of points on an ellipsoid, construct an arc (possibly elliptical) which
      *   * passes through both points
      *   * is completely within the ellipsoid surface
-     *   * has its centerEvaluate a point on the ellipsoid at angles give in radians.
+     *   * has its center
      * * If the ellipsoid is a sphere, this is the shortest great-circle arc between the two points.
      * * If the ellipsoid is not a sphere, this is close to but not precisely the shortest path.
      * @param thetaARadians longitude, in radians, for pointA
@@ -356,12 +356,12 @@ export class Ellipsoid {
      * @param thetaBRadians longitude, in radians, for pointB
      * @param phiBRadians latitude, in radians, for pointB
      * @param result optional preallocated result
+     * @param returns great arc, or undefined if input points are identical or diametrically opposite.
      */
     radiansPairToGreatArc(thetaARadians, phiARadians, thetaBRadians, phiBRadians, result) {
         SphereImplicit.radiansToUnitSphereXYZ(thetaARadians, phiARadians, this._workUnitVectorA);
         SphereImplicit.radiansToUnitSphereXYZ(thetaBRadians, phiBRadians, this._workUnitVectorB);
         const sweepAngle = this._workUnitVectorA.angleTo(this._workUnitVectorB);
-        // the unit vectors (on unit sphere) are never 0, so this cannot fail.
         const matrix = Matrix3d.createRigidFromColumns(this._workUnitVectorA, this._workUnitVectorB, AxisOrder.XYZ);
         if (matrix !== undefined) {
             const matrix1 = this._transform.matrix.multiplyMatrixMatrix(matrix);
@@ -369,9 +369,7 @@ export class Ellipsoid {
         }
         return undefined;
     }
-    /**
-     * See radiansPairToGreatArc, which does this computation with positions from `angleA` and `angleB` directly as radians
-     */
+    /** See [[radiansPairToGreatArc]], which does this computation with positions from `angleA` and `angleB` directly as radians. */
     anglePairToGreatArc(angleA, angleB, result) {
         return this.radiansPairToGreatArc(angleA.longitudeRadians, angleA.latitudeRadians, angleB.longitudeRadians, angleB.latitudeRadians, result);
     }
@@ -496,18 +494,35 @@ export class Ellipsoid {
         return Arc3d.create(center, vector0, vector90, longitudeSweep, result);
     }
     /**
-     * * create a section arc with and end at positions A and B, and in plane with the normal at a fractional
-     *    interpolation between.
+     * Create a section arc with start and end at positions A and B, and in plane with the normal at a fractional
+     * interpolation between.
      * @param angleA start point of arc (given as angles on this ellipsoid)
-     * @param intermediateNormalFraction
+     * @param intermediateNormalFraction fraction at which to interpolate normals at A and B to define arc plane
      * @param angleB end point of arc (given as angles on this ellipsoid)
+     * @param result optional pre-allocated object to populate and return
+     * @returns arc in the plane defined by the normal at the intermediate point, or `undefined` if calculation fails.
      */
-    sectionArcWithIntermediateNormal(angleA, intermediateNormalFraction, angleB) {
+    sectionArcInPlaneOfInterpolatedNormal(angleA, intermediateNormalFraction, angleB, result) {
         const normalA = this.radiansToUnitNormalRay(angleA.longitudeRadians, angleA.latitudeRadians);
         const normalB = this.radiansToUnitNormalRay(angleB.longitudeRadians, angleB.latitudeRadians);
+        if (!normalA || !normalB)
+            return undefined;
         const normal = normalA.direction.interpolate(intermediateNormalFraction, normalB.direction);
-        const arc = this.createSectionArcPointPointVectorInPlane(angleA, angleB, normal);
-        return arc;
+        return this.createSectionArcPointPointVectorInPlane(angleA, angleB, normal, result);
+    }
+    /**
+     * Create a section arc with start and end at positions A and B, and in plane with the normal at a fractional
+     * interpolation between.
+     * @param angleA start point of arc (given as angles on this ellipsoid)
+     * @param intermediateNormalFraction fraction at which to interpolate normals at A and B to define arc plane
+     * @param angleB end point of arc (given as angles on this ellipsoid)
+     * @returns arc in the plane defined by the normal at the intermediate point. If calculation fails, return an
+     * arc with zero matrix.
+     * @deprecated Prefer [[sectionArcInPlaneOfInterpolatedNormal]], which has expanded return type.
+     */
+    sectionArcWithIntermediateNormal(angleA, intermediateNormalFraction, angleB) {
+        const arc = this.sectionArcInPlaneOfInterpolatedNormal(angleA, intermediateNormalFraction, angleB);
+        return arc ?? Arc3d.createXYZXYZXYZ(0, 0, 0, 0, 0, 0, 0, 0, 0);
     }
     /**
      * Evaluate a point and derivatives with respect to angle on the ellipsoid at angles give in radians.