@itwin/core-geometry 5.3.0-dev.16 → 5.3.0-dev.18

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

@@ -11,21 +11,17 @@ import { SimpleNewton } from "../../numerics/Newton";
 /**
  * Czech cubic.
  * This is y= m*x^3 with
- * * x any point on the x axis
- * * `fraction` along the spiral goes to `x = fraction * L`
- * * m is gamma / (6RL)
+ * * x any point on the x axis.
+ * * `fraction` along the spiral goes to `x = fraction * L`.
+ * * m is gamma / (6RL).
  *    * 1/(6RL) is the leading term of the sine series.
  *    * `gamma = 2R/sqrt (4RR-LL)` pushes y down a little bit to simulate the lost series terms.
- * @param localToWorld
- * @param nominalL1
- * @param nominalR1
- * @param activeInterval
  * @internal
  */
 export class CzechSpiralEvaluator extends CubicEvaluator {
     nominalLength1;
     nominalRadius1;
-    /** Constructor is private.  Caller responsible for cubicM validity. */
+    /** Constructor is private. Caller responsible for cubicM validity. */
     constructor(length1, radius1, cubicM) {
         super(length1, cubicM);
         this.nominalLength1 = length1;
@@ -35,8 +31,6 @@ export class CzechSpiralEvaluator extends CubicEvaluator {
      * Return the scale factor between simple x^3 / (6RL) cubic and the czech correction.
      * * For typical case with l1/R1 smallish, this is just less than 1.0:
      *   (0.25==>0.99215), (0.15==>0.997184), (0.10==>0.998749), (0.05==>999687)
-     * @param length1
-     * @param radius1
      */
     static gammaConstant(length1, radius1) {
         return 2.0 * radius1 / Math.sqrt(4.0 * radius1 * radius1 - length1 * length1);
@@ -44,10 +38,10 @@ export class CzechSpiralEvaluator extends CubicEvaluator {
     /** Compute the czech cubic constant. */
     static computeCubicM(length1, radius1) {
         const gamma = CzechSpiralEvaluator.gammaConstant(length1, radius1);
-        // In the private update method, the LR values should have been vetted.
+        // in the private update method, the LR values should have been vetted
         if (gamma === undefined)
             return undefined;
-        // If radius is negative, it shows up in gamma.  But the a signed denominator undoes it.  So take abs of denominator.
+        // if radius is negative, it shows up in gamma; but the a signed denominator undoes it so take abs of denominator
         return gamma / Math.abs((6.0 * radius1 * length1));
     }
     static create(length1, radius1) {
@@ -61,9 +55,11 @@ export class CzechSpiralEvaluator extends CubicEvaluator {
         this.nominalRadius1 *= scaleFactor;
         super.scaleInPlace(scaleFactor);
     }
-    /** return a deep copy of the evaluator */
-    clone() { return new CzechSpiralEvaluator(this.nominalLength1, this.nominalRadius1, this.cubicM); }
-    /** Member by member matchup ... */
+    /** Return a deep copy of the evaluator. */
+    clone() {
+        return new CzechSpiralEvaluator(this.nominalLength1, this.nominalRadius1, this.cubicM);
+    }
+    /** Member by member matchup. */
     isAlmostEqual(other) {
         if (other instanceof CzechSpiralEvaluator) {
             return Geometry.isSameCoordinate(this.nominalLength1, other.nominalLength1)
@@ -74,7 +70,8 @@ export class CzechSpiralEvaluator extends CubicEvaluator {
     /**
      * Return a (fast but mediocre) approximation of spiral length as a function of x axis position.
      * * This x-to-distance relation is not as precise as the CurvePrimitive method moveSignedDistanceFromFraction.
-     * * It is supported here for users interested in replicating the Czech distance mapping rather than the more accurate CurvePrimitive measurements.
+     * * It is supported here for users interested in replicating the Czech distance mapping rather than the more accurate
+     * CurvePrimitive measurements.
      * @param x distance along the x axis.
      */
     xToCzechApproximateDistance(x) {
@@ -83,16 +80,15 @@ export class CzechSpiralEvaluator extends CubicEvaluator {
     /**
      * Return the inverse of the `forwardL2R2Map` function.
      * * The undefined result can only occur for distances outside the usual spirals.
-     * @param s (approximate) distance along the spiral.
-     *
+     * @param d (approximate) distance along the spiral.
      */
     czechApproximateDistanceToX(d) {
         return CzechSpiralEvaluator.inverseL2R2Map(d, 1.0, this.nominalLength1, this.nominalRadius1);
     }
     /**
-       * evaluate a series expansion that is used with varying signs (plus or minus 1) in czech and italian spirals.
-       * @param x distance along the x axis.
-       */
+     * Evaluate a series expansion that is used with varying signs (plus or minus 1) in czech and italian spirals.
+     * @param x distance along the x axis.
+     */
     static forwardL2R2Map(x, sign, length, radius) {
         const l2 = length * length;
         const r2 = radius * radius;
@@ -103,17 +99,16 @@ export class CzechSpiralEvaluator extends CubicEvaluator {
     /**
      * Return the inverse of the `forwardL2R2Map` function.
      * * The undefined result can only occur for distances outside the usual spirals.
-     * @param s (approximate) distance along the spiral.
-     *
+     * @param d (approximate) distance along the spiral.
      */
-    static inverseL2R2Map(b, sign, length, radius) {
+    static inverseL2R2Map(d, sign, length, radius) {
         const l2 = length * length;
         const r2 = radius * radius;
         const Q = 4.0 * r2 - l2;
         const a = sign / (10.0 * Q * l2);
-        return SimpleNewton.runNewton1D(b, (x) => {
+        return SimpleNewton.runNewton1D(d, (x) => {
             const xx = x * x;
-            return x * (1.0 + xx * xx * a) - b;
+            return x * (1.0 + xx * xx * a) - d;
         }, (x) => {
             const xx = x * x;
             return 1.0 + 5 * xx * xx * a;
@@ -123,39 +118,36 @@ export class CzechSpiralEvaluator extends CubicEvaluator {
 /**
  * Italian cubic.
  * This is y= m*x^3 with
- * * x any point on the x axis
- * * `fraction` along the spiral goes to `x = fraction * L`
- * * m is gamma / (6RL)
+ * * x any point on the x axis.
+ * * `fraction` along the spiral goes to `x = fraction * L`.
+ * * m is gamma / (6RL).
  *    * 1/(6RL) is the leading term of the sine series.
  *    * `gamma = 2R/sqrt (4RR-LL)` pushes y down a little bit to simulate the lost series terms.
- * @param localToWorld
- * @param nominalL1
- * @param nominalR1
- * @param activeInterval
  * @internal
  */
 export class ItalianSpiralEvaluator extends CubicEvaluator {
     nominalLength1;
     nominalRadius1;
-    /** Compute the czech cubic constant.
-     * ** funky mixture of lengths ....
+    /**
+     * Compute the czech cubic constant.
+     * * Funky mixture of lengths.
      */
     static computeCubicM(lengthXByForward, radius1) {
         const gamma = CzechSpiralEvaluator.gammaConstant(lengthXByForward, radius1);
-        // In the private update method, the LR values should have been vetted.
+        // in the private update method, the LR values should have been vetted
         if (gamma === undefined)
             return undefined;
-        // If radius is negative, it shows up in gamma.  But the a signed denominator undoes it.  So take abs of denominator.
+        // if radius is negative, it shows up in gamma; but the a signed denominator undoes it so take abs of denominator
         return gamma / Math.abs((6.0 * radius1 * lengthXByForward));
     }
-    /** Constructor is private.  Caller responsible for cubicM validity. */
+    /** Constructor is private. Caller responsible for cubicM validity. */
     constructor(length1, radius1, lengthX, cubicM) {
         super(lengthX, cubicM);
         this.nominalLength1 = length1;
         this.nominalRadius1 = radius1;
     }
     static create(length1, radius1) {
-        // um ... this seems goofy.  lengthX from forward, then invert for another .   But that's what the native code does ...
+        // this seems goofy.; lengthX from forward, then invert for another but that's what the native code does too
         const lengthX = CzechSpiralEvaluator.forwardL2R2Map(length1, -1.0, length1, radius1);
         const lengthX1 = CzechSpiralEvaluator.inverseL2R2Map(length1, 1.0, lengthX, radius1);
         if (lengthX1 === undefined)
@@ -170,9 +162,11 @@ export class ItalianSpiralEvaluator extends CubicEvaluator {
         this.nominalRadius1 *= scaleFactor;
         super.scaleInPlace(scaleFactor);
     }
-    /** return a deep copy of the evaluator */
-    clone() { return new ItalianSpiralEvaluator(this.nominalLength1, this.nominalRadius1, super.axisLength, this.cubicM); }
-    /** Member by member matchup ... */
+    /** Return a deep copy of the evaluator. */
+    clone() {
+        return new ItalianSpiralEvaluator(this.nominalLength1, this.nominalRadius1, super.axisLength, this.cubicM);
+    }
+    /** Member by member matchup. */
     isAlmostEqual(other) {
         if (other instanceof ItalianSpiralEvaluator) {
             return Geometry.isSameCoordinate(this.nominalLength1, other.nominalLength1)
@@ -181,19 +175,19 @@ export class ItalianSpiralEvaluator extends CubicEvaluator {
         return false;
     }
     /**
-       * Return a (fast but mediocre) approximation of spiral length as a function of x axis position.
-       * * This x-to-distance relation is not as precise as the CurvePrimitive method moveSignedDistanceFromFraction.
-       * * It is supported here for users interested in replicating the Czech distance mapping rather than the more accurate CurvePrimitive measurements.
-       * @param x distance along the x axis.
-       */
+     * Return a (fast but mediocre) approximation of spiral length as a function of x axis position.
+     * * This x-to-distance relation is not as precise as the CurvePrimitive method moveSignedDistanceFromFraction.
+     * * It is supported here for users interested in replicating the Czech distance mapping rather than the more accurate
+     * CurvePrimitive measurements.
+     * @param x distance along the x axis.
+     */
     distanceToItalianApproximateX(x) {
         return CzechSpiralEvaluator.forwardL2R2Map(x, -1.0, this.nominalLength1, this.nominalRadius1);
     }
     /**
      * Return the inverse of the `forwardL2R2Map` function.
      * * The undefined result can only occur for distances outside the usual spirals.
-     * @param s (approximate) distance along the spiral.
-     *
+     * @param d (approximate) distance along the spiral.
      */
     xToItalianApproximateDistance(d) {
         return CzechSpiralEvaluator.inverseL2R2Map(d, -1.0, this.nominalLength1, this.nominalRadius1);

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

@@ -1 +1 @@
-{"version":3,"file":"CzechSpiralEvaluator.js","sourceRoot":"","sources":["../../../../src/curve/spiral/CzechSpiralEvaluator.ts"],"names":[],"mappings":"AAAA;;;+FAG+F;AAE/F;;GAEG;AACH,OAAO,EAAE,QAAQ,EAAE,MAAM,gBAAgB,CAAC;AAC1C,OAAO,EAAE,cAAc,EAAE,MAAM,kBAAkB,CAAC;AAClD,OAAO,EAAE,YAAY,EAAE,MAAM,uBAAuB,CAAC;AACrD;;;;;;;;;;;;;GAaG;AACH,MAAM,OAAO,oBAAqB,SAAQ,cAAc;IAC/C,cAAc,CAAS;IACvB,cAAc,CAAS;IAE9B,uEAAuE;IACvE,YAAoB,OAAe,EAAE,OAAe,EAAE,MAAc;QAClE,KAAK,CAAC,OAAO,EAAE,MAAM,CAAC,CAAC;QACvB,IAAI,CAAC,cAAc,GAAG,OAAO,CAAC;QAC9B,IAAI,CAAC,cAAc,GAAG,OAAO,CAAC;IAChC,CAAC;IACD;;;;;;OAMG;IACI,MAAM,CAAC,aAAa,CAAC,OAAe,EAAE,OAAe;QAC1D,OAAO,GAAG,GAAG,OAAO,GAAG,IAAI,CAAC,IAAI,CAAC,GAAG,GAAG,OAAO,GAAG,OAAO,GAAG,OAAO,GAAG,OAAO,CAAC,CAAC;IAChF,CAAC;IACD,wCAAwC;IACjC,MAAM,CAAC,aAAa,CAAC,OAAe,EAAE,OAAe;QAC1D,MAAM,KAAK,GAAG,oBAAoB,CAAC,aAAa,CAAC,OAAO,EAAE,OAAO,CAAC,CAAC;QACnE,uEAAuE;QACvE,IAAI,KAAK,KAAK,SAAS;YACrB,OAAO,SAAS,CAAC;QACnB,qHAAqH;QACrH,OAAO,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,OAAO,GAAG,OAAO,CAAC,CAAC,CAAC;IACrD,CAAC;IAEM,MAAM,CAAC,MAAM,CAAC,OAAe,EAAE,OAAe;QACnD,MAAM,CAAC,GAAG,IAAI,CAAC,aAAa,CAAC,OAAO,EAAE,OAAO,CAAC,CAAC;QAC/C,IAAI,CAAC,KAAK,SAAS;YACjB,OAAO,SAAS,CAAC;QACnB,OAAO,IAAI,oBAAoB,CAAC,OAAO,EAAE,OAAO,EAAE,CAAC,CAAC,CAAC;IACvD,CAAC;IAEe,YAAY,CAAC,WAAmB;QAC9C,IAAI,CAAC,cAAc,IAAI,WAAW,CAAC;QACnC,IAAI,CAAC,cAAc,IAAI,WAAW,CAAC;QACnC,KAAK,CAAC,YAAY,CAAC,WAAW,CAAC,CAAC;IAClC,CAAC;IACD,0CAA0C;IACnC,KAAK,KAA2B,OAAO,IAAI,oBAAoB,CAAC,IAAI,CAAC,cAAc,EAAE,IAAI,CAAC,cAAc,EAAE,IAAI,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC;IAChI,mCAAmC;IAC5B,aAAa,CAAC,KAAU;QAC7B,IAAI,KAAK,YAAY,oBAAoB,EAAE,CAAC;YAC1C,OAAO,QAAQ,CAAC,gBAAgB,CAAC,IAAI,CAAC,cAAc,EAAE,KAAK,CAAC,cAAc,CAAC;mBACtE,QAAQ,CAAC,gBAAgB,CAAC,IAAI,CAAC,cAAc,EAAE,KAAK,CAAC,cAAc,CAAC,CAAC;QAC5E,CAAC;QACD,OAAO,KAAK,CAAC;IACf,CAAC;IAED;;;;;OAKG;IACI,2BAA2B,CAAC,CAAS;QAC1C,OAAO,oBAAoB,CAAC,cAAc,CAAC,CAAC,EAAE,GAAG,EAAE,IAAI,CAAC,cAAc,EAAE,IAAI,CAAC,cAAc,CAAC,CAAC;IAC/F,CAAC;IACD;;;;;OAKG;IACI,2BAA2B,CAAC,CAAS;QAC1C,OAAO,oBAAoB,CAAC,cAAc,CAAC,CAAC,EAAE,GAAG,EAAE,IAAI,CAAC,cAAc,EAAE,IAAI,CAAC,cAAc,CAAC,CAAC;IAC/F,CAAC;IAED;;;SAGK;IACE,MAAM,CAAC,cAAc,CAAC,CAAS,EAAE,IAAY,EAAE,MAAc,EAAE,MAAc;QAClF,MAAM,EAAE,GAAG,MAAM,GAAG,MAAM,CAAC;QAC3B,MAAM,EAAE,GAAG,MAAM,GAAG,MAAM,CAAC;QAC3B,MAAM,CAAC,GAAG,GAAG,GAAG,EAAE,GAAG,EAAE,CAAC;QACxB,MAAM,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC;QACjB,OAAO,CAAC,GAAG,CAAC,GAAG,GAAG,IAAI,GAAG,EAAE,GAAG,EAAE,GAAG,CAAC,IAAI,GAAG,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC;IACtD,CAAC;IACD;;;;;OAKG;IACI,MAAM,CAAC,cAAc,CAAC,CAAS,EAAE,IAAY,EAAE,MAAc,EAAE,MAAc;QAClF,MAAM,EAAE,GAAG,MAAM,GAAG,MAAM,CAAC;QAC3B,MAAM,EAAE,GAAG,MAAM,GAAG,MAAM,CAAC;QAC3B,MAAM,CAAC,GAAG,GAAG,GAAG,EAAE,GAAG,EAAE,CAAC;QACxB,MAAM,CAAC,GAAG,IAAI,GAAG,CAAC,IAAI,GAAG,CAAC,GAAG,EAAE,CAAC,CAAC;QACjC,OAAO,YAAY,CAAC,WAAW,CAAC,CAAC,EAC/B,CAAC,CAAS,EAAE,EAAE;YACZ,MAAM,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC;YACjB,OAAO,CAAC,GAAG,CAAC,GAAG,GAAG,EAAE,GAAG,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC;QACrC,CAAC,EACD,CAAC,CAAS,EAAE,EAAE;YACZ,MAAM,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC;YACjB,OAAO,GAAG,GAAG,CAAC,GAAG,EAAE,GAAG,EAAE,GAAG,CAAC,CAAC;QAC/B,CAAC,CAAC,CAAC;IACP,CAAC;CACF;AACD;;;;;;;;;;;;;GAaG;AACH,MAAM,OAAO,sBAAuB,SAAQ,cAAc;IACjD,cAAc,CAAS;IACvB,cAAc,CAAS;IAE9B;;OAEG;IACK,MAAM,CAAC,aAAa,CAAC,gBAAwB,EAAE,OAAe;QACpE,MAAM,KAAK,GAAG,oBAAoB,CAAC,aAAa,CAAC,gBAAgB,EAAE,OAAO,CAAC,CAAC;QAC5E,uEAAuE;QACvE,IAAI,KAAK,KAAK,SAAS;YACrB,OAAO,SAAS,CAAC;QACnB,qHAAqH;QACrH,OAAO,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,OAAO,GAAG,gBAAgB,CAAC,CAAC,CAAC;IAC9D,CAAC;IAED,uEAAuE;IACvE,YAAoB,OAAe,EAAE,OAAe,EAAE,OAAe,EAAE,MAAc;QACnF,KAAK,CAAC,OAAO,EAAE,MAAM,CAAC,CAAC;QACvB,IAAI,CAAC,cAAc,GAAG,OAAO,CAAC;QAC9B,IAAI,CAAC,cAAc,GAAG,OAAO,CAAC;IAChC,CAAC;IACM,MAAM,CAAC,MAAM,CAAC,OAAe,EAAE,OAAe;QACnD,uHAAuH;QACvH,MAAM,OAAO,GAAG,oBAAoB,CAAC,cAAc,CAAC,OAAO,EAAE,CAAC,GAAG,EAAE,OAAO,EAAE,OAAO,CAAC,CAAC;QACrF,MAAM,QAAQ,GAAG,oBAAoB,CAAC,cAAc,CAAC,OAAO,EAAE,GAAG,EAAE,OAAO,EAAE,OAAO,CAAC,CAAC;QACrF,IAAI,QAAQ,KAAK,SAAS;YACxB,OAAO,SAAS,CAAC;QACnB,MAAM,CAAC,GAAG,sBAAsB,CAAC,aAAa,CAAC,OAAO,EAAE,OAAO,CAAC,CAAC;QACjE,IAAI,CAAC,KAAK,SAAS;YACjB,OAAO,SAAS,CAAC;QACnB,OAAO,IAAI,sBAAsB,CAAC,OAAO,EAAE,OAAO,EAAE,QAAQ,EAAE,CAAC,CAAC,CAAC;IACnE,CAAC;IAEe,YAAY,CAAC,WAAmB;QAC9C,IAAI,CAAC,cAAc,IAAI,WAAW,CAAC;QACnC,IAAI,CAAC,cAAc,IAAI,WAAW,CAAC;QACnC,KAAK,CAAC,YAAY,CAAC,WAAW,CAAC,CAAC;IAClC,CAAC;IACD,0CAA0C;IACnC,KAAK,KAA6B,OAAO,IAAI,sBAAsB,CAAC,IAAI,CAAC,cAAc,EAAE,IAAI,CAAC,cAAc,EAAE,KAAK,CAAC,UAAU,EAAE,IAAI,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC;IACtJ,mCAAmC;IAC5B,aAAa,CAAC,KAAU;QAC7B,IAAI,KAAK,YAAY,sBAAsB,EAAE,CAAC;YAC5C,OAAO,QAAQ,CAAC,gBAAgB,CAAC,IAAI,CAAC,cAAc,EAAE,KAAK,CAAC,cAAc,CAAC;mBACtE,QAAQ,CAAC,gBAAgB,CAAC,IAAI,CAAC,cAAc,EAAE,KAAK,CAAC,cAAc,CAAC,CAAC;QAC5E,CAAC;QACD,OAAO,KAAK,CAAC;IACf,CAAC;IACD;;;;;SAKK;IACE,6BAA6B,CAAC,CAAS;QAC5C,OAAO,oBAAoB,CAAC,cAAc,CAAC,CAAC,EAAE,CAAC,GAAG,EAAE,IAAI,CAAC,cAAc,EAAE,IAAI,CAAC,cAAc,CAAC,CAAC;IAChG,CAAC;IACD;;;;;OAKG;IACI,6BAA6B,CAAC,CAAS;QAC5C,OAAO,oBAAoB,CAAC,cAAc,CAAC,CAAC,EAAE,CAAC,GAAG,EAAE,IAAI,CAAC,cAAc,EAAE,IAAI,CAAC,cAAc,CAAC,CAAC;IAChG,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 { CubicEvaluator } from \"./CubicEvaluator\";\r\nimport { SimpleNewton } from \"../../numerics/Newton\";\r\n/**\r\n * Czech cubic.\r\n * This is y= m*x^3 with\r\n * * x any point on the x axis\r\n * * `fraction` along the spiral goes to `x = fraction * L`\r\n * * m is gamma / (6RL)\r\n *    * 1/(6RL) is the leading term of the sine series.\r\n *    * `gamma = 2R/sqrt (4RR-LL)` pushes y down a little bit to simulate the lost series terms.\r\n * @param localToWorld\r\n * @param nominalL1\r\n * @param nominalR1\r\n * @param activeInterval\r\n * @internal\r\n */\r\nexport class CzechSpiralEvaluator extends CubicEvaluator {\r\n  public nominalLength1: number;\r\n  public nominalRadius1: number;\r\n\r\n  /** Constructor is private.  Caller responsible for cubicM validity. */\r\n  private constructor(length1: number, radius1: number, cubicM: number) {\r\n    super(length1, cubicM);\r\n    this.nominalLength1 = length1;\r\n    this.nominalRadius1 = radius1;\r\n  }\r\n  /**\r\n   * Return the scale factor between simple x^3 / (6RL) cubic and the czech correction.\r\n   * * For typical case with l1/R1 smallish, this is just less than 1.0:\r\n   *   (0.25==>0.99215), (0.15==>0.997184), (0.10==>0.998749), (0.05==>999687)\r\n   * @param length1\r\n   * @param radius1\r\n   */\r\n  public static gammaConstant(length1: number, radius1: number): number | undefined {\r\n    return 2.0 * radius1 / Math.sqrt(4.0 * radius1 * radius1 - length1 * length1);\r\n  }\r\n  /** Compute the czech cubic constant. */\r\n  public static computeCubicM(length1: number, radius1: number): number | undefined {\r\n    const gamma = CzechSpiralEvaluator.gammaConstant(length1, radius1);\r\n    // In the private update method, the LR values should have been vetted.\r\n    if (gamma === undefined)\r\n      return undefined;\r\n    // If radius is negative, it shows up in gamma.  But the a signed denominator undoes it.  So take abs of denominator.\r\n    return gamma / Math.abs((6.0 * radius1 * length1));\r\n  }\r\n\r\n  public static create(length1: number, radius1: number): CzechSpiralEvaluator | undefined {\r\n    const m = this.computeCubicM(length1, radius1);\r\n    if (m === undefined)\r\n      return undefined;\r\n    return new CzechSpiralEvaluator(length1, radius1, m);\r\n  }\r\n\r\n  public override scaleInPlace(scaleFactor: number) {\r\n    this.nominalLength1 *= scaleFactor;\r\n    this.nominalRadius1 *= scaleFactor;\r\n    super.scaleInPlace(scaleFactor);\r\n  }\r\n  /** return a deep copy of the evaluator */\r\n  public clone(): CzechSpiralEvaluator { return new CzechSpiralEvaluator(this.nominalLength1, this.nominalRadius1, this.cubicM); }\r\n  /** Member by member matchup ... */\r\n  public isAlmostEqual(other: any): boolean {\r\n    if (other instanceof CzechSpiralEvaluator) {\r\n      return Geometry.isSameCoordinate(this.nominalLength1, other.nominalLength1)\r\n        && Geometry.isSameCoordinate(this.nominalRadius1, other.nominalRadius1);\r\n    }\r\n    return false;\r\n  }\r\n\r\n  /**\r\n   * Return a (fast but mediocre) approximation of spiral length as a function of x axis position.\r\n   * * This x-to-distance relation is not as precise as the CurvePrimitive method moveSignedDistanceFromFraction.\r\n   * * It is supported here for users interested in replicating the Czech distance mapping rather than the more accurate CurvePrimitive measurements.\r\n   * @param x distance along the x axis.\r\n   */\r\n  public xToCzechApproximateDistance(x: number): number {\r\n    return CzechSpiralEvaluator.forwardL2R2Map(x, 1.0, this.nominalLength1, this.nominalRadius1);\r\n  }\r\n  /**\r\n   * Return the inverse of the `forwardL2R2Map` function.\r\n   * * The undefined result can only occur for distances outside the usual spirals.\r\n   * @param s (approximate) distance along the spiral.\r\n   *\r\n   */\r\n  public czechApproximateDistanceToX(d: number): number | undefined {\r\n    return CzechSpiralEvaluator.inverseL2R2Map(d, 1.0, this.nominalLength1, this.nominalRadius1);\r\n  }\r\n\r\n  /**\r\n     * evaluate a series expansion that is used with varying signs (plus or minus 1) in czech and italian spirals.\r\n     * @param x distance along the x axis.\r\n     */\r\n  public static forwardL2R2Map(x: number, sign: number, length: number, radius: number): number {\r\n    const l2 = length * length;\r\n    const r2 = radius * radius;\r\n    const Q = 4.0 * r2 - l2;\r\n    const xx = x * x;\r\n    return x * (1.0 + sign * xx * xx / (10.0 * Q * l2));\r\n  }\r\n  /**\r\n   * Return the inverse of the `forwardL2R2Map` function.\r\n   * * The undefined result can only occur for distances outside the usual spirals.\r\n   * @param s (approximate) distance along the spiral.\r\n   *\r\n   */\r\n  public static inverseL2R2Map(b: number, sign: number, length: number, radius: number): number | undefined {\r\n    const l2 = length * length;\r\n    const r2 = radius * radius;\r\n    const Q = 4.0 * r2 - l2;\r\n    const a = sign / (10.0 * Q * l2);\r\n    return SimpleNewton.runNewton1D(b,\r\n      (x: number) => {\r\n        const xx = x * x;\r\n        return x * (1.0 + xx * xx * a) - b;\r\n      },\r\n      (x: number) => {\r\n        const xx = x * x;\r\n        return 1.0 + 5 * xx * xx * a;\r\n      });\r\n  }\r\n}\r\n/**\r\n * Italian cubic.\r\n * This is y= m*x^3 with\r\n * * x any point on the x axis\r\n * * `fraction` along the spiral goes to `x = fraction * L`\r\n * * m is gamma / (6RL)\r\n *    * 1/(6RL) is the leading term of the sine series.\r\n *    * `gamma = 2R/sqrt (4RR-LL)` pushes y down a little bit to simulate the lost series terms.\r\n * @param localToWorld\r\n * @param nominalL1\r\n * @param nominalR1\r\n * @param activeInterval\r\n * @internal\r\n */\r\nexport class ItalianSpiralEvaluator extends CubicEvaluator {\r\n  public nominalLength1: number;\r\n  public nominalRadius1: number;\r\n\r\n  /** Compute the czech cubic constant.\r\n   * ** funky mixture of lengths ....\r\n   */\r\n  private static computeCubicM(lengthXByForward: number, radius1: number): number | undefined {\r\n    const gamma = CzechSpiralEvaluator.gammaConstant(lengthXByForward, radius1);\r\n    // In the private update method, the LR values should have been vetted.\r\n    if (gamma === undefined)\r\n      return undefined;\r\n    // If radius is negative, it shows up in gamma.  But the a signed denominator undoes it.  So take abs of denominator.\r\n    return gamma / Math.abs((6.0 * radius1 * lengthXByForward));\r\n  }\r\n\r\n  /** Constructor is private.  Caller responsible for cubicM validity. */\r\n  private constructor(length1: number, radius1: number, lengthX: number, cubicM: number) {\r\n    super(lengthX, cubicM);\r\n    this.nominalLength1 = length1;\r\n    this.nominalRadius1 = radius1;\r\n  }\r\n  public static create(length1: number, radius1: number): ItalianSpiralEvaluator | undefined {\r\n    // um ... this seems goofy.  lengthX from forward, then invert for another .   But that's what the native code does ...\r\n    const lengthX = CzechSpiralEvaluator.forwardL2R2Map(length1, -1.0, length1, radius1);\r\n    const lengthX1 = CzechSpiralEvaluator.inverseL2R2Map(length1, 1.0, lengthX, radius1);\r\n    if (lengthX1 === undefined)\r\n      return undefined;\r\n    const m = ItalianSpiralEvaluator.computeCubicM(lengthX, radius1);\r\n    if (m === undefined)\r\n      return undefined;\r\n    return new ItalianSpiralEvaluator(length1, radius1, lengthX1, m);\r\n  }\r\n\r\n  public override scaleInPlace(scaleFactor: number) {\r\n    this.nominalLength1 *= scaleFactor;\r\n    this.nominalRadius1 *= scaleFactor;\r\n    super.scaleInPlace(scaleFactor);\r\n  }\r\n  /** return a deep copy of the evaluator */\r\n  public clone(): ItalianSpiralEvaluator { return new ItalianSpiralEvaluator(this.nominalLength1, this.nominalRadius1, super.axisLength, this.cubicM); }\r\n  /** Member by member matchup ... */\r\n  public isAlmostEqual(other: any): boolean {\r\n    if (other instanceof ItalianSpiralEvaluator) {\r\n      return Geometry.isSameCoordinate(this.nominalLength1, other.nominalLength1)\r\n        && Geometry.isSameCoordinate(this.nominalRadius1, other.nominalRadius1);\r\n    }\r\n    return false;\r\n  }\r\n  /**\r\n     * Return a (fast but mediocre) approximation of spiral length as a function of x axis position.\r\n     * * This x-to-distance relation is not as precise as the CurvePrimitive method moveSignedDistanceFromFraction.\r\n     * * It is supported here for users interested in replicating the Czech distance mapping rather than the more accurate CurvePrimitive measurements.\r\n     * @param x distance along the x axis.\r\n     */\r\n  public distanceToItalianApproximateX(x: number): number {\r\n    return CzechSpiralEvaluator.forwardL2R2Map(x, -1.0, this.nominalLength1, this.nominalRadius1);\r\n  }\r\n  /**\r\n   * Return the inverse of the `forwardL2R2Map` function.\r\n   * * The undefined result can only occur for distances outside the usual spirals.\r\n   * @param s (approximate) distance along the spiral.\r\n   *\r\n   */\r\n  public xToItalianApproximateDistance(d: number): number | undefined {\r\n    return CzechSpiralEvaluator.inverseL2R2Map(d, -1.0, this.nominalLength1, this.nominalRadius1);\r\n  }\r\n}\r\n"]}
+{"version":3,"file":"CzechSpiralEvaluator.js","sourceRoot":"","sources":["../../../../src/curve/spiral/CzechSpiralEvaluator.ts"],"names":[],"mappings":"AAAA;;;+FAG+F;AAE/F;;GAEG;AACH,OAAO,EAAE,QAAQ,EAAE,MAAM,gBAAgB,CAAC;AAC1C,OAAO,EAAE,cAAc,EAAE,MAAM,kBAAkB,CAAC;AAClD,OAAO,EAAE,YAAY,EAAE,MAAM,uBAAuB,CAAC;AAErD;;;;;;;;;GASG;AACH,MAAM,OAAO,oBAAqB,SAAQ,cAAc;IAC/C,cAAc,CAAS;IACvB,cAAc,CAAS;IAC9B,sEAAsE;IACtE,YAAoB,OAAe,EAAE,OAAe,EAAE,MAAc;QAClE,KAAK,CAAC,OAAO,EAAE,MAAM,CAAC,CAAC;QACvB,IAAI,CAAC,cAAc,GAAG,OAAO,CAAC;QAC9B,IAAI,CAAC,cAAc,GAAG,OAAO,CAAC;IAChC,CAAC;IACD;;;;OAIG;IACI,MAAM,CAAC,aAAa,CAAC,OAAe,EAAE,OAAe;QAC1D,OAAO,GAAG,GAAG,OAAO,GAAG,IAAI,CAAC,IAAI,CAAC,GAAG,GAAG,OAAO,GAAG,OAAO,GAAG,OAAO,GAAG,OAAO,CAAC,CAAC;IAChF,CAAC;IACD,wCAAwC;IACjC,MAAM,CAAC,aAAa,CAAC,OAAe,EAAE,OAAe;QAC1D,MAAM,KAAK,GAAG,oBAAoB,CAAC,aAAa,CAAC,OAAO,EAAE,OAAO,CAAC,CAAC;QACnE,sEAAsE;QACtE,IAAI,KAAK,KAAK,SAAS;YACrB,OAAO,SAAS,CAAC;QACnB,iHAAiH;QACjH,OAAO,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,OAAO,GAAG,OAAO,CAAC,CAAC,CAAC;IACrD,CAAC;IACM,MAAM,CAAC,MAAM,CAAC,OAAe,EAAE,OAAe;QACnD,MAAM,CAAC,GAAG,IAAI,CAAC,aAAa,CAAC,OAAO,EAAE,OAAO,CAAC,CAAC;QAC/C,IAAI,CAAC,KAAK,SAAS;YACjB,OAAO,SAAS,CAAC;QACnB,OAAO,IAAI,oBAAoB,CAAC,OAAO,EAAE,OAAO,EAAE,CAAC,CAAC,CAAC;IACvD,CAAC;IACe,YAAY,CAAC,WAAmB;QAC9C,IAAI,CAAC,cAAc,IAAI,WAAW,CAAC;QACnC,IAAI,CAAC,cAAc,IAAI,WAAW,CAAC;QACnC,KAAK,CAAC,YAAY,CAAC,WAAW,CAAC,CAAC;IAClC,CAAC;IACD,2CAA2C;IACpC,KAAK;QACV,OAAO,IAAI,oBAAoB,CAAC,IAAI,CAAC,cAAc,EAAE,IAAI,CAAC,cAAc,EAAE,IAAI,CAAC,MAAM,CAAC,CAAC;IACzF,CAAC;IACD,gCAAgC;IACzB,aAAa,CAAC,KAAU;QAC7B,IAAI,KAAK,YAAY,oBAAoB,EAAE,CAAC;YAC1C,OAAO,QAAQ,CAAC,gBAAgB,CAAC,IAAI,CAAC,cAAc,EAAE,KAAK,CAAC,cAAc,CAAC;mBACtE,QAAQ,CAAC,gBAAgB,CAAC,IAAI,CAAC,cAAc,EAAE,KAAK,CAAC,cAAc,CAAC,CAAC;QAC5E,CAAC;QACD,OAAO,KAAK,CAAC;IACf,CAAC;IACD;;;;;;OAMG;IACI,2BAA2B,CAAC,CAAS;QAC1C,OAAO,oBAAoB,CAAC,cAAc,CAAC,CAAC,EAAE,GAAG,EAAE,IAAI,CAAC,cAAc,EAAE,IAAI,CAAC,cAAc,CAAC,CAAC;IAC/F,CAAC;IACD;;;;OAIG;IACI,2BAA2B,CAAC,CAAS;QAC1C,OAAO,oBAAoB,CAAC,cAAc,CAAC,CAAC,EAAE,GAAG,EAAE,IAAI,CAAC,cAAc,EAAE,IAAI,CAAC,cAAc,CAAC,CAAC;IAC/F,CAAC;IACD;;;OAGG;IACI,MAAM,CAAC,cAAc,CAAC,CAAS,EAAE,IAAY,EAAE,MAAc,EAAE,MAAc;QAClF,MAAM,EAAE,GAAG,MAAM,GAAG,MAAM,CAAC;QAC3B,MAAM,EAAE,GAAG,MAAM,GAAG,MAAM,CAAC;QAC3B,MAAM,CAAC,GAAG,GAAG,GAAG,EAAE,GAAG,EAAE,CAAC;QACxB,MAAM,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC;QACjB,OAAO,CAAC,GAAG,CAAC,GAAG,GAAG,IAAI,GAAG,EAAE,GAAG,EAAE,GAAG,CAAC,IAAI,GAAG,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC;IACtD,CAAC;IACD;;;;OAIG;IACI,MAAM,CAAC,cAAc,CAAC,CAAS,EAAE,IAAY,EAAE,MAAc,EAAE,MAAc;QAClF,MAAM,EAAE,GAAG,MAAM,GAAG,MAAM,CAAC;QAC3B,MAAM,EAAE,GAAG,MAAM,GAAG,MAAM,CAAC;QAC3B,MAAM,CAAC,GAAG,GAAG,GAAG,EAAE,GAAG,EAAE,CAAC;QACxB,MAAM,CAAC,GAAG,IAAI,GAAG,CAAC,IAAI,GAAG,CAAC,GAAG,EAAE,CAAC,CAAC;QACjC,OAAO,YAAY,CAAC,WAAW,CAC7B,CAAC,EACD,CAAC,CAAS,EAAE,EAAE;YACZ,MAAM,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC;YACjB,OAAO,CAAC,GAAG,CAAC,GAAG,GAAG,EAAE,GAAG,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC;QACrC,CAAC,EACD,CAAC,CAAS,EAAE,EAAE;YACZ,MAAM,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC;YACjB,OAAO,GAAG,GAAG,CAAC,GAAG,EAAE,GAAG,EAAE,GAAG,CAAC,CAAC;QAC/B,CAAC,CACF,CAAC;IACJ,CAAC;CACF;AAED;;;;;;;;;GASG;AACH,MAAM,OAAO,sBAAuB,SAAQ,cAAc;IACjD,cAAc,CAAS;IACvB,cAAc,CAAS;IAE9B;;;OAGG;IACK,MAAM,CAAC,aAAa,CAAC,gBAAwB,EAAE,OAAe;QACpE,MAAM,KAAK,GAAG,oBAAoB,CAAC,aAAa,CAAC,gBAAgB,EAAE,OAAO,CAAC,CAAC;QAC5E,sEAAsE;QACtE,IAAI,KAAK,KAAK,SAAS;YACrB,OAAO,SAAS,CAAC;QACnB,iHAAiH;QACjH,OAAO,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,OAAO,GAAG,gBAAgB,CAAC,CAAC,CAAC;IAC9D,CAAC;IACD,sEAAsE;IACtE,YAAoB,OAAe,EAAE,OAAe,EAAE,OAAe,EAAE,MAAc;QACnF,KAAK,CAAC,OAAO,EAAE,MAAM,CAAC,CAAC;QACvB,IAAI,CAAC,cAAc,GAAG,OAAO,CAAC;QAC9B,IAAI,CAAC,cAAc,GAAG,OAAO,CAAC;IAChC,CAAC;IACM,MAAM,CAAC,MAAM,CAAC,OAAe,EAAE,OAAe;QACnD,4GAA4G;QAC5G,MAAM,OAAO,GAAG,oBAAoB,CAAC,cAAc,CAAC,OAAO,EAAE,CAAC,GAAG,EAAE,OAAO,EAAE,OAAO,CAAC,CAAC;QACrF,MAAM,QAAQ,GAAG,oBAAoB,CAAC,cAAc,CAAC,OAAO,EAAE,GAAG,EAAE,OAAO,EAAE,OAAO,CAAC,CAAC;QACrF,IAAI,QAAQ,KAAK,SAAS;YACxB,OAAO,SAAS,CAAC;QACnB,MAAM,CAAC,GAAG,sBAAsB,CAAC,aAAa,CAAC,OAAO,EAAE,OAAO,CAAC,CAAC;QACjE,IAAI,CAAC,KAAK,SAAS;YACjB,OAAO,SAAS,CAAC;QACnB,OAAO,IAAI,sBAAsB,CAAC,OAAO,EAAE,OAAO,EAAE,QAAQ,EAAE,CAAC,CAAC,CAAC;IACnE,CAAC;IACe,YAAY,CAAC,WAAmB;QAC9C,IAAI,CAAC,cAAc,IAAI,WAAW,CAAC;QACnC,IAAI,CAAC,cAAc,IAAI,WAAW,CAAC;QACnC,KAAK,CAAC,YAAY,CAAC,WAAW,CAAC,CAAC;IAClC,CAAC;IACD,2CAA2C;IACpC,KAAK;QACV,OAAO,IAAI,sBAAsB,CAAC,IAAI,CAAC,cAAc,EAAE,IAAI,CAAC,cAAc,EAAE,KAAK,CAAC,UAAU,EAAE,IAAI,CAAC,MAAM,CAAC,CAAC;IAC7G,CAAC;IACD,gCAAgC;IACzB,aAAa,CAAC,KAAU;QAC7B,IAAI,KAAK,YAAY,sBAAsB,EAAE,CAAC;YAC5C,OAAO,QAAQ,CAAC,gBAAgB,CAAC,IAAI,CAAC,cAAc,EAAE,KAAK,CAAC,cAAc,CAAC;mBACtE,QAAQ,CAAC,gBAAgB,CAAC,IAAI,CAAC,cAAc,EAAE,KAAK,CAAC,cAAc,CAAC,CAAC;QAC5E,CAAC;QACD,OAAO,KAAK,CAAC;IACf,CAAC;IACD;;;;;;OAMG;IACI,6BAA6B,CAAC,CAAS;QAC5C,OAAO,oBAAoB,CAAC,cAAc,CAAC,CAAC,EAAE,CAAC,GAAG,EAAE,IAAI,CAAC,cAAc,EAAE,IAAI,CAAC,cAAc,CAAC,CAAC;IAChG,CAAC;IACD;;;;OAIG;IACI,6BAA6B,CAAC,CAAS;QAC5C,OAAO,oBAAoB,CAAC,cAAc,CAAC,CAAC,EAAE,CAAC,GAAG,EAAE,IAAI,CAAC,cAAc,EAAE,IAAI,CAAC,cAAc,CAAC,CAAC;IAChG,CAAC;CACF","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 { CubicEvaluator } from \"./CubicEvaluator\";\nimport { SimpleNewton } from \"../../numerics/Newton\";\n\n/**\n * Czech cubic.\n * This is y= m*x^3 with\n * * x any point on the x axis.\n * * `fraction` along the spiral goes to `x = fraction * L`.\n * * m is gamma / (6RL).\n *    * 1/(6RL) is the leading term of the sine series.\n *    * `gamma = 2R/sqrt (4RR-LL)` pushes y down a little bit to simulate the lost series terms.\n * @internal\n */\nexport class CzechSpiralEvaluator extends CubicEvaluator {\n  public nominalLength1: number;\n  public nominalRadius1: number;\n  /** Constructor is private. Caller responsible for cubicM validity. */\n  private constructor(length1: number, radius1: number, cubicM: number) {\n    super(length1, cubicM);\n    this.nominalLength1 = length1;\n    this.nominalRadius1 = radius1;\n  }\n  /**\n   * Return the scale factor between simple x^3 / (6RL) cubic and the czech correction.\n   * * For typical case with l1/R1 smallish, this is just less than 1.0:\n   *   (0.25==>0.99215), (0.15==>0.997184), (0.10==>0.998749), (0.05==>999687)\n   */\n  public static gammaConstant(length1: number, radius1: number): number | undefined {\n    return 2.0 * radius1 / Math.sqrt(4.0 * radius1 * radius1 - length1 * length1);\n  }\n  /** Compute the czech cubic constant. */\n  public static computeCubicM(length1: number, radius1: number): number | undefined {\n    const gamma = CzechSpiralEvaluator.gammaConstant(length1, radius1);\n    // in the private update method, the LR values should have been vetted\n    if (gamma === undefined)\n      return undefined;\n    // if radius is negative, it shows up in gamma; but the a signed denominator undoes it so take abs of denominator\n    return gamma / Math.abs((6.0 * radius1 * length1));\n  }\n  public static create(length1: number, radius1: number): CzechSpiralEvaluator | undefined {\n    const m = this.computeCubicM(length1, radius1);\n    if (m === undefined)\n      return undefined;\n    return new CzechSpiralEvaluator(length1, radius1, m);\n  }\n  public override scaleInPlace(scaleFactor: number) {\n    this.nominalLength1 *= scaleFactor;\n    this.nominalRadius1 *= scaleFactor;\n    super.scaleInPlace(scaleFactor);\n  }\n  /** Return a deep copy of the evaluator. */\n  public clone(): CzechSpiralEvaluator {\n    return new CzechSpiralEvaluator(this.nominalLength1, this.nominalRadius1, this.cubicM);\n  }\n  /** Member by member matchup. */\n  public isAlmostEqual(other: any): boolean {\n    if (other instanceof CzechSpiralEvaluator) {\n      return Geometry.isSameCoordinate(this.nominalLength1, other.nominalLength1)\n        && Geometry.isSameCoordinate(this.nominalRadius1, other.nominalRadius1);\n    }\n    return false;\n  }\n  /**\n   * Return a (fast but mediocre) approximation of spiral length as a function of x axis position.\n   * * This x-to-distance relation is not as precise as the CurvePrimitive method moveSignedDistanceFromFraction.\n   * * It is supported here for users interested in replicating the Czech distance mapping rather than the more accurate\n   * CurvePrimitive measurements.\n   * @param x distance along the x axis.\n   */\n  public xToCzechApproximateDistance(x: number): number {\n    return CzechSpiralEvaluator.forwardL2R2Map(x, 1.0, this.nominalLength1, this.nominalRadius1);\n  }\n  /**\n   * Return the inverse of the `forwardL2R2Map` function.\n   * * The undefined result can only occur for distances outside the usual spirals.\n   * @param d (approximate) distance along the spiral.\n   */\n  public czechApproximateDistanceToX(d: number): number | undefined {\n    return CzechSpiralEvaluator.inverseL2R2Map(d, 1.0, this.nominalLength1, this.nominalRadius1);\n  }\n  /**\n   * Evaluate a series expansion that is used with varying signs (plus or minus 1) in czech and italian spirals.\n   * @param x distance along the x axis.\n   */\n  public static forwardL2R2Map(x: number, sign: number, length: number, radius: number): number {\n    const l2 = length * length;\n    const r2 = radius * radius;\n    const Q = 4.0 * r2 - l2;\n    const xx = x * x;\n    return x * (1.0 + sign * xx * xx / (10.0 * Q * l2));\n  }\n  /**\n   * Return the inverse of the `forwardL2R2Map` function.\n   * * The undefined result can only occur for distances outside the usual spirals.\n   * @param d (approximate) distance along the spiral.\n   */\n  public static inverseL2R2Map(d: number, sign: number, length: number, radius: number): number | undefined {\n    const l2 = length * length;\n    const r2 = radius * radius;\n    const Q = 4.0 * r2 - l2;\n    const a = sign / (10.0 * Q * l2);\n    return SimpleNewton.runNewton1D(\n      d,\n      (x: number) => {\n        const xx = x * x;\n        return x * (1.0 + xx * xx * a) - d;\n      },\n      (x: number) => {\n        const xx = x * x;\n        return 1.0 + 5 * xx * xx * a;\n      },\n    );\n  }\n}\n\n/**\n * Italian cubic.\n * This is y= m*x^3 with\n * * x any point on the x axis.\n * * `fraction` along the spiral goes to `x = fraction * L`.\n * * m is gamma / (6RL).\n *    * 1/(6RL) is the leading term of the sine series.\n *    * `gamma = 2R/sqrt (4RR-LL)` pushes y down a little bit to simulate the lost series terms.\n * @internal\n */\nexport class ItalianSpiralEvaluator extends CubicEvaluator {\n  public nominalLength1: number;\n  public nominalRadius1: number;\n\n  /**\n   * Compute the czech cubic constant.\n   * * Funky mixture of lengths.\n   */\n  private static computeCubicM(lengthXByForward: number, radius1: number): number | undefined {\n    const gamma = CzechSpiralEvaluator.gammaConstant(lengthXByForward, radius1);\n    // in the private update method, the LR values should have been vetted\n    if (gamma === undefined)\n      return undefined;\n    // if radius is negative, it shows up in gamma; but the a signed denominator undoes it so take abs of denominator\n    return gamma / Math.abs((6.0 * radius1 * lengthXByForward));\n  }\n  /** Constructor is private. Caller responsible for cubicM validity. */\n  private constructor(length1: number, radius1: number, lengthX: number, cubicM: number) {\n    super(lengthX, cubicM);\n    this.nominalLength1 = length1;\n    this.nominalRadius1 = radius1;\n  }\n  public static create(length1: number, radius1: number): ItalianSpiralEvaluator | undefined {\n    // this seems goofy.; lengthX from forward, then invert for another but that's what the native code does too\n    const lengthX = CzechSpiralEvaluator.forwardL2R2Map(length1, -1.0, length1, radius1);\n    const lengthX1 = CzechSpiralEvaluator.inverseL2R2Map(length1, 1.0, lengthX, radius1);\n    if (lengthX1 === undefined)\n      return undefined;\n    const m = ItalianSpiralEvaluator.computeCubicM(lengthX, radius1);\n    if (m === undefined)\n      return undefined;\n    return new ItalianSpiralEvaluator(length1, radius1, lengthX1, m);\n  }\n  public override scaleInPlace(scaleFactor: number) {\n    this.nominalLength1 *= scaleFactor;\n    this.nominalRadius1 *= scaleFactor;\n    super.scaleInPlace(scaleFactor);\n  }\n  /** Return a deep copy of the evaluator. */\n  public clone(): ItalianSpiralEvaluator {\n    return new ItalianSpiralEvaluator(this.nominalLength1, this.nominalRadius1, super.axisLength, this.cubicM);\n  }\n  /** Member by member matchup. */\n  public isAlmostEqual(other: any): boolean {\n    if (other instanceof ItalianSpiralEvaluator) {\n      return Geometry.isSameCoordinate(this.nominalLength1, other.nominalLength1)\n        && Geometry.isSameCoordinate(this.nominalRadius1, other.nominalRadius1);\n    }\n    return false;\n  }\n  /**\n   * Return a (fast but mediocre) approximation of spiral length as a function of x axis position.\n   * * This x-to-distance relation is not as precise as the CurvePrimitive method moveSignedDistanceFromFraction.\n   * * It is supported here for users interested in replicating the Czech distance mapping rather than the more accurate\n   * CurvePrimitive measurements.\n   * @param x distance along the x axis.\n   */\n  public distanceToItalianApproximateX(x: number): number {\n    return CzechSpiralEvaluator.forwardL2R2Map(x, -1.0, this.nominalLength1, this.nominalRadius1);\n  }\n  /**\n   * Return the inverse of the `forwardL2R2Map` function.\n   * * The undefined result can only occur for distances outside the usual spirals.\n   * @param d (approximate) distance along the spiral.\n   */\n  public xToItalianApproximateDistance(d: number): number | undefined {\n    return CzechSpiralEvaluator.inverseL2R2Map(d, -1.0, this.nominalLength1, this.nominalRadius1);\n  }\n}\n"]}

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

@@ -11,9 +11,9 @@ export declare class DirectHalfCosineSpiralEvaluator extends XYCurveEvaluator {
     constructor(length1: number, radius1: number);
     private updateConstants;
     scaleInPlace(scaleFactor: number): void;
-    /** return a deep copy of the evaluator */
+    /** Return a deep copy of the evaluator. */
     clone(): DirectHalfCosineSpiralEvaluator;
-    /** Member by member matchup ... */
+    /** Member by member matchup. */
     isAlmostEqual(other: any): boolean;
     /** Evaluate X at fractional position. */
     fractionToX(fraction: number): number;

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

@@ -1 +1 @@
-{"version":3,"file":"DirectHalfCosineSpiralEvaluator.d.ts","sourceRoot":"","sources":["../../../../src/curve/spiral/DirectHalfCosineSpiralEvaluator.ts"],"names":[],"mappings":"AASA,OAAO,EAAE,gBAAgB,EAAE,MAAM,oBAAoB,CAAC;AAEtD;;GAEG;AACH,qBAAa,+BAAgC,SAAQ,gBAAgB;IAC5D,cAAc,EAAE,MAAM,CAAC;IACvB,cAAc,EAAE,MAAM,CAAC;IAC9B,OAAO,CAAC,EAAE,CAAS;IACnB,OAAO,CAAC,GAAG,CAAS;IACpB,OAAO,CAAC,GAAG,CAAS;gBACD,OAAO,EAAE,MAAM,EAAE,OAAO,EAAE,MAAM;IAWnD,OAAO,CAAC,eAAe;IAGhB,YAAY,CAAC,WAAW,EAAE,MAAM;IAKvC,0CAA0C;IACnC,KAAK,IAAI,+BAA+B;IAC/C,mCAAmC;IAC5B,aAAa,CAAC,KAAK,EAAE,GAAG,GAAG,OAAO;IAQzC,yCAAyC;IAClC,WAAW,CAAC,QAAQ,EAAE,MAAM,GAAG,MAAM;IAC5C,yCAAyC;IAClC,WAAW,CAAC,QAAQ,EAAE,MAAM,GAAG,MAAM;IAI5C,gFAAgF;IACzE,YAAY,CAAC,SAAS,EAAE,MAAM,GAAG,MAAM;IAE9C,gFAAgF;IAEzE,YAAY,CAAC,QAAQ,EAAE,MAAM,GAAG,MAAM;IAK7C,uFAAuF;IAChF,aAAa,CAAC,SAAS,EAAE,MAAM,GAAG,MAAM;IAE/C,sFAAsF;IAE/E,aAAa,CAAC,QAAQ,EAAE,MAAM,GAAG,MAAM;IAK9C,uFAAuF;IAChF,aAAa,CAAC,SAAS,EAAE,MAAM,GAAG,MAAM;IAE/C,sFAAsF;IAE/E,aAAa,CAAC,QAAQ,EAAE,MAAM,GAAG,MAAM;IAM9C,yEAAyE;IAEzD,0BAA0B,CAAC,QAAQ,EAAE,MAAM,GAAG,MAAM;IAIpE,mDAAmD;IAC5C,WAAW,CAAC,CAAC,EAAE,MAAM,GAAG,MAAM,GAAG,SAAS;CAUlD"}
+{"version":3,"file":"DirectHalfCosineSpiralEvaluator.d.ts","sourceRoot":"","sources":["../../../../src/curve/spiral/DirectHalfCosineSpiralEvaluator.ts"],"names":[],"mappings":"AASA,OAAO,EAAE,gBAAgB,EAAE,MAAM,oBAAoB,CAAC;AAGtD;;GAEG;AACH,qBAAa,+BAAgC,SAAQ,gBAAgB;IAC5D,cAAc,EAAE,MAAM,CAAC;IACvB,cAAc,EAAE,MAAM,CAAC;IAC9B,OAAO,CAAC,EAAE,CAAS;IACnB,OAAO,CAAC,GAAG,CAAS;IACpB,OAAO,CAAC,GAAG,CAAS;gBACD,OAAO,EAAE,MAAM,EAAE,OAAO,EAAE,MAAM;IAUnD,OAAO,CAAC,eAAe;IAGhB,YAAY,CAAC,WAAW,EAAE,MAAM;IAKvC,2CAA2C;IACpC,KAAK,IAAI,+BAA+B;IAG/C,gCAAgC;IACzB,aAAa,CAAC,KAAK,EAAE,GAAG,GAAG,OAAO;IAQzC,yCAAyC;IAClC,WAAW,CAAC,QAAQ,EAAE,MAAM,GAAG,MAAM;IAG5C,yCAAyC;IAClC,WAAW,CAAC,QAAQ,EAAE,MAAM,GAAG,MAAM;IAI5C,gFAAgF;IACzE,YAAY,CAAC,SAAS,EAAE,MAAM,GAAG,MAAM;IAG9C,gFAAgF;IACzE,YAAY,CAAC,QAAQ,EAAE,MAAM,GAAG,MAAM;IAK7C,uFAAuF;IAChF,aAAa,CAAC,SAAS,EAAE,MAAM,GAAG,MAAM;IAG/C,sFAAsF;IAC/E,aAAa,CAAC,QAAQ,EAAE,MAAM,GAAG,MAAM;IAK9C,uFAAuF;IAChF,aAAa,CAAC,SAAS,EAAE,MAAM,GAAG,MAAM;IAG/C,sFAAsF;IAC/E,aAAa,CAAC,QAAQ,EAAE,MAAM,GAAG,MAAM;IAK9C,yEAAyE;IACzD,0BAA0B,CAAC,QAAQ,EAAE,MAAM,GAAG,MAAM;IAGpE,mDAAmD;IAC5C,WAAW,CAAC,CAAC,EAAE,MAAM,GAAG,MAAM,GAAG,SAAS;CAWlD"}

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

@@ -35,9 +35,11 @@ export class DirectHalfCosineSpiralEvaluator extends XYCurveEvaluator {
         this.nominalRadius1 *= scaleFactor;
         this.updateConstants();
     }
-    /** return a deep copy of the evaluator */
-    clone() { return new DirectHalfCosineSpiralEvaluator(this.nominalLength1, this.nominalRadius1); }
-    /** Member by member matchup ... */
+    /** Return a deep copy of the evaluator. */
+    clone() {
+        return new DirectHalfCosineSpiralEvaluator(this.nominalLength1, this.nominalRadius1);
+    }
+    /** Member by member matchup. */
     isAlmostEqual(other) {
         if (other instanceof DirectHalfCosineSpiralEvaluator) {
             return Geometry.isSameCoordinate(this.nominalLength1, other.nominalLength1)
@@ -47,14 +49,18 @@ export class DirectHalfCosineSpiralEvaluator extends XYCurveEvaluator {
         return false;
     }
     /** Evaluate X at fractional position. */
-    fractionToX(fraction) { return fraction * this.nominalLength1; }
+    fractionToX(fraction) {
+        return fraction * this.nominalLength1;
+    }
     /** Evaluate Y at fractional position. */
     fractionToY(fraction) {
         const theta = fraction * Math.PI;
         return this._c * (this._c2 * fraction * fraction - this._c1 * (1.0 - Math.cos(theta)));
     }
     /** Evaluate derivative of X with respect to fraction at fractional position. */
-    fractionToDX(_fraction) { return this.nominalLength1; }
+    fractionToDX(_fraction) {
+        return this.nominalLength1;
+    }
     /** Evaluate derivative of Y with respect to fraction at fractional position. */
     fractionToDY(fraction) {
         const pi = Math.PI;
@@ -62,7 +68,9 @@ export class DirectHalfCosineSpiralEvaluator extends XYCurveEvaluator {
         return this._c * (2.0 * this._c2 * fraction - this._c1 * pi * Math.sin(theta));
     }
     /** Evaluate second derivative of X with respect to fraction at fractional position. */
-    fractionToDDX(_fraction) { return 0.0; }
+    fractionToDDX(_fraction) {
+        return 0.0;
+    }
     /** Evaluate third derivative of Y with respect to fraction at fractional position. */
     fractionToDDY(fraction) {
         const pi = Math.PI;
@@ -70,7 +78,9 @@ export class DirectHalfCosineSpiralEvaluator extends XYCurveEvaluator {
         return this._c * (2.0 * this._c2 - this._c1 * pi * pi * Math.cos(theta));
     }
     /** Evaluate second derivative of X with respect to fraction at fractional position. */
-    fractionToD3X(_fraction) { return 0.0; }
+    fractionToD3X(_fraction) {
+        return 0.0;
+    }
     /** Evaluate third derivative of Y with respect to fraction at fractional position. */
     fractionToD3Y(fraction) {
         const pi = Math.PI;

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

@@ -1 +1 @@
-{"version":3,"file":"DirectHalfCosineSpiralEvaluator.js","sourceRoot":"","sources":["../../../../src/curve/spiral/DirectHalfCosineSpiralEvaluator.ts"],"names":[],"mappings":"AAAA;;;+FAG+F;AAE/F;;GAEG;AACH,OAAO,EAAE,QAAQ,EAAE,MAAM,gBAAgB,CAAC;AAC1C,OAAO,EAAE,gBAAgB,EAAE,MAAM,oBAAoB,CAAC;AACtD,OAAO,EAAE,YAAY,EAAE,MAAM,uBAAuB,CAAC;AACrD;;GAEG;AACH,MAAM,OAAO,+BAAgC,SAAQ,gBAAgB;IAC5D,cAAc,CAAS;IACvB,cAAc,CAAS;IACtB,EAAE,CAAS;IACX,GAAG,CAAS;IACZ,GAAG,CAAS;IACpB,YAAmB,OAAe,EAAE,OAAe;QACjD,KAAK,EAAE,CAAC;QACR,IAAI,CAAC,cAAc,GAAG,OAAO,CAAC;QAC9B,IAAI,CAAC,cAAc,GAAG,OAAO,CAAC;QAC9B,MAAM,EAAE,GAAG,IAAI,CAAC,EAAE,CAAC;QACnB,IAAI,CAAC,GAAG,GAAG,GAAG,GAAG,CAAC,GAAG,GAAG,EAAE,GAAG,EAAE,CAAC,CAAC;QACjC,IAAI,CAAC,GAAG,GAAG,IAAI,CAAC;QAChB,IAAI,CAAC,EAAE,GAAG,GAAG,CAAC,CAAC,sBAAsB;QACrC,IAAI,CAAC,eAAe,EAAE,CAAC;IACzB,CAAC;IAEO,eAAe;QACrB,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,cAAc,GAAG,IAAI,CAAC,cAAc,GAAG,IAAI,CAAC,cAAc,CAAC;IAC5E,CAAC;IACM,YAAY,CAAC,WAAmB;QACrC,IAAI,CAAC,cAAc,IAAI,WAAW,CAAC;QACnC,IAAI,CAAC,cAAc,IAAI,WAAW,CAAC;QACnC,IAAI,CAAC,eAAe,EAAE,CAAC;IACzB,CAAC;IACD,0CAA0C;IACnC,KAAK,KAAsC,OAAO,IAAI,+BAA+B,CAAC,IAAI,CAAC,cAAc,EAAE,IAAI,CAAC,cAAc,CAAC,CAAC,CAAC,CAAC;IACzI,mCAAmC;IAC5B,aAAa,CAAC,KAAU;QAC7B,IAAI,KAAK,YAAY,+BAA+B,EAAE,CAAC;YACrD,OAAO,QAAQ,CAAC,gBAAgB,CAAC,IAAI,CAAC,cAAc,EAAE,KAAK,CAAC,cAAc,CAAC;mBACtE,QAAQ,CAAC,gBAAgB,CAAC,IAAI,CAAC,cAAc,EAAE,KAAK,CAAC,cAAc,CAAC,CAAC;YAC1E,oDAAoD;QACtD,CAAC;QACD,OAAO,KAAK,CAAC;IACf,CAAC;IACD,yCAAyC;IAClC,WAAW,CAAC,QAAgB,IAAY,OAAO,QAAQ,GAAG,IAAI,CAAC,cAAc,CAAC,CAAC,CAAC;IACvF,yCAAyC;IAClC,WAAW,CAAC,QAAgB;QACjC,MAAM,KAAK,GAAG,QAAQ,GAAG,IAAI,CAAC,EAAE,CAAC;QACjC,OAAO,IAAI,CAAC,EAAE,GAAG,CAAC,IAAI,CAAC,GAAG,GAAG,QAAQ,GAAG,QAAQ,GAAG,IAAI,CAAC,GAAG,GAAG,CAAC,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC;IACzF,CAAC;IACD,gFAAgF;IACzE,YAAY,CAAC,SAAiB,IAAY,OAAO,IAAI,CAAC,cAAc,CAAC,CAAC,CAAC;IAE9E,gFAAgF;IAEzE,YAAY,CAAC,QAAgB;QAClC,MAAM,EAAE,GAAG,IAAI,CAAC,EAAE,CAAC;QACnB,MAAM,KAAK,GAAG,QAAQ,GAAG,EAAE,CAAC;QAC5B,OAAO,IAAI,CAAC,EAAE,GAAG,CAAC,GAAG,GAAG,IAAI,CAAC,GAAG,GAAG,QAAQ,GAAG,IAAI,CAAC,GAAG,GAAG,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC;IACjF,CAAC;IACD,uFAAuF;IAChF,aAAa,CAAC,SAAiB,IAAY,OAAO,GAAG,CAAC,CAAC,CAAC;IAE/D,sFAAsF;IAE/E,aAAa,CAAC,QAAgB;QACnC,MAAM,EAAE,GAAG,IAAI,CAAC,EAAE,CAAC;QACnB,MAAM,KAAK,GAAG,QAAQ,GAAG,EAAE,CAAC;QAC5B,OAAO,IAAI,CAAC,EAAE,GAAG,CAAC,GAAG,GAAG,IAAI,CAAC,GAAG,GAAG,IAAI,CAAC,GAAG,GAAG,EAAE,GAAG,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC;IAC3E,CAAC;IACD,uFAAuF;IAChF,aAAa,CAAC,SAAiB,IAAY,OAAO,GAAG,CAAC,CAAC,CAAC;IAE/D,sFAAsF;IAE/E,aAAa,CAAC,QAAgB;QACnC,MAAM,EAAE,GAAG,IAAI,CAAC,EAAE,CAAC;QACnB,MAAM,KAAK,GAAG,QAAQ,GAAG,EAAE,CAAC;QAC5B,OAAO,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,GAAG,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC;IAC7D,CAAC;IAED,yEAAyE;IAEzD,0BAA0B,CAAC,QAAgB;QACzD,OAAO,QAAQ,CAAC,YAAY,CAAC,IAAI,CAAC,YAAY,CAAC,QAAQ,CAAC,EAAE,IAAI,CAAC,YAAY,CAAC,QAAQ,CAAC,CAAC,CAAC;IAEzF,CAAC;IACD,mDAAmD;IAC5C,WAAW,CAAC,CAAS;QAC1B,MAAM,SAAS,GAAG,CAAC,GAAG,IAAI,CAAC,cAAc,CAAC;QAC1C,MAAM,SAAS,GAAG,YAAY,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;IAEnB,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 { XYCurveEvaluator } from \"./XYCurveEvaluator\";\r\nimport { SimpleNewton } from \"../../numerics/Newton\";\r\n/**\r\n * @internal\r\n */\r\nexport class DirectHalfCosineSpiralEvaluator extends XYCurveEvaluator {\r\n  public nominalLength1: number;\r\n  public nominalRadius1: number;\r\n  private _c: number;\r\n  private _c1: number;\r\n  private _c2: number;\r\n  public constructor(length1: number, radius1: number) {\r\n    super();\r\n    this.nominalLength1 = length1;\r\n    this.nominalRadius1 = radius1;\r\n    const pi = Math.PI;\r\n    this._c1 = 1.0 / (2.0 * pi * pi);\r\n    this._c2 = 0.25;\r\n    this._c = 0.0; // TO BE UPDATED BELOW\r\n    this.updateConstants();\r\n  }\r\n\r\n  private updateConstants() {\r\n    this._c = this.nominalLength1 * this.nominalLength1 / this.nominalRadius1;\r\n  }\r\n  public scaleInPlace(scaleFactor: number) {\r\n    this.nominalLength1 *= scaleFactor;\r\n    this.nominalRadius1 *= scaleFactor;\r\n    this.updateConstants();\r\n  }\r\n  /** return a deep copy of the evaluator */\r\n  public clone(): DirectHalfCosineSpiralEvaluator { return new DirectHalfCosineSpiralEvaluator(this.nominalLength1, this.nominalRadius1); }\r\n  /** Member by member matchup ... */\r\n  public isAlmostEqual(other: any): boolean {\r\n    if (other instanceof DirectHalfCosineSpiralEvaluator) {\r\n      return Geometry.isSameCoordinate(this.nominalLength1, other.nominalLength1)\r\n        && Geometry.isSameCoordinate(this.nominalRadius1, other.nominalRadius1);\r\n      // remark: c,c1,c2 are computed, need not be tested.\r\n    }\r\n    return false;\r\n  }\r\n  /** Evaluate X at fractional position. */\r\n  public fractionToX(fraction: number): number { return fraction * this.nominalLength1; }\r\n  /** Evaluate Y at fractional position. */\r\n  public fractionToY(fraction: number): number {\r\n    const theta = fraction * Math.PI;\r\n    return this._c * (this._c2 * fraction * fraction - this._c1 * (1.0 - Math.cos(theta)));\r\n  }\r\n  /** Evaluate derivative of X with respect to fraction at fractional position. */\r\n  public fractionToDX(_fraction: number): number { return this.nominalLength1; }\r\n\r\n  /** Evaluate derivative of Y with respect to fraction at fractional position. */\r\n\r\n  public fractionToDY(fraction: number): number {\r\n    const pi = Math.PI;\r\n    const theta = fraction * pi;\r\n    return this._c * (2.0 * this._c2 * fraction - this._c1 * pi * Math.sin(theta));\r\n  }\r\n  /** Evaluate second derivative of X with respect to fraction at fractional position. */\r\n  public fractionToDDX(_fraction: number): number { return 0.0; }\r\n\r\n  /** Evaluate third derivative of Y with respect to fraction at fractional position. */\r\n\r\n  public fractionToDDY(fraction: number): number {\r\n    const pi = Math.PI;\r\n    const theta = fraction * pi;\r\n    return this._c * (2.0 * this._c2 - this._c1 * pi * pi * Math.cos(theta));\r\n  }\r\n  /** Evaluate second derivative of X with respect to fraction at fractional position. */\r\n  public fractionToD3X(_fraction: number): number { return 0.0; }\r\n\r\n  /** Evaluate third derivative of Y with respect to fraction at fractional position. */\r\n\r\n  public fractionToD3Y(fraction: number): number {\r\n    const pi = Math.PI;\r\n    const theta = fraction * pi;\r\n    return this._c * this._c1 * pi * pi * pi * Math.sin(theta);\r\n  }\r\n\r\n  /** Return the magnitude of the first vector at fractional coordinate. */\r\n\r\n  public override fractionToTangentMagnitude(fraction: number): number {\r\n    return Geometry.hypotenuseXY(this.fractionToDX(fraction), this.fractionToDY(fraction));\r\n\r\n  }\r\n  /** Invert the fractionToX function for given X. */\r\n  public xToFraction(x: number): number | undefined {\r\n    const fraction0 = x / this.nominalLength1;\r\n    const fraction1 = SimpleNewton.runNewton1D(fraction0,\r\n      (f: number) => (this.fractionToX(f) - x),\r\n      (f: number) => this.fractionToDX(f));\r\n    if (fraction1 === undefined)\r\n      return undefined;\r\n    return fraction1;\r\n\r\n  }\r\n}\r\n"]}
+{"version":3,"file":"DirectHalfCosineSpiralEvaluator.js","sourceRoot":"","sources":["../../../../src/curve/spiral/DirectHalfCosineSpiralEvaluator.ts"],"names":[],"mappings":"AAAA;;;+FAG+F;AAE/F;;GAEG;AACH,OAAO,EAAE,QAAQ,EAAE,MAAM,gBAAgB,CAAC;AAC1C,OAAO,EAAE,gBAAgB,EAAE,MAAM,oBAAoB,CAAC;AACtD,OAAO,EAAE,YAAY,EAAE,MAAM,uBAAuB,CAAC;AAErD;;GAEG;AACH,MAAM,OAAO,+BAAgC,SAAQ,gBAAgB;IAC5D,cAAc,CAAS;IACvB,cAAc,CAAS;IACtB,EAAE,CAAS;IACX,GAAG,CAAS;IACZ,GAAG,CAAS;IACpB,YAAmB,OAAe,EAAE,OAAe;QACjD,KAAK,EAAE,CAAC;QACR,IAAI,CAAC,cAAc,GAAG,OAAO,CAAC;QAC9B,IAAI,CAAC,cAAc,GAAG,OAAO,CAAC;QAC9B,MAAM,EAAE,GAAG,IAAI,CAAC,EAAE,CAAC;QACnB,IAAI,CAAC,GAAG,GAAG,GAAG,GAAG,CAAC,GAAG,GAAG,EAAE,GAAG,EAAE,CAAC,CAAC;QACjC,IAAI,CAAC,GAAG,GAAG,IAAI,CAAC;QAChB,IAAI,CAAC,EAAE,GAAG,GAAG,CAAC,CAAC,sBAAsB;QACrC,IAAI,CAAC,eAAe,EAAE,CAAC;IACzB,CAAC;IACO,eAAe;QACrB,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,cAAc,GAAG,IAAI,CAAC,cAAc,GAAG,IAAI,CAAC,cAAc,CAAC;IAC5E,CAAC;IACM,YAAY,CAAC,WAAmB;QACrC,IAAI,CAAC,cAAc,IAAI,WAAW,CAAC;QACnC,IAAI,CAAC,cAAc,IAAI,WAAW,CAAC;QACnC,IAAI,CAAC,eAAe,EAAE,CAAC;IACzB,CAAC;IACD,2CAA2C;IACpC,KAAK;QACV,OAAO,IAAI,+BAA+B,CAAC,IAAI,CAAC,cAAc,EAAE,IAAI,CAAC,cAAc,CAAC,CAAC;IACvF,CAAC;IACD,gCAAgC;IACzB,aAAa,CAAC,KAAU;QAC7B,IAAI,KAAK,YAAY,+BAA+B,EAAE,CAAC;YACrD,OAAO,QAAQ,CAAC,gBAAgB,CAAC,IAAI,CAAC,cAAc,EAAE,KAAK,CAAC,cAAc,CAAC;mBACtE,QAAQ,CAAC,gBAAgB,CAAC,IAAI,CAAC,cAAc,EAAE,KAAK,CAAC,cAAc,CAAC,CAAC;YAC1E,oDAAoD;QACtD,CAAC;QACD,OAAO,KAAK,CAAC;IACf,CAAC;IACD,yCAAyC;IAClC,WAAW,CAAC,QAAgB;QACjC,OAAO,QAAQ,GAAG,IAAI,CAAC,cAAc,CAAC;IACxC,CAAC;IACD,yCAAyC;IAClC,WAAW,CAAC,QAAgB;QACjC,MAAM,KAAK,GAAG,QAAQ,GAAG,IAAI,CAAC,EAAE,CAAC;QACjC,OAAO,IAAI,CAAC,EAAE,GAAG,CAAC,IAAI,CAAC,GAAG,GAAG,QAAQ,GAAG,QAAQ,GAAG,IAAI,CAAC,GAAG,GAAG,CAAC,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC;IACzF,CAAC;IACD,gFAAgF;IACzE,YAAY,CAAC,SAAiB;QACnC,OAAO,IAAI,CAAC,cAAc,CAAC;IAC7B,CAAC;IACD,gFAAgF;IACzE,YAAY,CAAC,QAAgB;QAClC,MAAM,EAAE,GAAG,IAAI,CAAC,EAAE,CAAC;QACnB,MAAM,KAAK,GAAG,QAAQ,GAAG,EAAE,CAAC;QAC5B,OAAO,IAAI,CAAC,EAAE,GAAG,CAAC,GAAG,GAAG,IAAI,CAAC,GAAG,GAAG,QAAQ,GAAG,IAAI,CAAC,GAAG,GAAG,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC;IACjF,CAAC;IACD,uFAAuF;IAChF,aAAa,CAAC,SAAiB;QACpC,OAAO,GAAG,CAAC;IACb,CAAC;IACD,sFAAsF;IAC/E,aAAa,CAAC,QAAgB;QACnC,MAAM,EAAE,GAAG,IAAI,CAAC,EAAE,CAAC;QACnB,MAAM,KAAK,GAAG,QAAQ,GAAG,EAAE,CAAC;QAC5B,OAAO,IAAI,CAAC,EAAE,GAAG,CAAC,GAAG,GAAG,IAAI,CAAC,GAAG,GAAG,IAAI,CAAC,GAAG,GAAG,EAAE,GAAG,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC;IAC3E,CAAC;IACD,uFAAuF;IAChF,aAAa,CAAC,SAAiB;QACpC,OAAO,GAAG,CAAC;IACb,CAAC;IACD,sFAAsF;IAC/E,aAAa,CAAC,QAAgB;QACnC,MAAM,EAAE,GAAG,IAAI,CAAC,EAAE,CAAC;QACnB,MAAM,KAAK,GAAG,QAAQ,GAAG,EAAE,CAAC;QAC5B,OAAO,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,GAAG,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC;IAC7D,CAAC;IACD,yEAAyE;IACzD,0BAA0B,CAAC,QAAgB;QACzD,OAAO,QAAQ,CAAC,YAAY,CAAC,IAAI,CAAC,YAAY,CAAC,QAAQ,CAAC,EAAE,IAAI,CAAC,YAAY,CAAC,QAAQ,CAAC,CAAC,CAAC;IACzF,CAAC;IACD,mDAAmD;IAC5C,WAAW,CAAC,CAAS;QAC1B,MAAM,SAAS,GAAG,CAAC,GAAG,IAAI,CAAC,cAAc,CAAC;QAC1C,MAAM,SAAS,GAAG,YAAY,CAAC,WAAW,CACxC,SAAS,EACT,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,CACpC,CAAC;QACF,IAAI,SAAS,KAAK,SAAS;YACzB,OAAO,SAAS,CAAC;QACnB,OAAO,SAAS,CAAC;IACnB,CAAC;CACF","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 { XYCurveEvaluator } from \"./XYCurveEvaluator\";\nimport { SimpleNewton } from \"../../numerics/Newton\";\n\n/**\n * @internal\n */\nexport class DirectHalfCosineSpiralEvaluator extends XYCurveEvaluator {\n  public nominalLength1: number;\n  public nominalRadius1: number;\n  private _c: number;\n  private _c1: number;\n  private _c2: number;\n  public constructor(length1: number, radius1: number) {\n    super();\n    this.nominalLength1 = length1;\n    this.nominalRadius1 = radius1;\n    const pi = Math.PI;\n    this._c1 = 1.0 / (2.0 * pi * pi);\n    this._c2 = 0.25;\n    this._c = 0.0; // TO BE UPDATED BELOW\n    this.updateConstants();\n  }\n  private updateConstants() {\n    this._c = this.nominalLength1 * this.nominalLength1 / this.nominalRadius1;\n  }\n  public scaleInPlace(scaleFactor: number) {\n    this.nominalLength1 *= scaleFactor;\n    this.nominalRadius1 *= scaleFactor;\n    this.updateConstants();\n  }\n  /** Return a deep copy of the evaluator. */\n  public clone(): DirectHalfCosineSpiralEvaluator {\n    return new DirectHalfCosineSpiralEvaluator(this.nominalLength1, this.nominalRadius1);\n  }\n  /** Member by member matchup. */\n  public isAlmostEqual(other: any): boolean {\n    if (other instanceof DirectHalfCosineSpiralEvaluator) {\n      return Geometry.isSameCoordinate(this.nominalLength1, other.nominalLength1)\n        && Geometry.isSameCoordinate(this.nominalRadius1, other.nominalRadius1);\n      // remark: c,c1,c2 are computed, need not be tested.\n    }\n    return false;\n  }\n  /** Evaluate X at fractional position. */\n  public fractionToX(fraction: number): number {\n    return fraction * this.nominalLength1;\n  }\n  /** Evaluate Y at fractional position. */\n  public fractionToY(fraction: number): number {\n    const theta = fraction * Math.PI;\n    return this._c * (this._c2 * fraction * fraction - this._c1 * (1.0 - Math.cos(theta)));\n  }\n  /** Evaluate derivative of X with respect to fraction at fractional position. */\n  public fractionToDX(_fraction: number): number {\n    return this.nominalLength1;\n  }\n  /** Evaluate derivative of Y with respect to fraction at fractional position. */\n  public fractionToDY(fraction: number): number {\n    const pi = Math.PI;\n    const theta = fraction * pi;\n    return this._c * (2.0 * this._c2 * fraction - this._c1 * pi * Math.sin(theta));\n  }\n  /** Evaluate second derivative of X with respect to fraction at fractional position. */\n  public fractionToDDX(_fraction: number): number {\n    return 0.0;\n  }\n  /** Evaluate third derivative of Y with respect to fraction at fractional position. */\n  public fractionToDDY(fraction: number): number {\n    const pi = Math.PI;\n    const theta = fraction * pi;\n    return this._c * (2.0 * this._c2 - this._c1 * pi * pi * Math.cos(theta));\n  }\n  /** Evaluate second derivative of X with respect to fraction at fractional position. */\n  public fractionToD3X(_fraction: number): number {\n    return 0.0;\n  }\n  /** Evaluate third derivative of Y with respect to fraction at fractional position. */\n  public fractionToD3Y(fraction: number): number {\n    const pi = Math.PI;\n    const theta = fraction * pi;\n    return this._c * this._c1 * pi * pi * pi * Math.sin(theta);\n  }\n  /** Return the magnitude of the first vector at fractional coordinate. */\n  public override fractionToTangentMagnitude(fraction: number): number {\n    return Geometry.hypotenuseXY(this.fractionToDX(fraction), this.fractionToDY(fraction));\n  }\n  /** Invert the fractionToX function for given X. */\n  public xToFraction(x: number): number | undefined {\n    const fraction0 = x / this.nominalLength1;\n    const fraction1 = SimpleNewton.runNewton1D(\n      fraction0,\n      (f: number) => (this.fractionToX(f) - x),\n      (f: number) => this.fractionToDX(f),\n    );\n    if (fraction1 === undefined)\n      return undefined;\n    return fraction1;\n  }\n}\n"]}