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

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

@@ -0,0 +1 @@
+{"version":3,"file":"PolishCubicSpiralEvaluator.js","sourceRoot":"","sources":["../../../../src/curve/spiral/PolishCubicSpiralEvaluator.ts"],"names":[],"mappings":"AAAA;;;+FAG+F;AAE/F;;GAEG;AACH,OAAO,EAAE,QAAQ,EAAE,MAAM,gBAAgB,CAAC;AAC1C,OAAO,EAAE,YAAY,EAAE,MAAM,uBAAuB,CAAC;AACrD,OAAO,EAAE,cAAc,EAAE,MAAM,kBAAkB,CAAC;AAClD;;;;;;;;;;;;;GAaG;AACH,MAAM,OAAO,oBAAqB,SAAQ,cAAc;IAItD,uEAAuE;IACvE,YAAoB,OAAe,EAAE,OAAe,EAAE,UAAkB,EAAE,MAAc;QACtF,KAAK,CAAC,UAAU,EAAE,MAAM,CAAC,CAAC;QAC1B,IAAI,CAAC,cAAc,GAAG,OAAO,CAAC;QAC9B,IAAI,CAAC,cAAc,GAAG,OAAO,CAAC;IAChC,CAAC;IAED,wCAAwC;IACjC,MAAM,CAAC,aAAa,CAAC,OAAe,EAAE,OAAe;QAC1D,OAAO,GAAG,GAAG,CAAC,GAAG,GAAG,OAAO,GAAG,OAAO,CAAC,CAAC;IACzC,CAAC;IAEM,MAAM,CAAC,MAAM,CAAC,OAAe,EAAE,OAAe;QACnD,MAAM,CAAC,GAAG,IAAI,CAAC,aAAa,CAAC,OAAO,EAAE,OAAO,CAAC,CAAC;QAC/C,IAAI,CAAC,KAAK,SAAS;YACjB,OAAO,SAAS,CAAC;QACnB,MAAM,IAAI,GAAG,oBAAoB,CAAC,2BAA2B,CAAC,OAAO,EAAE,OAAO,EAAE,OAAO,CAAC,CAAC;QACzF,IAAI,IAAI,KAAK,SAAS;YACpB,OAAO,SAAS,CAAC;QACnB,OAAO,IAAI,oBAAoB,CAAC,OAAO,EAAE,OAAO,EAAE,IAAI,EAAE,CAAC,CAAC,CAAC;IAC7D,CAAC;IAEe,YAAY,CAAC,WAAmB;QAC9C,IAAI,CAAC,cAAc,IAAI,WAAW,CAAC;QACnC,IAAI,CAAC,cAAc,IAAI,WAAW,CAAC;QACnC,KAAK,CAAC,YAAY,CAAC,WAAW,CAAC,CAAC;IAClC,CAAC;IACD,0CAA0C;IACnC,KAAK,KAA2B,OAAO,IAAI,oBAAoB,CAAC,IAAI,CAAC,cAAc,EAAE,IAAI,CAAC,cAAc,EAAE,KAAK,CAAC,WAAW,EAAE,IAAI,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC;IACnJ,mCAAmC;IAC5B,aAAa,CAAC,KAAU;QAC7B,IAAI,KAAK,YAAY,oBAAoB,EAAE;YACzC,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;SAC3E;QACD,OAAO,KAAK,CAAC;IACf,CAAC;IAED,2EAA2E;IACpE,MAAM,CAAC,0BAA0B,CAAC,OAAe,EAAE,OAAe;QACvE,OAAO,GAAG,GAAG,CAAC,GAAG,GAAG,OAAO,GAAG,OAAO,GAAG,OAAO,GAAG,OAAO,CAAC,CAAC;IAC7D,CAAC;IACD;;;;;;OAMG;IACI,MAAM,CAAC,sBAAsB,CAAC,CAAS,EAAE,OAAe,EAAE,OAAe;QAC9E,mFAAmF;QACnF,MAAM,EAAE,GAAG,IAAI,CAAC,0BAA0B,CAAC,OAAO,EAAE,OAAO,CAAC,CAAC;QAC7D,MAAM,GAAG,GAAG,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC;QACvB,MAAM,GAAG,GAAG,GAAG,GAAG,CAAC,CAAC;QACpB,MAAM,GAAG,GAAG,GAAG,GAAG,CAAC,CAAC;QACpB,MAAM,EAAE,GAAG,CAAC,GAAG,CAAC,GAAG,GAAG,GAAG,GAAG,CAAC,GAAG,GAAG,GAAG,GAAG,CAAC,CAAC,GAAG,GAAG,IAAI,GAAG,GAAG,GAAG,CAAC,GAAG,GAAG,KAAK,GAAG,GAAG,GAAG,GAAG,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC;QACpG,OAAO,EAAE,CAAC;IACZ,CAAC;IAED;;;;;;OAMG;IACI,MAAM,CAAC,gCAAgC,CAAC,CAAS,EAAE,OAAe,EAAE,OAAe;QACxF,mFAAmF;QACnF,MAAM,EAAE,GAAG,IAAI,CAAC,0BAA0B,CAAC,OAAO,EAAE,OAAO,CAAC,CAAC;QAC7D,MAAM,GAAG,GAAG,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC;QACvB,MAAM,GAAG,GAAG,GAAG,GAAG,CAAC,CAAC;QACpB,MAAM,GAAG,GAAG,GAAG,GAAG,CAAC,CAAC;QACpB,sBAAsB;QACtB,sDAAsD;QACtD,mDAAmD;QACnD,aAAa;QACb,qBAAqB;QACrB,sBAAsB;QACtB,0BAA0B;QAC1B,MAAM,EAAE,GAAG,CAAC,GAAG,GAAG,GAAG,GAAG,CAAC,GAAG,GAAG,GAAG,GAAG,CAAC,CAAC,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,CAAC,GAAG,GAAG,IAAI,GAAG,GAAG,GAAG,GAAG,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC;QAC7F,OAAO,EAAE,CAAC;IACZ,CAAC;IAED,kDAAkD;IAC3C,MAAM,CAAC,2BAA2B,CAAC,CAAS,EAAE,OAAe,EAAE,OAAe;QACnF,MAAM,IAAI,GAAG,YAAY,CAAC,WAAW,CAAC,CAAC,EACrC,CAAC,CAAS,EAAE,EAAE,CAAC,CAAC,IAAI,CAAC,sBAAsB,CAAC,CAAC,EAAE,OAAO,EAAE,OAAO,CAAC,GAAG,CAAC,CAAC,EACrE,CAAC,CAAS,EAAE,EAAE,CAAC,IAAI,CAAC,gCAAgC,CAAC,CAAC,EAAE,OAAO,EAAE,OAAO,CAAC,CAAC,CAAC;QAC7E,OAAO,IAAI,CAAC;IACd,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 { SimpleNewton } from \"../../numerics/Newton\";\r\nimport { CubicEvaluator } from \"./CubicEvaluator\";\r\n/**\r\n * Polish Cubic.\r\n * * Construction takes nominal length and end radius.\r\n * curve is  is y= m*x^3 with\r\n * * x any point on the x axis\r\n * * m is (1/6RL)\r\n * * Lx = x length is along the axis, determined by inversion of a distance series at nominal length\r\n * *\r\n * @param localToWorld\r\n * @param nominalL1\r\n * @param nominalR1\r\n * @param activeInterval\r\n * @internal\r\n */\r\nexport class PolishCubicEvaluator extends CubicEvaluator {\r\n  public nominalLength1: number;\r\n  public nominalRadius1: number;\r\n\r\n  /** Constructor is private.  Caller responsible for cubicM validity. */\r\n  private constructor(length1: number, radius1: number, axisLength: number, cubicM: number) {\r\n    super(axisLength, cubicM);\r\n    this.nominalLength1 = length1;\r\n    this.nominalRadius1 = radius1;\r\n  }\r\n\r\n  /** Compute the czech cubic constant. */\r\n  public static computeCubicM(length1: number, radius1: number): number {\r\n    return 1.0 / (6.0 * length1 * radius1);\r\n  }\r\n\r\n  public static create(length1: number, radius1: number): PolishCubicEvaluator | undefined {\r\n    const m = this.computeCubicM(length1, radius1);\r\n    if (m === undefined)\r\n      return undefined;\r\n    const xMax = PolishCubicEvaluator.approximateDistanceAlongToX(length1, radius1, length1);\r\n    if (xMax === undefined)\r\n      return undefined;\r\n    return new PolishCubicEvaluator(length1, radius1, xMax, m);\r\n  }\r\n\r\n  public override scaleInPlace(scaleFactor: number) {\r\n    this.nominalLength1 *= scaleFactor;\r\n    this.nominalRadius1 *= scaleFactor;\r\n    super.scaleInPlace(scaleFactor);\r\n  }\r\n  /** return a deep copy of the evaluator */\r\n  public clone(): PolishCubicEvaluator { return new PolishCubicEvaluator(this.nominalLength1, this.nominalRadius1, super._axisLength, this.cubicM); }\r\n  /** Member by member matchup ... */\r\n  public isAlmostEqual(other: any): boolean {\r\n    if (other instanceof PolishCubicEvaluator) {\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  /** Compute the coefficient of x^4 in the x-to-distance series expansion */\r\n  public static computeX4SeriesCoefficient(length1: number, radius1: number): number {\r\n    return 1.0 / (4.0 * length1 * length1 * radius1 * radius1);\r\n  }\r\n  /**\r\n   * Evaluate a series approximation of distance along the true curve.\r\n   * @param x distance along x axis\r\n   * @param radius1 nominal end radius\r\n   * @param length1 nominal length along curve\r\n   * @returns\r\n   */\r\n  public static xToApproximateDistance(x: number, radius1: number, length1: number): number {\r\n    // C31 * ( 1 + 1 / 10 * E31 - 1 / 72 * E31^2 + 1 / 208 * E31^3 - 5 / 2176 * E31^4 )\r\n    const a4 = this.computeX4SeriesCoefficient(length1, radius1);\r\n    const ax2 = a4 * x * x;\r\n    const ax3 = ax2 * x;\r\n    const ax4 = ax3 * x;\r\n    const s0 = x * (1.0 + ax4 * (0.1 + ax4 * (-1.0 / 72.0 + ax4 * (1.0 / 208.0 - 5.0 * ax4 / 2176.0))));\r\n    return s0;\r\n  }\r\n\r\n  /**\r\n   * Evaluate the derivative of the x-to-distance series.\r\n   * @param x distance along x axis\r\n   * @param radius1 nominal end radius\r\n   * @param length1 nominal length along curve\r\n   * @returns\r\n   */\r\n  public static xToApproximateDistanceDerivative(x: number, radius1: number, length1: number): number {\r\n    // C31 * ( 1 + 1 / 10 * E31 - 1 / 72 * E31^2 + 1 / 208 * E31^3 - 5 / 2176 * E31^4 )\r\n    const a4 = this.computeX4SeriesCoefficient(length1, radius1);\r\n    const ax2 = a4 * x * x;\r\n    const ax3 = ax2 * x;\r\n    const ax4 = ax3 * x;\r\n    // derivative notes ..\r\n    // take away leading x -- this reduces each power by 1\r\n    // multiply each coefficient by its original power:\r\n    //  0.1==>0.5\r\n    // 1/72==> 9/72 = 1/8\r\n    // 1/208==>13/208=1/16\r\n    // 1/2176==>17/2176= 1/128\r\n    const ds = (1.0 + ax4 * (0.5 + ax4 * (-1.0 / 8.0 + ax4 * (1.0 / 16.0 - 5.0 * ax4 / 128.0))));\r\n    return ds;\r\n  }\r\n\r\n  /** Invert the xToApproximateDistance function. */\r\n  public static approximateDistanceAlongToX(s: number, radius1: number, length1: number): number | undefined {\r\n    const root = SimpleNewton.runNewton1D(s,\r\n      (x: number) => (this.xToApproximateDistance(x, radius1, length1) - s),\r\n      (x: number) => this.xToApproximateDistanceDerivative(x, radius1, length1));\r\n    return root;\r\n  }\r\n}\r\n"]}

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

@@ -0,0 +1,53 @@
+import { Angle } from "../../geometry3d/Angle";
+/** A true transition spiral is a curve defined by its curvature, with the curvature function symmetric about midpoint.
+ * * The symmetry condition creates a relationship among the following 4 quantities:
+ * ** curvature0 = curvature (i.e. 1/radius) at start
+ * ** curvature1 = curvature (i.e. 1/radius) at end
+ * ** sweepRadians = signed turning angle from start to end
+ * ** arcLength = length of curve
+ * * The relationship is the equation
+ * ** `sweepRadians = arcLength * average Curvature = arcLength * 0.5 * (curvature0 + curvature1)`
+ * * That is, regardless of any curvature properties other than symmetry, specifying any 3 of the quantities fully determines the remaining one.
+ * @public
+ */
+export declare class TransitionConditionalProperties {
+    /** radius (or 0 at start) */
+    radius0: number | undefined;
+    /** radius (or 0) at end */
+    radius1: number | undefined;
+    /** bearing at start, measured from x towards y */
+    bearing0: Angle | undefined;
+    /** bearing at end, measured from x towards y */
+    bearing1: Angle | undefined;
+    /** curve length */
+    curveLength: number | undefined;
+    /**
+     * capture numeric or undefined values
+     * @param radius0 start radius or undefined
+     * @param radius1 end radius or undefined
+     * @param bearing0 start bearing or undefined
+     * @param bearing1 end bearing or undefined
+     * @param arcLength arc length or undefined
+     */
+    constructor(radius0: number | undefined, radius1: number | undefined, bearing0: Angle | undefined, bearing1: Angle | undefined, arcLength: number | undefined);
+    /** return the number of defined values among the 5 properties. */
+    numDefinedProperties(): number;
+    /** clone with all properties (i.e. preserve undefined states) */
+    clone(): TransitionConditionalProperties;
+    /** Return true if all components are defined and agree equationally. */
+    getIsValidCompleteSet(): boolean;
+    /** Examine which properties are defined and compute the (single) undefined.
+     * @returns Return true if the input state had precisely one undefined member.
+     */
+    tryResolveAnySingleUnknown(): boolean;
+    private almostEqualCoordinate;
+    private almostEqualBearing;
+    /**
+     * Test if this and other have matching numeric and undefined members.
+     */
+    isAlmostEqual(other?: TransitionConditionalProperties): boolean;
+    /** Apply a NONZERO scale factor to all distances. */
+    applyScaleFactor(a: number): void;
+    static areAlmostEqual(a: TransitionConditionalProperties | undefined, b: TransitionConditionalProperties | undefined): boolean;
+}
+//# sourceMappingURL=TransitionConditionalProperties.d.ts.map

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

@@ -0,0 +1 @@
+{"version":3,"file":"TransitionConditionalProperties.d.ts","sourceRoot":"","sources":["../../../../src/curve/spiral/TransitionConditionalProperties.ts"],"names":[],"mappings":"AASA,OAAO,EAAE,KAAK,EAAE,MAAM,wBAAwB,CAAC;AAG/C;;;;;;;;;;GAUG;AACH,qBAAa,+BAA+B;IAC1C,6BAA6B;IACtB,OAAO,EAAE,MAAM,GAAG,SAAS,CAAC;IACnC,2BAA2B;IACpB,OAAO,EAAE,MAAM,GAAG,SAAS,CAAC;IACnC,kDAAkD;IAC3C,QAAQ,EAAE,KAAK,GAAG,SAAS,CAAC;IACnC,gDAAgD;IACzC,QAAQ,EAAE,KAAK,GAAG,SAAS,CAAC;IACnC,mBAAmB;IACZ,WAAW,EAAE,MAAM,GAAG,SAAS,CAAC;IACvC;;;;;;;OAOG;gBAED,OAAO,EAAE,MAAM,GAAG,SAAS,EAC3B,OAAO,EAAE,MAAM,GAAG,SAAS,EAC3B,QAAQ,EAAE,KAAK,GAAG,SAAS,EAC3B,QAAQ,EAAE,KAAK,GAAG,SAAS,EAC3B,SAAS,EAAE,MAAM,GAAG,SAAS;IAO/B,kEAAkE;IAC3D,oBAAoB;IAO3B,iEAAiE;IAC1D,KAAK,IAAI,+BAA+B;IAQ/C,wEAAwE;IACjE,qBAAqB;IAS5B;;OAEG;IACI,0BAA0B,IAAI,OAAO;IAkC5C,OAAO,CAAC,qBAAqB;IAO7B,OAAO,CAAC,kBAAkB;IAQ1B;;OAEG;IAEI,aAAa,CAAC,KAAK,CAAC,EAAE,+BAA+B;IAe5D,qDAAqD;IAC9C,gBAAgB,CAAC,CAAC,EAAE,MAAM;WAQnB,cAAc,CAAC,CAAC,EAAE,+BAA+B,GAAG,SAAS,EAAE,CAAC,EAAE,+BAA+B,GAAG,SAAS,GAAG,OAAO;CAKtI"}

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

@@ -0,0 +1,142 @@
+/*---------------------------------------------------------------------------------------------
+* Copyright (c) Bentley Systems, Incorporated. All rights reserved.
+* See LICENSE.md in the project root for license terms and full copyright notice.
+*--------------------------------------------------------------------------------------------*/
+/** @packageDocumentation
+ * @module Curve
+ */
+import { Geometry } from "../../Geometry";
+import { Angle } from "../../geometry3d/Angle";
+import { TransitionSpiral3d } from "./TransitionSpiral3d";
+// import {} from "./";
+/** A true transition spiral is a curve defined by its curvature, with the curvature function symmetric about midpoint.
+ * * The symmetry condition creates a relationship among the following 4 quantities:
+ * ** curvature0 = curvature (i.e. 1/radius) at start
+ * ** curvature1 = curvature (i.e. 1/radius) at end
+ * ** sweepRadians = signed turning angle from start to end
+ * ** arcLength = length of curve
+ * * The relationship is the equation
+ * ** `sweepRadians = arcLength * average Curvature = arcLength * 0.5 * (curvature0 + curvature1)`
+ * * That is, regardless of any curvature properties other than symmetry, specifying any 3 of the quantities fully determines the remaining one.
+ * @public
+ */
+export class TransitionConditionalProperties {
+    /**
+     * capture numeric or undefined values
+     * @param radius0 start radius or undefined
+     * @param radius1 end radius or undefined
+     * @param bearing0 start bearing or undefined
+     * @param bearing1 end bearing or undefined
+     * @param arcLength arc length or undefined
+     */
+    constructor(radius0, radius1, bearing0, bearing1, arcLength) {
+        this.radius0 = radius0;
+        this.radius1 = radius1;
+        this.bearing0 = bearing0;
+        this.bearing1 = bearing1;
+        this.curveLength = arcLength;
+    }
+    /** return the number of defined values among the 5 properties. */
+    numDefinedProperties() {
+        return Geometry.defined01(this.radius0)
+            + Geometry.defined01(this.radius1)
+            + Geometry.defined01(this.bearing0)
+            + Geometry.defined01(this.bearing1)
+            + Geometry.defined01(this.curveLength);
+    }
+    /** clone with all properties (i.e. preserve undefined states) */
+    clone() {
+        return new TransitionConditionalProperties(this.radius0, this.radius1, this.bearing0 === undefined ? undefined : this.bearing0.clone(), this.bearing1 === undefined ? undefined : this.bearing1.clone(), this.curveLength);
+    }
+    /** Return true if all components are defined and agree equationally. */
+    getIsValidCompleteSet() {
+        if (this.curveLength !== undefined && this.bearing0 !== undefined && this.bearing1 !== undefined
+            && this.radius0 !== undefined && this.radius1 !== undefined) {
+            const length1 = TransitionSpiral3d.radiusRadiusSweepRadiansToArcLength(this.radius0, this.radius1, this.bearing1.radians - this.bearing0.radians);
+            return Geometry.isSameCoordinate(this.curveLength, length1);
+        }
+        return false;
+    }
+    /** Examine which properties are defined and compute the (single) undefined.
+     * @returns Return true if the input state had precisely one undefined member.
+     */
+    tryResolveAnySingleUnknown() {
+        if (this.getIsValidCompleteSet())
+            return true;
+        if (this.bearing0 && this.bearing1) {
+            const sweepRadians = this.bearing1.radians - this.bearing0.radians;
+            if (this.curveLength === undefined && this.radius0 !== undefined && this.radius1 !== undefined) {
+                this.curveLength = TransitionSpiral3d.radiusRadiusSweepRadiansToArcLength(this.radius0, this.radius1, sweepRadians);
+                return true;
+            }
+            if (this.curveLength !== undefined && this.radius0 === undefined && this.radius1 !== undefined) {
+                this.radius0 = TransitionSpiral3d.radius1LengthSweepRadiansToRadius0(this.radius1, this.curveLength, sweepRadians);
+                return true;
+            }
+            if (this.curveLength !== undefined && this.radius0 !== undefined && this.radius1 === undefined) {
+                this.radius1 = TransitionSpiral3d.radius0LengthSweepRadiansToRadius1(this.radius0, this.curveLength, sweepRadians);
+                return true;
+            }
+            return false;
+        }
+        // at least one bearing is undefined ...
+        if (this.curveLength === undefined || this.radius0 === undefined || this.radius1 === undefined)
+            return false;
+        if (this.bearing0) { // bearing 1 is undefined
+            this.bearing1 = Angle.createRadians(this.bearing0.radians + TransitionSpiral3d.radiusRadiusLengthToSweepRadians(this.radius0, this.radius1, this.curveLength));
+            return true;
+        }
+        if (this.bearing1) { // bearing 0 is undefined
+            this.bearing0 = Angle.createRadians(this.bearing1.radians - TransitionSpiral3d.radiusRadiusLengthToSweepRadians(this.radius0, this.radius1, this.curveLength));
+            return true;
+        }
+        return false;
+    }
+    almostEqualCoordinate(a, b) {
+        if (a === undefined && b === undefined)
+            return true;
+        if (a !== undefined && b !== undefined)
+            return Geometry.isSameCoordinate(a, b);
+        return false;
+    }
+    almostEqualBearing(a, b) {
+        if (a === undefined && b === undefined)
+            return true;
+        if (a !== undefined && b !== undefined)
+            return a.isAlmostEqualNoPeriodShift(b);
+        return false;
+    }
+    /**
+     * Test if this and other have matching numeric and undefined members.
+     */
+    isAlmostEqual(other) {
+        if (!other)
+            return false;
+        if (!this.almostEqualCoordinate(this.radius0, other.radius0))
+            return false;
+        if (!this.almostEqualCoordinate(this.radius1, other.radius1))
+            return false;
+        if (!this.almostEqualBearing(this.bearing0, other.bearing0))
+            return false;
+        if (!this.almostEqualBearing(this.bearing1, other.bearing1))
+            return false;
+        if (!this.almostEqualCoordinate(this.curveLength, other.curveLength))
+            return false;
+        return true;
+    }
+    /** Apply a NONZERO scale factor to all distances. */
+    applyScaleFactor(a) {
+        if (this.radius0 !== undefined)
+            this.radius0 *= a;
+        if (this.radius1 !== undefined)
+            this.radius1 *= a;
+        if (this.curveLength !== undefined)
+            this.curveLength *= a;
+    }
+    static areAlmostEqual(a, b) {
+        if (a === undefined)
+            return b === undefined;
+        return a.isAlmostEqual(b);
+    }
+}
+//# sourceMappingURL=TransitionConditionalProperties.js.map

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

@@ -0,0 +1 @@
+{"version":3,"file":"TransitionConditionalProperties.js","sourceRoot":"","sources":["../../../../src/curve/spiral/TransitionConditionalProperties.ts"],"names":[],"mappings":"AAAA;;;+FAG+F;AAE/F;;GAEG;AACH,OAAO,EAAE,QAAQ,EAAE,MAAM,gBAAgB,CAAC;AAC1C,OAAO,EAAE,KAAK,EAAE,MAAM,wBAAwB,CAAC;AAC/C,OAAO,EAAE,kBAAkB,EAAE,MAAM,sBAAsB,CAAC;AAC1D,uBAAuB;AACvB;;;;;;;;;;GAUG;AACH,MAAM,OAAO,+BAA+B;IAW1C;;;;;;;OAOG;IACH,YACE,OAA2B,EAC3B,OAA2B,EAC3B,QAA2B,EAC3B,QAA2B,EAC3B,SAA6B;QAC7B,IAAI,CAAC,OAAO,GAAG,OAAO,CAAC;QACvB,IAAI,CAAC,OAAO,GAAG,OAAO,CAAC;QACvB,IAAI,CAAC,QAAQ,GAAG,QAAQ,CAAC;QACzB,IAAI,CAAC,QAAQ,GAAG,QAAQ,CAAC;QACzB,IAAI,CAAC,WAAW,GAAG,SAAS,CAAC;IAC/B,CAAC;IACD,kEAAkE;IAC3D,oBAAoB;QACzB,OAAO,QAAQ,CAAC,SAAS,CAAC,IAAI,CAAC,OAAO,CAAC;cACnC,QAAQ,CAAC,SAAS,CAAC,IAAI,CAAC,OAAO,CAAC;cAChC,QAAQ,CAAC,SAAS,CAAC,IAAI,CAAC,QAAQ,CAAC;cACjC,QAAQ,CAAC,SAAS,CAAC,IAAI,CAAC,QAAQ,CAAC;cACjC,QAAQ,CAAC,SAAS,CAAC,IAAI,CAAC,WAAW,CAAC,CAAC;IAC3C,CAAC;IACD,iEAAiE;IAC1D,KAAK;QACV,OAAO,IAAI,+BAA+B,CACxC,IAAI,CAAC,OAAO,EACZ,IAAI,CAAC,OAAO,EACZ,IAAI,CAAC,QAAQ,KAAK,SAAS,CAAC,CAAC,CAAC,SAAS,CAAC,CAAC,CAAC,IAAI,CAAC,QAAQ,CAAC,KAAK,EAAE,EAC/D,IAAI,CAAC,QAAQ,KAAK,SAAS,CAAC,CAAC,CAAC,SAAS,CAAC,CAAC,CAAC,IAAI,CAAC,QAAQ,CAAC,KAAK,EAAE,EAC/D,IAAI,CAAC,WAAW,CAAC,CAAC;IACtB,CAAC;IACD,wEAAwE;IACjE,qBAAqB;QAC1B,IAAI,IAAI,CAAC,WAAW,KAAK,SAAS,IAAI,IAAI,CAAC,QAAQ,KAAK,SAAS,IAAI,IAAI,CAAC,QAAQ,KAAK,SAAS;eAC3F,IAAI,CAAC,OAAO,KAAK,SAAS,IAAI,IAAI,CAAC,OAAO,KAAK,SAAS,EAAE;YAC7D,MAAM,OAAO,GAAG,kBAAkB,CAAC,mCAAmC,CAAC,IAAI,CAAC,OAAO,EAAE,IAAI,CAAC,OAAO,EAC/F,IAAI,CAAC,QAAQ,CAAC,OAAO,GAAG,IAAI,CAAC,QAAQ,CAAC,OAAO,CAAC,CAAC;YACjD,OAAO,QAAQ,CAAC,gBAAgB,CAAC,IAAI,CAAC,WAAW,EAAE,OAAO,CAAC,CAAC;SAC7D;QACD,OAAO,KAAK,CAAC;IACf,CAAC;IACD;;OAEG;IACI,0BAA0B;QAC/B,IAAI,IAAI,CAAC,qBAAqB,EAAE;YAC9B,OAAO,IAAI,CAAC;QACd,IAAI,IAAI,CAAC,QAAQ,IAAI,IAAI,CAAC,QAAQ,EAAE;YAClC,MAAM,YAAY,GAAG,IAAI,CAAC,QAAQ,CAAC,OAAO,GAAG,IAAI,CAAC,QAAQ,CAAC,OAAO,CAAC;YACnE,IAAI,IAAI,CAAC,WAAW,KAAK,SAAS,IAAI,IAAI,CAAC,OAAO,KAAK,SAAS,IAAI,IAAI,CAAC,OAAO,KAAK,SAAS,EAAE;gBAC9F,IAAI,CAAC,WAAW,GAAG,kBAAkB,CAAC,mCAAmC,CAAC,IAAI,CAAC,OAAO,EAAE,IAAI,CAAC,OAAO,EAAE,YAAY,CAAC,CAAC;gBACpH,OAAO,IAAI,CAAC;aACb;YACD,IAAI,IAAI,CAAC,WAAW,KAAK,SAAS,IAAI,IAAI,CAAC,OAAO,KAAK,SAAS,IAAI,IAAI,CAAC,OAAO,KAAK,SAAS,EAAE;gBAC9F,IAAI,CAAC,OAAO,GAAG,kBAAkB,CAAC,kCAAkC,CAAC,IAAI,CAAC,OAAO,EAAE,IAAI,CAAC,WAAW,EAAE,YAAY,CAAC,CAAC;gBACnH,OAAO,IAAI,CAAC;aACb;YACD,IAAI,IAAI,CAAC,WAAW,KAAK,SAAS,IAAI,IAAI,CAAC,OAAO,KAAK,SAAS,IAAI,IAAI,CAAC,OAAO,KAAK,SAAS,EAAE;gBAC9F,IAAI,CAAC,OAAO,GAAG,kBAAkB,CAAC,kCAAkC,CAAC,IAAI,CAAC,OAAO,EAAE,IAAI,CAAC,WAAW,EAAE,YAAY,CAAC,CAAC;gBACnH,OAAO,IAAI,CAAC;aACb;YACD,OAAO,KAAK,CAAC;SACd;QACD,wCAAwC;QACxC,IAAI,IAAI,CAAC,WAAW,KAAK,SAAS,IAAI,IAAI,CAAC,OAAO,KAAK,SAAS,IAAI,IAAI,CAAC,OAAO,KAAK,SAAS;YAC5F,OAAO,KAAK,CAAC;QAEf,IAAI,IAAI,CAAC,QAAQ,EAAE,EAAE,yBAAyB;YAC5C,IAAI,CAAC,QAAQ,GAAG,KAAK,CAAC,aAAa,CAAC,IAAI,CAAC,QAAQ,CAAC,OAAO,GAAG,kBAAkB,CAAC,gCAAgC,CAAC,IAAI,CAAC,OAAO,EAAE,IAAI,CAAC,OAAO,EAAE,IAAI,CAAC,WAAW,CAAC,CAAC,CAAC;YAC/J,OAAO,IAAI,CAAC;SACb;QAED,IAAI,IAAI,CAAC,QAAQ,EAAE,EAAE,yBAAyB;YAC5C,IAAI,CAAC,QAAQ,GAAG,KAAK,CAAC,aAAa,CAAC,IAAI,CAAC,QAAQ,CAAC,OAAO,GAAG,kBAAkB,CAAC,gCAAgC,CAAC,IAAI,CAAC,OAAO,EAAE,IAAI,CAAC,OAAO,EAAE,IAAI,CAAC,WAAW,CAAC,CAAC,CAAC;YAC/J,OAAO,IAAI,CAAC;SACb;QACD,OAAO,KAAK,CAAC;IACf,CAAC;IACO,qBAAqB,CAAC,CAAqB,EAAE,CAAqB;QACxE,IAAI,CAAC,KAAK,SAAS,IAAI,CAAC,KAAK,SAAS;YACpC,OAAO,IAAI,CAAC;QACd,IAAI,CAAC,KAAK,SAAS,IAAI,CAAC,KAAK,SAAS;YACpC,OAAO,QAAQ,CAAC,gBAAgB,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;QACzC,OAAO,KAAK,CAAC;IACf,CAAC;IACO,kBAAkB,CAAC,CAAoB,EAAE,CAAoB;QACnE,IAAI,CAAC,KAAK,SAAS,IAAI,CAAC,KAAK,SAAS;YACpC,OAAO,IAAI,CAAC;QACd,IAAI,CAAC,KAAK,SAAS,IAAI,CAAC,KAAK,SAAS;YACpC,OAAO,CAAC,CAAC,0BAA0B,CAAC,CAAC,CAAC,CAAC;QACzC,OAAO,KAAK,CAAC;IACf,CAAC;IAED;;OAEG;IAEI,aAAa,CAAC,KAAuC;QAC1D,IAAI,CAAC,KAAK;YACR,OAAO,KAAK,CAAC;QACf,IAAI,CAAC,IAAI,CAAC,qBAAqB,CAAC,IAAI,CAAC,OAAO,EAAE,KAAK,CAAC,OAAO,CAAC;YAC1D,OAAO,KAAK,CAAC;QACf,IAAI,CAAC,IAAI,CAAC,qBAAqB,CAAC,IAAI,CAAC,OAAO,EAAE,KAAK,CAAC,OAAO,CAAC;YAC1D,OAAO,KAAK,CAAC;QACf,IAAI,CAAC,IAAI,CAAC,kBAAkB,CAAC,IAAI,CAAC,QAAQ,EAAE,KAAK,CAAC,QAAQ,CAAC;YACzD,OAAO,KAAK,CAAC;QACf,IAAI,CAAC,IAAI,CAAC,kBAAkB,CAAC,IAAI,CAAC,QAAQ,EAAE,KAAK,CAAC,QAAQ,CAAC;YACzD,OAAO,KAAK,CAAC;QACf,IAAI,CAAC,IAAI,CAAC,qBAAqB,CAAC,IAAI,CAAC,WAAW,EAAE,KAAK,CAAC,WAAW,CAAC;YAClE,OAAO,KAAK,CAAC;QACf,OAAO,IAAI,CAAC;IACd,CAAC;IACD,qDAAqD;IAC9C,gBAAgB,CAAC,CAAS;QAC/B,IAAI,IAAI,CAAC,OAAO,KAAK,SAAS;YAC5B,IAAI,CAAC,OAAO,IAAI,CAAC,CAAC;QACpB,IAAI,IAAI,CAAC,OAAO,KAAK,SAAS;YAC5B,IAAI,CAAC,OAAO,IAAI,CAAC,CAAC;QACpB,IAAI,IAAI,CAAC,WAAW,KAAK,SAAS;YAChC,IAAI,CAAC,WAAW,IAAI,CAAC,CAAC;IAC1B,CAAC;IACM,MAAM,CAAC,cAAc,CAAC,CAA8C,EAAE,CAA8C;QACzH,IAAI,CAAC,KAAK,SAAS;YACjB,OAAO,CAAC,KAAK,SAAS,CAAC;QACzB,OAAO,CAAC,CAAC,aAAa,CAAC,CAAC,CAAC,CAAC;IAC5B,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 { Angle } from \"../../geometry3d/Angle\";\r\nimport { TransitionSpiral3d } from \"./TransitionSpiral3d\";\r\n// import {} from \"./\";\r\n/** A true transition spiral is a curve defined by its curvature, with the curvature function symmetric about midpoint.\r\n * * The symmetry condition creates a relationship among the following 4 quantities:\r\n * ** curvature0 = curvature (i.e. 1/radius) at start\r\n * ** curvature1 = curvature (i.e. 1/radius) at end\r\n * ** sweepRadians = signed turning angle from start to end\r\n * ** arcLength = length of curve\r\n * * The relationship is the equation\r\n * ** `sweepRadians = arcLength * average Curvature = arcLength * 0.5 * (curvature0 + curvature1)`\r\n * * That is, regardless of any curvature properties other than symmetry, specifying any 3 of the quantities fully determines the remaining one.\r\n * @public\r\n */\r\nexport class TransitionConditionalProperties {\r\n  /** radius (or 0 at start) */\r\n  public radius0: number | undefined;\r\n  /** radius (or 0) at end */\r\n  public radius1: number | undefined;\r\n  /** bearing at start, measured from x towards y */\r\n  public bearing0: Angle | undefined;\r\n  /** bearing at end, measured from x towards y */\r\n  public bearing1: Angle | undefined;\r\n  /** curve length */\r\n  public curveLength: number | undefined;\r\n  /**\r\n   * capture numeric or undefined values\r\n   * @param radius0 start radius or undefined\r\n   * @param radius1 end radius or undefined\r\n   * @param bearing0 start bearing or undefined\r\n   * @param bearing1 end bearing or undefined\r\n   * @param arcLength arc length or undefined\r\n   */\r\n  public constructor(\r\n    radius0: number | undefined,\r\n    radius1: number | undefined,\r\n    bearing0: Angle | undefined,\r\n    bearing1: Angle | undefined,\r\n    arcLength: number | undefined) {\r\n    this.radius0 = radius0;\r\n    this.radius1 = radius1;\r\n    this.bearing0 = bearing0;\r\n    this.bearing1 = bearing1;\r\n    this.curveLength = arcLength;\r\n  }\r\n  /** return the number of defined values among the 5 properties. */\r\n  public numDefinedProperties() {\r\n    return Geometry.defined01(this.radius0)\r\n      + Geometry.defined01(this.radius1)\r\n      + Geometry.defined01(this.bearing0)\r\n      + Geometry.defined01(this.bearing1)\r\n      + Geometry.defined01(this.curveLength);\r\n  }\r\n  /** clone with all properties (i.e. preserve undefined states) */\r\n  public clone(): TransitionConditionalProperties {\r\n    return new TransitionConditionalProperties(\r\n      this.radius0,\r\n      this.radius1,\r\n      this.bearing0 === undefined ? undefined : this.bearing0.clone(),\r\n      this.bearing1 === undefined ? undefined : this.bearing1.clone(),\r\n      this.curveLength);\r\n  }\r\n  /** Return true if all components are defined and agree equationally. */\r\n  public getIsValidCompleteSet() {\r\n    if (this.curveLength !== undefined && this.bearing0 !== undefined && this.bearing1 !== undefined\r\n      && this.radius0 !== undefined && this.radius1 !== undefined) {\r\n      const length1 = TransitionSpiral3d.radiusRadiusSweepRadiansToArcLength(this.radius0, this.radius1,\r\n        this.bearing1.radians - this.bearing0.radians);\r\n      return Geometry.isSameCoordinate(this.curveLength, length1);\r\n    }\r\n    return false;\r\n  }\r\n  /** Examine which properties are defined and compute the (single) undefined.\r\n   * @returns Return true if the input state had precisely one undefined member.\r\n   */\r\n  public tryResolveAnySingleUnknown(): boolean {\r\n    if (this.getIsValidCompleteSet())\r\n      return true;\r\n    if (this.bearing0 && this.bearing1) {\r\n      const sweepRadians = this.bearing1.radians - this.bearing0.radians;\r\n      if (this.curveLength === undefined && this.radius0 !== undefined && this.radius1 !== undefined) {\r\n        this.curveLength = TransitionSpiral3d.radiusRadiusSweepRadiansToArcLength(this.radius0, this.radius1, sweepRadians);\r\n        return true;\r\n      }\r\n      if (this.curveLength !== undefined && this.radius0 === undefined && this.radius1 !== undefined) {\r\n        this.radius0 = TransitionSpiral3d.radius1LengthSweepRadiansToRadius0(this.radius1, this.curveLength, sweepRadians);\r\n        return true;\r\n      }\r\n      if (this.curveLength !== undefined && this.radius0 !== undefined && this.radius1 === undefined) {\r\n        this.radius1 = TransitionSpiral3d.radius0LengthSweepRadiansToRadius1(this.radius0, this.curveLength, sweepRadians);\r\n        return true;\r\n      }\r\n      return false;\r\n    }\r\n    // at least one bearing is undefined ...\r\n    if (this.curveLength === undefined || this.radius0 === undefined || this.radius1 === undefined)\r\n      return false;\r\n\r\n    if (this.bearing0) { // bearing 1 is undefined\r\n      this.bearing1 = Angle.createRadians(this.bearing0.radians + TransitionSpiral3d.radiusRadiusLengthToSweepRadians(this.radius0, this.radius1, this.curveLength));\r\n      return true;\r\n    }\r\n\r\n    if (this.bearing1) { // bearing 0 is undefined\r\n      this.bearing0 = Angle.createRadians(this.bearing1.radians - TransitionSpiral3d.radiusRadiusLengthToSweepRadians(this.radius0, this.radius1, this.curveLength));\r\n      return true;\r\n    }\r\n    return false;\r\n  }\r\n  private almostEqualCoordinate(a: number | undefined, b: number | undefined): boolean {\r\n    if (a === undefined && b === undefined)\r\n      return true;\r\n    if (a !== undefined && b !== undefined)\r\n      return Geometry.isSameCoordinate(a, b);\r\n    return false;\r\n  }\r\n  private almostEqualBearing(a: Angle | undefined, b: Angle | undefined): boolean {\r\n    if (a === undefined && b === undefined)\r\n      return true;\r\n    if (a !== undefined && b !== undefined)\r\n      return a.isAlmostEqualNoPeriodShift(b);\r\n    return false;\r\n  }\r\n\r\n  /**\r\n   * Test if this and other have matching numeric and undefined members.\r\n   */\r\n\r\n  public isAlmostEqual(other?: TransitionConditionalProperties) {\r\n    if (!other)\r\n      return false;\r\n    if (!this.almostEqualCoordinate(this.radius0, other.radius0))\r\n      return false;\r\n    if (!this.almostEqualCoordinate(this.radius1, other.radius1))\r\n      return false;\r\n    if (!this.almostEqualBearing(this.bearing0, other.bearing0))\r\n      return false;\r\n    if (!this.almostEqualBearing(this.bearing1, other.bearing1))\r\n      return false;\r\n    if (!this.almostEqualCoordinate(this.curveLength, other.curveLength))\r\n      return false;\r\n    return true;\r\n  }\r\n  /** Apply a NONZERO scale factor to all distances. */\r\n  public applyScaleFactor(a: number) {\r\n    if (this.radius0 !== undefined)\r\n      this.radius0 *= a;\r\n    if (this.radius1 !== undefined)\r\n      this.radius1 *= a;\r\n    if (this.curveLength !== undefined)\r\n      this.curveLength *= a;\r\n  }\r\n  public static areAlmostEqual(a: TransitionConditionalProperties | undefined, b: TransitionConditionalProperties | undefined): boolean {\r\n    if (a === undefined)\r\n      return b === undefined;\r\n    return a.isAlmostEqual(b);\r\n  }\r\n}\r\n"]}

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

@@ -0,0 +1,101 @@
+/** @packageDocumentation
+ * @module Curve
+ */
+import { CurvePrimitive } from "../CurvePrimitive";
+import { Matrix3d } from "../../geometry3d/Matrix3d";
+import { Segment1d } from "../../geometry3d/Segment1d";
+import { Transform } from "../../geometry3d/Transform";
+import { TransitionConditionalProperties } from "./TransitionConditionalProperties";
+import { LineString3d } from "../LineString3d";
+/**
+ * This is the set of valid type names for "integrated" spirals
+ * * Behavior is expressed by a `NormalizedTransition` snap function.
+ * * The snap function varies smoothly from f(0)=0 to f(1)=1
+ * * The various snap functions are:
+ *   * clothoid: linear
+ *   * biquadratic: 2 quadratics pieced together, joining with 1st derivative continuity at f(0.) = 0.5, with zero slope f'(0)=0 and f'(1)= 0
+ *   * bloss: A single cubic with zero slope at 0 and 1
+ *   * cosine: half of a cosine wave, centered around 0.5
+ *   * sine: full period of a sine wave added to the line f(u)=u
+ * *
+ * @public
+ */
+export declare type IntegratedSpiralTypeName = "clothoid" | "bloss" | "biquadratic" | "cosine" | "sine";
+/**
+ * This is the set of valid type names for "direct" spirals.
+ * "Direct" spirals can evaluate fractionToPoint by direct equations, i.e. not requiring the numeric integrations in "Integrated" spiral types.
+ * @public
+ */
+export declare type DirectSpiralTypeName = "JapaneseCubic" | "Arema" | "ChineseCubic" | "HalfCosine" | "AustralianRailCorp" | "WesternAustralian" | "Czech" | "MXCubicAlongArc" | "Polish" | "Italian";
+/**
+ * TransitionSpiral3d is a base class for multiple variants of spirals.
+ * * The menagerie of spiral types have 2 broad categories:
+ *   * IntegratedSpiral3d -- a spiral whose direct function for curvature versus distance must be integrated to determine x,y
+ *     * The IntegratedSpiral3d types are enumerated in `IntegratedSpiralTypes`
+ *   * DirectSpiral3d -- a spiral implemented with direct calculation of x,y from fractional position along the spiral.
+ *     * The direct spiral types are enumerated in the `DirectSpiralType`
+ * * The method set for CurvePrimitive support includes a `handleTransitionSpiral(g: TransitionSpiral3d)` which receives all the spiral types.
+ * * The spiral class may impose expectations that its inflection is at the origin, with tangent along the x axis.
+ *   * This is generally necessary for direct spirals.
+ *   * This is not necessary for integrated spirals.
+ * @public
+ */
+export declare abstract class TransitionSpiral3d extends CurvePrimitive {
+    /** string name of spiral type */
+    protected _spiralType: string;
+    /** Original defining properties. */
+    protected _designProperties: TransitionConditionalProperties | undefined;
+    /** Fractional interval for the "active" part of a containing spiral.
+     * (The radius, angle, and length conditions define a complete spiral, and some portion of it is "active")
+     */
+    protected _activeFractionInterval: Segment1d;
+    /** Return (reference to) the active portion of the reference spiral. */
+    get activeFractionInterval(): Segment1d;
+    /** strokes in the active portion */
+    abstract get activeStrokes(): LineString3d;
+    /** Placement transform */
+    protected _localToWorld: Transform;
+    /** (reference to) placement transform. */
+    get localToWorld(): Transform;
+    protected constructor(spiralType: string | undefined, localToWorld: Transform, activeFractionInterval: Segment1d | undefined, designProperties: TransitionConditionalProperties | undefined);
+    get spiralType(): string;
+    /** Return 1/r with convention that if true zero is given as radius it represents infinite radius (0 curvature, straight line) */
+    static radiusToCurvature(radius: number): number;
+    /** Return 1/k with convention that if near-zero is given as curvature, its infinite radius is returned as 0 */
+    static curvatureToRadius(curvature: number): number;
+    /** Return the average of the start and end curvatures. */
+    static averageCurvature(radiusLimits: Segment1d): number;
+    /**
+     * Given two radii (or zeros for 0 curvature) return the average curvature
+     * @param r0 start radius, or 0 for line
+     * @param r1 end radius, or 0 for line
+     */
+    static averageCurvatureR0R1(r0: number, r1: number): number;
+    /**
+     * Given two radii (or zeros for 0 curvature) return the average curvature
+     * @param r0 start radius, or 0 for line
+     * @param r1 end radius, or 0 for line
+     */
+    static interpolateCurvatureR0R1(r0: number, fraction: number, r1: number): number;
+    /** Return the arc length of a transition spiral with given sweep and radius pair. */
+    static radiusRadiusSweepRadiansToArcLength(radius0: number, radius1: number, sweepRadians: number): number;
+    /** Return the turn angle for spiral of given length between two radii */
+    static radiusRadiusLengthToSweepRadians(radius0: number, radius1: number, arcLength: number): number;
+    /** Return the end radius for spiral of given start radius, length, and turn angle. */
+    static radius0LengthSweepRadiansToRadius1(radius0: number, arcLength: number, sweepRadians: number): number;
+    /** Return the start radius for spiral of given end radius, length, and turn angle. */
+    static radius1LengthSweepRadiansToRadius0(radius1: number, arcLength: number, sweepRadians: number): number;
+    /** Return the original defining properties (if any) saved by the constructor. */
+    get designProperties(): TransitionConditionalProperties | undefined;
+    /**
+     * * If transformA is rigid with uniform scale, apply the rigid part of transformA to the localToWorld transform and return the scale and rigid separation.
+     * * If not rigid, do nothing and return undefined.
+     * * Also apply the scale factor to the designProperties.
+     * @param transformA
+     */
+    protected applyRigidPartOfTransform(transformA: Transform): {
+        rigidAxes: Matrix3d;
+        scale: number;
+    } | undefined;
+}
+//# sourceMappingURL=TransitionSpiral3d.d.ts.map

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

@@ -0,0 +1 @@
+{"version":3,"file":"TransitionSpiral3d.d.ts","sourceRoot":"","sources":["../../../../src/curve/spiral/TransitionSpiral3d.ts"],"names":[],"mappings":"AAKA;;GAEG;AAEH,OAAO,EAAE,cAAc,EAAE,MAAM,mBAAmB,CAAC;AACnD,OAAO,EAAE,QAAQ,EAAE,MAAM,2BAA2B,CAAC;AACrD,OAAO,EAAE,SAAS,EAAE,MAAM,4BAA4B,CAAC;AACvD,OAAO,EAAE,SAAS,EAAE,MAAM,4BAA4B,CAAC;AACvD,OAAO,EAAE,+BAA+B,EAAE,MAAM,mCAAmC,CAAC;AAEpF,OAAO,EAAE,YAAY,EAAE,MAAM,iBAAiB,CAAC;AAC/C;;;;;;;;;;;;GAYG;AACH,oBAAY,wBAAwB,GAAG,UAAU,GAAG,OAAO,GAAG,aAAa,GAAG,QAAQ,GAAG,MAAM,CAAC;AAEhG;;;;GAIG;AACH,oBAAY,oBAAoB,GAC7B,eAAe,GACd,OAAO,GACP,cAAc,GACd,YAAY,GACZ,oBAAoB,GACpB,mBAAmB,GACnB,OAAO,GACP,iBAAiB,GACjB,QAAQ,GACR,SAAS,CACV;AAEH;;;;;;;;;;;;GAYG;AACH,8BAAsB,kBAAmB,SAAQ,cAAc;IAC7D,iCAAiC;IACjC,SAAS,CAAC,WAAW,EAAE,MAAM,CAAC;IAC9B,oCAAoC;IACpC,SAAS,CAAC,iBAAiB,EAAE,+BAA+B,GAAG,SAAS,CAAC;IAEzE;;OAEG;IACH,SAAS,CAAC,uBAAuB,EAAE,SAAS,CAAC;IAC7C,wEAAwE;IACxE,IAAW,sBAAsB,IAAI,SAAS,CAAyC;IACvF,oCAAoC;IACpC,aAAoB,aAAa,IAAI,YAAY,CAAC;IAClD,0BAA0B;IAC1B,SAAS,CAAC,aAAa,EAAE,SAAS,CAAC;IACnC,0CAA0C;IAC1C,IAAW,YAAY,IAAI,SAAS,CAA+B;IAEnE,SAAS,aAAa,UAAU,EAAE,MAAM,GAAG,SAAS,EAAE,YAAY,EAAE,SAAS,EAAE,sBAAsB,EAAE,SAAS,GAAG,SAAS,EAAE,gBAAgB,EAAE,+BAA+B,GAAG,SAAS;IAQ3L,IAAW,UAAU,IAAI,MAAM,CAA6B;IAC5D,iIAAiI;WACnH,iBAAiB,CAAC,MAAM,EAAE,MAAM,GAAG,MAAM;IAEvD,+GAA+G;WACjG,iBAAiB,CAAC,SAAS,EAAE,MAAM,GAAG,MAAM;IAM1D,0DAA0D;WAC5C,gBAAgB,CAAC,YAAY,EAAE,SAAS,GAAG,MAAM;IAG/D;;;;OAIG;WACW,oBAAoB,CAAC,EAAE,EAAE,MAAM,EAAE,EAAE,EAAE,MAAM,GAAG,MAAM;IAGlE;;;;OAIG;WACW,wBAAwB,CAAC,EAAE,EAAE,MAAM,EAAE,QAAQ,EAAE,MAAM,EAAE,EAAE,EAAE,MAAM,GAAG,MAAM;IAIxF,qFAAqF;WACvE,mCAAmC,CAAC,OAAO,EAAE,MAAM,EAAE,OAAO,EAAE,MAAM,EAAE,YAAY,EAAE,MAAM,GAAG,MAAM;IAIjH,yEAAyE;WAC3D,gCAAgC,CAAC,OAAO,EAAE,MAAM,EAAE,OAAO,EAAE,MAAM,EAAE,SAAS,EAAE,MAAM,GAAG,MAAM;IAI3G,sFAAsF;WACxE,kCAAkC,CAAC,OAAO,EAAE,MAAM,EAAE,SAAS,EAAE,MAAM,EAAE,YAAY,EAAE,MAAM;IAGzG,sFAAsF;WACxE,kCAAkC,CAAC,OAAO,EAAE,MAAM,EAAE,SAAS,EAAE,MAAM,EAAE,YAAY,EAAE,MAAM;IAGzG,iFAAiF;IACjF,IAAW,gBAAgB,IAAI,+BAA+B,GAAG,SAAS,CAAmC;IAC7G;;;;;OAKG;IACH,SAAS,CAAC,yBAAyB,CAAC,UAAU,EAAE,SAAS,GAAG;QAAE,SAAS,EAAE,QAAQ,CAAC;QAAC,KAAK,EAAE,MAAM,CAAA;KAAE,GAAG,SAAS;CAgB/G"}

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

@@ -0,0 +1,106 @@
+/*---------------------------------------------------------------------------------------------
+* Copyright (c) Bentley Systems, Incorporated. All rights reserved.
+* See LICENSE.md in the project root for license terms and full copyright notice.
+*--------------------------------------------------------------------------------------------*/
+/** @packageDocumentation
+ * @module Curve
+ */
+import { CurvePrimitive } from "../CurvePrimitive";
+import { Segment1d } from "../../geometry3d/Segment1d";
+import { Transform } from "../../geometry3d/Transform";
+import { Geometry } from "../../Geometry";
+/**
+ * TransitionSpiral3d is a base class for multiple variants of spirals.
+ * * The menagerie of spiral types have 2 broad categories:
+ *   * IntegratedSpiral3d -- a spiral whose direct function for curvature versus distance must be integrated to determine x,y
+ *     * The IntegratedSpiral3d types are enumerated in `IntegratedSpiralTypes`
+ *   * DirectSpiral3d -- a spiral implemented with direct calculation of x,y from fractional position along the spiral.
+ *     * The direct spiral types are enumerated in the `DirectSpiralType`
+ * * The method set for CurvePrimitive support includes a `handleTransitionSpiral(g: TransitionSpiral3d)` which receives all the spiral types.
+ * * The spiral class may impose expectations that its inflection is at the origin, with tangent along the x axis.
+ *   * This is generally necessary for direct spirals.
+ *   * This is not necessary for integrated spirals.
+ * @public
+ */
+export class TransitionSpiral3d extends CurvePrimitive {
+    constructor(spiralType, localToWorld, activeFractionInterval, designProperties) {
+        super();
+        this._spiralType = spiralType ? spiralType : "unknownSpiralType";
+        this._designProperties = designProperties;
+        this._localToWorld = localToWorld;
+        this._activeFractionInterval = activeFractionInterval ? activeFractionInterval : Segment1d.create(0, 1);
+    }
+    /** Return (reference to) the active portion of the reference spiral. */
+    get activeFractionInterval() { return this._activeFractionInterval; }
+    /** (reference to) placement transform. */
+    get localToWorld() { return this._localToWorld; }
+    get spiralType() { return this._spiralType; }
+    /** Return 1/r with convention that if true zero is given as radius it represents infinite radius (0 curvature, straight line) */
+    static radiusToCurvature(radius) { return (radius === 0.0) ? 0.0 : 1.0 / radius; }
+    /** Return 1/k with convention that if near-zero is given as curvature, its infinite radius is returned as 0 */
+    static curvatureToRadius(curvature) {
+        if (Math.abs(curvature) < Geometry.smallAngleRadians)
+            return 0.0;
+        return 1.0 / curvature;
+    }
+    /** Return the average of the start and end curvatures. */
+    static averageCurvature(radiusLimits) {
+        return 0.5 * (TransitionSpiral3d.radiusToCurvature(radiusLimits.x0) + TransitionSpiral3d.radiusToCurvature(radiusLimits.x1));
+    }
+    /**
+     * Given two radii (or zeros for 0 curvature) return the average curvature
+     * @param r0 start radius, or 0 for line
+     * @param r1 end radius, or 0 for line
+     */
+    static averageCurvatureR0R1(r0, r1) {
+        return 0.5 * (TransitionSpiral3d.radiusToCurvature(r0) + TransitionSpiral3d.radiusToCurvature(r1));
+    }
+    /**
+     * Given two radii (or zeros for 0 curvature) return the average curvature
+     * @param r0 start radius, or 0 for line
+     * @param r1 end radius, or 0 for line
+     */
+    static interpolateCurvatureR0R1(r0, fraction, r1) {
+        return Geometry.interpolate(TransitionSpiral3d.radiusToCurvature(r0), fraction, TransitionSpiral3d.radiusToCurvature(r1));
+    }
+    /** Return the arc length of a transition spiral with given sweep and radius pair. */
+    static radiusRadiusSweepRadiansToArcLength(radius0, radius1, sweepRadians) {
+        return Math.abs(sweepRadians / TransitionSpiral3d.averageCurvatureR0R1(radius0, radius1));
+    }
+    /** Return the turn angle for spiral of given length between two radii */
+    static radiusRadiusLengthToSweepRadians(radius0, radius1, arcLength) {
+        return TransitionSpiral3d.averageCurvatureR0R1(radius0, radius1) * arcLength;
+    }
+    /** Return the end radius for spiral of given start radius, length, and turn angle. */
+    static radius0LengthSweepRadiansToRadius1(radius0, arcLength, sweepRadians) {
+        return TransitionSpiral3d.curvatureToRadius((2.0 * sweepRadians / arcLength) - TransitionSpiral3d.radiusToCurvature(radius0));
+    }
+    /** Return the start radius for spiral of given end radius, length, and turn angle. */
+    static radius1LengthSweepRadiansToRadius0(radius1, arcLength, sweepRadians) {
+        return TransitionSpiral3d.curvatureToRadius((2.0 * sweepRadians / arcLength) - TransitionSpiral3d.radiusToCurvature(radius1));
+    }
+    /** Return the original defining properties (if any) saved by the constructor. */
+    get designProperties() { return this._designProperties; }
+    /**
+     * * If transformA is rigid with uniform scale, apply the rigid part of transformA to the localToWorld transform and return the scale and rigid separation.
+     * * If not rigid, do nothing and return undefined.
+     * * Also apply the scale factor to the designProperties.
+     * @param transformA
+     */
+    applyRigidPartOfTransform(transformA) {
+        const rigidData = transformA.matrix.factorRigidWithSignedScale();
+        if (rigidData !== undefined) {
+            // [sQ a][R b] = [sQ*R sQb+a]
+            // but we save it as [Q*R sQb+a] with spiral data scaled by s.
+            const transformC0 = transformA.multiplyTransformTransform(this.localToWorld);
+            // BUT pull the scale part out of the matrix ...
+            const matrixC = rigidData.rigidAxes.multiplyMatrixMatrix(this.localToWorld.matrix);
+            this._localToWorld = Transform.createOriginAndMatrix(transformC0.origin, matrixC);
+            if (this.designProperties)
+                this.designProperties.applyScaleFactor(rigidData.scale);
+            return rigidData;
+        }
+        return undefined;
+    }
+}
+//# sourceMappingURL=TransitionSpiral3d.js.map

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

@@ -0,0 +1 @@
+{"version":3,"file":"TransitionSpiral3d.js","sourceRoot":"","sources":["../../../../src/curve/spiral/TransitionSpiral3d.ts"],"names":[],"mappings":"AAAA;;;+FAG+F;AAE/F;;GAEG;AAEH,OAAO,EAAE,cAAc,EAAE,MAAM,mBAAmB,CAAC;AAEnD,OAAO,EAAE,SAAS,EAAE,MAAM,4BAA4B,CAAC;AACvD,OAAO,EAAE,SAAS,EAAE,MAAM,4BAA4B,CAAC;AAEvD,OAAO,EAAE,QAAQ,EAAE,MAAM,gBAAgB,CAAC;AAmC1C;;;;;;;;;;;;GAYG;AACH,MAAM,OAAgB,kBAAmB,SAAQ,cAAc;IAmB7D,YAAsB,UAA8B,EAAE,YAAuB,EAAE,sBAA6C,EAAE,gBAA6D;QACzL,KAAK,EAAE,CAAC;QACR,IAAI,CAAC,WAAW,GAAG,UAAU,CAAC,CAAC,CAAC,UAAU,CAAC,CAAC,CAAC,mBAAmB,CAAC;QACjE,IAAI,CAAC,iBAAiB,GAAG,gBAAgB,CAAC;QAC1C,IAAI,CAAC,aAAa,GAAG,YAAY,CAAC;QAClC,IAAI,CAAC,uBAAuB,GAAG,sBAAsB,CAAC,CAAC,CAAC,sBAAsB,CAAC,CAAC,CAAC,SAAS,CAAC,MAAM,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;IAC1G,CAAC;IAfD,wEAAwE;IACxE,IAAW,sBAAsB,KAAgB,OAAO,IAAI,CAAC,uBAAuB,CAAC,CAAC,CAAC;IAKvF,0CAA0C;IAC1C,IAAW,YAAY,KAAgB,OAAO,IAAI,CAAC,aAAa,CAAC,CAAC,CAAC;IAUnE,IAAW,UAAU,KAAa,OAAO,IAAI,CAAC,WAAW,CAAC,CAAC,CAAC;IAC5D,iIAAiI;IAC1H,MAAM,CAAC,iBAAiB,CAAC,MAAc,IAAY,OAAO,CAAC,MAAM,KAAK,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,GAAG,MAAM,CAAC,CAAC,CAAC;IAEzG,+GAA+G;IACxG,MAAM,CAAC,iBAAiB,CAAC,SAAiB;QAC/C,IAAI,IAAI,CAAC,GAAG,CAAC,SAAS,CAAC,GAAG,QAAQ,CAAC,iBAAiB;YAClD,OAAO,GAAG,CAAC;QACb,OAAO,GAAG,GAAG,SAAS,CAAC;IACzB,CAAC;IAED,0DAA0D;IACnD,MAAM,CAAC,gBAAgB,CAAC,YAAuB;QACpD,OAAO,GAAG,GAAG,CAAC,kBAAkB,CAAC,iBAAiB,CAAC,YAAY,CAAC,EAAE,CAAC,GAAG,kBAAkB,CAAC,iBAAiB,CAAC,YAAY,CAAC,EAAE,CAAC,CAAC,CAAC;IAC/H,CAAC;IACD;;;;OAIG;IACI,MAAM,CAAC,oBAAoB,CAAC,EAAU,EAAE,EAAU;QACvD,OAAO,GAAG,GAAG,CAAC,kBAAkB,CAAC,iBAAiB,CAAC,EAAE,CAAC,GAAG,kBAAkB,CAAC,iBAAiB,CAAC,EAAE,CAAC,CAAC,CAAC;IACrG,CAAC;IACD;;;;OAIG;IACI,MAAM,CAAC,wBAAwB,CAAC,EAAU,EAAE,QAAgB,EAAE,EAAU;QAC7E,OAAO,QAAQ,CAAC,WAAW,CAAC,kBAAkB,CAAC,iBAAiB,CAAC,EAAE,CAAC,EAAE,QAAQ,EAAE,kBAAkB,CAAC,iBAAiB,CAAC,EAAE,CAAC,CAAC,CAAC;IAC5H,CAAC;IAED,qFAAqF;IAC9E,MAAM,CAAC,mCAAmC,CAAC,OAAe,EAAE,OAAe,EAAE,YAAoB;QACtG,OAAO,IAAI,CAAC,GAAG,CAAC,YAAY,GAAG,kBAAkB,CAAC,oBAAoB,CAAC,OAAO,EAAE,OAAO,CAAC,CAAC,CAAC;IAC5F,CAAC;IAED,yEAAyE;IAClE,MAAM,CAAC,gCAAgC,CAAC,OAAe,EAAE,OAAe,EAAE,SAAiB;QAChG,OAAO,kBAAkB,CAAC,oBAAoB,CAAC,OAAO,EAAE,OAAO,CAAC,GAAG,SAAS,CAAC;IAC/E,CAAC;IAED,sFAAsF;IAC/E,MAAM,CAAC,kCAAkC,CAAC,OAAe,EAAE,SAAiB,EAAE,YAAoB;QACvG,OAAO,kBAAkB,CAAC,iBAAiB,CAAC,CAAC,GAAG,GAAG,YAAY,GAAG,SAAS,CAAC,GAAG,kBAAkB,CAAC,iBAAiB,CAAC,OAAO,CAAC,CAAC,CAAC;IAChI,CAAC;IACD,sFAAsF;IAC/E,MAAM,CAAC,kCAAkC,CAAC,OAAe,EAAE,SAAiB,EAAE,YAAoB;QACvG,OAAO,kBAAkB,CAAC,iBAAiB,CAAC,CAAC,GAAG,GAAG,YAAY,GAAG,SAAS,CAAC,GAAG,kBAAkB,CAAC,iBAAiB,CAAC,OAAO,CAAC,CAAC,CAAC;IAChI,CAAC;IACD,iFAAiF;IACjF,IAAW,gBAAgB,KAAkD,OAAO,IAAI,CAAC,iBAAiB,CAAC,CAAC,CAAC;IAC7G;;;;;OAKG;IACO,yBAAyB,CAAC,UAAqB;QACvD,MAAM,SAAS,GAAG,UAAU,CAAC,MAAM,CAAC,0BAA0B,EAAE,CAAC;QACjE,IAAI,SAAS,KAAK,SAAS,EAAE;YAC3B,6BAA6B;YAC7B,8DAA8D;YAC9D,MAAM,WAAW,GAAG,UAAU,CAAC,0BAA0B,CAAC,IAAI,CAAC,YAAY,CAAC,CAAC;YAC7E,gDAAgD;YAChD,MAAM,OAAO,GAAG,SAAS,CAAC,SAAS,CAAC,oBAAoB,CAAC,IAAI,CAAC,YAAY,CAAC,MAAM,CAAC,CAAC;YACnF,IAAI,CAAC,aAAa,GAAG,SAAS,CAAC,qBAAqB,CAAC,WAAW,CAAC,MAAM,EAAE,OAAO,CAAC,CAAC;YAClF,IAAI,IAAI,CAAC,gBAAgB;gBACvB,IAAI,CAAC,gBAAgB,CAAC,gBAAgB,CAAC,SAAS,CAAC,KAAK,CAAC,CAAC;YAE1D,OAAO,SAAS,CAAC;SAClB;QACD,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 { CurvePrimitive } from \"../CurvePrimitive\";\r\nimport { Matrix3d } from \"../../geometry3d/Matrix3d\";\r\nimport { Segment1d } from \"../../geometry3d/Segment1d\";\r\nimport { Transform } from \"../../geometry3d/Transform\";\r\nimport { TransitionConditionalProperties } from \"./TransitionConditionalProperties\";\r\nimport { Geometry } from \"../../Geometry\";\r\nimport { LineString3d } from \"../LineString3d\";\r\n/**\r\n * This is the set of valid type names for \"integrated\" spirals\r\n * * Behavior is expressed by a `NormalizedTransition` snap function.\r\n * * The snap function varies smoothly from f(0)=0 to f(1)=1\r\n * * The various snap functions are:\r\n *   * clothoid: linear\r\n *   * biquadratic: 2 quadratics pieced together, joining with 1st derivative continuity at f(0.) = 0.5, with zero slope f'(0)=0 and f'(1)= 0\r\n *   * bloss: A single cubic with zero slope at 0 and 1\r\n *   * cosine: half of a cosine wave, centered around 0.5\r\n *   * sine: full period of a sine wave added to the line f(u)=u\r\n * *\r\n * @public\r\n */\r\nexport type IntegratedSpiralTypeName = \"clothoid\" | \"bloss\" | \"biquadratic\" | \"cosine\" | \"sine\";\r\n\r\n/**\r\n * This is the set of valid type names for \"direct\" spirals.\r\n * \"Direct\" spirals can evaluate fractionToPoint by direct equations, i.e. not requiring the numeric integrations in \"Integrated\" spiral types.\r\n * @public\r\n */\r\nexport type DirectSpiralTypeName =\r\n   \"JapaneseCubic\"  // 1 term from each of the X,Y clothoid series expansions:  y = x^3 / (6RL)\r\n  | \"Arema\"       // 2 terms from each of the X,Y clothoid series expansions.  Identical to ChineseCubic!\r\n  | \"ChineseCubic\"  // Identical to Arema!\r\n  | \"HalfCosine\"  // high continuity cosine variation from quadratic.\r\n  | \"AustralianRailCorp\" // cubic with high accuracy distance series\r\n  | \"WesternAustralian\"  // simple cubic -- 2 terms of x series, 1 term of y series.\r\n  | \"Czech\"  // simple cubic with two term distance approximation\r\n  | \"MXCubicAlongArc\"  // x obtained from fractional distance via 2-terms from series, y = x^3/ (6RL)\r\n  | \"Polish\"\r\n  | \"Italian\"\r\n  ;\r\n\r\n/**\r\n * TransitionSpiral3d is a base class for multiple variants of spirals.\r\n * * The menagerie of spiral types have 2 broad categories:\r\n *   * IntegratedSpiral3d -- a spiral whose direct function for curvature versus distance must be integrated to determine x,y\r\n *     * The IntegratedSpiral3d types are enumerated in `IntegratedSpiralTypes`\r\n *   * DirectSpiral3d -- a spiral implemented with direct calculation of x,y from fractional position along the spiral.\r\n *     * The direct spiral types are enumerated in the `DirectSpiralType`\r\n * * The method set for CurvePrimitive support includes a `handleTransitionSpiral(g: TransitionSpiral3d)` which receives all the spiral types.\r\n * * The spiral class may impose expectations that its inflection is at the origin, with tangent along the x axis.\r\n *   * This is generally necessary for direct spirals.\r\n *   * This is not necessary for integrated spirals.\r\n * @public\r\n */\r\nexport abstract class TransitionSpiral3d extends CurvePrimitive {\r\n  /** string name of spiral type */\r\n  protected _spiralType: string;\r\n  /** Original defining properties. */\r\n  protected _designProperties: TransitionConditionalProperties | undefined;\r\n\r\n  /** Fractional interval for the \"active\" part of a containing spiral.\r\n   * (The radius, angle, and length conditions define a complete spiral, and some portion of it is \"active\")\r\n   */\r\n  protected _activeFractionInterval: Segment1d;\r\n  /** Return (reference to) the active portion of the reference spiral. */\r\n  public get activeFractionInterval(): Segment1d { return this._activeFractionInterval; }\r\n  /** strokes in the active portion */\r\n  public abstract get activeStrokes(): LineString3d;\r\n  /** Placement transform */\r\n  protected _localToWorld: Transform;\r\n  /** (reference to) placement transform. */\r\n  public get localToWorld(): Transform { return this._localToWorld; }\r\n\r\n  protected constructor(spiralType: string | undefined, localToWorld: Transform, activeFractionInterval: Segment1d | undefined, designProperties: TransitionConditionalProperties | undefined) {\r\n    super();\r\n    this._spiralType = spiralType ? spiralType : \"unknownSpiralType\";\r\n    this._designProperties = designProperties;\r\n    this._localToWorld = localToWorld;\r\n    this._activeFractionInterval = activeFractionInterval ? activeFractionInterval : Segment1d.create(0, 1);\r\n  }\r\n\r\n  public get spiralType(): string { return this._spiralType; }\r\n  /** Return 1/r with convention that if true zero is given as radius it represents infinite radius (0 curvature, straight line) */\r\n  public static radiusToCurvature(radius: number): number { return (radius === 0.0) ? 0.0 : 1.0 / radius; }\r\n\r\n  /** Return 1/k with convention that if near-zero is given as curvature, its infinite radius is returned as 0 */\r\n  public static curvatureToRadius(curvature: number): number {\r\n    if (Math.abs(curvature) < Geometry.smallAngleRadians)\r\n      return 0.0;\r\n    return 1.0 / curvature;\r\n  }\r\n\r\n  /** Return the average of the start and end curvatures. */\r\n  public static averageCurvature(radiusLimits: Segment1d): number {\r\n    return 0.5 * (TransitionSpiral3d.radiusToCurvature(radiusLimits.x0) + TransitionSpiral3d.radiusToCurvature(radiusLimits.x1));\r\n  }\r\n  /**\r\n   * Given two radii (or zeros for 0 curvature) return the average curvature\r\n   * @param r0 start radius, or 0 for line\r\n   * @param r1 end radius, or 0 for line\r\n   */\r\n  public static averageCurvatureR0R1(r0: number, r1: number): number {\r\n    return 0.5 * (TransitionSpiral3d.radiusToCurvature(r0) + TransitionSpiral3d.radiusToCurvature(r1));\r\n  }\r\n  /**\r\n   * Given two radii (or zeros for 0 curvature) return the average curvature\r\n   * @param r0 start radius, or 0 for line\r\n   * @param r1 end radius, or 0 for line\r\n   */\r\n  public static interpolateCurvatureR0R1(r0: number, fraction: number, r1: number): number {\r\n    return Geometry.interpolate(TransitionSpiral3d.radiusToCurvature(r0), fraction, TransitionSpiral3d.radiusToCurvature(r1));\r\n  }\r\n\r\n  /** Return the arc length of a transition spiral with given sweep and radius pair. */\r\n  public static radiusRadiusSweepRadiansToArcLength(radius0: number, radius1: number, sweepRadians: number): number {\r\n    return Math.abs(sweepRadians / TransitionSpiral3d.averageCurvatureR0R1(radius0, radius1));\r\n  }\r\n\r\n  /** Return the turn angle for spiral of given length between two radii */\r\n  public static radiusRadiusLengthToSweepRadians(radius0: number, radius1: number, arcLength: number): number {\r\n    return TransitionSpiral3d.averageCurvatureR0R1(radius0, radius1) * arcLength;\r\n  }\r\n\r\n  /** Return the end radius for spiral of given start radius, length, and turn angle. */\r\n  public static radius0LengthSweepRadiansToRadius1(radius0: number, arcLength: number, sweepRadians: number) {\r\n    return TransitionSpiral3d.curvatureToRadius((2.0 * sweepRadians / arcLength) - TransitionSpiral3d.radiusToCurvature(radius0));\r\n  }\r\n  /** Return the start radius for spiral of given end radius, length, and turn angle. */\r\n  public static radius1LengthSweepRadiansToRadius0(radius1: number, arcLength: number, sweepRadians: number) {\r\n    return TransitionSpiral3d.curvatureToRadius((2.0 * sweepRadians / arcLength) - TransitionSpiral3d.radiusToCurvature(radius1));\r\n  }\r\n  /** Return the original defining properties (if any) saved by the constructor. */\r\n  public get designProperties(): TransitionConditionalProperties | undefined { return this._designProperties; }\r\n  /**\r\n   * * If transformA is rigid with uniform scale, apply the rigid part of transformA to the localToWorld transform and return the scale and rigid separation.\r\n   * * If not rigid, do nothing and return undefined.\r\n   * * Also apply the scale factor to the designProperties.\r\n   * @param transformA\r\n   */\r\n  protected applyRigidPartOfTransform(transformA: Transform): { rigidAxes: Matrix3d, scale: number } | undefined {\r\n    const rigidData = transformA.matrix.factorRigidWithSignedScale();\r\n    if (rigidData !== undefined) {\r\n      // [sQ a][R b] = [sQ*R sQb+a]\r\n      // but we save it as [Q*R sQb+a] with spiral data scaled by s.\r\n      const transformC0 = transformA.multiplyTransformTransform(this.localToWorld);\r\n      // BUT pull the scale part out of the matrix ...\r\n      const matrixC = rigidData.rigidAxes.multiplyMatrixMatrix(this.localToWorld.matrix);\r\n      this._localToWorld = Transform.createOriginAndMatrix(transformC0.origin, matrixC);\r\n      if (this.designProperties)\r\n        this.designProperties.applyScaleFactor(rigidData.scale);\r\n\r\n      return rigidData;\r\n    }\r\n    return undefined;\r\n  }\r\n}\r\n"]}

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

@@ -0,0 +1,80 @@
+import { Plane3dByOriginAndVectors } from "../../geometry3d/Plane3dByOriginAndVectors";
+import { Point3d, Vector3d } from "../../geometry3d/Point3dVector3d";
+import { Ray3d } from "../../geometry3d/Ray3d";
+/**
+ * XYCurveEvaluator is an abstract with methods for evaluating X and Y parts of a curve parameterized by a fraction.
+ * * The required methods call for independent X and Y evaluation.
+ * * Base class methods package those (multiple) calls into point, ray, and plane structures.
+ *    * A implementation that has evaluation substantial cost that can be shared among x,y parts or between
+ *       primary functions and derivatives might choose to implement the point and derivative methods directly.
+ * @internal
+ */
+export declare abstract class XYCurveEvaluator {
+    /** return a deep copy of the evaluator */
+    abstract clone(): XYCurveEvaluator;
+    /** test for near identical evaluator */
+    abstract isAlmostEqual(other: any): boolean;
+    /** Evaluate X at fractional position. */
+    abstract fractionToX(fraction: number): number;
+    /** Evaluate Y at fractional position. */
+    abstract fractionToY(fraction: number): number;
+    /** Evaluate derivative of X with respect to fraction at fractional position. */
+    abstract fractionToDX(fraction: number): number;
+    /** Evaluate derivative of Y with respect to fraction at fractional position. */
+    abstract fractionToDY(fraction: number): number;
+    /** Evaluate second derivative of X with respect to fraction at fractional position. */
+    abstract fractionToDDX(fraction: number): number;
+    /** Evaluate second derivative of Y with respect to fraction at fractional position. */
+    abstract fractionToDDY(fraction: number): number;
+    /** Evaluate both X and Y at fractional coordinate, return bundled as a point. */
+    /** Evaluate second derivative of X with respect to fraction at fractional position. */
+    abstract fractionToD3X(fraction: number): number;
+    /** Evaluate second derivative of Y with respect to fraction at fractional position. */
+    abstract fractionToD3Y(fraction: number): number;
+    /** Evaluate both X and Y at fractional coordinate, return bundled as a point. */
+    fractionToPoint(fraction: number, result?: Point3d): Point3d;
+    /** Evaluate both X and Y and their first derivatives at fractional coordinate, return bundled as origin and (non-unit) direction vector. */
+    fractionToPointAndDerivative(fraction: number, result?: Ray3d): Ray3d;
+    /** Evaluate both X and Y and their second derivatives at fractional coordinate, return bundled as origin and (non-unit) vectorU an vectorV. */
+    fractionToPointAnd2Derivatives(fraction: number, result?: Plane3dByOriginAndVectors): Plane3dByOriginAndVectors;
+    /**
+     * Return the magnitude of the tangent vector at fraction.
+     * @param fraction fractional position along the curve
+     */
+    fractionToTangentMagnitude(fraction: number): number;
+    /** Invert the fractionToX function for given X. */
+    abstract xToFraction(x: number): number | undefined;
+    /** Initialize class level work arrays for 5 point Gaussian Quadrature. */
+    protected static _gaussX: Float64Array;
+    protected static _gaussWeight: Float64Array;
+    protected static _gaussMapper: (xA: number, xB: number, arrayX: Float64Array, arrayW: Float64Array) => number;
+    static initWorkSpace(): void;
+    /**
+     * Integrate between nominal fractions with default gauss rule.
+     * * The caller is expected to choose nearby fractions so that the single gauss interval accuracy is good.
+     * @param fraction0
+     * @param fraction1
+     */
+    integrateDistanceBetweenFractions(fraction0: number, fraction1: number): number;
+    /**
+     * Inverse integrated distance
+     * @param fraction0 start of fraction interval
+     * @param fraction1 end of fraction interval
+     * @param distance0 distance at start
+     * @param distance1 distance at end
+     * @param targetDistance intermediate distance.
+     */
+    inverseDistanceFraction(fraction0: number, fraction1: number, distance0: number, distance1: number, targetDistance: number): number | undefined;
+    /**
+     *
+     * @param fraction fractional position along x axis
+     * @param xy xy coordinates of point on the curve
+     * @param d1xy
+     * @param d2xy
+     * @param d3xy
+     */
+    fractionToPointAnd3Derivatives(fraction: number, xy: Point3d, d1xy?: Vector3d, d2xy?: Vector3d, d3xy?: Vector3d): void;
+    /** Apply a uniform scale around the origin. */
+    abstract scaleInPlace(scaleFactor: number): void;
+}
+//# sourceMappingURL=XYCurveEvaluator.d.ts.map

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

@@ -0,0 +1 @@
+{"version":3,"file":"XYCurveEvaluator.d.ts","sourceRoot":"","sources":["../../../../src/curve/spiral/XYCurveEvaluator.ts"],"names":[],"mappings":"AASA,OAAO,EAAE,yBAAyB,EAAE,MAAM,4CAA4C,CAAC;AACvF,OAAO,EAAE,OAAO,EAAE,QAAQ,EAAE,MAAM,kCAAkC,CAAC;AACrE,OAAO,EAAE,KAAK,EAAE,MAAM,wBAAwB,CAAC;AAG/C;;;;;;;GAOG;AACH,8BAAsB,gBAAgB;IACpC,0CAA0C;aAC1B,KAAK,IAAI,gBAAgB;IACzC,wCAAwC;aACxB,aAAa,CAAC,KAAK,EAAE,GAAG,GAAG,OAAO;IAClD,yCAAyC;aACzB,WAAW,CAAC,QAAQ,EAAE,MAAM,GAAG,MAAM;IACrD,yCAAyC;aACzB,WAAW,CAAC,QAAQ,EAAE,MAAM,GAAG,MAAM;IACrD,gFAAgF;aAChE,YAAY,CAAC,QAAQ,EAAE,MAAM,GAAG,MAAM;IACtD,gFAAgF;aAChE,YAAY,CAAC,QAAQ,EAAE,MAAM,GAAG,MAAM;IACtD,uFAAuF;aACvE,aAAa,CAAC,QAAQ,EAAE,MAAM,GAAG,MAAM;IACvD,uFAAuF;aACvE,aAAa,CAAC,QAAQ,EAAE,MAAM,GAAG,MAAM;IACvD,iFAAiF;IACjF,uFAAuF;aACvE,aAAa,CAAC,QAAQ,EAAE,MAAM,GAAG,MAAM;IACvD,uFAAuF;aACvE,aAAa,CAAC,QAAQ,EAAE,MAAM,GAAG,MAAM;IACvD,iFAAiF;IAC1E,eAAe,CAAC,QAAQ,EAAE,MAAM,EAAE,MAAM,CAAC,EAAE,OAAO,GAAG,OAAO;IAGnE,4IAA4I;IACrI,4BAA4B,CAAC,QAAQ,EAAE,MAAM,EAAE,MAAM,CAAC,EAAE,KAAK,GAAG,KAAK;IAK5E,+IAA+I;IACxI,8BAA8B,CAAC,QAAQ,EAAE,MAAM,EAAE,MAAM,CAAC,EAAE,yBAAyB,GAAG,yBAAyB;IAOtH;;;OAGG;IACI,0BAA0B,CAAC,QAAQ,EAAE,MAAM,GAAG,MAAM;IAM3D,mDAAmD;aAEnC,WAAW,CAAC,CAAC,EAAE,MAAM,GAAG,MAAM,GAAG,SAAS;IAC1D,0EAA0E;IAG1E,SAAS,CAAC,MAAM,CAAC,OAAO,EAAE,YAAY,CAAC;IACvC,SAAS,CAAC,MAAM,CAAC,YAAY,EAAE,YAAY,CAAC;IAC5C,SAAS,CAAC,MAAM,CAAC,YAAY,EAAE,CAAC,EAAE,EAAE,MAAM,EAAE,EAAE,EAAE,MAAM,EAAE,MAAM,EAAE,YAAY,EAAE,MAAM,EAAE,YAAY,KAAK,MAAM,CAAC;WAChG,aAAa;IAK3B;;;;;OAKG;IACI,iCAAiC,CAAC,SAAS,EAAE,MAAM,EAAE,SAAS,EAAE,MAAM,GAAG,MAAM;IAUtF;;;;;;;OAOG;IACI,uBAAuB,CAAC,SAAS,EAAE,MAAM,EAAE,SAAS,EAAE,MAAM,EAAE,SAAS,EAAE,MAAM,EAAE,SAAS,EAAE,MAAM,EAAE,cAAc,EAAE,MAAM,GAAG,MAAM,GAAG,SAAS;IAatJ;;;;;;;OAOG;IACI,8BAA8B,CAAC,QAAQ,EAAE,MAAM,EAAE,EAAE,EAAE,OAAO,EAAE,IAAI,CAAC,EAAE,QAAQ,EAAE,IAAI,CAAC,EAAE,QAAQ,EAAE,IAAI,CAAC,EAAE,QAAQ;IAStH,+CAA+C;aAC/B,YAAY,CAAC,WAAW,EAAE,MAAM,GAAG,IAAI;CACxD"}