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

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

@@ -13,10 +13,10 @@ 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.
+ * 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 {
@@ -39,7 +39,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 +53,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 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(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,7 +120,7 @@ class ClothoidSeriesRLEvaluator extends XYCurveEvaluator_1.XYCurveEvaluator {
         return result;
     }
     fractionToYGo(fraction, numTerms) {
-        // Write the series for sin (theta)
+        // write the 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) - ...)
@@ -130,7 +142,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 +197,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":["/*---------------------------------------------------------------------------------------------\r\n* Copyright (c) Bentley Systems, Incorporated. All rights reserved.\r\n* See LICENSE.md in the project root for license terms and full copyright notice.\r\n*--------------------------------------------------------------------------------------------*/\r\n\r\n/** @packageDocumentation\r\n * @module Curve\r\n */\r\n\r\nimport { Geometry } from \"../../Geometry\";\r\nimport { XYCurveEvaluator } from \"./XYCurveEvaluator\";\r\nimport { SimpleNewton } from \"../../numerics/Newton\";\r\n/**\r\n * Methods to evaluate caller-specified number of terms of the x and y series for a clothoid.\r\n * Each instance has\r\n * * Number of x and y terms to use.\r\n * * constant for theta=c * x * x\r\n *    * This value is c=1/(2 R L)  for curve length L measured from inflection to point with radius R.\r\n * @internal\r\n */\r\nexport class ClothoidSeriesRLEvaluator extends XYCurveEvaluator {\r\n  public numXTerms: number;\r\n  public numYTerms: number;\r\n  public constantDiv2LR: number;\r\n  public nominalLength1: number;\r\n  public constructor(nominalLength1: number, constantDiv2LR: number, numXTerms: number = 4, numYTerms: number = 4) {\r\n    super();\r\n    this.nominalLength1 = nominalLength1;\r\n    this.constantDiv2LR = constantDiv2LR;\r\n    this.numXTerms = numXTerms;\r\n    this.numYTerms = numYTerms;\r\n  }\r\n  /** Return a deep clone. */\r\n  public clone(): ClothoidSeriesRLEvaluator {\r\n    return new ClothoidSeriesRLEvaluator(this.nominalLength1, this.constantDiv2LR, this.numXTerms, this.numYTerms);\r\n  }\r\n  public scaleInPlace(scaleFactor: number) {\r\n    this.nominalLength1 *= scaleFactor;\r\n    this.constantDiv2LR /= (scaleFactor * scaleFactor);\r\n  }\r\n  /** Member by member matchup ... */\r\n  public isAlmostEqual(other: any): boolean {\r\n    if (other instanceof ClothoidSeriesRLEvaluator) {\r\n      return this.numXTerms === other.numXTerms\r\n        && this.numYTerms === other.numYTerms\r\n        && Geometry.isAlmostEqualNumber(this.constantDiv2LR, other.constantDiv2LR)\r\n        && Geometry.isSameCoordinate(this.nominalLength1, other.nominalLength1);\r\n    }\r\n    return false;\r\n  }\r\n  /**\r\n   * Evaluate the X series at a nominal distance along the curve.\r\n   * @param fraction fractional position along the curve.\r\n   */\r\n  public fractionToX(fraction: number): number { return this.fractionToXGo(fraction, this.numXTerms); }\r\n  /**\r\n   * Evaluate the Y series at a nominal distance along the curve.\r\n   * @param fraction fractional position along the curve.\r\n   */\r\n  public fractionToY(fraction: number): number { return this.fractionToYGo(fraction, this.numYTerms); }\r\n\r\n  /**\r\n   * Evaluate the derivative of the X series at a nominal distance along the curve.\r\n   * @param fraction fractional position along the curve.\r\n   */\r\n  public fractionToDX(fraction: number): number { return this.fractionToDXGo(fraction, this.numXTerms); }\r\n  /**\r\n   * Evaluate the derivative of the Y series at a nominal distance along the curve.\r\n   * @param fraction fractional position along the curve.\r\n   */\r\n  public fractionToDY(fraction: number): number { return this.fractionToDYGo(fraction, this.numYTerms); }\r\n\r\n  /**\r\n   * Evaluate the derivative of the X series at a nominal distance along the curve.\r\n   * @param fraction fractional position along the curve.\r\n   */\r\n  public fractionToDDX(fraction: number): number { return this.fractionToDDXGo(fraction, this.numXTerms); }\r\n  /**\r\n   * Evaluate the derivative of the Y series at a nominal distance along the curve.\r\n   * @param fraction fractional position along the curve.\r\n   */\r\n  public fractionToDDY(fraction: number): number { return this.fractionToDDYGo(fraction, this.numYTerms); }\r\n\r\n  /**\r\n   * Evaluate the X series at a nominal distance along the curve.\r\n   * @param fraction fractional position along the curve.\r\n   * @param numTerms number of terms to use.\r\n   */\r\n  public fractionToXGo(fraction: number, numTerms: number): number {\r\n    // Write the series for cos (theta)\r\n    // replace theta by s*s*c\r\n    // integrate wrt s\r\n    //  x = s - s^5 c^4/ 2 + s^9 c^8/(4!) - s^13 c^12 / 6!\r\n    //  x = s(1 - (s^4 c^2/2) ( 1/5 -s^4 c^2 / (3*4)  ( 1/9 - ....) ) )\r\n    const s = fraction * this.nominalLength1;\r\n    let result = s;\r\n    if (numTerms < 2)\r\n      return result;\r\n    const q1 = s * s * this.constantDiv2LR;\r\n    const beta = - q1 * q1;\r\n    let alpha = s;\r\n    let m = 1;\r\n    let n = 5;\r\n    for (let i = 1; i < numTerms; i++) {\r\n      alpha *= beta / (m * (m + 1));\r\n      result += alpha / n;\r\n      m += 2;\r\n      n += 4;\r\n    }\r\n    return result;\r\n  }\r\n  public fractionToYGo(fraction: number, numTerms: number): number {\r\n    // Write the series for sin (theta)\r\n    // replace theta by s*s*c\r\n    // integrate wrt s\r\n    //  x = s^3 c^2/ 3( (1/3)) - s^7 c^6/(3!) ((1/7)) - s^11 c^10 / 5! ((1/9) - ...)\r\n    const s = fraction * this.nominalLength1;\r\n    const q1 = s * s * this.constantDiv2LR;\r\n    let result = q1 * s / 3;\r\n    if (numTerms < 2)\r\n      return result;\r\n    const beta = - q1 * q1;\r\n    let alpha = q1 * s;\r\n    let m = 2;\r\n    let n = 7;\r\n    for (let i = 1; i < numTerms; i++) {\r\n      alpha *= beta / (m * (m + 1));\r\n      result += alpha / n;\r\n      m += 2;\r\n      n += 4;\r\n    }\r\n    return result;\r\n  }\r\n  public fractionToDXGo(fraction: number, numTerms: number): number {\r\n    // Yes -- this does happen during derivatives of cosines with more than 0 terms !!\r\n    if (numTerms <= 0)\r\n      return 0;\r\n    // dX = 1 - s^4c^2/2 + s^8 c^4 / 4! -\r\n    // new Term = old Term * beta / (m(m+1))\r\n    const s = fraction * this.nominalLength1;\r\n    let result = 1;\r\n    if (numTerms < 2) {\r\n      return result * this.nominalLength1;\r\n    }\r\n    const q1 = s * s * this.constantDiv2LR;\r\n    const beta = - q1 * q1;\r\n    let alpha = 1.0;\r\n    let m = 1;\r\n    for (let i = 1; i < numTerms; i++) {\r\n      alpha *= beta / (m * (m + 1));\r\n      result += alpha;\r\n      m += 2;\r\n    }\r\n    return result * this.nominalLength1;\r\n  }\r\n  public fractionToDYGo(fraction: number, numTerms: number): number {\r\n    if (numTerms <= 0)\r\n      return 0;\r\n    // dY = q - q^3/3!\r\n    // q = s^2 c\r\n    // dY = s^2 c - s^6 c^3/3! + s^10 c^5/ 5!\r\n    // recurrence  advancing m by 2  alpha *= -(s^4 c^2) / (m(m+1))\r\n    const s = fraction * this.nominalLength1;\r\n    const q1 = s * s * this.constantDiv2LR;\r\n    let result = q1;\r\n    if (numTerms < 2)\r\n      return result * this.nominalLength1;\r\n    const beta = - q1 * q1;\r\n    let alpha = q1;\r\n    let m = 2;\r\n    for (let i = 1; i < numTerms; i++) {\r\n      alpha *= beta / (m * (m + 1));\r\n      result += alpha;\r\n      m += 2;\r\n    }\r\n    return result * this.nominalLength1;\r\n  }\r\n\r\n  public fractionToDDXGo(fraction: number, numTerms: number): number {\r\n    // DX is \"cosine\"\r\n    // DDX is \"- sine\" series times chain rule dTheta/ds = 2 * s * this.constantDivLR\r\n    const s = fraction * this.nominalLength1;\r\n\r\n    const dTheta = 2 * this.constantDiv2LR * s;\r\n    const sine = this.fractionToDYGo(fraction, numTerms - 1);\r\n    const resultA = (- dTheta * sine * this.nominalLength1);\r\n    return resultA;\r\n  }\r\n  public fractionToDDYGo(fraction: number, numTerms: number): number {\r\n    // DY is \"sine\"\r\n    // DDY is \"cosine\" series times chain rule dTheta/ds = 2 * s * this.constantDivLR\r\n    // BUT  .... derivative of the cosine series leading term is zero ... use one less term!\r\n    const s = fraction * this.nominalLength1;\r\n    const dTheta = 2 * this.constantDiv2LR * s;\r\n    const cosine = this.fractionToDXGo(fraction, numTerms);\r\n    return cosine * dTheta * this.nominalLength1;\r\n  }\r\n\r\n  public fractionToD3X(fraction: number): number {\r\n    if (this.numXTerms <= 1)\r\n      return 0.0;\r\n    // DX is \"cosine\"\r\n    // DDX is \"- sine\" series times chain rule dTheta/ds = 2 * s * this.constantDivLR\r\n    const s = fraction * this.nominalLength1;\r\n    const dTheta = 2.0 * this.constantDiv2LR * s;\r\n    const d2Theta = 2.0 * this.constantDiv2LR;\r\n    const sine = this.fractionToDYGo(fraction, this.numXTerms - 1);\r\n    const cosine = this.fractionToDXGo(fraction, this.numXTerms - 1);\r\n    return (- cosine * dTheta * dTheta - sine * d2Theta) * this.nominalLength1 * this.nominalLength1;\r\n  }\r\n  public fractionToD3Y(fraction: number): number {\r\n    // DY is \"sine\"\r\n    // DDY is \"cosine\" series times chain rule dTheta/ds = 2 * s * this.constantDivLR\r\n    const s = fraction * this.nominalLength1;\r\n    const dTheta = 2.0 * this.constantDiv2LR * s;\r\n    const d2Theta = 2.0 * this.constantDiv2LR;\r\n    // dY is sine series with numYTerms.\r\n    // ddY is cosine series.  Leading term of sine series is non-constant, so numYTerms here also\r\n    // d3Y is sine series. Derivative of preceding cosine killed first term.\r\n    const cosine = this.fractionToDXGo(fraction, this.numYTerms);\r\n    const sine = this.fractionToDYGo(fraction, this.numYTerms - 1);\r\n    return (-sine * dTheta * dTheta + cosine * d2Theta) * this.nominalLength1 * this.nominalLength1;\r\n  }\r\n\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"]}
+{"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;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,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,kCAAkC;QAClC,yBAAyB;QACzB,kBAAkB;QAClB,sDAAsD;QACtD,gEAAgE;QAChE,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,kCAAkC;QAClC,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,+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;AA1ND,8DA0NC","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 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): 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 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  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":["/*---------------------------------------------------------------------------------------------\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 { XYCurveEvaluator } from \"./XYCurveEvaluator\";\r\n/** Intermediate class for evaluation of bare cubic spirals `y = m ^x^3` with x axis from [0..xLength]\r\n * * This implements all the computations among fraction, x, and y.\r\n * * Derived classes implement specialized logic such as (a) precomputing `m` and (b) domain-specific fraction-to-distance approximations.\r\n * @internal\r\n */\r\nexport abstract class CubicEvaluator extends XYCurveEvaluator {\r\n  protected _cubicM: number;\r\n  protected _axisLength: number;\r\n  protected constructor(axisLength: number, cubicM: number) {\r\n    super();\r\n    this._cubicM = cubicM;\r\n    this._axisLength = axisLength;\r\n  }\r\n  /** Update both constants. */\r\n  public setConstants(axisLength: number, cubicM: number) {\r\n    this._axisLength = axisLength;\r\n    this._cubicM = cubicM;\r\n  }\r\n  public get axisLength() { return this._axisLength; }\r\n  public get cubicM() { return this._cubicM; }\r\n  /**\r\n   * Apply `scaleFactor` to the xLength and cubicM.\r\n   * * Derived classes commonly call this as `super.scaleInPlace()`, and additionally apply the scale to their members.\r\n   * @param scaleFactor\r\n   */\r\n  public scaleInPlace(scaleFactor: number) {\r\n    this._axisLength *= scaleFactor;\r\n    // \"x\" arriving at \"m * x^3\" will be scaled. \"m\" has to be divided by the scale to cancel 2 of the 3 . .\r\n    this._cubicM /= (scaleFactor * scaleFactor);\r\n  }\r\n  /** Evaluate X at fraction. */\r\n  public fractionToX(fraction: number): number { return fraction * this._axisLength; }\r\n  /** Evaluate derivative of X with respect to fraction */\r\n  public fractionToDX(_fraction: number): number { return this._axisLength; }\r\n  /** Evaluate second derivative of X with respect to fraction */\r\n  public fractionToDDX(_fraction: number): number { return 0.0; }\r\n  /** Evaluate third derivative of X with respect to fraction */\r\n  public fractionToD3X(_fraction: number): number { return 0.0; }\r\n  /** Evaluate Y at fraction. */\r\n  public fractionToY(fraction: number): number {\r\n    const x = fraction * this._axisLength;\r\n    return this._cubicM * x * x * x;\r\n  }\r\n  /** Evaluate derivative of Y with respect to fraction. */\r\n  public fractionToDY(fraction: number): number {\r\n    const x = fraction * this._axisLength;\r\n    return 3.0 * this._cubicM * x * x * this._axisLength;\r\n  }\r\n  /** Evaluate second derivative of Y with respect to fraction. */\r\n  public fractionToDDY(fraction: number): number {\r\n    const x = fraction * this._axisLength;\r\n    return 6.0 * this._cubicM * x * this._axisLength * this._axisLength;\r\n  }\r\n  /** Evaluate third derivative of Y with respect to fraction. */\r\n  public fractionToD3Y(_fraction: number): number {\r\n    return 6.0 * this._cubicM * this._axisLength * this._axisLength * this._axisLength;\r\n  }\r\n  /** Evaluate fraction at x. */\r\n  public xToFraction(x: number): number { return x / this._axisLength; }\r\n}\r\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"]}

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

@@ -2,106 +2,96 @@ import { CubicEvaluator } from "./CubicEvaluator";
 /**
  * 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 declare class CzechSpiralEvaluator extends CubicEvaluator {
     nominalLength1: number;
     nominalRadius1: number;
-    /** Constructor is private.  Caller responsible for cubicM validity. */
+    /** Constructor is private. Caller responsible for cubicM validity. */
     private constructor();
     /**
      * 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: number, radius1: number): number | undefined;
     /** Compute the czech cubic constant. */
     static computeCubicM(length1: number, radius1: number): number | undefined;
     static create(length1: number, radius1: number): CzechSpiralEvaluator | undefined;
     scaleInPlace(scaleFactor: number): void;
-    /** return a deep copy of the evaluator */
+    /** Return a deep copy of the evaluator. */
     clone(): CzechSpiralEvaluator;
-    /** Member by member matchup ... */
+    /** Member by member matchup. */
     isAlmostEqual(other: any): boolean;
     /**
      * 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: number): number;
     /**
      * 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: number): number | undefined;
     /**
-       * 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: number, sign: number, length: number, radius: number): number;
     /**
      * 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: number, sign: number, length: number, radius: number): number | undefined;
+    static inverseL2R2Map(d: number, sign: number, length: number, radius: number): number | undefined;
 }
 /**
  * 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 declare class ItalianSpiralEvaluator extends CubicEvaluator {
     nominalLength1: number;
     nominalRadius1: number;
-    /** Compute the czech cubic constant.
-     * ** funky mixture of lengths ....
+    /**
+     * Compute the czech cubic constant.
+     * * Funky mixture of lengths.
      */
     private static computeCubicM;
-    /** Constructor is private.  Caller responsible for cubicM validity. */
+    /** Constructor is private. Caller responsible for cubicM validity. */
     private constructor();
     static create(length1: number, radius1: number): ItalianSpiralEvaluator | undefined;
     scaleInPlace(scaleFactor: number): void;
-    /** return a deep copy of the evaluator */
+    /** Return a deep copy of the evaluator. */
     clone(): ItalianSpiralEvaluator;
-    /** Member by member matchup ... */
+    /** Member by member matchup. */
     isAlmostEqual(other: any): boolean;
     /**
-       * 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: number): number;
     /**
      * 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: number): number | undefined;
 }

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

@@ -1 +1 @@
-{"version":3,"file":"CzechSpiralEvaluator.d.ts","sourceRoot":"","sources":["../../../../src/curve/spiral/CzechSpiralEvaluator.ts"],"names":[],"mappings":"AASA,OAAO,EAAE,cAAc,EAAE,MAAM,kBAAkB,CAAC;AAElD;;;;;;;;;;;;;GAaG;AACH,qBAAa,oBAAqB,SAAQ,cAAc;IAC/C,cAAc,EAAE,MAAM,CAAC;IACvB,cAAc,EAAE,MAAM,CAAC;IAE9B,uEAAuE;IACvE,OAAO;IAKP;;;;;;OAMG;WACW,aAAa,CAAC,OAAO,EAAE,MAAM,EAAE,OAAO,EAAE,MAAM,GAAG,MAAM,GAAG,SAAS;IAGjF,wCAAwC;WAC1B,aAAa,CAAC,OAAO,EAAE,MAAM,EAAE,OAAO,EAAE,MAAM,GAAG,MAAM,GAAG,SAAS;WASnE,MAAM,CAAC,OAAO,EAAE,MAAM,EAAE,OAAO,EAAE,MAAM,GAAG,oBAAoB,GAAG,SAAS;IAOxE,YAAY,CAAC,WAAW,EAAE,MAAM;IAKhD,0CAA0C;IACnC,KAAK,IAAI,oBAAoB;IACpC,mCAAmC;IAC5B,aAAa,CAAC,KAAK,EAAE,GAAG,GAAG,OAAO;IAQzC;;;;;OAKG;IACI,2BAA2B,CAAC,CAAC,EAAE,MAAM,GAAG,MAAM;IAGrD;;;;;OAKG;IACI,2BAA2B,CAAC,CAAC,EAAE,MAAM,GAAG,MAAM,GAAG,SAAS;IAIjE;;;SAGK;WACS,cAAc,CAAC,CAAC,EAAE,MAAM,EAAE,IAAI,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,GAAG,MAAM;IAO7F;;;;;OAKG;WACW,cAAc,CAAC,CAAC,EAAE,MAAM,EAAE,IAAI,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,GAAG,MAAM,GAAG,SAAS;CAe1G;AACD;;;;;;;;;;;;;GAaG;AACH,qBAAa,sBAAuB,SAAQ,cAAc;IACjD,cAAc,EAAE,MAAM,CAAC;IACvB,cAAc,EAAE,MAAM,CAAC;IAE9B;;OAEG;IACH,OAAO,CAAC,MAAM,CAAC,aAAa;IAS5B,uEAAuE;IACvE,OAAO;WAKO,MAAM,CAAC,OAAO,EAAE,MAAM,EAAE,OAAO,EAAE,MAAM,GAAG,sBAAsB,GAAG,SAAS;IAY1E,YAAY,CAAC,WAAW,EAAE,MAAM;IAKhD,0CAA0C;IACnC,KAAK,IAAI,sBAAsB;IACtC,mCAAmC;IAC5B,aAAa,CAAC,KAAK,EAAE,GAAG,GAAG,OAAO;IAOzC;;;;;SAKK;IACE,6BAA6B,CAAC,CAAC,EAAE,MAAM,GAAG,MAAM;IAGvD;;;;;OAKG;IACI,6BAA6B,CAAC,CAAC,EAAE,MAAM,GAAG,MAAM,GAAG,SAAS;CAGpE"}
+{"version":3,"file":"CzechSpiralEvaluator.d.ts","sourceRoot":"","sources":["../../../../src/curve/spiral/CzechSpiralEvaluator.ts"],"names":[],"mappings":"AASA,OAAO,EAAE,cAAc,EAAE,MAAM,kBAAkB,CAAC;AAGlD;;;;;;;;;GASG;AACH,qBAAa,oBAAqB,SAAQ,cAAc;IAC/C,cAAc,EAAE,MAAM,CAAC;IACvB,cAAc,EAAE,MAAM,CAAC;IAC9B,sEAAsE;IACtE,OAAO;IAKP;;;;OAIG;WACW,aAAa,CAAC,OAAO,EAAE,MAAM,EAAE,OAAO,EAAE,MAAM,GAAG,MAAM,GAAG,SAAS;IAGjF,wCAAwC;WAC1B,aAAa,CAAC,OAAO,EAAE,MAAM,EAAE,OAAO,EAAE,MAAM,GAAG,MAAM,GAAG,SAAS;WAQnE,MAAM,CAAC,OAAO,EAAE,MAAM,EAAE,OAAO,EAAE,MAAM,GAAG,oBAAoB,GAAG,SAAS;IAMxE,YAAY,CAAC,WAAW,EAAE,MAAM;IAKhD,2CAA2C;IACpC,KAAK,IAAI,oBAAoB;IAGpC,gCAAgC;IACzB,aAAa,CAAC,KAAK,EAAE,GAAG,GAAG,OAAO;IAOzC;;;;;;OAMG;IACI,2BAA2B,CAAC,CAAC,EAAE,MAAM,GAAG,MAAM;IAGrD;;;;OAIG;IACI,2BAA2B,CAAC,CAAC,EAAE,MAAM,GAAG,MAAM,GAAG,SAAS;IAGjE;;;OAGG;WACW,cAAc,CAAC,CAAC,EAAE,MAAM,EAAE,IAAI,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,GAAG,MAAM;IAO7F;;;;OAIG;WACW,cAAc,CAAC,CAAC,EAAE,MAAM,EAAE,IAAI,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,GAAG,MAAM,GAAG,SAAS;CAiB1G;AAED;;;;;;;;;GASG;AACH,qBAAa,sBAAuB,SAAQ,cAAc;IACjD,cAAc,EAAE,MAAM,CAAC;IACvB,cAAc,EAAE,MAAM,CAAC;IAE9B;;;OAGG;IACH,OAAO,CAAC,MAAM,CAAC,aAAa;IAQ5B,sEAAsE;IACtE,OAAO;WAKO,MAAM,CAAC,OAAO,EAAE,MAAM,EAAE,OAAO,EAAE,MAAM,GAAG,sBAAsB,GAAG,SAAS;IAW1E,YAAY,CAAC,WAAW,EAAE,MAAM;IAKhD,2CAA2C;IACpC,KAAK,IAAI,sBAAsB;IAGtC,gCAAgC;IACzB,aAAa,CAAC,KAAK,EAAE,GAAG,GAAG,OAAO;IAOzC;;;;;;OAMG;IACI,6BAA6B,CAAC,CAAC,EAAE,MAAM,GAAG,MAAM;IAGvD;;;;OAIG;IACI,6BAA6B,CAAC,CAAC,EAAE,MAAM,GAAG,MAAM,GAAG,SAAS;CAGpE"}