@itwin/core-geometry 5.3.2 → 5.4.0-dev.11

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

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

@@ -20,7 +20,7 @@ export declare class CzechSpiralEvaluator extends CubicEvaluator {
      *   (0.25==>0.99215), (0.15==>0.997184), (0.10==>0.998749), (0.05==>999687)
      */
     static gammaConstant(length1: number, radius1: number): number | undefined;
-    /** Compute the czech cubic constant. */
+    /** 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;
@@ -67,10 +67,7 @@ export declare class CzechSpiralEvaluator extends CubicEvaluator {
 export declare class ItalianSpiralEvaluator extends CubicEvaluator {
     nominalLength1: number;
     nominalRadius1: number;
-    /**
-     * Compute the czech cubic constant.
-     * * Funky mixture of lengths.
-     */
+    /** Compute the Italian cubic constant. */
     private static computeCubicM;
     /** Constructor is private. Caller responsible for cubicM validity. */
     private constructor();

package/lib/esm/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;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"}
+{"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;IAC9B,0CAA0C;IAC1C,OAAO,CAAC,MAAM,CAAC,aAAa;IAG5B,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"}

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

@@ -35,7 +35,7 @@ export class CzechSpiralEvaluator extends CubicEvaluator {
     static gammaConstant(length1, radius1) {
         return 2.0 * radius1 / Math.sqrt(4.0 * radius1 * radius1 - length1 * length1);
     }
-    /** Compute the czech cubic constant. */
+    /** Compute the Czech cubic constant. */
     static computeCubicM(length1, radius1) {
         const gamma = CzechSpiralEvaluator.gammaConstant(length1, radius1);
         // in the private update method, the LR values should have been vetted
@@ -128,17 +128,9 @@ export class CzechSpiralEvaluator extends CubicEvaluator {
 export class ItalianSpiralEvaluator extends CubicEvaluator {
     nominalLength1;
     nominalRadius1;
-    /**
-     * Compute the czech cubic constant.
-     * * Funky mixture of lengths.
-     */
+    /** Compute the Italian cubic constant. */
     static computeCubicM(lengthXByForward, radius1) {
-        const gamma = CzechSpiralEvaluator.gammaConstant(lengthXByForward, radius1);
-        // in the private update method, the LR values should have been vetted
-        if (gamma === undefined)
-            return undefined;
-        // if radius is negative, it shows up in gamma; but the a signed denominator undoes it so take abs of denominator
-        return gamma / Math.abs((6.0 * radius1 * lengthXByForward));
+        return CzechSpiralEvaluator.computeCubicM(lengthXByForward, radius1);
     }
     /** Constructor is private. Caller responsible for cubicM validity. */
     constructor(length1, radius1, lengthX, cubicM) {
@@ -147,7 +139,7 @@ export class ItalianSpiralEvaluator extends CubicEvaluator {
         this.nominalRadius1 = radius1;
     }
     static create(length1, radius1) {
-        // this seems goofy.; lengthX from forward, then invert for another but that's what the native code does too
+        // this seems goofy; lengthX from forward, then invert for another but that's what the native code does too
         const lengthX = CzechSpiralEvaluator.forwardL2R2Map(length1, -1.0, length1, radius1);
         const lengthX1 = CzechSpiralEvaluator.inverseL2R2Map(length1, 1.0, lengthX, radius1);
         if (lengthX1 === undefined)

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

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

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

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

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

@@ -15,16 +15,17 @@ import { XYCurveEvaluator } from "./XYCurveEvaluator";
 * DirectSpiral3d acts like a TransitionSpiral3d for serialization purposes, but implements spiral types that have
 * "direct" xy calculations without the integrations required for IntegratedSpiral3d.
 * * Each DirectSpiral3d carries an XYCurveEvaluator to give it specialized behavior.
-* * Direct spirals that flow through serialization to native imodel02 are created with these static methods:
+* * Direct spirals are created with these static methods:
 *   * createArema
 *   * createJapaneseCubic
-*   * createAustralianRail
-*   * createDirectHalfCosine
 *   * createChineseCubic
+*   * createWesternAustralian
+*   * createDirectHalfCosine
+*   * createAustralianRail
 *   * createCzechCubic
 *   * createPolishCubic
+*   * createMXCubicAlongArc
 *   * createItalian
-*   * createWesternAustralian
 * @public
 */
 export declare class DirectSpiral3d extends TransitionSpiral3d {
@@ -110,11 +111,11 @@ export declare class DirectSpiral3d extends TransitionSpiral3d {
     /**
      * Create an MX Cubic whose nominal length is close to along the curve.
      * This is y = m*x^3 with
-     * * m is 1/(6RL1).
+     * * m is 1/(6RL).
      *    * 1/(6RL) is the leading term of the sine series.
      * * L1 is an along-the-x-axis distance that is slightly LESS THAN the nominal length.
      * * x is axis position that is slightly LESS than nominal distance along.
-     * * L1, x use the approximation `x = s * ( 1 - s^4/ (40 R R L L))
+     * * L1, x use the approximation `x = s * ( 1 - s^4/ (40 RR LL))
      * @param localToWorld
      * @param nominalL1
      * @param nominalR1

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

@@ -1 +1 @@
-{"version":3,"file":"DirectSpiral3d.d.ts","sourceRoot":"","sources":["../../../../src/curve/spiral/DirectSpiral3d.ts"],"names":[],"mappings":"AAQA,OAAO,EAAE,KAAK,EAAE,MAAM,wBAAwB,CAAC;AAC/C,OAAO,EAAE,eAAe,EAAE,cAAc,EAAE,MAAM,kCAAkC,CAAC;AACnF,OAAO,EAAE,4BAA4B,EAAE,MAAM,+CAA+C,CAAC;AAC7F,OAAO,EAAE,yBAAyB,EAAE,MAAM,4CAA4C,CAAC;AACvF,OAAO,EAAE,OAAO,EAAE,MAAM,kCAAkC,CAAC;AAC3D,OAAO,EAAE,KAAK,EAAE,MAAM,wBAAwB,CAAC;AAC/C,OAAO,EAAE,SAAS,EAAE,MAAM,4BAA4B,CAAC;AACvD,OAAO,EAAE,SAAS,EAAE,MAAM,4BAA4B,CAAC;AAEvD,OAAO,EAAE,YAAY,EAAE,MAAM,iBAAiB,CAAC;AAC/C,OAAO,EAAE,aAAa,EAAE,MAAM,kBAAkB,CAAC;AAOjD,OAAO,EAAE,+BAA+B,EAAE,MAAM,mCAAmC,CAAC;AACpF,OAAO,EAAE,kBAAkB,EAAE,MAAM,sBAAsB,CAAC;AAC1D,OAAO,EAAE,gBAAgB,EAAE,MAAM,oBAAoB,CAAC;AAEtD;;;;;;;;;;;;;;;EAeE;AACF,qBAAa,cAAe,SAAQ,kBAAkB;IACpD,yCAAyC;IACzC,SAAgB,kBAAkB,sBAAsB;IACxD,wFAAwF;IACxF,OAAO,CAAC,cAAc,CAAe;IACrC;;;;;OAKG;IACH,OAAO,CAAC,cAAc,CAAC,CAAe;IACtC,yEAAyE;IACzE,IAAW,aAAa,IAAI,YAAY,CAEvC;IACD,OAAO,CAAC,UAAU,CAAS;IAC3B,OAAO,CAAC,UAAU,CAAS;IAC3B,OAAO,CAAC,UAAU,CAAmB;IACrC,qCAAqC;IACrC,IAAW,SAAS,IAAI,MAAM,CAE7B;IACD,+DAA+D;IAC/D,IAAW,SAAS,IAAI,MAAM,CAE7B;IACD,uCAAuC;IACvC,IAAW,iBAAiB,IAAI,MAAM,CAErC;IACD;;;OAGG;IACH,IAAW,SAAS,IAAI,gBAAgB,CAEvC;gBAIC,YAAY,EAAE,SAAS,EACvB,UAAU,EAAE,MAAM,GAAG,SAAS,EAC9B,kBAAkB,EAAE,+BAA+B,GAAG,SAAS,EAC/D,SAAS,EAAE,MAAM,EACjB,SAAS,EAAE,MAAM,EACjB,sBAAsB,EAAE,SAAS,GAAG,SAAS,EAC7C,SAAS,EAAE,gBAAgB;IAW7B;;;;;;;OAOG;IACH,OAAO,CAAC,cAAc;IA6BtB,yBAAyB;IACT,yBAAyB;IAWzC;;;;;;;;;;OAUG;WACW,uBAAuB,CACnC,UAAU,EAAE,MAAM,EAClB,YAAY,EAAE,SAAS,EACvB,QAAQ,EAAE,MAAM,EAChB,QAAQ,EAAE,MAAM,EAChB,kBAAkB,EAAE,+BAA+B,GAAG,SAAS,EAC/D,SAAS,EAAE,MAAM,EACjB,SAAS,EAAE,MAAM,EACjB,cAAc,EAAE,SAAS,GAAG,SAAS,GACpC,cAAc,GAAG,SAAS;IAgB7B;;;;;;;;;;;;OAYG;WACW,gBAAgB,CAC5B,YAAY,EAAE,SAAS,EAAE,SAAS,EAAE,MAAM,EAAE,SAAS,EAAE,MAAM,EAAE,cAAc,CAAC,EAAE,SAAS,GACxF,cAAc,GAAG,SAAS;IAc7B;;;;;;;;;;;;OAYG;WACW,aAAa,CACzB,YAAY,EAAE,SAAS,EAAE,SAAS,EAAE,MAAM,EAAE,SAAS,EAAE,MAAM,EAAE,cAAc,CAAC,EAAE,SAAS,GACxF,cAAc,GAAG,SAAS;IAc7B;;;;;;;;;;;;OAYG;WACW,qBAAqB,CACjC,YAAY,EAAE,SAAS,EAAE,SAAS,EAAE,MAAM,EAAE,SAAS,EAAE,MAAM,EAAE,cAAc,CAAC,EAAE,SAAS,GACxF,cAAc,GAAG,SAAS;IAc7B;;;;;;;;;OASG;WACW,iBAAiB,CAC7B,YAAY,EAAE,SAAS,EAAE,SAAS,EAAE,MAAM,EAAE,SAAS,EAAE,MAAM,EAAE,cAAc,CAAC,EAAE,SAAS,GACxF,cAAc,GAAG,SAAS;IAc7B;;;;;;;;;;;;OAYG;WACW,oBAAoB,CAChC,YAAY,EAAE,SAAS,EAAE,SAAS,EAAE,MAAM,EAAE,SAAS,EAAE,MAAM,EAAE,cAAc,CAAC,EAAE,SAAS,GACxF,cAAc,GAAG,SAAS;WAcf,sBAAsB,CAClC,YAAY,EAAE,SAAS,EAAE,SAAS,EAAE,MAAM,EAAE,SAAS,EAAE,MAAM,EAAE,cAAc,CAAC,EAAE,SAAS,GACxF,cAAc,GAAG,SAAS;IAW7B;;;;;;;;;OASG;WACW,WAAW,CACvB,YAAY,EAAE,SAAS,EAAE,SAAS,EAAE,MAAM,EAAE,SAAS,EAAE,MAAM,EAAE,cAAc,CAAC,EAAE,SAAS,GACxF,cAAc,GAAG,SAAS;IAG7B;;;;;;;;;OASG;WACW,mBAAmB,CAC/B,YAAY,EAAE,SAAS,EAAE,SAAS,EAAE,MAAM,EAAE,SAAS,EAAE,MAAM,EAAE,cAAc,CAAC,EAAE,SAAS,GACxF,cAAc,GAAG,SAAS;IAK7B;;;;;;;;;OASG;WACW,kBAAkB,CAC9B,YAAY,EAAE,SAAS,EAAE,SAAS,EAAE,MAAM,EAAE,SAAS,EAAE,MAAM,EAAE,cAAc,CAAC,EAAE,SAAS,GACxF,cAAc,GAAG,SAAS;IAG7B;;;;;;;;OAQG;WACW,uBAAuB,CACnC,YAAY,EAAE,SAAS,EAAE,SAAS,EAAE,MAAM,EAAE,SAAS,EAAE,MAAM,EAAE,cAAc,CAAC,EAAE,SAAS,GACxF,cAAc,GAAG,SAAS;IAG7B;;;;;;;;;;;;;;;;;;;OAmBG;WACW,yBAAyB,CACrC,UAAU,EAAE,MAAM,EAClB,OAAO,EAAE,MAAM,GAAG,SAAS,EAC3B,OAAO,EAAE,MAAM,GAAG,SAAS,EAC3B,QAAQ,EAAE,KAAK,GAAG,SAAS,EAC3B,SAAS,EAAE,KAAK,GAAG,SAAS,EAC5B,SAAS,EAAE,MAAM,GAAG,SAAS,EAC7B,cAAc,EAAE,SAAS,GAAG,SAAS,EACrC,YAAY,EAAE,SAAS,GACtB,kBAAkB,GAAG,SAAS;IA+BjC,iCAAiC;IACjB,KAAK,IAAI,cAAc;IAWvC,mEAAmE;IAC5D,mBAAmB,CAAC,UAAU,EAAE,SAAS,GAAG,OAAO;IAU1D,qCAAqC;IACrB,UAAU,CAAC,MAAM,CAAC,EAAE,OAAO,GAAG,OAAO;IAGrD,mCAAmC;IACnB,QAAQ,CAAC,MAAM,CAAC,EAAE,OAAO,GAAG,OAAO;IAGnD,oFAAoF;IAC7E,SAAS,CAAC,KAAK,EAAE,4BAA4B,GAAG,OAAO;IAK9D;;;OAGG;IACI,WAAW;IAUlB,6DAA6D;IACtD,mBAAmB,CAAC,KAAK,EAAE,GAAG,GAAG,OAAO;IAG/C;;;;OAIG;IACI,WAAW,CAAC,IAAI,EAAE,YAAY,EAAE,OAAO,CAAC,EAAE,aAAa,GAAG,IAAI;IAGrE,+CAA+C;IACxC,kBAAkB,CAAC,IAAI,EAAE,cAAc,EAAE,OAAO,CAAC,EAAE,aAAa,GAAG,IAAI;IAmB9E;;;OAGG;IACI,4BAA4B,CAAC,OAAO,CAAC,EAAE,aAAa,GAAG,MAAM;IAcpE;;;OAGG;IACI,cAAc,IAAI,IAAI;IAM7B,qDAAqD;IAC9C,eAAe,CAAC,cAAc,EAAE,MAAM,EAAE,MAAM,CAAC,EAAE,OAAO,GAAG,OAAO;IAMzE,oEAAoE;IAC7D,4BAA4B,CAAC,cAAc,EAAE,MAAM,EAAE,MAAM,CAAC,EAAE,KAAK,GAAG,KAAK;IAOlF;;;;;;OAMG;IACI,8BAA8B,CACnC,cAAc,EAAE,MAAM,EAAE,MAAM,CAAC,EAAE,yBAAyB,GACzD,yBAAyB,GAAG,SAAS;IASxC,mFAAmF;IAC5E,yBAAyB,CAAC,OAAO,EAAE,eAAe,GAAG,GAAG;IAG/D,6CAA6C;IAC7B,aAAa,CAAC,KAAK,EAAE,GAAG,GAAG,OAAO;CAUnD"}
+{"version":3,"file":"DirectSpiral3d.d.ts","sourceRoot":"","sources":["../../../../src/curve/spiral/DirectSpiral3d.ts"],"names":[],"mappings":"AAQA,OAAO,EAAE,KAAK,EAAE,MAAM,wBAAwB,CAAC;AAC/C,OAAO,EAAE,eAAe,EAAE,cAAc,EAAE,MAAM,kCAAkC,CAAC;AACnF,OAAO,EAAE,4BAA4B,EAAE,MAAM,+CAA+C,CAAC;AAC7F,OAAO,EAAE,yBAAyB,EAAE,MAAM,4CAA4C,CAAC;AACvF,OAAO,EAAE,OAAO,EAAE,MAAM,kCAAkC,CAAC;AAC3D,OAAO,EAAE,KAAK,EAAE,MAAM,wBAAwB,CAAC;AAC/C,OAAO,EAAE,SAAS,EAAE,MAAM,4BAA4B,CAAC;AACvD,OAAO,EAAE,SAAS,EAAE,MAAM,4BAA4B,CAAC;AAEvD,OAAO,EAAE,YAAY,EAAE,MAAM,iBAAiB,CAAC;AAC/C,OAAO,EAAE,aAAa,EAAE,MAAM,kBAAkB,CAAC;AAOjD,OAAO,EAAE,+BAA+B,EAAE,MAAM,mCAAmC,CAAC;AACpF,OAAO,EAAE,kBAAkB,EAAE,MAAM,sBAAsB,CAAC;AAC1D,OAAO,EAAE,gBAAgB,EAAE,MAAM,oBAAoB,CAAC;AAEtD;;;;;;;;;;;;;;;;EAgBE;AAEF,qBAAa,cAAe,SAAQ,kBAAkB;IACpD,yCAAyC;IACzC,SAAgB,kBAAkB,sBAAsB;IACxD,wFAAwF;IACxF,OAAO,CAAC,cAAc,CAAe;IACrC;;;;;OAKG;IACH,OAAO,CAAC,cAAc,CAAC,CAAe;IACtC,yEAAyE;IACzE,IAAW,aAAa,IAAI,YAAY,CAEvC;IACD,OAAO,CAAC,UAAU,CAAS;IAC3B,OAAO,CAAC,UAAU,CAAS;IAC3B,OAAO,CAAC,UAAU,CAAmB;IACrC,qCAAqC;IACrC,IAAW,SAAS,IAAI,MAAM,CAE7B;IACD,+DAA+D;IAC/D,IAAW,SAAS,IAAI,MAAM,CAE7B;IACD,uCAAuC;IACvC,IAAW,iBAAiB,IAAI,MAAM,CAErC;IACD;;;OAGG;IACH,IAAW,SAAS,IAAI,gBAAgB,CAEvC;gBAIC,YAAY,EAAE,SAAS,EACvB,UAAU,EAAE,MAAM,GAAG,SAAS,EAC9B,kBAAkB,EAAE,+BAA+B,GAAG,SAAS,EAC/D,SAAS,EAAE,MAAM,EACjB,SAAS,EAAE,MAAM,EACjB,sBAAsB,EAAE,SAAS,GAAG,SAAS,EAC7C,SAAS,EAAE,gBAAgB;IAW7B;;;;;;;OAOG;IACH,OAAO,CAAC,cAAc;IA6BtB,yBAAyB;IACT,yBAAyB;IAWzC;;;;;;;;;;OAUG;WACW,uBAAuB,CACnC,UAAU,EAAE,MAAM,EAClB,YAAY,EAAE,SAAS,EACvB,QAAQ,EAAE,MAAM,EAChB,QAAQ,EAAE,MAAM,EAChB,kBAAkB,EAAE,+BAA+B,GAAG,SAAS,EAC/D,SAAS,EAAE,MAAM,EACjB,SAAS,EAAE,MAAM,EACjB,cAAc,EAAE,SAAS,GAAG,SAAS,GACpC,cAAc,GAAG,SAAS;IAgB7B;;;;;;;;;;;;OAYG;WACW,gBAAgB,CAC5B,YAAY,EAAE,SAAS,EAAE,SAAS,EAAE,MAAM,EAAE,SAAS,EAAE,MAAM,EAAE,cAAc,CAAC,EAAE,SAAS,GACxF,cAAc,GAAG,SAAS;IAc7B;;;;;;;;;;;;OAYG;WACW,aAAa,CACzB,YAAY,EAAE,SAAS,EAAE,SAAS,EAAE,MAAM,EAAE,SAAS,EAAE,MAAM,EAAE,cAAc,CAAC,EAAE,SAAS,GACxF,cAAc,GAAG,SAAS;IAc7B;;;;;;;;;;;;OAYG;WACW,qBAAqB,CACjC,YAAY,EAAE,SAAS,EAAE,SAAS,EAAE,MAAM,EAAE,SAAS,EAAE,MAAM,EAAE,cAAc,CAAC,EAAE,SAAS,GACxF,cAAc,GAAG,SAAS;IAc7B;;;;;;;;;OASG;WACW,iBAAiB,CAC7B,YAAY,EAAE,SAAS,EAAE,SAAS,EAAE,MAAM,EAAE,SAAS,EAAE,MAAM,EAAE,cAAc,CAAC,EAAE,SAAS,GACxF,cAAc,GAAG,SAAS;IAc7B;;;;;;;;;;;;OAYG;WACW,oBAAoB,CAChC,YAAY,EAAE,SAAS,EAAE,SAAS,EAAE,MAAM,EAAE,SAAS,EAAE,MAAM,EAAE,cAAc,CAAC,EAAE,SAAS,GACxF,cAAc,GAAG,SAAS;WAcf,sBAAsB,CAClC,YAAY,EAAE,SAAS,EAAE,SAAS,EAAE,MAAM,EAAE,SAAS,EAAE,MAAM,EAAE,cAAc,CAAC,EAAE,SAAS,GACxF,cAAc,GAAG,SAAS;IAW7B;;;;;;;;;OASG;WACW,WAAW,CACvB,YAAY,EAAE,SAAS,EAAE,SAAS,EAAE,MAAM,EAAE,SAAS,EAAE,MAAM,EAAE,cAAc,CAAC,EAAE,SAAS,GACxF,cAAc,GAAG,SAAS;IAG7B;;;;;;;;;OASG;WACW,mBAAmB,CAC/B,YAAY,EAAE,SAAS,EAAE,SAAS,EAAE,MAAM,EAAE,SAAS,EAAE,MAAM,EAAE,cAAc,CAAC,EAAE,SAAS,GACxF,cAAc,GAAG,SAAS;IAK7B;;;;;;;;;OASG;WACW,kBAAkB,CAC9B,YAAY,EAAE,SAAS,EAAE,SAAS,EAAE,MAAM,EAAE,SAAS,EAAE,MAAM,EAAE,cAAc,CAAC,EAAE,SAAS,GACxF,cAAc,GAAG,SAAS;IAG7B;;;;;;;;OAQG;WACW,uBAAuB,CACnC,YAAY,EAAE,SAAS,EAAE,SAAS,EAAE,MAAM,EAAE,SAAS,EAAE,MAAM,EAAE,cAAc,CAAC,EAAE,SAAS,GACxF,cAAc,GAAG,SAAS;IAG7B;;;;;;;;;;;;;;;;;;;OAmBG;WACW,yBAAyB,CACrC,UAAU,EAAE,MAAM,EAClB,OAAO,EAAE,MAAM,GAAG,SAAS,EAC3B,OAAO,EAAE,MAAM,GAAG,SAAS,EAC3B,QAAQ,EAAE,KAAK,GAAG,SAAS,EAC3B,SAAS,EAAE,KAAK,GAAG,SAAS,EAC5B,SAAS,EAAE,MAAM,GAAG,SAAS,EAC7B,cAAc,EAAE,SAAS,GAAG,SAAS,EACrC,YAAY,EAAE,SAAS,GACtB,kBAAkB,GAAG,SAAS;IA+BjC,iCAAiC;IACjB,KAAK,IAAI,cAAc;IAWvC,mEAAmE;IAC5D,mBAAmB,CAAC,UAAU,EAAE,SAAS,GAAG,OAAO;IAU1D,qCAAqC;IACrB,UAAU,CAAC,MAAM,CAAC,EAAE,OAAO,GAAG,OAAO;IAGrD,mCAAmC;IACnB,QAAQ,CAAC,MAAM,CAAC,EAAE,OAAO,GAAG,OAAO;IAGnD,oFAAoF;IAC7E,SAAS,CAAC,KAAK,EAAE,4BAA4B,GAAG,OAAO;IAK9D;;;OAGG;IACI,WAAW;IAUlB,6DAA6D;IACtD,mBAAmB,CAAC,KAAK,EAAE,GAAG,GAAG,OAAO;IAG/C;;;;OAIG;IACI,WAAW,CAAC,IAAI,EAAE,YAAY,EAAE,OAAO,CAAC,EAAE,aAAa,GAAG,IAAI;IAGrE,+CAA+C;IACxC,kBAAkB,CAAC,IAAI,EAAE,cAAc,EAAE,OAAO,CAAC,EAAE,aAAa,GAAG,IAAI;IAmB9E;;;OAGG;IACI,4BAA4B,CAAC,OAAO,CAAC,EAAE,aAAa,GAAG,MAAM;IAcpE;;;OAGG;IACI,cAAc,IAAI,IAAI;IAM7B,qDAAqD;IAC9C,eAAe,CAAC,cAAc,EAAE,MAAM,EAAE,MAAM,CAAC,EAAE,OAAO,GAAG,OAAO;IAMzE,oEAAoE;IAC7D,4BAA4B,CAAC,cAAc,EAAE,MAAM,EAAE,MAAM,CAAC,EAAE,KAAK,GAAG,KAAK;IAOlF;;;;;;OAMG;IACI,8BAA8B,CACnC,cAAc,EAAE,MAAM,EAAE,MAAM,CAAC,EAAE,yBAAyB,GACzD,yBAAyB,GAAG,SAAS;IASxC,mFAAmF;IAC5E,yBAAyB,CAAC,OAAO,EAAE,eAAe,GAAG,GAAG;IAG/D,6CAA6C;IAC7B,aAAa,CAAC,KAAK,EAAE,GAAG,GAAG,OAAO;CAUnD"}

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

@@ -21,18 +21,20 @@ import { TransitionSpiral3d } from "./TransitionSpiral3d";
 * DirectSpiral3d acts like a TransitionSpiral3d for serialization purposes, but implements spiral types that have
 * "direct" xy calculations without the integrations required for IntegratedSpiral3d.
 * * Each DirectSpiral3d carries an XYCurveEvaluator to give it specialized behavior.
-* * Direct spirals that flow through serialization to native imodel02 are created with these static methods:
+* * Direct spirals are created with these static methods:
 *   * createArema
 *   * createJapaneseCubic
-*   * createAustralianRail
-*   * createDirectHalfCosine
 *   * createChineseCubic
+*   * createWesternAustralian
+*   * createDirectHalfCosine
+*   * createAustralianRail
 *   * createCzechCubic
 *   * createPolishCubic
+*   * createMXCubicAlongArc
 *   * createItalian
-*   * createWesternAustralian
 * @public
 */
+// see internaldocs/Spiral.md for more info
 export class DirectSpiral3d extends TransitionSpiral3d {
     /** String name for schema properties. */
     curvePrimitiveType = "transitionSpiral";
@@ -188,11 +190,11 @@ export class DirectSpiral3d extends TransitionSpiral3d {
     /**
      * Create an MX Cubic whose nominal length is close to along the curve.
      * This is y = m*x^3 with
-     * * m is 1/(6RL1).
+     * * m is 1/(6RL).
      *    * 1/(6RL) is the leading term of the sine series.
      * * L1 is an along-the-x-axis distance that is slightly LESS THAN the nominal length.
      * * x is axis position that is slightly LESS than nominal distance along.
-     * * L1, x use the approximation `x = s * ( 1 - s^4/ (40 R R L L))
+     * * L1, x use the approximation `x = s * ( 1 - s^4/ (40 RR LL))
      * @param localToWorld
      * @param nominalL1
      * @param nominalR1