@itwin/core-geometry 5.3.0-dev.9 → 5.3.0

package/lib/cjs/curve/spiral/AustralianRailCorpXYEvaluator.js

@@ -26,9 +26,15 @@ class AustralianRailCorpXYEvaluator extends CubicEvaluator_1.CubicEvaluator {
         this._nominalLength1 = nominalLength1;
         this._nominalRadius1 = nominalRadius1;
     }
-    get nominalLength1() { return this._nominalLength1; }
-    get nominalRadius1() { return this._nominalRadius1; }
-    clone() { return new AustralianRailCorpXYEvaluator(this._nominalLength1, this._nominalRadius1, this._axisLength, this._cubicM); }
+    get nominalLength1() {
+        return this._nominalLength1;
+    }
+    get nominalRadius1() {
+        return this._nominalRadius1;
+    }
+    clone() {
+        return new AustralianRailCorpXYEvaluator(this._nominalLength1, this._nominalRadius1, this._axisLength, this._cubicM);
+    }
     static create(nominalLength1, nominalRadius1) {
         const axisLength = AustralianRailCorpXYEvaluator.radiusAndNominalLengthToAxisLength(nominalRadius1, nominalLength1);
         const phi = this.radiusAndAxisLengthToPhi(nominalRadius1, axisLength);
@@ -36,11 +42,7 @@ class AustralianRailCorpXYEvaluator extends CubicEvaluator_1.CubicEvaluator {
         const cubicM = Math.tan(phi) / (3.0 * xc2);
         return new AustralianRailCorpXYEvaluator(nominalLength1, nominalRadius1, axisLength, cubicM);
     }
-    /**
-     * Compute the phi constant for AustralianRail spiral with given end radius and length along axis.
-     * @param nominalRadius1
-     * @param axisLength
-     */
+    /** Compute the phi constant for AustralianRail spiral with given end radius and length along axis. */
     static radiusAndAxisLengthToPhi(nominalRadius1, axisLength) {
         const xc = axisLength;
         const expr1 = (2. / Math.sqrt(3.));
@@ -53,26 +55,23 @@ class AustralianRailCorpXYEvaluator extends CubicEvaluator_1.CubicEvaluator {
         return Math.asin(expr1 * Math.cos(Math.acos(expr2) / 3. + expr3));
     }
     scaleInPlace(scaleFactor) {
-        // apply the scale factor to all contents.
-        // all distances scale directly . . .
+        // apply the scale factor to all contents; all distances scale directly
         this._nominalLength1 *= scaleFactor;
         this._nominalRadius1 *= scaleFactor;
         super.scaleInPlace(scaleFactor);
     }
-    /** Compute length along axis for AustralianRail spiral nominal radius and length.
-     *
-     */
+    /** Compute length along axis for AustralianRail spiral nominal radius and length. */
     static radiusAndNominalLengthToAxisLength(nominalRadius1, nominalLength1, tolerance = 1.0e-5, requiredConvergenceCount = 2) {
         const R = nominalRadius1;
         let idx = 0;
         let m, phi, xc2;
         let xc = .7 * nominalLength1;
         let convergenceCount = 0;
-        // remark: This converges quickly --
+        // remark: This converges quickly
         // for L=100, R=400
         //   ** full precision at 7th iteration.
         //   ** classic tolerance 1.0e-5 (7 digits from L) with requiredConvergenceCount = 1 gives 11 digits after 3 iterations
-        //   ** each iteration adds about 2 digits.   This is quite good for a successive replacement without derivative !!!
+        //   ** each iteration adds about 2 digits. This is quite good for a successive replacement without derivative
         //   ** Unanswerable question: If this is only done once and reused over all evaluations, do you want:
         //       * run the 7 iterations to get full precision
         //       * stop with the classic tolerance to get compatibility?
@@ -81,8 +80,7 @@ class AustralianRailCorpXYEvaluator extends CubicEvaluator_1.CubicEvaluator {
             xc2 = xc * xc;
             m = Math.tan(phi) / (3.0 * xc2);
             const m2x4 = m * m * xc2 * xc2;
-            const correction = xc * m2x4 * ((9. / 10) + m2x4 * (-(9. / 8.) + m2x4 * (+(729. / 208.) + m2x4 *
-                -(32805. / 2176.))));
+            const correction = xc * m2x4 * ((9. / 10) + m2x4 * (-(9. / 8.) + m2x4 * (+(729. / 208.) + m2x4 * -(32805. / 2176.))));
             const correctedLength = xc + correction;
             xc = (nominalLength1 / correctedLength) * xc;
             if (Math.abs(nominalLength1 - correctedLength) < tolerance) {
@@ -107,12 +105,15 @@ class AustralianRailCorpXYEvaluator extends CubicEvaluator_1.CubicEvaluator {
     }
     /**
      * Return a (quite good approximation) of fraction along x axis for given distance along spiral.
-     * * The AustralianRailSpiral has a supporting power series to approximately map distance along the spiral to an x coordinate.
-     * * The `xToFraction(x)` method quickly (with a single divide) converts this x to fraction used fro this.fractionToX (fraction), this.fractionToY(fraction) etc to get coordinates and derivatives.
+     * * The AustralianRailSpiral has a supporting power series to approximately map distance along the spiral to
+     * an x coordinate.
+     * * The `xToFraction(x)` method quickly (with a single divide) converts this x to fraction used for
+     * this.fractionToX (fraction), this.fractionToY(fraction) etc to get coordinates and derivatives.
      * * The x-to-distance relation is not as precise as the CurvePrimitive method moveSignedDistanceFromFraction.
-     * * It is supported here for users interested in replicating the AustralianRail distance mapping rather than the more accurate CurvePrimitive measurements.
-     * * Round tripping distance through (a) distanceAlongSpiralToAustralianApproximateX, (b) xToFraction, and (c) curveLengthBetweenFractions has
-     *   * 10 digit accuracy for L/R = 4, 12 digit accuracy for L/R = 10
+     * * It is supported here for users interested in replicating the AustralianRail distance mapping rather than the
+     * more accurate CurvePrimitive measurements.
+     * * Round tripping distance through (a) distanceAlongSpiralToAustralianApproximateX, (b) xToFraction, and
+     * (c) curveLengthBetweenFractions has 10 digit accuracy for L/R = 4, 12 digit accuracy for L/R = 10.
      * @param s distance along the axis.
      */
     distanceAlongSpiralToAustralianApproximateX(s) {

package/lib/cjs/curve/spiral/AustralianRailCorpXYEvaluator.js.map

@@ -1 +1 @@
-{"version":3,"file":"AustralianRailCorpXYEvaluator.js","sourceRoot":"","sources":["../../../../src/curve/spiral/AustralianRailCorpXYEvaluator.ts"],"names":[],"mappings":";AAAA;;;+FAG+F;;;AAE/F;;GAEG;AACH,6CAA0C;AAC1C,kDAA+C;AAC/C,qDAAkD;AAClD;;;;;;GAMG;AACH,MAAa,6BAA8B,SAAQ,+BAAc;IACvD,eAAe,CAAS;IACxB,eAAe,CAAS;IAChC,YAAoB,cAAsB,EAAE,cAAsB,EAAE,UAAkB,EAAE,MAAc;QACpG,KAAK,CAAC,UAAU,EAAE,MAAM,CAAC,CAAC;QAC1B,IAAI,CAAC,eAAe,GAAG,cAAc,CAAC;QACtC,IAAI,CAAC,eAAe,GAAG,cAAc,CAAC;IACxC,CAAC;IACD,IAAW,cAAc,KAAK,OAAO,IAAI,CAAC,eAAe,CAAC,CAAC,CAAC;IAC5D,IAAW,cAAc,KAAK,OAAO,IAAI,CAAC,eAAe,CAAC,CAAC,CAAC;IACrD,KAAK,KAAoC,OAAO,IAAI,6BAA6B,CAAC,IAAI,CAAC,eAAe,EAAE,IAAI,CAAC,eAAe,EAAE,IAAI,CAAC,WAAW,EAAE,IAAI,CAAC,OAAO,CAAC,CAAC,CAAC,CAAC;IAChK,MAAM,CAAC,MAAM,CAAC,cAAsB,EAAE,cAAsB;QACjE,MAAM,UAAU,GAAG,6BAA6B,CAAC,kCAAkC,CAAC,cAAc,EAAE,cAAc,CAAC,CAAC;QACpH,MAAM,GAAG,GAAG,IAAI,CAAC,wBAAwB,CAAC,cAAc,EAAE,UAAU,CAAC,CAAC;QACtE,MAAM,GAAG,GAAG,UAAU,GAAG,UAAU,CAAC;QACpC,MAAM,MAAM,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,GAAG,GAAG,CAAC,CAAC;QAC3C,OAAO,IAAI,6BAA6B,CAAC,cAAc,EAAE,cAAc,EAAE,UAAU,EAAE,MAAM,CAAC,CAAC;IAC/F,CAAC;IACD;;;;OAIG;IACI,MAAM,CAAC,wBAAwB,CAAC,cAAsB,EAAE,UAAkB;QAC/E,MAAM,EAAE,GAAG,UAAU,CAAC;QACtB,MAAM,KAAK,GAAG,CAAC,EAAE,GAAG,IAAI,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC;QACnC,IAAI,KAAK,GAAG,CAAC,CAAC,CAAC,EAAE,GAAG,EAAE,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,EAAE,CAAC,GAAG,EAAE,GAAG,cAAc,CAAC,CAAC;QAC/D,IAAI,KAAK,GAAG,CAAC,GAAG;YACd,KAAK,GAAG,CAAC,GAAG,CAAC;QACf,IAAI,KAAK,GAAG,GAAG;YACb,KAAK,GAAG,GAAG,CAAC;QACd,MAAM,KAAK,GAAG,aAAK,CAAC,gBAAgB,CAAC,GAAG,CAAC,CAAC;QAE1C,OAAO,IAAI,CAAC,IAAI,CAAC,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,IAAI,CAAC,KAAK,CAAC,GAAG,EAAE,GAAG,KAAK,CAAC,CAAC,CAAC;IACpE,CAAC;IACe,YAAY,CAAC,WAAmB;QAC9C,0CAA0C;QAC1C,qCAAqC;QACrC,IAAI,CAAC,eAAe,IAAI,WAAW,CAAC;QACpC,IAAI,CAAC,eAAe,IAAI,WAAW,CAAC;QACpC,KAAK,CAAC,YAAY,CAAC,WAAW,CAAC,CAAC;IAClC,CAAC;IACD;;OAEG;IACI,MAAM,CAAC,kCAAkC,CAAC,cAAsB,EAAE,cAAsB,EAAE,YAAoB,MAAM,EACzH,2BAAmC,CAAC;QACpC,MAAM,CAAC,GAAG,cAAc,CAAC;QACzB,IAAI,GAAG,GAAG,CAAC,CAAC;QACZ,IAAI,CAAC,EAAE,GAAG,EAAE,GAAG,CAAC;QAEhB,IAAI,EAAE,GAAG,EAAE,GAAG,cAAc,CAAC;QAC7B,IAAI,gBAAgB,GAAG,CAAC,CAAC;QACzB,oCAAoC;QACpC,mBAAmB;QACnB,wCAAwC;QACxC,uHAAuH;QACvH,oHAAoH;QACpH,sGAAsG;QACtG,qDAAqD;QACrD,gEAAgE;QAChE,KAAK,GAAG,GAAG,CAAC,EAAE,GAAG,GAAG,GAAG,EAAE,EAAE,GAAG,EAAE,CAAC;YAC/B,GAAG,GAAG,IAAI,CAAC,wBAAwB,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC;YAC3C,GAAG,GAAG,EAAE,GAAG,EAAE,CAAC;YACd,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,GAAG,GAAG,CAAC,CAAC;YAChC,MAAM,IAAI,GAAG,CAAC,GAAG,CAAC,GAAG,GAAG,GAAG,GAAG,CAAC;YAC/B,MAAM,UAAU,GAAG,EAAE,GAAG,IAAI,GAAG,CAC7B,CAAC,EAAE,GAAG,EAAE,CAAC,GAAG,IAAI,GAAG,CACjB,CAAC,CAAC,EAAE,GAAG,EAAE,CAAC,GAAG,IAAI,GAAG,CAClB,CAAC,CAAC,IAAI,GAAG,IAAI,CAAC,GAAG,IAAI;gBACrB,CAAC,CAAC,MAAM,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC;YAC3B,MAAM,eAAe,GAAG,EAAE,GAAG,UAAU,CAAC;YACxC,EAAE,GAAG,CAAC,cAAc,GAAG,eAAe,CAAC,GAAG,EAAE,CAAC;YAE7C,IAAI,IAAI,CAAC,GAAG,CAAC,cAAc,GAAG,eAAe,CAAC,GAAG,SAAS,EAAE,CAAC;gBAC3D,gBAAgB,EAAE,CAAC;gBACnB,IAAI,gBAAgB,IAAI,wBAAwB;oBAC9C,MAAM;YACV,CAAC;iBAAM,CAAC;gBACN,gBAAgB,GAAG,CAAC,CAAC;YACvB,CAAC;QACH,CAAC;QACD,OAAO,EAAE,CAAC;IACZ,CAAC;IACM,aAAa,CAAC,KAAU;QAC7B,IAAI,KAAK,YAAY,6BAA6B,EAAE,CAAC;YACnD,OAAO,mBAAQ,CAAC,mBAAmB,CAAC,IAAI,CAAC,OAAO,EAAE,KAAK,CAAC,OAAO,CAAC;mBAC3D,mBAAQ,CAAC,mBAAmB,CAAC,IAAI,CAAC,WAAW,EAAE,KAAK,CAAC,WAAW,CAAC;mBACjE,mBAAQ,CAAC,mBAAmB,CAAC,IAAI,CAAC,eAAe,EAAE,KAAK,CAAC,eAAe,CAAC;mBACzE,mBAAQ,CAAC,mBAAmB,CAAC,IAAI,CAAC,eAAe,EAAE,KAAK,CAAC,eAAe,CAAC,CAAC;QACjF,CAAC;QACD,OAAO,KAAK,CAAC;IACf,CAAC;IACD;;;;;;;;;OASG;IACI,2CAA2C,CAAC,CAAS;QAC1D,MAAM,EAAE,GAAG,MAAM,CAAC;QAClB,MAAM,EAAE,GAAG,MAAM,CAAC;QAClB,MAAM,EAAE,GAAG,OAAO,CAAC;QACnB,MAAM,EAAE,GAAG,QAAQ,CAAC;QACpB,MAAM,CAAC,GAAG,IAAI,CAAC,OAAO,CAAC;QACvB,MAAM,IAAI,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;QACnC,MAAM,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,GAAG,IAAI,GAAG,CAAC,EAAE,GAAG,IAAI,GAAG,CAAC,EAAE,GAAG,IAAI,GAAG,CAAC,EAAE,GAAG,IAAI,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;QAC1E,OAAO,CAAC,CAAC;IACX,CAAC;CACF;AAjHD,sEAiHC","sourcesContent":["/*---------------------------------------------------------------------------------------------\n* Copyright (c) Bentley Systems, Incorporated. All rights reserved.\n* See LICENSE.md in the project root for license terms and full copyright notice.\n*--------------------------------------------------------------------------------------------*/\n\n/** @packageDocumentation\n * @module Curve\n */\nimport { Geometry } from \"../../Geometry\";\nimport { Angle } from \"../../geometry3d/Angle\";\nimport { CubicEvaluator } from \"./CubicEvaluator\";\n/**\n * AustralianRailCorp spiral (also known as New South Wales spiral)\n * * The ultimate curve is a cubic `y = m * x^3`.\n * * `m` is a constant throughout the curve.\n * * Computation of m from the R and L is an complicated sequence, but is only done at construction time.\n * @internal\n */\nexport class AustralianRailCorpXYEvaluator extends CubicEvaluator {\n  private _nominalLength1: number;\n  private _nominalRadius1: number;\n  private constructor(nominalLength1: number, nominalRadius1: number, axisLength: number, cubicM: number) {\n    super(axisLength, cubicM);\n    this._nominalLength1 = nominalLength1;\n    this._nominalRadius1 = nominalRadius1;\n  }\n  public get nominalLength1() { return this._nominalLength1; }\n  public get nominalRadius1() { return this._nominalRadius1; }\n  public clone(): AustralianRailCorpXYEvaluator { return new AustralianRailCorpXYEvaluator(this._nominalLength1, this._nominalRadius1, this._axisLength, this._cubicM); }\n  public static create(nominalLength1: number, nominalRadius1: number): AustralianRailCorpXYEvaluator | undefined {\n    const axisLength = AustralianRailCorpXYEvaluator.radiusAndNominalLengthToAxisLength(nominalRadius1, nominalLength1);\n    const phi = this.radiusAndAxisLengthToPhi(nominalRadius1, axisLength);\n    const xc2 = axisLength * axisLength;\n    const cubicM = Math.tan(phi) / (3.0 * xc2);\n    return new AustralianRailCorpXYEvaluator(nominalLength1, nominalRadius1, axisLength, cubicM);\n  }\n  /**\n   * Compute the phi constant for AustralianRail spiral with given end radius and length along axis.\n   * @param nominalRadius1\n   * @param axisLength\n   */\n  public static radiusAndAxisLengthToPhi(nominalRadius1: number, axisLength: number): number {\n    const xc = axisLength;\n    const expr1 = (2. / Math.sqrt(3.));\n    let expr2 = (-(3. / 4.) * Math.sqrt(3.) * xc / nominalRadius1);\n    if (expr2 < -1.0)\n      expr2 = -1.0;\n    if (expr2 > 1.0)\n      expr2 = 1.0;\n    const expr3 = Angle.degreesToRadians(240);\n\n    return Math.asin(expr1 * Math.cos(Math.acos(expr2) / 3. + expr3));\n  }\n  public override scaleInPlace(scaleFactor: number) {\n    // apply the scale factor to all contents.\n    // all distances scale directly . . .\n    this._nominalLength1 *= scaleFactor;\n    this._nominalRadius1 *= scaleFactor;\n    super.scaleInPlace(scaleFactor);\n  }\n  /** Compute length along axis for AustralianRail spiral nominal radius and length.\n   *\n   */\n  public static radiusAndNominalLengthToAxisLength(nominalRadius1: number, nominalLength1: number, tolerance: number = 1.0e-5,\n    requiredConvergenceCount: number = 2) {\n    const R = nominalRadius1;\n    let idx = 0;\n    let m, phi, xc2;\n\n    let xc = .7 * nominalLength1;\n    let convergenceCount = 0;\n    // remark: This converges quickly --\n    // for L=100, R=400\n    //   ** full precision at 7th iteration.\n    //   ** classic tolerance 1.0e-5 (7 digits from L) with requiredConvergenceCount = 1 gives 11 digits after 3 iterations\n    //   ** each iteration adds about 2 digits.   This is quite good for a successive replacement without derivative !!!\n    //   ** Unanswerable question: If this is only done once and reused over all evaluations, do you want:\n    //       * run the 7 iterations to get full precision\n    //       * stop with the classic tolerance to get compatibility?\n    for (idx = 0; idx < 100; ++idx) {\n      phi = this.radiusAndAxisLengthToPhi(R, xc);\n      xc2 = xc * xc;\n      m = Math.tan(phi) / (3.0 * xc2);\n      const m2x4 = m * m * xc2 * xc2;\n      const correction = xc * m2x4 * (\n        (9. / 10) + m2x4 * (\n          -(9. / 8.) + m2x4 * (\n            +(729. / 208.) + m2x4 *\n            -(32805. / 2176.))));\n      const correctedLength = xc + correction;\n      xc = (nominalLength1 / correctedLength) * xc;\n\n      if (Math.abs(nominalLength1 - correctedLength) < tolerance) {\n        convergenceCount++;\n        if (convergenceCount >= requiredConvergenceCount)\n          break;\n      } else {\n        convergenceCount = 0;\n      }\n    }\n    return xc;\n  }\n  public isAlmostEqual(other: any): boolean {\n    if (other instanceof AustralianRailCorpXYEvaluator) {\n      return Geometry.isAlmostEqualNumber(this._cubicM, other._cubicM)\n        && Geometry.isAlmostEqualNumber(this._axisLength, other._axisLength)\n        && Geometry.isAlmostEqualNumber(this._nominalLength1, other._nominalLength1)\n        && Geometry.isAlmostEqualNumber(this._nominalRadius1, other._nominalRadius1);\n    }\n    return false;\n  }\n  /**\n   * Return a (quite good approximation) of fraction along x axis for given distance along spiral.\n   * * The AustralianRailSpiral has a supporting power series to approximately map distance along the spiral to an x coordinate.\n   * * The `xToFraction(x)` method quickly (with a single divide) converts this x to fraction used fro this.fractionToX (fraction), this.fractionToY(fraction) etc to get coordinates and derivatives.\n   * * The x-to-distance relation is not as precise as the CurvePrimitive method moveSignedDistanceFromFraction.\n   * * It is supported here for users interested in replicating the AustralianRail distance mapping rather than the more accurate CurvePrimitive measurements.\n   * * Round tripping distance through (a) distanceAlongSpiralToAustralianApproximateX, (b) xToFraction, and (c) curveLengthBetweenFractions has\n   *   * 10 digit accuracy for L/R = 4, 12 digit accuracy for L/R = 10\n   * @param s distance along the axis.\n   */\n  public distanceAlongSpiralToAustralianApproximateX(s: number): number {\n    const a1 = 0.9000;\n    const a2 = 5.1750;\n    const a3 = 43.1948;\n    const a4 = 426.0564;\n    const m = this._cubicM;\n    const m2s4 = m * m * s * s * s * s;\n    const x = s * (1.0 - m2s4 * (a1 - m2s4 * (a2 - m2s4 * (a3 - m2s4 * a4))));\n    return x;\n  }\n}\n"]}
+{"version":3,"file":"AustralianRailCorpXYEvaluator.js","sourceRoot":"","sources":["../../../../src/curve/spiral/AustralianRailCorpXYEvaluator.ts"],"names":[],"mappings":";AAAA;;;+FAG+F;;;AAE/F;;GAEG;AACH,6CAA0C;AAC1C,kDAA+C;AAC/C,qDAAkD;AAElD;;;;;;GAMG;AACH,MAAa,6BAA8B,SAAQ,+BAAc;IACvD,eAAe,CAAS;IACxB,eAAe,CAAS;IAChC,YAAoB,cAAsB,EAAE,cAAsB,EAAE,UAAkB,EAAE,MAAc;QACpG,KAAK,CAAC,UAAU,EAAE,MAAM,CAAC,CAAC;QAC1B,IAAI,CAAC,eAAe,GAAG,cAAc,CAAC;QACtC,IAAI,CAAC,eAAe,GAAG,cAAc,CAAC;IACxC,CAAC;IACD,IAAW,cAAc;QACvB,OAAO,IAAI,CAAC,eAAe,CAAC;IAC9B,CAAC;IACD,IAAW,cAAc;QACvB,OAAO,IAAI,CAAC,eAAe,CAAC;IAC9B,CAAC;IACM,KAAK;QACV,OAAO,IAAI,6BAA6B,CAAC,IAAI,CAAC,eAAe,EAAE,IAAI,CAAC,eAAe,EAAE,IAAI,CAAC,WAAW,EAAE,IAAI,CAAC,OAAO,CAAC,CAAC;IACvH,CAAC;IACM,MAAM,CAAC,MAAM,CAAC,cAAsB,EAAE,cAAsB;QACjE,MAAM,UAAU,GAAG,6BAA6B,CAAC,kCAAkC,CAAC,cAAc,EAAE,cAAc,CAAC,CAAC;QACpH,MAAM,GAAG,GAAG,IAAI,CAAC,wBAAwB,CAAC,cAAc,EAAE,UAAU,CAAC,CAAC;QACtE,MAAM,GAAG,GAAG,UAAU,GAAG,UAAU,CAAC;QACpC,MAAM,MAAM,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,GAAG,GAAG,CAAC,CAAC;QAC3C,OAAO,IAAI,6BAA6B,CAAC,cAAc,EAAE,cAAc,EAAE,UAAU,EAAE,MAAM,CAAC,CAAC;IAC/F,CAAC;IACD,sGAAsG;IAC/F,MAAM,CAAC,wBAAwB,CAAC,cAAsB,EAAE,UAAkB;QAC/E,MAAM,EAAE,GAAG,UAAU,CAAC;QACtB,MAAM,KAAK,GAAG,CAAC,EAAE,GAAG,IAAI,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC;QACnC,IAAI,KAAK,GAAG,CAAC,CAAC,CAAC,EAAE,GAAG,EAAE,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,EAAE,CAAC,GAAG,EAAE,GAAG,cAAc,CAAC,CAAC;QAC/D,IAAI,KAAK,GAAG,CAAC,GAAG;YACd,KAAK,GAAG,CAAC,GAAG,CAAC;QACf,IAAI,KAAK,GAAG,GAAG;YACb,KAAK,GAAG,GAAG,CAAC;QACd,MAAM,KAAK,GAAG,aAAK,CAAC,gBAAgB,CAAC,GAAG,CAAC,CAAC;QAC1C,OAAO,IAAI,CAAC,IAAI,CAAC,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,IAAI,CAAC,KAAK,CAAC,GAAG,EAAE,GAAG,KAAK,CAAC,CAAC,CAAC;IACpE,CAAC;IACe,YAAY,CAAC,WAAmB;QAC9C,uEAAuE;QACvE,IAAI,CAAC,eAAe,IAAI,WAAW,CAAC;QACpC,IAAI,CAAC,eAAe,IAAI,WAAW,CAAC;QACpC,KAAK,CAAC,YAAY,CAAC,WAAW,CAAC,CAAC;IAClC,CAAC;IACD,qFAAqF;IAC9E,MAAM,CAAC,kCAAkC,CAC9C,cAAsB,EAAE,cAAsB,EAAE,YAAoB,MAAM,EAAE,2BAAmC,CAAC;QAEhH,MAAM,CAAC,GAAG,cAAc,CAAC;QACzB,IAAI,GAAG,GAAG,CAAC,CAAC;QACZ,IAAI,CAAC,EAAE,GAAG,EAAE,GAAG,CAAC;QAChB,IAAI,EAAE,GAAG,EAAE,GAAG,cAAc,CAAC;QAC7B,IAAI,gBAAgB,GAAG,CAAC,CAAC;QACzB,iCAAiC;QACjC,mBAAmB;QACnB,wCAAwC;QACxC,uHAAuH;QACvH,8GAA8G;QAC9G,sGAAsG;QACtG,qDAAqD;QACrD,gEAAgE;QAChE,KAAK,GAAG,GAAG,CAAC,EAAE,GAAG,GAAG,GAAG,EAAE,EAAE,GAAG,EAAE,CAAC;YAC/B,GAAG,GAAG,IAAI,CAAC,wBAAwB,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC;YAC3C,GAAG,GAAG,EAAE,GAAG,EAAE,CAAC;YACd,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,GAAG,GAAG,CAAC,CAAC;YAChC,MAAM,IAAI,GAAG,CAAC,GAAG,CAAC,GAAG,GAAG,GAAG,GAAG,CAAC;YAC/B,MAAM,UAAU,GAAG,EAAE,GAAG,IAAI,GAAG,CAAC,CAAC,EAAE,GAAG,EAAE,CAAC,GAAG,IAAI,GAAG,CAAC,CAAC,CAAC,EAAE,GAAG,EAAE,CAAC,GAAG,IAAI,GAAG,CAAC,CAAC,CAAC,IAAI,GAAG,IAAI,CAAC,GAAG,IAAI,GAAG,CAAC,CAAC,MAAM,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC;YACtH,MAAM,eAAe,GAAG,EAAE,GAAG,UAAU,CAAC;YACxC,EAAE,GAAG,CAAC,cAAc,GAAG,eAAe,CAAC,GAAG,EAAE,CAAC;YAC7C,IAAI,IAAI,CAAC,GAAG,CAAC,cAAc,GAAG,eAAe,CAAC,GAAG,SAAS,EAAE,CAAC;gBAC3D,gBAAgB,EAAE,CAAC;gBACnB,IAAI,gBAAgB,IAAI,wBAAwB;oBAC9C,MAAM;YACV,CAAC;iBAAM,CAAC;gBACN,gBAAgB,GAAG,CAAC,CAAC;YACvB,CAAC;QACH,CAAC;QACD,OAAO,EAAE,CAAC;IACZ,CAAC;IACM,aAAa,CAAC,KAAU;QAC7B,IAAI,KAAK,YAAY,6BAA6B,EAAE,CAAC;YACnD,OAAO,mBAAQ,CAAC,mBAAmB,CAAC,IAAI,CAAC,OAAO,EAAE,KAAK,CAAC,OAAO,CAAC;mBAC3D,mBAAQ,CAAC,mBAAmB,CAAC,IAAI,CAAC,WAAW,EAAE,KAAK,CAAC,WAAW,CAAC;mBACjE,mBAAQ,CAAC,mBAAmB,CAAC,IAAI,CAAC,eAAe,EAAE,KAAK,CAAC,eAAe,CAAC;mBACzE,mBAAQ,CAAC,mBAAmB,CAAC,IAAI,CAAC,eAAe,EAAE,KAAK,CAAC,eAAe,CAAC,CAAC;QACjF,CAAC;QACD,OAAO,KAAK,CAAC;IACf,CAAC;IACD;;;;;;;;;;;;OAYG;IACI,2CAA2C,CAAC,CAAS;QAC1D,MAAM,EAAE,GAAG,MAAM,CAAC;QAClB,MAAM,EAAE,GAAG,MAAM,CAAC;QAClB,MAAM,EAAE,GAAG,OAAO,CAAC;QACnB,MAAM,EAAE,GAAG,QAAQ,CAAC;QACpB,MAAM,CAAC,GAAG,IAAI,CAAC,OAAO,CAAC;QACvB,MAAM,IAAI,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;QACnC,MAAM,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,GAAG,IAAI,GAAG,CAAC,EAAE,GAAG,IAAI,GAAG,CAAC,EAAE,GAAG,IAAI,GAAG,CAAC,EAAE,GAAG,IAAI,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;QAC1E,OAAO,CAAC,CAAC;IACX,CAAC;CACF;AA7GD,sEA6GC","sourcesContent":["/*---------------------------------------------------------------------------------------------\n* Copyright (c) Bentley Systems, Incorporated. All rights reserved.\n* See LICENSE.md in the project root for license terms and full copyright notice.\n*--------------------------------------------------------------------------------------------*/\n\n/** @packageDocumentation\n * @module Curve\n */\nimport { Geometry } from \"../../Geometry\";\nimport { Angle } from \"../../geometry3d/Angle\";\nimport { CubicEvaluator } from \"./CubicEvaluator\";\n\n/**\n * AustralianRailCorp spiral (also known as New South Wales spiral)\n * * The ultimate curve is a cubic `y = m * x^3`.\n * * `m` is a constant throughout the curve.\n * * Computation of m from the R and L is an complicated sequence, but is only done at construction time.\n * @internal\n */\nexport class AustralianRailCorpXYEvaluator extends CubicEvaluator {\n  private _nominalLength1: number;\n  private _nominalRadius1: number;\n  private constructor(nominalLength1: number, nominalRadius1: number, axisLength: number, cubicM: number) {\n    super(axisLength, cubicM);\n    this._nominalLength1 = nominalLength1;\n    this._nominalRadius1 = nominalRadius1;\n  }\n  public get nominalLength1() {\n    return this._nominalLength1;\n  }\n  public get nominalRadius1() {\n    return this._nominalRadius1;\n  }\n  public clone(): AustralianRailCorpXYEvaluator {\n    return new AustralianRailCorpXYEvaluator(this._nominalLength1, this._nominalRadius1, this._axisLength, this._cubicM);\n  }\n  public static create(nominalLength1: number, nominalRadius1: number): AustralianRailCorpXYEvaluator | undefined {\n    const axisLength = AustralianRailCorpXYEvaluator.radiusAndNominalLengthToAxisLength(nominalRadius1, nominalLength1);\n    const phi = this.radiusAndAxisLengthToPhi(nominalRadius1, axisLength);\n    const xc2 = axisLength * axisLength;\n    const cubicM = Math.tan(phi) / (3.0 * xc2);\n    return new AustralianRailCorpXYEvaluator(nominalLength1, nominalRadius1, axisLength, cubicM);\n  }\n  /** Compute the phi constant for AustralianRail spiral with given end radius and length along axis. */\n  public static radiusAndAxisLengthToPhi(nominalRadius1: number, axisLength: number): number {\n    const xc = axisLength;\n    const expr1 = (2. / Math.sqrt(3.));\n    let expr2 = (-(3. / 4.) * Math.sqrt(3.) * xc / nominalRadius1);\n    if (expr2 < -1.0)\n      expr2 = -1.0;\n    if (expr2 > 1.0)\n      expr2 = 1.0;\n    const expr3 = Angle.degreesToRadians(240);\n    return Math.asin(expr1 * Math.cos(Math.acos(expr2) / 3. + expr3));\n  }\n  public override scaleInPlace(scaleFactor: number): void {\n    // apply the scale factor to all contents; all distances scale directly\n    this._nominalLength1 *= scaleFactor;\n    this._nominalRadius1 *= scaleFactor;\n    super.scaleInPlace(scaleFactor);\n  }\n  /** Compute length along axis for AustralianRail spiral nominal radius and length. */\n  public static radiusAndNominalLengthToAxisLength(\n    nominalRadius1: number, nominalLength1: number, tolerance: number = 1.0e-5, requiredConvergenceCount: number = 2,\n  ): number {\n    const R = nominalRadius1;\n    let idx = 0;\n    let m, phi, xc2;\n    let xc = .7 * nominalLength1;\n    let convergenceCount = 0;\n    // remark: This converges quickly\n    // for L=100, R=400\n    //   ** full precision at 7th iteration.\n    //   ** classic tolerance 1.0e-5 (7 digits from L) with requiredConvergenceCount = 1 gives 11 digits after 3 iterations\n    //   ** each iteration adds about 2 digits. This is quite good for a successive replacement without derivative\n    //   ** Unanswerable question: If this is only done once and reused over all evaluations, do you want:\n    //       * run the 7 iterations to get full precision\n    //       * stop with the classic tolerance to get compatibility?\n    for (idx = 0; idx < 100; ++idx) {\n      phi = this.radiusAndAxisLengthToPhi(R, xc);\n      xc2 = xc * xc;\n      m = Math.tan(phi) / (3.0 * xc2);\n      const m2x4 = m * m * xc2 * xc2;\n      const correction = xc * m2x4 * ((9. / 10) + m2x4 * (-(9. / 8.) + m2x4 * (+(729. / 208.) + m2x4 * -(32805. / 2176.))));\n      const correctedLength = xc + correction;\n      xc = (nominalLength1 / correctedLength) * xc;\n      if (Math.abs(nominalLength1 - correctedLength) < tolerance) {\n        convergenceCount++;\n        if (convergenceCount >= requiredConvergenceCount)\n          break;\n      } else {\n        convergenceCount = 0;\n      }\n    }\n    return xc;\n  }\n  public isAlmostEqual(other: any): boolean {\n    if (other instanceof AustralianRailCorpXYEvaluator) {\n      return Geometry.isAlmostEqualNumber(this._cubicM, other._cubicM)\n        && Geometry.isAlmostEqualNumber(this._axisLength, other._axisLength)\n        && Geometry.isAlmostEqualNumber(this._nominalLength1, other._nominalLength1)\n        && Geometry.isAlmostEqualNumber(this._nominalRadius1, other._nominalRadius1);\n    }\n    return false;\n  }\n  /**\n   * Return a (quite good approximation) of fraction along x axis for given distance along spiral.\n   * * The AustralianRailSpiral has a supporting power series to approximately map distance along the spiral to\n   * an x coordinate.\n   * * The `xToFraction(x)` method quickly (with a single divide) converts this x to fraction used for\n   * this.fractionToX (fraction), this.fractionToY(fraction) etc to get coordinates and derivatives.\n   * * The x-to-distance relation is not as precise as the CurvePrimitive method moveSignedDistanceFromFraction.\n   * * It is supported here for users interested in replicating the AustralianRail distance mapping rather than the\n   * more accurate CurvePrimitive measurements.\n   * * Round tripping distance through (a) distanceAlongSpiralToAustralianApproximateX, (b) xToFraction, and\n   * (c) curveLengthBetweenFractions has 10 digit accuracy for L/R = 4, 12 digit accuracy for L/R = 10.\n   * @param s distance along the axis.\n   */\n  public distanceAlongSpiralToAustralianApproximateX(s: number): number {\n    const a1 = 0.9000;\n    const a2 = 5.1750;\n    const a3 = 43.1948;\n    const a4 = 426.0564;\n    const m = this._cubicM;\n    const m2s4 = m * m * s * s * s * s;\n    const x = s * (1.0 - m2s4 * (a1 - m2s4 * (a2 - m2s4 * (a3 - m2s4 * a4))));\n    return x;\n  }\n}\n"]}

package/lib/cjs/curve/spiral/ClothoidSeries.d.ts

@@ -1,22 +1,26 @@
 import { XYCurveEvaluator } from "./XYCurveEvaluator";
 /**
- * Methods to evaluate caller-specified number of terms of the x and y series for a clothoid.
- * Each instance has
- * * Number of x and y terms to use.
- * * constant for theta=c * x * x
- *    * This value is c=1/(2 R L)  for curve length L measured from inflection to point with radius R.
+ * Methods to evaluate caller-specified number of terms of the x and y Taylor series for a clothoid.
+ * Each instance has:
+ * * number of x and y terms to use.
+ * * constant for theta = c*x*x.
+ *    * This value is c = 1/(2*R*L) for curve length L measured from inflection to point with radius R.
  * @internal
  */
 export declare class ClothoidSeriesRLEvaluator extends XYCurveEvaluator {
+    /** Number of terms to use in x Taylor series. */
     numXTerms: number;
+    /** Number of terms to use in y Taylor series. */
     numYTerms: number;
+    /** Constant c = 1/(2*R*L) in theta = c*s*s. */
     constantDiv2LR: number;
+    /** The nominal curve length. */
     nominalLength1: number;
     constructor(nominalLength1: number, constantDiv2LR: number, numXTerms?: number, numYTerms?: number);
     /** Return a deep clone. */
     clone(): ClothoidSeriesRLEvaluator;
     scaleInPlace(scaleFactor: number): void;
-    /** Member by member matchup ... */
+    /** Member by member matchup. */
     isAlmostEqual(other: any): boolean;
     /**
      * Evaluate the X series at a nominal distance along the curve.

package/lib/cjs/curve/spiral/ClothoidSeries.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"ClothoidSeries.d.ts","sourceRoot":"","sources":["../../../../src/curve/spiral/ClothoidSeries.ts"],"names":[],"mappings":"AAUA,OAAO,EAAE,gBAAgB,EAAE,MAAM,oBAAoB,CAAC;AAEtD;;;;;;;GAOG;AACH,qBAAa,yBAA0B,SAAQ,gBAAgB;IACtD,SAAS,EAAE,MAAM,CAAC;IAClB,SAAS,EAAE,MAAM,CAAC;IAClB,cAAc,EAAE,MAAM,CAAC;IACvB,cAAc,EAAE,MAAM,CAAC;gBACX,cAAc,EAAE,MAAM,EAAE,cAAc,EAAE,MAAM,EAAE,SAAS,GAAE,MAAU,EAAE,SAAS,GAAE,MAAU;IAO/G,2BAA2B;IACpB,KAAK,IAAI,yBAAyB;IAGlC,YAAY,CAAC,WAAW,EAAE,MAAM;IAIvC,mCAAmC;IAC5B,aAAa,CAAC,KAAK,EAAE,GAAG,GAAG,OAAO;IASzC;;;OAGG;IACI,WAAW,CAAC,QAAQ,EAAE,MAAM,GAAG,MAAM;IAC5C;;;OAGG;IACI,WAAW,CAAC,QAAQ,EAAE,MAAM,GAAG,MAAM;IAE5C;;;OAGG;IACI,YAAY,CAAC,QAAQ,EAAE,MAAM,GAAG,MAAM;IAC7C;;;OAGG;IACI,YAAY,CAAC,QAAQ,EAAE,MAAM,GAAG,MAAM;IAE7C;;;OAGG;IACI,aAAa,CAAC,QAAQ,EAAE,MAAM,GAAG,MAAM;IAC9C;;;OAGG;IACI,aAAa,CAAC,QAAQ,EAAE,MAAM,GAAG,MAAM;IAE9C;;;;OAIG;IACI,aAAa,CAAC,QAAQ,EAAE,MAAM,EAAE,QAAQ,EAAE,MAAM,GAAG,MAAM;IAuBzD,aAAa,CAAC,QAAQ,EAAE,MAAM,EAAE,QAAQ,EAAE,MAAM,GAAG,MAAM;IAsBzD,cAAc,CAAC,QAAQ,EAAE,MAAM,EAAE,QAAQ,EAAE,MAAM,GAAG,MAAM;IAsB1D,cAAc,CAAC,QAAQ,EAAE,MAAM,EAAE,QAAQ,EAAE,MAAM,GAAG,MAAM;IAuB1D,eAAe,CAAC,QAAQ,EAAE,MAAM,EAAE,QAAQ,EAAE,MAAM,GAAG,MAAM;IAU3D,eAAe,CAAC,QAAQ,EAAE,MAAM,EAAE,QAAQ,EAAE,MAAM,GAAG,MAAM;IAU3D,aAAa,CAAC,QAAQ,EAAE,MAAM,GAAG,MAAM;IAYvC,aAAa,CAAC,QAAQ,EAAE,MAAM,GAAG,MAAM;IAcvC,WAAW,CAAC,CAAC,EAAE,MAAM,GAAG,MAAM,GAAG,SAAS;CASlD"}
+{"version":3,"file":"ClothoidSeries.d.ts","sourceRoot":"","sources":["../../../../src/curve/spiral/ClothoidSeries.ts"],"names":[],"mappings":"AAUA,OAAO,EAAE,gBAAgB,EAAE,MAAM,oBAAoB,CAAC;AAGtD;;;;;;;GAOG;AACH,qBAAa,yBAA0B,SAAQ,gBAAgB;IAC7D,iDAAiD;IAC1C,SAAS,EAAE,MAAM,CAAC;IACzB,iDAAiD;IAC1C,SAAS,EAAE,MAAM,CAAC;IACzB,+CAA+C;IACxC,cAAc,EAAE,MAAM,CAAC;IAC9B,gCAAgC;IACzB,cAAc,EAAE,MAAM,CAAC;gBACX,cAAc,EAAE,MAAM,EAAE,cAAc,EAAE,MAAM,EAAE,SAAS,GAAE,MAAU,EAAE,SAAS,GAAE,MAAU;IAO/G,2BAA2B;IACpB,KAAK,IAAI,yBAAyB;IAGlC,YAAY,CAAC,WAAW,EAAE,MAAM,GAAG,IAAI;IAI9C,gCAAgC;IACzB,aAAa,CAAC,KAAK,EAAE,GAAG,GAAG,OAAO;IASzC;;;OAGG;IACI,WAAW,CAAC,QAAQ,EAAE,MAAM,GAAG,MAAM;IAG5C;;;OAGG;IACI,WAAW,CAAC,QAAQ,EAAE,MAAM,GAAG,MAAM;IAG5C;;;OAGG;IACI,YAAY,CAAC,QAAQ,EAAE,MAAM,GAAG,MAAM;IAG7C;;;OAGG;IACI,YAAY,CAAC,QAAQ,EAAE,MAAM,GAAG,MAAM;IAG7C;;;OAGG;IACI,aAAa,CAAC,QAAQ,EAAE,MAAM,GAAG,MAAM;IAG9C;;;OAGG;IACI,aAAa,CAAC,QAAQ,EAAE,MAAM,GAAG,MAAM;IAG9C;;;;OAIG;IACI,aAAa,CAAC,QAAQ,EAAE,MAAM,EAAE,QAAQ,EAAE,MAAM,GAAG,MAAM;IAuBzD,aAAa,CAAC,QAAQ,EAAE,MAAM,EAAE,QAAQ,EAAE,MAAM,GAAG,MAAM;IAuBzD,cAAc,CAAC,QAAQ,EAAE,MAAM,EAAE,QAAQ,EAAE,MAAM,GAAG,MAAM;IAsB1D,cAAc,CAAC,QAAQ,EAAE,MAAM,EAAE,QAAQ,EAAE,MAAM,GAAG,MAAM;IAsB1D,eAAe,CAAC,QAAQ,EAAE,MAAM,EAAE,QAAQ,EAAE,MAAM,GAAG,MAAM;IAS3D,eAAe,CAAC,QAAQ,EAAE,MAAM,EAAE,QAAQ,EAAE,MAAM,GAAG,MAAM;IAS3D,aAAa,CAAC,QAAQ,EAAE,MAAM,GAAG,MAAM;IAYvC,aAAa,CAAC,QAAQ,EAAE,MAAM,GAAG,MAAM;IAavC,WAAW,CAAC,CAAC,EAAE,MAAM,GAAG,MAAM,GAAG,SAAS;CASlD"}

package/lib/cjs/curve/spiral/ClothoidSeries.js

@@ -12,17 +12,21 @@ const Geometry_1 = require("../../Geometry");
 const XYCurveEvaluator_1 = require("./XYCurveEvaluator");
 const Newton_1 = require("../../numerics/Newton");
 /**
- * Methods to evaluate caller-specified number of terms of the x and y series for a clothoid.
- * Each instance has
- * * Number of x and y terms to use.
- * * constant for theta=c * x * x
- *    * This value is c=1/(2 R L)  for curve length L measured from inflection to point with radius R.
+ * Methods to evaluate caller-specified number of terms of the x and y Taylor series for a clothoid.
+ * Each instance has:
+ * * number of x and y terms to use.
+ * * constant for theta = c*x*x.
+ *    * This value is c = 1/(2*R*L) for curve length L measured from inflection to point with radius R.
  * @internal
  */
 class ClothoidSeriesRLEvaluator extends XYCurveEvaluator_1.XYCurveEvaluator {
+    /** Number of terms to use in x Taylor series. */
     numXTerms;
+    /** Number of terms to use in y Taylor series. */
     numYTerms;
+    /** Constant c = 1/(2*R*L) in theta = c*s*s. */
     constantDiv2LR;
+    /** The nominal curve length. */
     nominalLength1;
     constructor(nominalLength1, constantDiv2LR, numXTerms = 4, numYTerms = 4) {
         super();
@@ -39,7 +43,7 @@ class ClothoidSeriesRLEvaluator extends XYCurveEvaluator_1.XYCurveEvaluator {
         this.nominalLength1 *= scaleFactor;
         this.constantDiv2LR /= (scaleFactor * scaleFactor);
     }
-    /** Member by member matchup ... */
+    /** Member by member matchup. */
     isAlmostEqual(other) {
         if (other instanceof ClothoidSeriesRLEvaluator) {
             return this.numXTerms === other.numXTerms
@@ -53,43 +57,55 @@ class ClothoidSeriesRLEvaluator extends XYCurveEvaluator_1.XYCurveEvaluator {
      * Evaluate the X series at a nominal distance along the curve.
      * @param fraction fractional position along the curve.
      */
-    fractionToX(fraction) { return this.fractionToXGo(fraction, this.numXTerms); }
+    fractionToX(fraction) {
+        return this.fractionToXGo(fraction, this.numXTerms);
+    }
     /**
      * Evaluate the Y series at a nominal distance along the curve.
      * @param fraction fractional position along the curve.
      */
-    fractionToY(fraction) { return this.fractionToYGo(fraction, this.numYTerms); }
+    fractionToY(fraction) {
+        return this.fractionToYGo(fraction, this.numYTerms);
+    }
     /**
      * Evaluate the derivative of the X series at a nominal distance along the curve.
      * @param fraction fractional position along the curve.
      */
-    fractionToDX(fraction) { return this.fractionToDXGo(fraction, this.numXTerms); }
+    fractionToDX(fraction) {
+        return this.fractionToDXGo(fraction, this.numXTerms);
+    }
     /**
      * Evaluate the derivative of the Y series at a nominal distance along the curve.
      * @param fraction fractional position along the curve.
      */
-    fractionToDY(fraction) { return this.fractionToDYGo(fraction, this.numYTerms); }
+    fractionToDY(fraction) {
+        return this.fractionToDYGo(fraction, this.numYTerms);
+    }
     /**
      * Evaluate the derivative of the X series at a nominal distance along the curve.
      * @param fraction fractional position along the curve.
      */
-    fractionToDDX(fraction) { return this.fractionToDDXGo(fraction, this.numXTerms); }
+    fractionToDDX(fraction) {
+        return this.fractionToDDXGo(fraction, this.numXTerms);
+    }
     /**
      * Evaluate the derivative of the Y series at a nominal distance along the curve.
      * @param fraction fractional position along the curve.
      */
-    fractionToDDY(fraction) { return this.fractionToDDYGo(fraction, this.numYTerms); }
+    fractionToDDY(fraction) {
+        return this.fractionToDDYGo(fraction, this.numYTerms);
+    }
     /**
      * Evaluate the X series at a nominal distance along the curve.
      * @param fraction fractional position along the curve.
      * @param numTerms number of terms to use.
      */
     fractionToXGo(fraction, numTerms) {
-        // Write the series for cos (theta)
+        // write the Taylor series for cos(theta)
         // replace theta by s*s*c
         // integrate wrt s
-        //  x = s - s^5 c^4/ 2 + s^9 c^8/(4!) - s^13 c^12 / 6!
-        //  x = s(1 - (s^4 c^2/2) ( 1/5 -s^4 c^2 / (3*4)  ( 1/9 - ....) ) )
+        //  x = s - s^5 c^2/ 5*2! + s^9 c^4 / 9*4! - s^13 c^6 / 13*6! + ...
+        //  x = s(1 - (s^4 c^2/2) ( 1/5 - (s^4 c^2 / 3*4) (1/9 - ...) ) )
         const s = fraction * this.nominalLength1;
         let result = s;
         if (numTerms < 2)
@@ -108,10 +124,11 @@ class ClothoidSeriesRLEvaluator extends XYCurveEvaluator_1.XYCurveEvaluator {
         return result;
     }
     fractionToYGo(fraction, numTerms) {
-        // Write the series for sin (theta)
+        // write the Taylor series for sin(theta)
         // replace theta by s*s*c
         // integrate wrt s
-        //  x = s^3 c^2/ 3( (1/3)) - s^7 c^6/(3!) ((1/7)) - s^11 c^10 / 5! ((1/9) - ...)
+        // y = s^3 c / 3 - s^7 c^3 / 7*3! + s^11 c^5 / 11*5! - s^15 c^7 / 15*7! + ...
+        // y = s^3 c ( 1/3 - s^4 c^2/ 3! ( (1/7) - (s^4 c^2 / 4*5) (1/11 - ...) ) )
         const s = fraction * this.nominalLength1;
         const q1 = s * s * this.constantDiv2LR;
         let result = q1 * s / 3;
@@ -130,7 +147,7 @@ class ClothoidSeriesRLEvaluator extends XYCurveEvaluator_1.XYCurveEvaluator {
         return result;
     }
     fractionToDXGo(fraction, numTerms) {
-        // Yes -- this does happen during derivatives of cosines with more than 0 terms !!
+        // yes -- this does happen during derivatives of cosines with more than 0 terms
         if (numTerms <= 0)
             return 0;
         // dX = 1 - s^4c^2/2 + s^8 c^4 / 4! -
@@ -185,7 +202,7 @@ class ClothoidSeriesRLEvaluator extends XYCurveEvaluator_1.XYCurveEvaluator {
     fractionToDDYGo(fraction, numTerms) {
         // DY is "sine"
         // DDY is "cosine" series times chain rule dTheta/ds = 2 * s * this.constantDivLR
-        // BUT  .... derivative of the cosine series leading term is zero ... use one less term!
+        // BUT  .... derivative of the cosine series leading term is zero ... use one less term
         const s = fraction * this.nominalLength1;
         const dTheta = 2 * this.constantDiv2LR * s;
         const cosine = this.fractionToDXGo(fraction, numTerms);

package/lib/cjs/curve/spiral/ClothoidSeries.js.map

@@ -1 +1 @@
-{"version":3,"file":"ClothoidSeries.js","sourceRoot":"","sources":["../../../../src/curve/spiral/ClothoidSeries.ts"],"names":[],"mappings":";AAAA;;;+FAG+F;;;AAE/F;;GAEG;AAEH,6CAA0C;AAC1C,yDAAsD;AACtD,kDAAqD;AACrD;;;;;;;GAOG;AACH,MAAa,yBAA0B,SAAQ,mCAAgB;IACtD,SAAS,CAAS;IAClB,SAAS,CAAS;IAClB,cAAc,CAAS;IACvB,cAAc,CAAS;IAC9B,YAAmB,cAAsB,EAAE,cAAsB,EAAE,YAAoB,CAAC,EAAE,YAAoB,CAAC;QAC7G,KAAK,EAAE,CAAC;QACR,IAAI,CAAC,cAAc,GAAG,cAAc,CAAC;QACrC,IAAI,CAAC,cAAc,GAAG,cAAc,CAAC;QACrC,IAAI,CAAC,SAAS,GAAG,SAAS,CAAC;QAC3B,IAAI,CAAC,SAAS,GAAG,SAAS,CAAC;IAC7B,CAAC;IACD,2BAA2B;IACpB,KAAK;QACV,OAAO,IAAI,yBAAyB,CAAC,IAAI,CAAC,cAAc,EAAE,IAAI,CAAC,cAAc,EAAE,IAAI,CAAC,SAAS,EAAE,IAAI,CAAC,SAAS,CAAC,CAAC;IACjH,CAAC;IACM,YAAY,CAAC,WAAmB;QACrC,IAAI,CAAC,cAAc,IAAI,WAAW,CAAC;QACnC,IAAI,CAAC,cAAc,IAAI,CAAC,WAAW,GAAG,WAAW,CAAC,CAAC;IACrD,CAAC;IACD,mCAAmC;IAC5B,aAAa,CAAC,KAAU;QAC7B,IAAI,KAAK,YAAY,yBAAyB,EAAE,CAAC;YAC/C,OAAO,IAAI,CAAC,SAAS,KAAK,KAAK,CAAC,SAAS;mBACpC,IAAI,CAAC,SAAS,KAAK,KAAK,CAAC,SAAS;mBAClC,mBAAQ,CAAC,mBAAmB,CAAC,IAAI,CAAC,cAAc,EAAE,KAAK,CAAC,cAAc,CAAC;mBACvE,mBAAQ,CAAC,gBAAgB,CAAC,IAAI,CAAC,cAAc,EAAE,KAAK,CAAC,cAAc,CAAC,CAAC;QAC5E,CAAC;QACD,OAAO,KAAK,CAAC;IACf,CAAC;IACD;;;OAGG;IACI,WAAW,CAAC,QAAgB,IAAY,OAAO,IAAI,CAAC,aAAa,CAAC,QAAQ,EAAE,IAAI,CAAC,SAAS,CAAC,CAAC,CAAC,CAAC;IACrG;;;OAGG;IACI,WAAW,CAAC,QAAgB,IAAY,OAAO,IAAI,CAAC,aAAa,CAAC,QAAQ,EAAE,IAAI,CAAC,SAAS,CAAC,CAAC,CAAC,CAAC;IAErG;;;OAGG;IACI,YAAY,CAAC,QAAgB,IAAY,OAAO,IAAI,CAAC,cAAc,CAAC,QAAQ,EAAE,IAAI,CAAC,SAAS,CAAC,CAAC,CAAC,CAAC;IACvG;;;OAGG;IACI,YAAY,CAAC,QAAgB,IAAY,OAAO,IAAI,CAAC,cAAc,CAAC,QAAQ,EAAE,IAAI,CAAC,SAAS,CAAC,CAAC,CAAC,CAAC;IAEvG;;;OAGG;IACI,aAAa,CAAC,QAAgB,IAAY,OAAO,IAAI,CAAC,eAAe,CAAC,QAAQ,EAAE,IAAI,CAAC,SAAS,CAAC,CAAC,CAAC,CAAC;IACzG;;;OAGG;IACI,aAAa,CAAC,QAAgB,IAAY,OAAO,IAAI,CAAC,eAAe,CAAC,QAAQ,EAAE,IAAI,CAAC,SAAS,CAAC,CAAC,CAAC,CAAC;IAEzG;;;;OAIG;IACI,aAAa,CAAC,QAAgB,EAAE,QAAgB;QACrD,mCAAmC;QACnC,yBAAyB;QACzB,kBAAkB;QAClB,sDAAsD;QACtD,mEAAmE;QACnE,MAAM,CAAC,GAAG,QAAQ,GAAG,IAAI,CAAC,cAAc,CAAC;QACzC,IAAI,MAAM,GAAG,CAAC,CAAC;QACf,IAAI,QAAQ,GAAG,CAAC;YACd,OAAO,MAAM,CAAC;QAChB,MAAM,EAAE,GAAG,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,cAAc,CAAC;QACvC,MAAM,IAAI,GAAG,CAAE,EAAE,GAAG,EAAE,CAAC;QACvB,IAAI,KAAK,GAAG,CAAC,CAAC;QACd,IAAI,CAAC,GAAG,CAAC,CAAC;QACV,IAAI,CAAC,GAAG,CAAC,CAAC;QACV,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,QAAQ,EAAE,CAAC,EAAE,EAAE,CAAC;YAClC,KAAK,IAAI,IAAI,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;YAC9B,MAAM,IAAI,KAAK,GAAG,CAAC,CAAC;YACpB,CAAC,IAAI,CAAC,CAAC;YACP,CAAC,IAAI,CAAC,CAAC;QACT,CAAC;QACD,OAAO,MAAM,CAAC;IAChB,CAAC;IACM,aAAa,CAAC,QAAgB,EAAE,QAAgB;QACrD,mCAAmC;QACnC,yBAAyB;QACzB,kBAAkB;QAClB,gFAAgF;QAChF,MAAM,CAAC,GAAG,QAAQ,GAAG,IAAI,CAAC,cAAc,CAAC;QACzC,MAAM,EAAE,GAAG,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,cAAc,CAAC;QACvC,IAAI,MAAM,GAAG,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC;QACxB,IAAI,QAAQ,GAAG,CAAC;YACd,OAAO,MAAM,CAAC;QAChB,MAAM,IAAI,GAAG,CAAE,EAAE,GAAG,EAAE,CAAC;QACvB,IAAI,KAAK,GAAG,EAAE,GAAG,CAAC,CAAC;QACnB,IAAI,CAAC,GAAG,CAAC,CAAC;QACV,IAAI,CAAC,GAAG,CAAC,CAAC;QACV,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,QAAQ,EAAE,CAAC,EAAE,EAAE,CAAC;YAClC,KAAK,IAAI,IAAI,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;YAC9B,MAAM,IAAI,KAAK,GAAG,CAAC,CAAC;YACpB,CAAC,IAAI,CAAC,CAAC;YACP,CAAC,IAAI,CAAC,CAAC;QACT,CAAC;QACD,OAAO,MAAM,CAAC;IAChB,CAAC;IACM,cAAc,CAAC,QAAgB,EAAE,QAAgB;QACtD,kFAAkF;QAClF,IAAI,QAAQ,IAAI,CAAC;YACf,OAAO,CAAC,CAAC;QACX,qCAAqC;QACrC,wCAAwC;QACxC,MAAM,CAAC,GAAG,QAAQ,GAAG,IAAI,CAAC,cAAc,CAAC;QACzC,IAAI,MAAM,GAAG,CAAC,CAAC;QACf,IAAI,QAAQ,GAAG,CAAC,EAAE,CAAC;YACjB,OAAO,MAAM,GAAG,IAAI,CAAC,cAAc,CAAC;QACtC,CAAC;QACD,MAAM,EAAE,GAAG,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,cAAc,CAAC;QACvC,MAAM,IAAI,GAAG,CAAE,EAAE,GAAG,EAAE,CAAC;QACvB,IAAI,KAAK,GAAG,GAAG,CAAC;QAChB,IAAI,CAAC,GAAG,CAAC,CAAC;QACV,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,QAAQ,EAAE,CAAC,EAAE,EAAE,CAAC;YAClC,KAAK,IAAI,IAAI,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;YAC9B,MAAM,IAAI,KAAK,CAAC;YAChB,CAAC,IAAI,CAAC,CAAC;QACT,CAAC;QACD,OAAO,MAAM,GAAG,IAAI,CAAC,cAAc,CAAC;IACtC,CAAC;IACM,cAAc,CAAC,QAAgB,EAAE,QAAgB;QACtD,IAAI,QAAQ,IAAI,CAAC;YACf,OAAO,CAAC,CAAC;QACX,kBAAkB;QAClB,YAAY;QACZ,yCAAyC;QACzC,+DAA+D;QAC/D,MAAM,CAAC,GAAG,QAAQ,GAAG,IAAI,CAAC,cAAc,CAAC;QACzC,MAAM,EAAE,GAAG,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,cAAc,CAAC;QACvC,IAAI,MAAM,GAAG,EAAE,CAAC;QAChB,IAAI,QAAQ,GAAG,CAAC;YACd,OAAO,MAAM,GAAG,IAAI,CAAC,cAAc,CAAC;QACtC,MAAM,IAAI,GAAG,CAAE,EAAE,GAAG,EAAE,CAAC;QACvB,IAAI,KAAK,GAAG,EAAE,CAAC;QACf,IAAI,CAAC,GAAG,CAAC,CAAC;QACV,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,QAAQ,EAAE,CAAC,EAAE,EAAE,CAAC;YAClC,KAAK,IAAI,IAAI,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;YAC9B,MAAM,IAAI,KAAK,CAAC;YAChB,CAAC,IAAI,CAAC,CAAC;QACT,CAAC;QACD,OAAO,MAAM,GAAG,IAAI,CAAC,cAAc,CAAC;IACtC,CAAC;IAEM,eAAe,CAAC,QAAgB,EAAE,QAAgB;QACvD,iBAAiB;QACjB,iFAAiF;QACjF,MAAM,CAAC,GAAG,QAAQ,GAAG,IAAI,CAAC,cAAc,CAAC;QAEzC,MAAM,MAAM,GAAG,CAAC,GAAG,IAAI,CAAC,cAAc,GAAG,CAAC,CAAC;QAC3C,MAAM,IAAI,GAAG,IAAI,CAAC,cAAc,CAAC,QAAQ,EAAE,QAAQ,GAAG,CAAC,CAAC,CAAC;QACzD,MAAM,OAAO,GAAG,CAAC,CAAE,MAAM,GAAG,IAAI,GAAG,IAAI,CAAC,cAAc,CAAC,CAAC;QACxD,OAAO,OAAO,CAAC;IACjB,CAAC;IACM,eAAe,CAAC,QAAgB,EAAE,QAAgB;QACvD,eAAe;QACf,iFAAiF;QACjF,wFAAwF;QACxF,MAAM,CAAC,GAAG,QAAQ,GAAG,IAAI,CAAC,cAAc,CAAC;QACzC,MAAM,MAAM,GAAG,CAAC,GAAG,IAAI,CAAC,cAAc,GAAG,CAAC,CAAC;QAC3C,MAAM,MAAM,GAAG,IAAI,CAAC,cAAc,CAAC,QAAQ,EAAE,QAAQ,CAAC,CAAC;QACvD,OAAO,MAAM,GAAG,MAAM,GAAG,IAAI,CAAC,cAAc,CAAC;IAC/C,CAAC;IAEM,aAAa,CAAC,QAAgB;QACnC,IAAI,IAAI,CAAC,SAAS,IAAI,CAAC;YACrB,OAAO,GAAG,CAAC;QACb,iBAAiB;QACjB,iFAAiF;QACjF,MAAM,CAAC,GAAG,QAAQ,GAAG,IAAI,CAAC,cAAc,CAAC;QACzC,MAAM,MAAM,GAAG,GAAG,GAAG,IAAI,CAAC,cAAc,GAAG,CAAC,CAAC;QAC7C,MAAM,OAAO,GAAG,GAAG,GAAG,IAAI,CAAC,cAAc,CAAC;QAC1C,MAAM,IAAI,GAAG,IAAI,CAAC,cAAc,CAAC,QAAQ,EAAE,IAAI,CAAC,SAAS,GAAG,CAAC,CAAC,CAAC;QAC/D,MAAM,MAAM,GAAG,IAAI,CAAC,cAAc,CAAC,QAAQ,EAAE,IAAI,CAAC,SAAS,GAAG,CAAC,CAAC,CAAC;QACjE,OAAO,CAAC,CAAE,MAAM,GAAG,MAAM,GAAG,MAAM,GAAG,IAAI,GAAG,OAAO,CAAC,GAAG,IAAI,CAAC,cAAc,GAAG,IAAI,CAAC,cAAc,CAAC;IACnG,CAAC;IACM,aAAa,CAAC,QAAgB;QACnC,eAAe;QACf,iFAAiF;QACjF,MAAM,CAAC,GAAG,QAAQ,GAAG,IAAI,CAAC,cAAc,CAAC;QACzC,MAAM,MAAM,GAAG,GAAG,GAAG,IAAI,CAAC,cAAc,GAAG,CAAC,CAAC;QAC7C,MAAM,OAAO,GAAG,GAAG,GAAG,IAAI,CAAC,cAAc,CAAC;QAC1C,oCAAoC;QACpC,6FAA6F;QAC7F,wEAAwE;QACxE,MAAM,MAAM,GAAG,IAAI,CAAC,cAAc,CAAC,QAAQ,EAAE,IAAI,CAAC,SAAS,CAAC,CAAC;QAC7D,MAAM,IAAI,GAAG,IAAI,CAAC,cAAc,CAAC,QAAQ,EAAE,IAAI,CAAC,SAAS,GAAG,CAAC,CAAC,CAAC;QAC/D,OAAO,CAAC,CAAC,IAAI,GAAG,MAAM,GAAG,MAAM,GAAG,MAAM,GAAG,OAAO,CAAC,GAAG,IAAI,CAAC,cAAc,GAAG,IAAI,CAAC,cAAc,CAAC;IAClG,CAAC;IAEM,WAAW,CAAC,CAAS;QAC1B,MAAM,SAAS,GAAG,CAAC,GAAG,IAAI,CAAC,cAAc,CAAC;QAC1C,MAAM,SAAS,GAAG,qBAAY,CAAC,WAAW,CAAC,SAAS,EAClD,CAAC,CAAS,EAAE,EAAE,CAAC,CAAC,IAAI,CAAC,WAAW,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,EACxC,CAAC,CAAS,EAAE,EAAE,CAAC,IAAI,CAAC,YAAY,CAAC,CAAC,CAAC,CAAC,CAAC;QACvC,IAAI,SAAS,KAAK,SAAS;YACzB,OAAO,SAAS,CAAC;QACnB,OAAO,SAAS,CAAC;IACnB,CAAC;CACF;AArND,8DAqNC","sourcesContent":["/*---------------------------------------------------------------------------------------------\n* Copyright (c) Bentley Systems, Incorporated. All rights reserved.\n* See LICENSE.md in the project root for license terms and full copyright notice.\n*--------------------------------------------------------------------------------------------*/\n\n/** @packageDocumentation\n * @module Curve\n */\n\nimport { Geometry } from \"../../Geometry\";\nimport { XYCurveEvaluator } from \"./XYCurveEvaluator\";\nimport { SimpleNewton } from \"../../numerics/Newton\";\n/**\n * Methods to evaluate caller-specified number of terms of the x and y series for a clothoid.\n * Each instance has\n * * Number of x and y terms to use.\n * * constant for theta=c * x * x\n *    * This value is c=1/(2 R L)  for curve length L measured from inflection to point with radius R.\n * @internal\n */\nexport class ClothoidSeriesRLEvaluator extends XYCurveEvaluator {\n  public numXTerms: number;\n  public numYTerms: number;\n  public constantDiv2LR: number;\n  public nominalLength1: number;\n  public constructor(nominalLength1: number, constantDiv2LR: number, numXTerms: number = 4, numYTerms: number = 4) {\n    super();\n    this.nominalLength1 = nominalLength1;\n    this.constantDiv2LR = constantDiv2LR;\n    this.numXTerms = numXTerms;\n    this.numYTerms = numYTerms;\n  }\n  /** Return a deep clone. */\n  public clone(): ClothoidSeriesRLEvaluator {\n    return new ClothoidSeriesRLEvaluator(this.nominalLength1, this.constantDiv2LR, this.numXTerms, this.numYTerms);\n  }\n  public scaleInPlace(scaleFactor: number) {\n    this.nominalLength1 *= scaleFactor;\n    this.constantDiv2LR /= (scaleFactor * scaleFactor);\n  }\n  /** Member by member matchup ... */\n  public isAlmostEqual(other: any): boolean {\n    if (other instanceof ClothoidSeriesRLEvaluator) {\n      return this.numXTerms === other.numXTerms\n        && this.numYTerms === other.numYTerms\n        && Geometry.isAlmostEqualNumber(this.constantDiv2LR, other.constantDiv2LR)\n        && Geometry.isSameCoordinate(this.nominalLength1, other.nominalLength1);\n    }\n    return false;\n  }\n  /**\n   * Evaluate the X series at a nominal distance along the curve.\n   * @param fraction fractional position along the curve.\n   */\n  public fractionToX(fraction: number): number { return this.fractionToXGo(fraction, this.numXTerms); }\n  /**\n   * Evaluate the Y series at a nominal distance along the curve.\n   * @param fraction fractional position along the curve.\n   */\n  public fractionToY(fraction: number): number { return this.fractionToYGo(fraction, this.numYTerms); }\n\n  /**\n   * Evaluate the derivative of the X series at a nominal distance along the curve.\n   * @param fraction fractional position along the curve.\n   */\n  public fractionToDX(fraction: number): number { return this.fractionToDXGo(fraction, this.numXTerms); }\n  /**\n   * Evaluate the derivative of the Y series at a nominal distance along the curve.\n   * @param fraction fractional position along the curve.\n   */\n  public fractionToDY(fraction: number): number { return this.fractionToDYGo(fraction, this.numYTerms); }\n\n  /**\n   * Evaluate the derivative of the X series at a nominal distance along the curve.\n   * @param fraction fractional position along the curve.\n   */\n  public fractionToDDX(fraction: number): number { return this.fractionToDDXGo(fraction, this.numXTerms); }\n  /**\n   * Evaluate the derivative of the Y series at a nominal distance along the curve.\n   * @param fraction fractional position along the curve.\n   */\n  public fractionToDDY(fraction: number): number { return this.fractionToDDYGo(fraction, this.numYTerms); }\n\n  /**\n   * Evaluate the X series at a nominal distance along the curve.\n   * @param fraction fractional position along the curve.\n   * @param numTerms number of terms to use.\n   */\n  public fractionToXGo(fraction: number, numTerms: number): number {\n    // Write the series for cos (theta)\n    // replace theta by s*s*c\n    // integrate wrt s\n    //  x = s - s^5 c^4/ 2 + s^9 c^8/(4!) - s^13 c^12 / 6!\n    //  x = s(1 - (s^4 c^2/2) ( 1/5 -s^4 c^2 / (3*4)  ( 1/9 - ....) ) )\n    const s = fraction * this.nominalLength1;\n    let result = s;\n    if (numTerms < 2)\n      return result;\n    const q1 = s * s * this.constantDiv2LR;\n    const beta = - q1 * q1;\n    let alpha = s;\n    let m = 1;\n    let n = 5;\n    for (let i = 1; i < numTerms; i++) {\n      alpha *= beta / (m * (m + 1));\n      result += alpha / n;\n      m += 2;\n      n += 4;\n    }\n    return result;\n  }\n  public fractionToYGo(fraction: number, numTerms: number): number {\n    // Write the series for sin (theta)\n    // replace theta by s*s*c\n    // integrate wrt s\n    //  x = s^3 c^2/ 3( (1/3)) - s^7 c^6/(3!) ((1/7)) - s^11 c^10 / 5! ((1/9) - ...)\n    const s = fraction * this.nominalLength1;\n    const q1 = s * s * this.constantDiv2LR;\n    let result = q1 * s / 3;\n    if (numTerms < 2)\n      return result;\n    const beta = - q1 * q1;\n    let alpha = q1 * s;\n    let m = 2;\n    let n = 7;\n    for (let i = 1; i < numTerms; i++) {\n      alpha *= beta / (m * (m + 1));\n      result += alpha / n;\n      m += 2;\n      n += 4;\n    }\n    return result;\n  }\n  public fractionToDXGo(fraction: number, numTerms: number): number {\n    // Yes -- this does happen during derivatives of cosines with more than 0 terms !!\n    if (numTerms <= 0)\n      return 0;\n    // dX = 1 - s^4c^2/2 + s^8 c^4 / 4! -\n    // new Term = old Term * beta / (m(m+1))\n    const s = fraction * this.nominalLength1;\n    let result = 1;\n    if (numTerms < 2) {\n      return result * this.nominalLength1;\n    }\n    const q1 = s * s * this.constantDiv2LR;\n    const beta = - q1 * q1;\n    let alpha = 1.0;\n    let m = 1;\n    for (let i = 1; i < numTerms; i++) {\n      alpha *= beta / (m * (m + 1));\n      result += alpha;\n      m += 2;\n    }\n    return result * this.nominalLength1;\n  }\n  public fractionToDYGo(fraction: number, numTerms: number): number {\n    if (numTerms <= 0)\n      return 0;\n    // dY = q - q^3/3!\n    // q = s^2 c\n    // dY = s^2 c - s^6 c^3/3! + s^10 c^5/ 5!\n    // recurrence  advancing m by 2  alpha *= -(s^4 c^2) / (m(m+1))\n    const s = fraction * this.nominalLength1;\n    const q1 = s * s * this.constantDiv2LR;\n    let result = q1;\n    if (numTerms < 2)\n      return result * this.nominalLength1;\n    const beta = - q1 * q1;\n    let alpha = q1;\n    let m = 2;\n    for (let i = 1; i < numTerms; i++) {\n      alpha *= beta / (m * (m + 1));\n      result += alpha;\n      m += 2;\n    }\n    return result * this.nominalLength1;\n  }\n\n  public fractionToDDXGo(fraction: number, numTerms: number): number {\n    // DX is \"cosine\"\n    // DDX is \"- sine\" series times chain rule dTheta/ds = 2 * s * this.constantDivLR\n    const s = fraction * this.nominalLength1;\n\n    const dTheta = 2 * this.constantDiv2LR * s;\n    const sine = this.fractionToDYGo(fraction, numTerms - 1);\n    const resultA = (- dTheta * sine * this.nominalLength1);\n    return resultA;\n  }\n  public fractionToDDYGo(fraction: number, numTerms: number): number {\n    // DY is \"sine\"\n    // DDY is \"cosine\" series times chain rule dTheta/ds = 2 * s * this.constantDivLR\n    // BUT  .... derivative of the cosine series leading term is zero ... use one less term!\n    const s = fraction * this.nominalLength1;\n    const dTheta = 2 * this.constantDiv2LR * s;\n    const cosine = this.fractionToDXGo(fraction, numTerms);\n    return cosine * dTheta * this.nominalLength1;\n  }\n\n  public fractionToD3X(fraction: number): number {\n    if (this.numXTerms <= 1)\n      return 0.0;\n    // DX is \"cosine\"\n    // DDX is \"- sine\" series times chain rule dTheta/ds = 2 * s * this.constantDivLR\n    const s = fraction * this.nominalLength1;\n    const dTheta = 2.0 * this.constantDiv2LR * s;\n    const d2Theta = 2.0 * this.constantDiv2LR;\n    const sine = this.fractionToDYGo(fraction, this.numXTerms - 1);\n    const cosine = this.fractionToDXGo(fraction, this.numXTerms - 1);\n    return (- cosine * dTheta * dTheta - sine * d2Theta) * this.nominalLength1 * this.nominalLength1;\n  }\n  public fractionToD3Y(fraction: number): number {\n    // DY is \"sine\"\n    // DDY is \"cosine\" series times chain rule dTheta/ds = 2 * s * this.constantDivLR\n    const s = fraction * this.nominalLength1;\n    const dTheta = 2.0 * this.constantDiv2LR * s;\n    const d2Theta = 2.0 * this.constantDiv2LR;\n    // dY is sine series with numYTerms.\n    // ddY is cosine series.  Leading term of sine series is non-constant, so numYTerms here also\n    // d3Y is sine series. Derivative of preceding cosine killed first term.\n    const cosine = this.fractionToDXGo(fraction, this.numYTerms);\n    const sine = this.fractionToDYGo(fraction, this.numYTerms - 1);\n    return (-sine * dTheta * dTheta + cosine * d2Theta) * this.nominalLength1 * this.nominalLength1;\n  }\n\n  public xToFraction(x: number): number | undefined {\n    const fraction0 = x / this.nominalLength1;\n    const fraction1 = SimpleNewton.runNewton1D(fraction0,\n      (f: number) => (this.fractionToX(f) - x),\n      (f: number) => this.fractionToDX(f));\n    if (fraction1 === undefined)\n      return undefined;\n    return fraction1;\n  }\n}\n"]}
+{"version":3,"file":"ClothoidSeries.js","sourceRoot":"","sources":["../../../../src/curve/spiral/ClothoidSeries.ts"],"names":[],"mappings":";AAAA;;;+FAG+F;;;AAE/F;;GAEG;AAEH,6CAA0C;AAC1C,yDAAsD;AACtD,kDAAqD;AAErD;;;;;;;GAOG;AACH,MAAa,yBAA0B,SAAQ,mCAAgB;IAC7D,iDAAiD;IAC1C,SAAS,CAAS;IACzB,iDAAiD;IAC1C,SAAS,CAAS;IACzB,+CAA+C;IACxC,cAAc,CAAS;IAC9B,gCAAgC;IACzB,cAAc,CAAS;IAC9B,YAAmB,cAAsB,EAAE,cAAsB,EAAE,YAAoB,CAAC,EAAE,YAAoB,CAAC;QAC7G,KAAK,EAAE,CAAC;QACR,IAAI,CAAC,cAAc,GAAG,cAAc,CAAC;QACrC,IAAI,CAAC,cAAc,GAAG,cAAc,CAAC;QACrC,IAAI,CAAC,SAAS,GAAG,SAAS,CAAC;QAC3B,IAAI,CAAC,SAAS,GAAG,SAAS,CAAC;IAC7B,CAAC;IACD,2BAA2B;IACpB,KAAK;QACV,OAAO,IAAI,yBAAyB,CAAC,IAAI,CAAC,cAAc,EAAE,IAAI,CAAC,cAAc,EAAE,IAAI,CAAC,SAAS,EAAE,IAAI,CAAC,SAAS,CAAC,CAAC;IACjH,CAAC;IACM,YAAY,CAAC,WAAmB;QACrC,IAAI,CAAC,cAAc,IAAI,WAAW,CAAC;QACnC,IAAI,CAAC,cAAc,IAAI,CAAC,WAAW,GAAG,WAAW,CAAC,CAAC;IACrD,CAAC;IACD,gCAAgC;IACzB,aAAa,CAAC,KAAU;QAC7B,IAAI,KAAK,YAAY,yBAAyB,EAAE,CAAC;YAC/C,OAAO,IAAI,CAAC,SAAS,KAAK,KAAK,CAAC,SAAS;mBACpC,IAAI,CAAC,SAAS,KAAK,KAAK,CAAC,SAAS;mBAClC,mBAAQ,CAAC,mBAAmB,CAAC,IAAI,CAAC,cAAc,EAAE,KAAK,CAAC,cAAc,CAAC;mBACvE,mBAAQ,CAAC,gBAAgB,CAAC,IAAI,CAAC,cAAc,EAAE,KAAK,CAAC,cAAc,CAAC,CAAC;QAC5E,CAAC;QACD,OAAO,KAAK,CAAC;IACf,CAAC;IACD;;;OAGG;IACI,WAAW,CAAC,QAAgB;QACjC,OAAO,IAAI,CAAC,aAAa,CAAC,QAAQ,EAAE,IAAI,CAAC,SAAS,CAAC,CAAC;IACtD,CAAC;IACD;;;OAGG;IACI,WAAW,CAAC,QAAgB;QACjC,OAAO,IAAI,CAAC,aAAa,CAAC,QAAQ,EAAE,IAAI,CAAC,SAAS,CAAC,CAAC;IACtD,CAAC;IACD;;;OAGG;IACI,YAAY,CAAC,QAAgB;QAClC,OAAO,IAAI,CAAC,cAAc,CAAC,QAAQ,EAAE,IAAI,CAAC,SAAS,CAAC,CAAC;IACvD,CAAC;IACD;;;OAGG;IACI,YAAY,CAAC,QAAgB;QAClC,OAAO,IAAI,CAAC,cAAc,CAAC,QAAQ,EAAE,IAAI,CAAC,SAAS,CAAC,CAAC;IACvD,CAAC;IACD;;;OAGG;IACI,aAAa,CAAC,QAAgB;QACnC,OAAO,IAAI,CAAC,eAAe,CAAC,QAAQ,EAAE,IAAI,CAAC,SAAS,CAAC,CAAC;IACxD,CAAC;IACD;;;OAGG;IACI,aAAa,CAAC,QAAgB;QACnC,OAAO,IAAI,CAAC,eAAe,CAAC,QAAQ,EAAE,IAAI,CAAC,SAAS,CAAC,CAAC;IACxD,CAAC;IACD;;;;OAIG;IACI,aAAa,CAAC,QAAgB,EAAE,QAAgB;QACrD,yCAAyC;QACzC,yBAAyB;QACzB,kBAAkB;QAClB,mEAAmE;QACnE,iEAAiE;QACjE,MAAM,CAAC,GAAG,QAAQ,GAAG,IAAI,CAAC,cAAc,CAAC;QACzC,IAAI,MAAM,GAAG,CAAC,CAAC;QACf,IAAI,QAAQ,GAAG,CAAC;YACd,OAAO,MAAM,CAAC;QAChB,MAAM,EAAE,GAAG,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,cAAc,CAAC;QACvC,MAAM,IAAI,GAAG,CAAE,EAAE,GAAG,EAAE,CAAC;QACvB,IAAI,KAAK,GAAG,CAAC,CAAC;QACd,IAAI,CAAC,GAAG,CAAC,CAAC;QACV,IAAI,CAAC,GAAG,CAAC,CAAC;QACV,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,QAAQ,EAAE,CAAC,EAAE,EAAE,CAAC;YAClC,KAAK,IAAI,IAAI,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;YAC9B,MAAM,IAAI,KAAK,GAAG,CAAC,CAAC;YACpB,CAAC,IAAI,CAAC,CAAC;YACP,CAAC,IAAI,CAAC,CAAC;QACT,CAAC;QACD,OAAO,MAAM,CAAC;IAChB,CAAC;IACM,aAAa,CAAC,QAAgB,EAAE,QAAgB;QACrD,yCAAyC;QACzC,yBAAyB;QACzB,kBAAkB;QAClB,6EAA6E;QAC7E,2EAA2E;QAC3E,MAAM,CAAC,GAAG,QAAQ,GAAG,IAAI,CAAC,cAAc,CAAC;QACzC,MAAM,EAAE,GAAG,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,cAAc,CAAC;QACvC,IAAI,MAAM,GAAG,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC;QACxB,IAAI,QAAQ,GAAG,CAAC;YACd,OAAO,MAAM,CAAC;QAChB,MAAM,IAAI,GAAG,CAAE,EAAE,GAAG,EAAE,CAAC;QACvB,IAAI,KAAK,GAAG,EAAE,GAAG,CAAC,CAAC;QACnB,IAAI,CAAC,GAAG,CAAC,CAAC;QACV,IAAI,CAAC,GAAG,CAAC,CAAC;QACV,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,QAAQ,EAAE,CAAC,EAAE,EAAE,CAAC;YAClC,KAAK,IAAI,IAAI,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;YAC9B,MAAM,IAAI,KAAK,GAAG,CAAC,CAAC;YACpB,CAAC,IAAI,CAAC,CAAC;YACP,CAAC,IAAI,CAAC,CAAC;QACT,CAAC;QACD,OAAO,MAAM,CAAC;IAChB,CAAC;IACM,cAAc,CAAC,QAAgB,EAAE,QAAgB;QACtD,+EAA+E;QAC/E,IAAI,QAAQ,IAAI,CAAC;YACf,OAAO,CAAC,CAAC;QACX,qCAAqC;QACrC,wCAAwC;QACxC,MAAM,CAAC,GAAG,QAAQ,GAAG,IAAI,CAAC,cAAc,CAAC;QACzC,IAAI,MAAM,GAAG,CAAC,CAAC;QACf,IAAI,QAAQ,GAAG,CAAC,EAAE,CAAC;YACjB,OAAO,MAAM,GAAG,IAAI,CAAC,cAAc,CAAC;QACtC,CAAC;QACD,MAAM,EAAE,GAAG,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,cAAc,CAAC;QACvC,MAAM,IAAI,GAAG,CAAE,EAAE,GAAG,EAAE,CAAC;QACvB,IAAI,KAAK,GAAG,GAAG,CAAC;QAChB,IAAI,CAAC,GAAG,CAAC,CAAC;QACV,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,QAAQ,EAAE,CAAC,EAAE,EAAE,CAAC;YAClC,KAAK,IAAI,IAAI,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;YAC9B,MAAM,IAAI,KAAK,CAAC;YAChB,CAAC,IAAI,CAAC,CAAC;QACT,CAAC;QACD,OAAO,MAAM,GAAG,IAAI,CAAC,cAAc,CAAC;IACtC,CAAC;IACM,cAAc,CAAC,QAAgB,EAAE,QAAgB;QACtD,IAAI,QAAQ,IAAI,CAAC;YACf,OAAO,CAAC,CAAC;QACX,kBAAkB;QAClB,YAAY;QACZ,yCAAyC;QACzC,+DAA+D;QAC/D,MAAM,CAAC,GAAG,QAAQ,GAAG,IAAI,CAAC,cAAc,CAAC;QACzC,MAAM,EAAE,GAAG,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,cAAc,CAAC;QACvC,IAAI,MAAM,GAAG,EAAE,CAAC;QAChB,IAAI,QAAQ,GAAG,CAAC;YACd,OAAO,MAAM,GAAG,IAAI,CAAC,cAAc,CAAC;QACtC,MAAM,IAAI,GAAG,CAAE,EAAE,GAAG,EAAE,CAAC;QACvB,IAAI,KAAK,GAAG,EAAE,CAAC;QACf,IAAI,CAAC,GAAG,CAAC,CAAC;QACV,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,QAAQ,EAAE,CAAC,EAAE,EAAE,CAAC;YAClC,KAAK,IAAI,IAAI,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;YAC9B,MAAM,IAAI,KAAK,CAAC;YAChB,CAAC,IAAI,CAAC,CAAC;QACT,CAAC;QACD,OAAO,MAAM,GAAG,IAAI,CAAC,cAAc,CAAC;IACtC,CAAC;IACM,eAAe,CAAC,QAAgB,EAAE,QAAgB;QACvD,iBAAiB;QACjB,iFAAiF;QACjF,MAAM,CAAC,GAAG,QAAQ,GAAG,IAAI,CAAC,cAAc,CAAC;QACzC,MAAM,MAAM,GAAG,CAAC,GAAG,IAAI,CAAC,cAAc,GAAG,CAAC,CAAC;QAC3C,MAAM,IAAI,GAAG,IAAI,CAAC,cAAc,CAAC,QAAQ,EAAE,QAAQ,GAAG,CAAC,CAAC,CAAC;QACzD,MAAM,OAAO,GAAG,CAAC,CAAE,MAAM,GAAG,IAAI,GAAG,IAAI,CAAC,cAAc,CAAC,CAAC;QACxD,OAAO,OAAO,CAAC;IACjB,CAAC;IACM,eAAe,CAAC,QAAgB,EAAE,QAAgB;QACvD,eAAe;QACf,iFAAiF;QACjF,uFAAuF;QACvF,MAAM,CAAC,GAAG,QAAQ,GAAG,IAAI,CAAC,cAAc,CAAC;QACzC,MAAM,MAAM,GAAG,CAAC,GAAG,IAAI,CAAC,cAAc,GAAG,CAAC,CAAC;QAC3C,MAAM,MAAM,GAAG,IAAI,CAAC,cAAc,CAAC,QAAQ,EAAE,QAAQ,CAAC,CAAC;QACvD,OAAO,MAAM,GAAG,MAAM,GAAG,IAAI,CAAC,cAAc,CAAC;IAC/C,CAAC;IACM,aAAa,CAAC,QAAgB;QACnC,IAAI,IAAI,CAAC,SAAS,IAAI,CAAC;YACrB,OAAO,GAAG,CAAC;QACb,iBAAiB;QACjB,iFAAiF;QACjF,MAAM,CAAC,GAAG,QAAQ,GAAG,IAAI,CAAC,cAAc,CAAC;QACzC,MAAM,MAAM,GAAG,GAAG,GAAG,IAAI,CAAC,cAAc,GAAG,CAAC,CAAC;QAC7C,MAAM,OAAO,GAAG,GAAG,GAAG,IAAI,CAAC,cAAc,CAAC;QAC1C,MAAM,IAAI,GAAG,IAAI,CAAC,cAAc,CAAC,QAAQ,EAAE,IAAI,CAAC,SAAS,GAAG,CAAC,CAAC,CAAC;QAC/D,MAAM,MAAM,GAAG,IAAI,CAAC,cAAc,CAAC,QAAQ,EAAE,IAAI,CAAC,SAAS,GAAG,CAAC,CAAC,CAAC;QACjE,OAAO,CAAC,CAAE,MAAM,GAAG,MAAM,GAAG,MAAM,GAAG,IAAI,GAAG,OAAO,CAAC,GAAG,IAAI,CAAC,cAAc,GAAG,IAAI,CAAC,cAAc,CAAC;IACnG,CAAC;IACM,aAAa,CAAC,QAAgB;QACnC,eAAe;QACf,iFAAiF;QACjF,MAAM,CAAC,GAAG,QAAQ,GAAG,IAAI,CAAC,cAAc,CAAC;QACzC,MAAM,MAAM,GAAG,GAAG,GAAG,IAAI,CAAC,cAAc,GAAG,CAAC,CAAC;QAC7C,MAAM,OAAO,GAAG,GAAG,GAAG,IAAI,CAAC,cAAc,CAAC;QAC1C,oCAAoC;QACpC,6FAA6F;QAC7F,wEAAwE;QACxE,MAAM,MAAM,GAAG,IAAI,CAAC,cAAc,CAAC,QAAQ,EAAE,IAAI,CAAC,SAAS,CAAC,CAAC;QAC7D,MAAM,IAAI,GAAG,IAAI,CAAC,cAAc,CAAC,QAAQ,EAAE,IAAI,CAAC,SAAS,GAAG,CAAC,CAAC,CAAC;QAC/D,OAAO,CAAC,CAAC,IAAI,GAAG,MAAM,GAAG,MAAM,GAAG,MAAM,GAAG,OAAO,CAAC,GAAG,IAAI,CAAC,cAAc,GAAG,IAAI,CAAC,cAAc,CAAC;IAClG,CAAC;IACM,WAAW,CAAC,CAAS;QAC1B,MAAM,SAAS,GAAG,CAAC,GAAG,IAAI,CAAC,cAAc,CAAC;QAC1C,MAAM,SAAS,GAAG,qBAAY,CAAC,WAAW,CAAC,SAAS,EAClD,CAAC,CAAS,EAAE,EAAE,CAAC,CAAC,IAAI,CAAC,WAAW,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,EACxC,CAAC,CAAS,EAAE,EAAE,CAAC,IAAI,CAAC,YAAY,CAAC,CAAC,CAAC,CAAC,CAAC;QACvC,IAAI,SAAS,KAAK,SAAS;YACzB,OAAO,SAAS,CAAC;QACnB,OAAO,SAAS,CAAC;IACnB,CAAC;CACF;AA/ND,8DA+NC","sourcesContent":["/*---------------------------------------------------------------------------------------------\n* Copyright (c) Bentley Systems, Incorporated. All rights reserved.\n* See LICENSE.md in the project root for license terms and full copyright notice.\n*--------------------------------------------------------------------------------------------*/\n\n/** @packageDocumentation\n * @module Curve\n */\n\nimport { Geometry } from \"../../Geometry\";\nimport { XYCurveEvaluator } from \"./XYCurveEvaluator\";\nimport { SimpleNewton } from \"../../numerics/Newton\";\n\n/**\n * Methods to evaluate caller-specified number of terms of the x and y Taylor series for a clothoid.\n * Each instance has:\n * * number of x and y terms to use.\n * * constant for theta = c*x*x.\n *    * This value is c = 1/(2*R*L) for curve length L measured from inflection to point with radius R.\n * @internal\n */\nexport class ClothoidSeriesRLEvaluator extends XYCurveEvaluator {\n  /** Number of terms to use in x Taylor series. */\n  public numXTerms: number;\n  /** Number of terms to use in y Taylor series. */\n  public numYTerms: number;\n  /** Constant c = 1/(2*R*L) in theta = c*s*s. */\n  public constantDiv2LR: number;\n  /** The nominal curve length. */\n  public nominalLength1: number;\n  public constructor(nominalLength1: number, constantDiv2LR: number, numXTerms: number = 4, numYTerms: number = 4) {\n    super();\n    this.nominalLength1 = nominalLength1;\n    this.constantDiv2LR = constantDiv2LR;\n    this.numXTerms = numXTerms;\n    this.numYTerms = numYTerms;\n  }\n  /** Return a deep clone. */\n  public clone(): ClothoidSeriesRLEvaluator {\n    return new ClothoidSeriesRLEvaluator(this.nominalLength1, this.constantDiv2LR, this.numXTerms, this.numYTerms);\n  }\n  public scaleInPlace(scaleFactor: number): void {\n    this.nominalLength1 *= scaleFactor;\n    this.constantDiv2LR /= (scaleFactor * scaleFactor);\n  }\n  /** Member by member matchup. */\n  public isAlmostEqual(other: any): boolean {\n    if (other instanceof ClothoidSeriesRLEvaluator) {\n      return this.numXTerms === other.numXTerms\n        && this.numYTerms === other.numYTerms\n        && Geometry.isAlmostEqualNumber(this.constantDiv2LR, other.constantDiv2LR)\n        && Geometry.isSameCoordinate(this.nominalLength1, other.nominalLength1);\n    }\n    return false;\n  }\n  /**\n   * Evaluate the X series at a nominal distance along the curve.\n   * @param fraction fractional position along the curve.\n   */\n  public fractionToX(fraction: number): number {\n    return this.fractionToXGo(fraction, this.numXTerms);\n  }\n  /**\n   * Evaluate the Y series at a nominal distance along the curve.\n   * @param fraction fractional position along the curve.\n   */\n  public fractionToY(fraction: number): number {\n    return this.fractionToYGo(fraction, this.numYTerms);\n  }\n  /**\n   * Evaluate the derivative of the X series at a nominal distance along the curve.\n   * @param fraction fractional position along the curve.\n   */\n  public fractionToDX(fraction: number): number {\n    return this.fractionToDXGo(fraction, this.numXTerms);\n  }\n  /**\n   * Evaluate the derivative of the Y series at a nominal distance along the curve.\n   * @param fraction fractional position along the curve.\n   */\n  public fractionToDY(fraction: number): number {\n    return this.fractionToDYGo(fraction, this.numYTerms);\n  }\n  /**\n   * Evaluate the derivative of the X series at a nominal distance along the curve.\n   * @param fraction fractional position along the curve.\n   */\n  public fractionToDDX(fraction: number): number {\n    return this.fractionToDDXGo(fraction, this.numXTerms);\n  }\n  /**\n   * Evaluate the derivative of the Y series at a nominal distance along the curve.\n   * @param fraction fractional position along the curve.\n   */\n  public fractionToDDY(fraction: number): number {\n    return this.fractionToDDYGo(fraction, this.numYTerms);\n  }\n  /**\n   * Evaluate the X series at a nominal distance along the curve.\n   * @param fraction fractional position along the curve.\n   * @param numTerms number of terms to use.\n   */\n  public fractionToXGo(fraction: number, numTerms: number): number {\n    // write the Taylor series for cos(theta)\n    // replace theta by s*s*c\n    // integrate wrt s\n    //  x = s - s^5 c^2/ 5*2! + s^9 c^4 / 9*4! - s^13 c^6 / 13*6! + ...\n    //  x = s(1 - (s^4 c^2/2) ( 1/5 - (s^4 c^2 / 3*4) (1/9 - ...) ) )\n    const s = fraction * this.nominalLength1;\n    let result = s;\n    if (numTerms < 2)\n      return result;\n    const q1 = s * s * this.constantDiv2LR;\n    const beta = - q1 * q1;\n    let alpha = s;\n    let m = 1;\n    let n = 5;\n    for (let i = 1; i < numTerms; i++) {\n      alpha *= beta / (m * (m + 1));\n      result += alpha / n;\n      m += 2;\n      n += 4;\n    }\n    return result;\n  }\n  public fractionToYGo(fraction: number, numTerms: number): number {\n    // write the Taylor series for sin(theta)\n    // replace theta by s*s*c\n    // integrate wrt s\n    // y = s^3 c / 3 - s^7 c^3 / 7*3! + s^11 c^5 / 11*5! - s^15 c^7 / 15*7! + ...\n    // y = s^3 c ( 1/3 - s^4 c^2/ 3! ( (1/7) - (s^4 c^2 / 4*5) (1/11 - ...) ) )\n    const s = fraction * this.nominalLength1;\n    const q1 = s * s * this.constantDiv2LR;\n    let result = q1 * s / 3;\n    if (numTerms < 2)\n      return result;\n    const beta = - q1 * q1;\n    let alpha = q1 * s;\n    let m = 2;\n    let n = 7;\n    for (let i = 1; i < numTerms; i++) {\n      alpha *= beta / (m * (m + 1));\n      result += alpha / n;\n      m += 2;\n      n += 4;\n    }\n    return result;\n  }\n  public fractionToDXGo(fraction: number, numTerms: number): number {\n    // yes -- this does happen during derivatives of cosines with more than 0 terms\n    if (numTerms <= 0)\n      return 0;\n    // dX = 1 - s^4c^2/2 + s^8 c^4 / 4! -\n    // new Term = old Term * beta / (m(m+1))\n    const s = fraction * this.nominalLength1;\n    let result = 1;\n    if (numTerms < 2) {\n      return result * this.nominalLength1;\n    }\n    const q1 = s * s * this.constantDiv2LR;\n    const beta = - q1 * q1;\n    let alpha = 1.0;\n    let m = 1;\n    for (let i = 1; i < numTerms; i++) {\n      alpha *= beta / (m * (m + 1));\n      result += alpha;\n      m += 2;\n    }\n    return result * this.nominalLength1;\n  }\n  public fractionToDYGo(fraction: number, numTerms: number): number {\n    if (numTerms <= 0)\n      return 0;\n    // dY = q - q^3/3!\n    // q = s^2 c\n    // dY = s^2 c - s^6 c^3/3! + s^10 c^5/ 5!\n    // recurrence  advancing m by 2  alpha *= -(s^4 c^2) / (m(m+1))\n    const s = fraction * this.nominalLength1;\n    const q1 = s * s * this.constantDiv2LR;\n    let result = q1;\n    if (numTerms < 2)\n      return result * this.nominalLength1;\n    const beta = - q1 * q1;\n    let alpha = q1;\n    let m = 2;\n    for (let i = 1; i < numTerms; i++) {\n      alpha *= beta / (m * (m + 1));\n      result += alpha;\n      m += 2;\n    }\n    return result * this.nominalLength1;\n  }\n  public fractionToDDXGo(fraction: number, numTerms: number): number {\n    // DX is \"cosine\"\n    // DDX is \"- sine\" series times chain rule dTheta/ds = 2 * s * this.constantDivLR\n    const s = fraction * this.nominalLength1;\n    const dTheta = 2 * this.constantDiv2LR * s;\n    const sine = this.fractionToDYGo(fraction, numTerms - 1);\n    const resultA = (- dTheta * sine * this.nominalLength1);\n    return resultA;\n  }\n  public fractionToDDYGo(fraction: number, numTerms: number): number {\n    // DY is \"sine\"\n    // DDY is \"cosine\" series times chain rule dTheta/ds = 2 * s * this.constantDivLR\n    // BUT  .... derivative of the cosine series leading term is zero ... use one less term\n    const s = fraction * this.nominalLength1;\n    const dTheta = 2 * this.constantDiv2LR * s;\n    const cosine = this.fractionToDXGo(fraction, numTerms);\n    return cosine * dTheta * this.nominalLength1;\n  }\n  public fractionToD3X(fraction: number): number {\n    if (this.numXTerms <= 1)\n      return 0.0;\n    // DX is \"cosine\"\n    // DDX is \"- sine\" series times chain rule dTheta/ds = 2 * s * this.constantDivLR\n    const s = fraction * this.nominalLength1;\n    const dTheta = 2.0 * this.constantDiv2LR * s;\n    const d2Theta = 2.0 * this.constantDiv2LR;\n    const sine = this.fractionToDYGo(fraction, this.numXTerms - 1);\n    const cosine = this.fractionToDXGo(fraction, this.numXTerms - 1);\n    return (- cosine * dTheta * dTheta - sine * d2Theta) * this.nominalLength1 * this.nominalLength1;\n  }\n  public fractionToD3Y(fraction: number): number {\n    // DY is \"sine\"\n    // DDY is \"cosine\" series times chain rule dTheta/ds = 2 * s * this.constantDivLR\n    const s = fraction * this.nominalLength1;\n    const dTheta = 2.0 * this.constantDiv2LR * s;\n    const d2Theta = 2.0 * this.constantDiv2LR;\n    // dY is sine series with numYTerms.\n    // ddY is cosine series.  Leading term of sine series is non-constant, so numYTerms here also\n    // d3Y is sine series. Derivative of preceding cosine killed first term.\n    const cosine = this.fractionToDXGo(fraction, this.numYTerms);\n    const sine = this.fractionToDYGo(fraction, this.numYTerms - 1);\n    return (-sine * dTheta * dTheta + cosine * d2Theta) * this.nominalLength1 * this.nominalLength1;\n  }\n  public xToFraction(x: number): number | undefined {\n    const fraction0 = x / this.nominalLength1;\n    const fraction1 = SimpleNewton.runNewton1D(fraction0,\n      (f: number) => (this.fractionToX(f) - x),\n      (f: number) => this.fractionToDX(f));\n    if (fraction1 === undefined)\n      return undefined;\n    return fraction1;\n  }\n}\n"]}

package/lib/cjs/curve/spiral/CubicEvaluator.d.ts

@@ -2,9 +2,11 @@
  * @module Curve
  */
 import { XYCurveEvaluator } from "./XYCurveEvaluator";
-/** Intermediate class for evaluation of bare cubic spirals `y = m ^x^3` with x axis from [0..xLength]
+/**
+ * Intermediate class for evaluation of bare cubic spirals `y = m*x^3` with x axis from [0..xLength]
  * * This implements all the computations among fraction, x, and y.
- * * Derived classes implement specialized logic such as (a) precomputing `m` and (b) domain-specific fraction-to-distance approximations.
+ * * Derived classes implement specialized logic such as (a) precomputing `m` and (b) domain-specific fraction-to-distance
+ * approximations.
  * @internal
  */
 export declare abstract class CubicEvaluator extends XYCurveEvaluator {
@@ -18,16 +20,15 @@ export declare abstract class CubicEvaluator extends XYCurveEvaluator {
     /**
      * Apply `scaleFactor` to the xLength and cubicM.
      * * Derived classes commonly call this as `super.scaleInPlace()`, and additionally apply the scale to their members.
-     * @param scaleFactor
      */
     scaleInPlace(scaleFactor: number): void;
     /** Evaluate X at fraction. */
     fractionToX(fraction: number): number;
-    /** Evaluate derivative of X with respect to fraction */
+    /** Evaluate derivative of X with respect to fraction. */
     fractionToDX(_fraction: number): number;
-    /** Evaluate second derivative of X with respect to fraction */
+    /** Evaluate second derivative of X with respect to fraction. */
     fractionToDDX(_fraction: number): number;
-    /** Evaluate third derivative of X with respect to fraction */
+    /** Evaluate third derivative of X with respect to fraction. */
     fractionToD3X(_fraction: number): number;
     /** Evaluate Y at fraction. */
     fractionToY(fraction: number): number;

package/lib/cjs/curve/spiral/CubicEvaluator.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"CubicEvaluator.d.ts","sourceRoot":"","sources":["../../../../src/curve/spiral/CubicEvaluator.ts"],"names":[],"mappings":"AAKA;;GAEG;AACH,OAAO,EAAE,gBAAgB,EAAE,MAAM,oBAAoB,CAAC;AACtD;;;;GAIG;AACH,8BAAsB,cAAe,SAAQ,gBAAgB;IAC3D,SAAS,CAAC,OAAO,EAAE,MAAM,CAAC;IAC1B,SAAS,CAAC,WAAW,EAAE,MAAM,CAAC;IAC9B,SAAS,aAAa,UAAU,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM;IAKxD,6BAA6B;IACtB,YAAY,CAAC,UAAU,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM;IAItD,IAAW,UAAU,WAA+B;IACpD,IAAW,MAAM,WAA2B;IAC5C;;;;OAIG;IACI,YAAY,CAAC,WAAW,EAAE,MAAM;IAKvC,8BAA8B;IACvB,WAAW,CAAC,QAAQ,EAAE,MAAM,GAAG,MAAM;IAC5C,wDAAwD;IACjD,YAAY,CAAC,SAAS,EAAE,MAAM,GAAG,MAAM;IAC9C,+DAA+D;IACxD,aAAa,CAAC,SAAS,EAAE,MAAM,GAAG,MAAM;IAC/C,8DAA8D;IACvD,aAAa,CAAC,SAAS,EAAE,MAAM,GAAG,MAAM;IAC/C,8BAA8B;IACvB,WAAW,CAAC,QAAQ,EAAE,MAAM,GAAG,MAAM;IAI5C,yDAAyD;IAClD,YAAY,CAAC,QAAQ,EAAE,MAAM,GAAG,MAAM;IAI7C,gEAAgE;IACzD,aAAa,CAAC,QAAQ,EAAE,MAAM,GAAG,MAAM;IAI9C,+DAA+D;IACxD,aAAa,CAAC,SAAS,EAAE,MAAM,GAAG,MAAM;IAG/C,8BAA8B;IACvB,WAAW,CAAC,CAAC,EAAE,MAAM,GAAG,MAAM;CACtC"}
+{"version":3,"file":"CubicEvaluator.d.ts","sourceRoot":"","sources":["../../../../src/curve/spiral/CubicEvaluator.ts"],"names":[],"mappings":"AAKA;;GAEG;AACH,OAAO,EAAE,gBAAgB,EAAE,MAAM,oBAAoB,CAAC;AAEtD;;;;;;GAMG;AACH,8BAAsB,cAAe,SAAQ,gBAAgB;IAC3D,SAAS,CAAC,OAAO,EAAE,MAAM,CAAC;IAC1B,SAAS,CAAC,WAAW,EAAE,MAAM,CAAC;IAC9B,SAAS,aAAa,UAAU,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM;IAKxD,6BAA6B;IACtB,YAAY,CAAC,UAAU,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM;IAItD,IAAW,UAAU,WAEpB;IACD,IAAW,MAAM,WAEhB;IACD;;;OAGG;IACI,YAAY,CAAC,WAAW,EAAE,MAAM;IAKvC,8BAA8B;IACvB,WAAW,CAAC,QAAQ,EAAE,MAAM,GAAG,MAAM;IAG5C,yDAAyD;IAClD,YAAY,CAAC,SAAS,EAAE,MAAM,GAAG,MAAM;IAG9C,gEAAgE;IACzD,aAAa,CAAC,SAAS,EAAE,MAAM,GAAG,MAAM;IAG/C,+DAA+D;IACxD,aAAa,CAAC,SAAS,EAAE,MAAM,GAAG,MAAM;IAG/C,8BAA8B;IACvB,WAAW,CAAC,QAAQ,EAAE,MAAM,GAAG,MAAM;IAI5C,yDAAyD;IAClD,YAAY,CAAC,QAAQ,EAAE,MAAM,GAAG,MAAM;IAI7C,gEAAgE;IACzD,aAAa,CAAC,QAAQ,EAAE,MAAM,GAAG,MAAM;IAI9C,+DAA+D;IACxD,aAAa,CAAC,SAAS,EAAE,MAAM,GAAG,MAAM;IAG/C,8BAA8B;IACvB,WAAW,CAAC,CAAC,EAAE,MAAM,GAAG,MAAM;CAGtC"}

package/lib/cjs/curve/spiral/CubicEvaluator.js

@@ -9,9 +9,11 @@ exports.CubicEvaluator = void 0;
  * @module Curve
  */
 const XYCurveEvaluator_1 = require("./XYCurveEvaluator");
-/** Intermediate class for evaluation of bare cubic spirals `y = m ^x^3` with x axis from [0..xLength]
+/**
+ * Intermediate class for evaluation of bare cubic spirals `y = m*x^3` with x axis from [0..xLength]
  * * This implements all the computations among fraction, x, and y.
- * * Derived classes implement specialized logic such as (a) precomputing `m` and (b) domain-specific fraction-to-distance approximations.
+ * * Derived classes implement specialized logic such as (a) precomputing `m` and (b) domain-specific fraction-to-distance
+ * approximations.
  * @internal
  */
 class CubicEvaluator extends XYCurveEvaluator_1.XYCurveEvaluator {
@@ -27,26 +29,37 @@ class CubicEvaluator extends XYCurveEvaluator_1.XYCurveEvaluator {
         this._axisLength = axisLength;
         this._cubicM = cubicM;
     }
-    get axisLength() { return this._axisLength; }
-    get cubicM() { return this._cubicM; }
+    get axisLength() {
+        return this._axisLength;
+    }
+    get cubicM() {
+        return this._cubicM;
+    }
     /**
      * Apply `scaleFactor` to the xLength and cubicM.
      * * Derived classes commonly call this as `super.scaleInPlace()`, and additionally apply the scale to their members.
-     * @param scaleFactor
      */
     scaleInPlace(scaleFactor) {
         this._axisLength *= scaleFactor;
-        // "x" arriving at "m * x^3" will be scaled. "m" has to be divided by the scale to cancel 2 of the 3 . .
+        // "x" arriving at "m * x^3" will be scaled. "m" has to be divided by the scale to cancel 2 of the 3
         this._cubicM /= (scaleFactor * scaleFactor);
     }
     /** Evaluate X at fraction. */
-    fractionToX(fraction) { return fraction * this._axisLength; }
-    /** Evaluate derivative of X with respect to fraction */
-    fractionToDX(_fraction) { return this._axisLength; }
-    /** Evaluate second derivative of X with respect to fraction */
-    fractionToDDX(_fraction) { return 0.0; }
-    /** Evaluate third derivative of X with respect to fraction */
-    fractionToD3X(_fraction) { return 0.0; }
+    fractionToX(fraction) {
+        return fraction * this._axisLength;
+    }
+    /** Evaluate derivative of X with respect to fraction. */
+    fractionToDX(_fraction) {
+        return this._axisLength;
+    }
+    /** Evaluate second derivative of X with respect to fraction. */
+    fractionToDDX(_fraction) {
+        return 0.0;
+    }
+    /** Evaluate third derivative of X with respect to fraction. */
+    fractionToD3X(_fraction) {
+        return 0.0;
+    }
     /** Evaluate Y at fraction. */
     fractionToY(fraction) {
         const x = fraction * this._axisLength;
@@ -67,7 +80,9 @@ class CubicEvaluator extends XYCurveEvaluator_1.XYCurveEvaluator {
         return 6.0 * this._cubicM * this._axisLength * this._axisLength * this._axisLength;
     }
     /** Evaluate fraction at x. */
-    xToFraction(x) { return x / this._axisLength; }
+    xToFraction(x) {
+        return x / this._axisLength;
+    }
 }
 exports.CubicEvaluator = CubicEvaluator;
 //# sourceMappingURL=CubicEvaluator.js.map

package/lib/cjs/curve/spiral/CubicEvaluator.js.map

@@ -1 +1 @@
-{"version":3,"file":"CubicEvaluator.js","sourceRoot":"","sources":["../../../../src/curve/spiral/CubicEvaluator.ts"],"names":[],"mappings":";AAAA;;;+FAG+F;;;AAE/F;;GAEG;AACH,yDAAsD;AACtD;;;;GAIG;AACH,MAAsB,cAAe,SAAQ,mCAAgB;IACjD,OAAO,CAAS;IAChB,WAAW,CAAS;IAC9B,YAAsB,UAAkB,EAAE,MAAc;QACtD,KAAK,EAAE,CAAC;QACR,IAAI,CAAC,OAAO,GAAG,MAAM,CAAC;QACtB,IAAI,CAAC,WAAW,GAAG,UAAU,CAAC;IAChC,CAAC;IACD,6BAA6B;IACtB,YAAY,CAAC,UAAkB,EAAE,MAAc;QACpD,IAAI,CAAC,WAAW,GAAG,UAAU,CAAC;QAC9B,IAAI,CAAC,OAAO,GAAG,MAAM,CAAC;IACxB,CAAC;IACD,IAAW,UAAU,KAAK,OAAO,IAAI,CAAC,WAAW,CAAC,CAAC,CAAC;IACpD,IAAW,MAAM,KAAK,OAAO,IAAI,CAAC,OAAO,CAAC,CAAC,CAAC;IAC5C;;;;OAIG;IACI,YAAY,CAAC,WAAmB;QACrC,IAAI,CAAC,WAAW,IAAI,WAAW,CAAC;QAChC,wGAAwG;QACxG,IAAI,CAAC,OAAO,IAAI,CAAC,WAAW,GAAG,WAAW,CAAC,CAAC;IAC9C,CAAC;IACD,8BAA8B;IACvB,WAAW,CAAC,QAAgB,IAAY,OAAO,QAAQ,GAAG,IAAI,CAAC,WAAW,CAAC,CAAC,CAAC;IACpF,wDAAwD;IACjD,YAAY,CAAC,SAAiB,IAAY,OAAO,IAAI,CAAC,WAAW,CAAC,CAAC,CAAC;IAC3E,+DAA+D;IACxD,aAAa,CAAC,SAAiB,IAAY,OAAO,GAAG,CAAC,CAAC,CAAC;IAC/D,8DAA8D;IACvD,aAAa,CAAC,SAAiB,IAAY,OAAO,GAAG,CAAC,CAAC,CAAC;IAC/D,8BAA8B;IACvB,WAAW,CAAC,QAAgB;QACjC,MAAM,CAAC,GAAG,QAAQ,GAAG,IAAI,CAAC,WAAW,CAAC;QACtC,OAAO,IAAI,CAAC,OAAO,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;IAClC,CAAC;IACD,yDAAyD;IAClD,YAAY,CAAC,QAAgB;QAClC,MAAM,CAAC,GAAG,QAAQ,GAAG,IAAI,CAAC,WAAW,CAAC;QACtC,OAAO,GAAG,GAAG,IAAI,CAAC,OAAO,GAAG,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,WAAW,CAAC;IACvD,CAAC;IACD,gEAAgE;IACzD,aAAa,CAAC,QAAgB;QACnC,MAAM,CAAC,GAAG,QAAQ,GAAG,IAAI,CAAC,WAAW,CAAC;QACtC,OAAO,GAAG,GAAG,IAAI,CAAC,OAAO,GAAG,CAAC,GAAG,IAAI,CAAC,WAAW,GAAG,IAAI,CAAC,WAAW,CAAC;IACtE,CAAC;IACD,+DAA+D;IACxD,aAAa,CAAC,SAAiB;QACpC,OAAO,GAAG,GAAG,IAAI,CAAC,OAAO,GAAG,IAAI,CAAC,WAAW,GAAG,IAAI,CAAC,WAAW,GAAG,IAAI,CAAC,WAAW,CAAC;IACrF,CAAC;IACD,8BAA8B;IACvB,WAAW,CAAC,CAAS,IAAY,OAAO,CAAC,GAAG,IAAI,CAAC,WAAW,CAAC,CAAC,CAAC;CACvE;AAtDD,wCAsDC","sourcesContent":["/*---------------------------------------------------------------------------------------------\n* Copyright (c) Bentley Systems, Incorporated. All rights reserved.\n* See LICENSE.md in the project root for license terms and full copyright notice.\n*--------------------------------------------------------------------------------------------*/\n\n/** @packageDocumentation\n * @module Curve\n */\nimport { XYCurveEvaluator } from \"./XYCurveEvaluator\";\n/** Intermediate class for evaluation of bare cubic spirals `y = m ^x^3` with x axis from [0..xLength]\n * * This implements all the computations among fraction, x, and y.\n * * Derived classes implement specialized logic such as (a) precomputing `m` and (b) domain-specific fraction-to-distance approximations.\n * @internal\n */\nexport abstract class CubicEvaluator extends XYCurveEvaluator {\n  protected _cubicM: number;\n  protected _axisLength: number;\n  protected constructor(axisLength: number, cubicM: number) {\n    super();\n    this._cubicM = cubicM;\n    this._axisLength = axisLength;\n  }\n  /** Update both constants. */\n  public setConstants(axisLength: number, cubicM: number) {\n    this._axisLength = axisLength;\n    this._cubicM = cubicM;\n  }\n  public get axisLength() { return this._axisLength; }\n  public get cubicM() { return this._cubicM; }\n  /**\n   * Apply `scaleFactor` to the xLength and cubicM.\n   * * Derived classes commonly call this as `super.scaleInPlace()`, and additionally apply the scale to their members.\n   * @param scaleFactor\n   */\n  public scaleInPlace(scaleFactor: number) {\n    this._axisLength *= scaleFactor;\n    // \"x\" arriving at \"m * x^3\" will be scaled. \"m\" has to be divided by the scale to cancel 2 of the 3 . .\n    this._cubicM /= (scaleFactor * scaleFactor);\n  }\n  /** Evaluate X at fraction. */\n  public fractionToX(fraction: number): number { return fraction * this._axisLength; }\n  /** Evaluate derivative of X with respect to fraction */\n  public fractionToDX(_fraction: number): number { return this._axisLength; }\n  /** Evaluate second derivative of X with respect to fraction */\n  public fractionToDDX(_fraction: number): number { return 0.0; }\n  /** Evaluate third derivative of X with respect to fraction */\n  public fractionToD3X(_fraction: number): number { return 0.0; }\n  /** Evaluate Y at fraction. */\n  public fractionToY(fraction: number): number {\n    const x = fraction * this._axisLength;\n    return this._cubicM * x * x * x;\n  }\n  /** Evaluate derivative of Y with respect to fraction. */\n  public fractionToDY(fraction: number): number {\n    const x = fraction * this._axisLength;\n    return 3.0 * this._cubicM * x * x * this._axisLength;\n  }\n  /** Evaluate second derivative of Y with respect to fraction. */\n  public fractionToDDY(fraction: number): number {\n    const x = fraction * this._axisLength;\n    return 6.0 * this._cubicM * x * this._axisLength * this._axisLength;\n  }\n  /** Evaluate third derivative of Y with respect to fraction. */\n  public fractionToD3Y(_fraction: number): number {\n    return 6.0 * this._cubicM * this._axisLength * this._axisLength * this._axisLength;\n  }\n  /** Evaluate fraction at x. */\n  public xToFraction(x: number): number { return x / this._axisLength; }\n}\n"]}
+{"version":3,"file":"CubicEvaluator.js","sourceRoot":"","sources":["../../../../src/curve/spiral/CubicEvaluator.ts"],"names":[],"mappings":";AAAA;;;+FAG+F;;;AAE/F;;GAEG;AACH,yDAAsD;AAEtD;;;;;;GAMG;AACH,MAAsB,cAAe,SAAQ,mCAAgB;IACjD,OAAO,CAAS;IAChB,WAAW,CAAS;IAC9B,YAAsB,UAAkB,EAAE,MAAc;QACtD,KAAK,EAAE,CAAC;QACR,IAAI,CAAC,OAAO,GAAG,MAAM,CAAC;QACtB,IAAI,CAAC,WAAW,GAAG,UAAU,CAAC;IAChC,CAAC;IACD,6BAA6B;IACtB,YAAY,CAAC,UAAkB,EAAE,MAAc;QACpD,IAAI,CAAC,WAAW,GAAG,UAAU,CAAC;QAC9B,IAAI,CAAC,OAAO,GAAG,MAAM,CAAC;IACxB,CAAC;IACD,IAAW,UAAU;QACnB,OAAO,IAAI,CAAC,WAAW,CAAC;IAC1B,CAAC;IACD,IAAW,MAAM;QACf,OAAO,IAAI,CAAC,OAAO,CAAC;IACtB,CAAC;IACD;;;OAGG;IACI,YAAY,CAAC,WAAmB;QACrC,IAAI,CAAC,WAAW,IAAI,WAAW,CAAC;QAChC,oGAAoG;QACpG,IAAI,CAAC,OAAO,IAAI,CAAC,WAAW,GAAG,WAAW,CAAC,CAAC;IAC9C,CAAC;IACD,8BAA8B;IACvB,WAAW,CAAC,QAAgB;QACjC,OAAO,QAAQ,GAAG,IAAI,CAAC,WAAW,CAAC;IACrC,CAAC;IACD,yDAAyD;IAClD,YAAY,CAAC,SAAiB;QACnC,OAAO,IAAI,CAAC,WAAW,CAAC;IAC1B,CAAC;IACD,gEAAgE;IACzD,aAAa,CAAC,SAAiB;QACpC,OAAO,GAAG,CAAC;IACb,CAAC;IACD,+DAA+D;IACxD,aAAa,CAAC,SAAiB;QACpC,OAAO,GAAG,CAAC;IACb,CAAC;IACD,8BAA8B;IACvB,WAAW,CAAC,QAAgB;QACjC,MAAM,CAAC,GAAG,QAAQ,GAAG,IAAI,CAAC,WAAW,CAAC;QACtC,OAAO,IAAI,CAAC,OAAO,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;IAClC,CAAC;IACD,yDAAyD;IAClD,YAAY,CAAC,QAAgB;QAClC,MAAM,CAAC,GAAG,QAAQ,GAAG,IAAI,CAAC,WAAW,CAAC;QACtC,OAAO,GAAG,GAAG,IAAI,CAAC,OAAO,GAAG,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,WAAW,CAAC;IACvD,CAAC;IACD,gEAAgE;IACzD,aAAa,CAAC,QAAgB;QACnC,MAAM,CAAC,GAAG,QAAQ,GAAG,IAAI,CAAC,WAAW,CAAC;QACtC,OAAO,GAAG,GAAG,IAAI,CAAC,OAAO,GAAG,CAAC,GAAG,IAAI,CAAC,WAAW,GAAG,IAAI,CAAC,WAAW,CAAC;IACtE,CAAC;IACD,+DAA+D;IACxD,aAAa,CAAC,SAAiB;QACpC,OAAO,GAAG,GAAG,IAAI,CAAC,OAAO,GAAG,IAAI,CAAC,WAAW,GAAG,IAAI,CAAC,WAAW,GAAG,IAAI,CAAC,WAAW,CAAC;IACrF,CAAC;IACD,8BAA8B;IACvB,WAAW,CAAC,CAAS;QAC1B,OAAO,CAAC,GAAG,IAAI,CAAC,WAAW,CAAC;IAC9B,CAAC;CACF;AAnED,wCAmEC","sourcesContent":["/*---------------------------------------------------------------------------------------------\n* Copyright (c) Bentley Systems, Incorporated. All rights reserved.\n* See LICENSE.md in the project root for license terms and full copyright notice.\n*--------------------------------------------------------------------------------------------*/\n\n/** @packageDocumentation\n * @module Curve\n */\nimport { XYCurveEvaluator } from \"./XYCurveEvaluator\";\n\n/**\n * Intermediate class for evaluation of bare cubic spirals `y = m*x^3` with x axis from [0..xLength]\n * * This implements all the computations among fraction, x, and y.\n * * Derived classes implement specialized logic such as (a) precomputing `m` and (b) domain-specific fraction-to-distance\n * approximations.\n * @internal\n */\nexport abstract class CubicEvaluator extends XYCurveEvaluator {\n  protected _cubicM: number;\n  protected _axisLength: number;\n  protected constructor(axisLength: number, cubicM: number) {\n    super();\n    this._cubicM = cubicM;\n    this._axisLength = axisLength;\n  }\n  /** Update both constants. */\n  public setConstants(axisLength: number, cubicM: number) {\n    this._axisLength = axisLength;\n    this._cubicM = cubicM;\n  }\n  public get axisLength() {\n    return this._axisLength;\n  }\n  public get cubicM() {\n    return this._cubicM;\n  }\n  /**\n   * Apply `scaleFactor` to the xLength and cubicM.\n   * * Derived classes commonly call this as `super.scaleInPlace()`, and additionally apply the scale to their members.\n   */\n  public scaleInPlace(scaleFactor: number) {\n    this._axisLength *= scaleFactor;\n    // \"x\" arriving at \"m * x^3\" will be scaled. \"m\" has to be divided by the scale to cancel 2 of the 3\n    this._cubicM /= (scaleFactor * scaleFactor);\n  }\n  /** Evaluate X at fraction. */\n  public fractionToX(fraction: number): number {\n    return fraction * this._axisLength;\n  }\n  /** Evaluate derivative of X with respect to fraction. */\n  public fractionToDX(_fraction: number): number {\n    return this._axisLength;\n  }\n  /** Evaluate second derivative of X with respect to fraction. */\n  public fractionToDDX(_fraction: number): number {\n    return 0.0;\n  }\n  /** Evaluate third derivative of X with respect to fraction. */\n  public fractionToD3X(_fraction: number): number {\n    return 0.0;\n  }\n  /** Evaluate Y at fraction. */\n  public fractionToY(fraction: number): number {\n    const x = fraction * this._axisLength;\n    return this._cubicM * x * x * x;\n  }\n  /** Evaluate derivative of Y with respect to fraction. */\n  public fractionToDY(fraction: number): number {\n    const x = fraction * this._axisLength;\n    return 3.0 * this._cubicM * x * x * this._axisLength;\n  }\n  /** Evaluate second derivative of Y with respect to fraction. */\n  public fractionToDDY(fraction: number): number {\n    const x = fraction * this._axisLength;\n    return 6.0 * this._cubicM * x * this._axisLength * this._axisLength;\n  }\n  /** Evaluate third derivative of Y with respect to fraction. */\n  public fractionToD3Y(_fraction: number): number {\n    return 6.0 * this._cubicM * this._axisLength * this._axisLength * this._axisLength;\n  }\n  /** Evaluate fraction at x. */\n  public xToFraction(x: number): number {\n    return x / this._axisLength;\n  }\n}\n"]}