@itwin/core-geometry 3.0.0-dev.76 → 3.0.0-dev.80

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

@@ -0,0 +1,220 @@
+/*---------------------------------------------------------------------------------------------
+* Copyright (c) Bentley Systems, Incorporated. All rights reserved.
+* See LICENSE.md in the project root for license terms and full copyright notice.
+*--------------------------------------------------------------------------------------------*/
+/** @packageDocumentation
+ * @module Curve
+ */
+import { Geometry } from "../../Geometry";
+import { XYCurveEvaluator } from "./XYCurveEvaluator";
+import { SimpleNewton } from "../../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.
+ * @internal
+ */
+export class ClothoidSeriesRLEvaluator extends XYCurveEvaluator {
+    constructor(nominalLength1, constantDiv2LR, numXTerms = 4, numYTerms = 4) {
+        super();
+        this.nominalLength1 = nominalLength1;
+        this.constantDiv2LR = constantDiv2LR;
+        this.numXTerms = numXTerms;
+        this.numYTerms = numYTerms;
+    }
+    /** Return a deep clone. */
+    clone() {
+        return new ClothoidSeriesRLEvaluator(this.nominalLength1, this.constantDiv2LR, this.numXTerms, this.numYTerms);
+    }
+    scaleInPlace(scaleFactor) {
+        this.nominalLength1 *= scaleFactor;
+        this.constantDiv2LR /= (scaleFactor * scaleFactor);
+    }
+    /** Member by member matchup ... */
+    isAlmostEqual(other) {
+        if (other instanceof ClothoidSeriesRLEvaluator) {
+            return this.numXTerms === other.numXTerms
+                && this.numYTerms === other.numYTerms
+                && Geometry.isAlmostEqualNumber(this.constantDiv2LR, other.constantDiv2LR)
+                && Geometry.isSameCoordinate(this.nominalLength1, other.nominalLength1);
+        }
+        return false;
+    }
+    /**
+     * 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); }
+    /**
+     * 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); }
+    /**
+     * 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); }
+    /**
+     * 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); }
+    /**
+     * 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); }
+    /**
+     * 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); }
+    /**
+     * 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)
+        // 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 - ....) ) )
+        const s = fraction * this.nominalLength1;
+        let result = s;
+        if (numTerms < 2)
+            return result;
+        const q1 = s * s * this.constantDiv2LR;
+        const beta = -q1 * q1;
+        let alpha = s;
+        let m = 1;
+        let n = 5;
+        for (let i = 1; i < numTerms; i++) {
+            alpha *= beta / (m * (m + 1));
+            result += alpha / n;
+            m += 2;
+            n += 4;
+        }
+        return result;
+    }
+    fractionToYGo(fraction, numTerms) {
+        // 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) - ...)
+        const s = fraction * this.nominalLength1;
+        const q1 = s * s * this.constantDiv2LR;
+        let result = q1 * s / 3;
+        if (numTerms < 2)
+            return result;
+        const beta = -q1 * q1;
+        let alpha = q1 * s;
+        let m = 2;
+        let n = 7;
+        for (let i = 1; i < numTerms; i++) {
+            alpha *= beta / (m * (m + 1));
+            result += alpha / n;
+            m += 2;
+            n += 4;
+        }
+        return result;
+    }
+    fractionToDXGo(fraction, numTerms) {
+        // 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! -
+        // new Term = old Term * beta / (m(m+1))
+        const s = fraction * this.nominalLength1;
+        let result = 1;
+        if (numTerms < 2) {
+            return result * this.nominalLength1;
+        }
+        const q1 = s * s * this.constantDiv2LR;
+        const beta = -q1 * q1;
+        let alpha = 1.0;
+        let m = 1;
+        for (let i = 1; i < numTerms; i++) {
+            alpha *= beta / (m * (m + 1));
+            result += alpha;
+            m += 2;
+        }
+        return result * this.nominalLength1;
+    }
+    fractionToDYGo(fraction, numTerms) {
+        if (numTerms <= 0)
+            return 0;
+        // dY = q - q^3/3!
+        // q = s^2 c
+        // dY = s^2 c - s^6 c^3/3! + s^10 c^5/ 5!
+        // recurrence  advancing m by 2  alpha *= -(s^4 c^2) / (m(m+1))
+        const s = fraction * this.nominalLength1;
+        const q1 = s * s * this.constantDiv2LR;
+        let result = q1;
+        if (numTerms < 2)
+            return result * this.nominalLength1;
+        const beta = -q1 * q1;
+        let alpha = q1;
+        let m = 2;
+        for (let i = 1; i < numTerms; i++) {
+            alpha *= beta / (m * (m + 1));
+            result += alpha;
+            m += 2;
+        }
+        return result * this.nominalLength1;
+    }
+    fractionToDDXGo(fraction, numTerms) {
+        // DX is "cosine"
+        // DDX is "- sine" series times chain rule dTheta/ds = 2 * s * this.constantDivLR
+        const s = fraction * this.nominalLength1;
+        const dTheta = 2 * this.constantDiv2LR * s;
+        const sine = this.fractionToDYGo(fraction, numTerms - 1);
+        const resultA = (-dTheta * sine * this.nominalLength1);
+        return resultA;
+    }
+    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!
+        const s = fraction * this.nominalLength1;
+        const dTheta = 2 * this.constantDiv2LR * s;
+        const cosine = this.fractionToDXGo(fraction, numTerms);
+        return cosine * dTheta * this.nominalLength1;
+    }
+    fractionToD3X(fraction) {
+        if (this.numXTerms <= 1)
+            return 0.0;
+        // DX is "cosine"
+        // DDX is "- sine" series times chain rule dTheta/ds = 2 * s * this.constantDivLR
+        const s = fraction * this.nominalLength1;
+        const dTheta = 2.0 * this.constantDiv2LR * s;
+        const d2Theta = 2.0 * this.constantDiv2LR;
+        const sine = this.fractionToDYGo(fraction, this.numXTerms - 1);
+        const cosine = this.fractionToDXGo(fraction, this.numXTerms - 1);
+        return (-cosine * dTheta * dTheta - sine * d2Theta) * this.nominalLength1 * this.nominalLength1;
+    }
+    fractionToD3Y(fraction) {
+        // DY is "sine"
+        // DDY is "cosine" series times chain rule dTheta/ds = 2 * s * this.constantDivLR
+        const s = fraction * this.nominalLength1;
+        const dTheta = 2.0 * this.constantDiv2LR * s;
+        const d2Theta = 2.0 * this.constantDiv2LR;
+        // dY is sine series with numYTerms.
+        // ddY is cosine series.  Leading term of sine series is non-constant, so numYTerms here also
+        // d3Y is sine series. Derivative of preceding cosine killed first term.
+        const cosine = this.fractionToDXGo(fraction, this.numYTerms);
+        const sine = this.fractionToDYGo(fraction, this.numYTerms - 1);
+        return (-sine * dTheta * dTheta + cosine * d2Theta) * this.nominalLength1 * this.nominalLength1;
+    }
+    xToFraction(x) {
+        const fraction0 = x / this.nominalLength1;
+        const fraction1 = SimpleNewton.runNewton1D(fraction0, (f) => (this.fractionToX(f) - x), (f) => this.fractionToDX(f));
+        if (fraction1 === undefined)
+            return undefined;
+        return fraction1;
+    }
+}
+//# sourceMappingURL=ClothoidSeries.js.map

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

@@ -0,0 +1 @@
+{"version":3,"file":"ClothoidSeries.js","sourceRoot":"","sources":["../../../../src/curve/spiral/ClothoidSeries.ts"],"names":[],"mappings":"AAAA;;;+FAG+F;AAE/F;;GAEG;AAEH,OAAO,EAAE,QAAQ,EAAE,MAAM,gBAAgB,CAAC;AAC1C,OAAO,EAAE,gBAAgB,EAAE,MAAM,oBAAoB,CAAC;AACtD,OAAO,EAAE,YAAY,EAAE,MAAM,uBAAuB,CAAC;AACrD;;;;;;;GAOG;AACH,MAAM,OAAO,yBAA0B,SAAQ,gBAAgB;IAK7D,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;YAC9C,OAAO,IAAI,CAAC,SAAS,KAAK,KAAK,CAAC,SAAS;mBACpC,IAAI,CAAC,SAAS,KAAK,KAAK,CAAC,SAAS;mBAClC,QAAQ,CAAC,mBAAmB,CAAC,IAAI,CAAC,cAAc,EAAE,KAAK,CAAC,cAAc,CAAC;mBACvE,QAAQ,CAAC,gBAAgB,CAAC,IAAI,CAAC,cAAc,EAAE,KAAK,CAAC,cAAc,CAAC,CAAC;SAC3E;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;YACjC,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;SACR;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;YACjC,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;SACR;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;YAChB,OAAO,MAAM,GAAG,IAAI,CAAC,cAAc,CAAC;SACrC;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;YACjC,KAAK,IAAI,IAAI,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;YAC9B,MAAM,IAAI,KAAK,CAAC;YAChB,CAAC,IAAI,CAAC,CAAC;SACR;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;YACjC,KAAK,IAAI,IAAI,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;YAC9B,MAAM,IAAI,KAAK,CAAC;YAChB,CAAC,IAAI,CAAC,CAAC;SACR;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,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;IACnB,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\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"]}

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

@@ -0,0 +1,43 @@
+/** @packageDocumentation
+ * @module Curve
+ */
+import { XYCurveEvaluator } from "./XYCurveEvaluator";
+/** 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.
+ * @internal
+ */
+export declare abstract class CubicEvaluator extends XYCurveEvaluator {
+    protected _cubicM: number;
+    protected _axisLength: number;
+    protected constructor(axisLength: number, cubicM: number);
+    /** Update both constants. */
+    setConstants(axisLength: number, cubicM: number): void;
+    get axisLength(): number;
+    get cubicM(): number;
+    /**
+     * 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 */
+    fractionToDX(_fraction: number): number;
+    /** Evaluate second derivative of X with respect to fraction */
+    fractionToDDX(_fraction: number): number;
+    /** Evaluate third derivative of X with respect to fraction */
+    fractionToD3X(_fraction: number): number;
+    /** Evaluate Y at fraction. */
+    fractionToY(fraction: number): number;
+    /** Evaluate derivative of Y with respect to fraction. */
+    fractionToDY(fraction: number): number;
+    /** Evaluate second derivative of Y with respect to fraction. */
+    fractionToDDY(fraction: number): number;
+    /** Evaluate third derivative of Y with respect to fraction. */
+    fractionToD3Y(_fraction: number): number;
+    /** Evaluate fraction at x. */
+    xToFraction(x: number): number;
+}
+//# sourceMappingURL=CubicEvaluator.d.ts.map

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

@@ -0,0 +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"}

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

@@ -0,0 +1,67 @@
+/*---------------------------------------------------------------------------------------------
+* Copyright (c) Bentley Systems, Incorporated. All rights reserved.
+* See LICENSE.md in the project root for license terms and full copyright notice.
+*--------------------------------------------------------------------------------------------*/
+/** @packageDocumentation
+ * @module Curve
+ */
+import { XYCurveEvaluator } from "./XYCurveEvaluator";
+/** 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.
+ * @internal
+ */
+export class CubicEvaluator extends XYCurveEvaluator {
+    constructor(axisLength, cubicM) {
+        super();
+        this._cubicM = cubicM;
+        this._axisLength = axisLength;
+    }
+    /** Update both constants. */
+    setConstants(axisLength, cubicM) {
+        this._axisLength = axisLength;
+        this._cubicM = 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 . .
+        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; }
+    /** Evaluate Y at fraction. */
+    fractionToY(fraction) {
+        const x = fraction * this._axisLength;
+        return this._cubicM * x * x * x;
+    }
+    /** Evaluate derivative of Y with respect to fraction. */
+    fractionToDY(fraction) {
+        const x = fraction * this._axisLength;
+        return 3.0 * this._cubicM * x * x * this._axisLength;
+    }
+    /** Evaluate second derivative of Y with respect to fraction. */
+    fractionToDDY(fraction) {
+        const x = fraction * this._axisLength;
+        return 6.0 * this._cubicM * x * this._axisLength * this._axisLength;
+    }
+    /** Evaluate third derivative of Y with respect to fraction. */
+    fractionToD3Y(_fraction) {
+        return 6.0 * this._cubicM * this._axisLength * this._axisLength * this._axisLength;
+    }
+    /** Evaluate fraction at x. */
+    xToFraction(x) { return x / this._axisLength; }
+}
+//# sourceMappingURL=CubicEvaluator.js.map

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

@@ -0,0 +1 @@
+{"version":3,"file":"CubicEvaluator.js","sourceRoot":"","sources":["../../../../src/curve/spiral/CubicEvaluator.ts"],"names":[],"mappings":"AAAA;;;+FAG+F;AAE/F;;GAEG;AACH,OAAO,EAAE,gBAAgB,EAAE,MAAM,oBAAoB,CAAC;AACtD;;;;GAIG;AACH,MAAM,OAAgB,cAAe,SAAQ,gBAAgB;IAG3D,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","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"]}

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

@@ -0,0 +1,108 @@
+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)
+ *    * 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. */
+    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 */
+    clone(): CzechSpiralEvaluator;
+    /** 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.
+     */
+    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.
+     *
+     */
+    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.
+       */
+    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.
+     *
+     */
+    static inverseL2R2Map(b: 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)
+ *    * 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 ....
+     */
+    private static computeCubicM;
+    /** 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 */
+    clone(): ItalianSpiralEvaluator;
+    /** 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.
+       */
+    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.
+     *
+     */
+    xToItalianApproximateDistance(d: number): number | undefined;
+}
+//# sourceMappingURL=CzechSpiralEvaluator.d.ts.map

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

@@ -0,0 +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"}