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

package/lib/curve/spiral/DirectSpiral3d.js

@@ -1,464 +0,0 @@
-"use strict";
-Object.defineProperty(exports, "__esModule", { value: true });
-exports.DirectSpiral3d = void 0;
-/*---------------------------------------------------------------------------------------------
-* 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
- */
-const Geometry_1 = require("../../Geometry");
-const Segment1d_1 = require("../../geometry3d/Segment1d");
-const LineString3d_1 = require("../LineString3d");
-const StrokeOptions_1 = require("../StrokeOptions");
-const ClothoidSeries_1 = require("./ClothoidSeries");
-const CzechSpiralEvaluator_1 = require("./CzechSpiralEvaluator");
-const DirectHalfCosineSpiralEvaluator_1 = require("./DirectHalfCosineSpiralEvaluator");
-const AustralianRailCorpXYEvaluator_1 = require("./AustralianRailCorpXYEvaluator");
-const TransitionSpiral3d_1 = require("./TransitionSpiral3d");
-const MXCubicAlongArcSpiralEvaluator_1 = require("./MXCubicAlongArcSpiralEvaluator");
-const PolishCubicSpiralEvaluator_1 = require("./PolishCubicSpiralEvaluator");
-/**
-* 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 create with these static methods:
-*   * createArema
-*   * createJapaneseCubic
-*   * createAustralianRail
-*   * createDirectHalfCosine
-*   * createChineseCubic
-*   * createCzechCubic
-*   * createPolishCubic
-*   * createItalian
-*   * createWesternAustralian
-* @public
-*/
-class DirectSpiral3d extends TransitionSpiral3d_1.TransitionSpiral3d {
-    // constructor demands radius1 and distance1 for nominal construction.
-    // caller is responsible for managing intervals of partial spiral
-    constructor(localToWorld, spiralType, originalProperties, nominalL1, nominalR1, activeFractionInterval, evaluator) {
-        super(spiralType, localToWorld, activeFractionInterval, originalProperties);
-        /** String name for schema properties */
-        this.curvePrimitiveType = "transitionSpiral";
-        this._nominalL1 = nominalL1;
-        this._nominalR1 = nominalR1;
-        this._evaluator = evaluator;
-        this._globalStrokes = LineString3d_1.LineString3d.create();
-        this._activeStrokes = LineString3d_1.LineString3d.create();
-        // initialize for compiler -- but this will be recomputed in refreshComputeProperties ...
-        this.refreshComputedProperties();
-    }
-    /** Return the internal stroked form of the (possibly partial) spiral   */
-    get activeStrokes() { return this._activeStrokes !== undefined ? this._activeStrokes : this._globalStrokes; }
-    /** Return the nominal end radius. */
-    get nominalR1() { return this._nominalR1; }
-    /** Return the nominal distance from inflection to endpoint. */
-    get nominalL1() { return this._nominalL1; }
-    /** Return the nominal end curvature */
-    get nominalCurvature1() { return TransitionSpiral3d_1.TransitionSpiral3d.radiusToCurvature(this._nominalR1); }
-    /** Return the low level evaluator
-     * @internal
-     */
-    get evaluator() { return this._evaluator; }
-    /**
-     * Compute stroke data in an interval.
-     * @param strokes strokes to clear and refill.
-     * @param fraction0 start fraction
-     * @param fraction1 end fraction
-     */
-    computeStrokes(strokes, fractionA, fractionB, numInterval) {
-        if (numInterval < 1)
-            numInterval = 1;
-        strokes.clear();
-        strokes.ensureEmptyUVParams();
-        strokes.ensureEmptyFractions();
-        const distances = strokes.packedUVParams;
-        const nominalIntervalLength = Math.abs(fractionB - fractionA) * this._nominalL1;
-        for (let i = 0; i <= numInterval; i++) {
-            const fraction = Geometry_1.Geometry.interpolate(fractionA, i / numInterval, fractionB);
-            const nominalDistanceAlong = fraction * nominalIntervalLength;
-            strokes.packedPoints.pushXYZ(this._evaluator.fractionToX(fraction), this._evaluator.fractionToY(fraction), 0);
-            distances.pushXY(fraction, nominalDistanceAlong); // the second distance will be updated below
-        }
-        let fraction0 = distances.getXAtUncheckedPointIndex(0);
-        let trueDistance0 = distances.getYAtUncheckedPointIndex(0); // whatever was assigned as start distance is fine
-        let trueDistance1, fraction1;
-        for (let i = 1; i <= numInterval; i++) {
-            fraction1 = distances.getXAtUncheckedPointIndex(i);
-            trueDistance1 = trueDistance0 + this._evaluator.integrateDistanceBetweenFractions(fraction0, fraction1);
-            distances.setXYZAtCheckedPointIndex(i, fraction1, trueDistance1);
-            fraction0 = fraction1;
-            trueDistance0 = trueDistance1;
-        }
-    }
-    /** Recompute strokes */
-    refreshComputedProperties() {
-        const sweepRadians = this.nominalL1 / (2.0 * this.nominalR1);
-        const radiansStep = 0.02;
-        const numInterval = StrokeOptions_1.StrokeOptions.applyAngleTol(undefined, 4, sweepRadians, radiansStep);
-        this.computeStrokes(this._globalStrokes, 0, 1, numInterval);
-        const numActiveInterval = Math.ceil(this._activeFractionInterval.absoluteDelta() * numInterval);
-        this._activeStrokes = LineString3d_1.LineString3d.create();
-        this.computeStrokes(this._activeStrokes, this._activeFractionInterval.x0, this._activeFractionInterval.x1, numActiveInterval);
-    }
-    /**
-     * Create a spiral object which uses numXTerm terms from the clothoid X series and numYTerm from the clothoid Y series.
-     * @param numXTerm  number of terms to use from X series
-     * @param numYTerm number of terms to use from Y series
-     * @param localToWorld placement frame.  Inflection point is at origin, initial direction is along x axis.
-     * @param nominalL1 design distance from inflection to end point.
-     * @param nominalR1 design radius at end point.
-     * @param activeInterval active interval (as fractions of nominalL1 !!!)
-     */
-    static createTruncatedClothoid(spiralType, localToWorld, numXTerm, numYTerm, originalProperties, nominalL1, nominalR1, activeInterval) {
-        if (numXTerm < 1)
-            numXTerm = 1;
-        if (numYTerm < 1)
-            numYTerm = 1;
-        const evaluator = new ClothoidSeries_1.ClothoidSeriesRLEvaluator(nominalL1, 1.0 / (2.0 * nominalL1 * nominalR1), numXTerm, numYTerm);
-        return new DirectSpiral3d(localToWorld.clone(), spiralType, originalProperties, nominalL1, nominalR1, activeInterval ? activeInterval.clone() : Segment1d_1.Segment1d.create(0, 1), evaluator);
-    }
-    /**
-     * Create an Japanese spiral clothoid approximation
-     *   * X is 1 terms of the clothoid series as a function of nominal distance along.
-     *   * Y is 1 terms f the clothoid series as a function of nominal distance along.
-     *   * Remark: This is identical to the ChineseCubic
-     * @param localToWorld axes with inflection at origin, tangent along x axis
-     * @param nominalL1 nominal length as used in series LR terms.
-     * @param nominalR1 nominal final radius as used in series LR terms
-     * @param activeInterval fractional interval with (0, nominalL1) range for nominal distance along
-     */
-    static createJapaneseCubic(localToWorld, nominalL1, nominalR1, activeInterval) {
-        return this.createTruncatedClothoid("JapaneseCubic", localToWorld, 1, 1, undefined, nominalL1, nominalR1, activeInterval);
-    }
-    /**
-     * Create a czech cubic.
-     * This is y= m*x^3 with
-     * * x any point on the x axis
-     * * `fraction` along the spiral goes to `x = fraction * L`
-     * * m is gamma / (6RL)
-     *    * 1/(6RL) is the leading term of the sine series.
-     *    * `gamma = 2R/sqrt (4RR-LL)` pushes y up a little bit to simulate the lost series terms.
-     * @param localToWorld
-     * @param nominalLx nominal length along x axis
-     * @param nominalR1
-     * @param activeInterval
-     */
-    static createCzechCubic(localToWorld, nominalLx, nominalR1, activeInterval) {
-        const evaluator = CzechSpiralEvaluator_1.CzechSpiralEvaluator.create(nominalLx, nominalR1);
-        if (evaluator === undefined)
-            return undefined;
-        return new DirectSpiral3d(localToWorld.clone(), "Czech", undefined, nominalLx, nominalR1, activeInterval ? activeInterval.clone() : Segment1d_1.Segment1d.create(0, 1), evaluator);
-    }
-    /**
-     * Create an italian spiral
-     * This is y= m*x^3 with
-     * * x any point on the x axis
-     * * `fraction` along the spiral goes to `x = fraction * L`
-     * * m is gamma / (6RL)
-     *    * 1/(6RL) is the leading term of the sine series.
-     *    * `gamma = 2R/sqrt (4RR-LL)` pushes y up a little bit to simulate the lost series terms.
-     * * L in gamma and m is the
-     * @param localToWorld
-     * @param nominalL1 nominal length along the spiral
-     * @param nominalR1
-     * @param activeInterval
-     */
-    static createItalian(localToWorld, nominalL1, nominalR1, activeInterval) {
-        const evaluator = CzechSpiralEvaluator_1.ItalianSpiralEvaluator.create(nominalL1, nominalR1);
-        if (evaluator === undefined)
-            return undefined;
-        return new DirectSpiral3d(localToWorld.clone(), "Italian", undefined, nominalL1, nominalR1, activeInterval ? activeInterval.clone() : Segment1d_1.Segment1d.create(0, 1), evaluator);
-    }
-    /**
-     * Create an MX Cubic whose nominal length is close to along the curve.
-     * This is y= m*x^3 with
-     * * m is 1/ (6RL1)
-     *    * 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))
-     * @param localToWorld
-     * @param nominalL1
-     * @param nominalR1
-     * @param activeInterval
-     */
-    static createMXCubicAlongArc(localToWorld, nominalL1, nominalR1, activeInterval) {
-        const evaluator = MXCubicAlongArcSpiralEvaluator_1.MXCubicAlongArcEvaluator.create(nominalL1, nominalR1);
-        if (evaluator === undefined)
-            return undefined;
-        return new DirectSpiral3d(localToWorld.clone(), "MXCubicAlongArc", undefined, nominalL1, nominalR1, activeInterval ? activeInterval.clone() : Segment1d_1.Segment1d.create(0, 1), evaluator);
-    }
-    /**
-     * Create a polish cubic
-     * This is y= m*x^3 with
-     * * m is 1/ (6RL)
-     *    * 1/(6RL) is the leading term of the sine series.
-     * * L is nominal length
-     * * R is nominal end radius.
-     * * x ranges up to the x axis distance for which the polish distance series produces f(x)=L
-     * * The support class PolishCubicEvaluator has static methods for the distance series and its inversion.
-     */
-    static createPolishCubic(localToWorld, nominalL1, nominalR1, activeInterval) {
-        const evaluator = PolishCubicSpiralEvaluator_1.PolishCubicEvaluator.create(nominalL1, nominalR1);
-        if (evaluator === undefined)
-            return undefined;
-        return new DirectSpiral3d(localToWorld.clone(), "PolishCubic", undefined, nominalL1, nominalR1, activeInterval ? activeInterval.clone() : Segment1d_1.Segment1d.create(0, 1), evaluator);
-    }
-    /**
-     * Create an AustralianRailCorp spiral
-     * This is y= m*x^3 with
-     * * x any point on the x axis
-     * * `fraction` along the spiral goes to `x = fraction * L`
-     * * m is gamma / (6RL)
-     *    * 1/(6RL) is the leading term of the sine series.
-     *    * `gamma = 2R/sqrt (4RR-LL)` pushes y up a little bit to simulate the lost series terms.
-     * @param localToWorld
-     * @param nominalL1
-     * @param nominalR1
-     * @param activeInterval
-     */
-    static createAustralianRail(localToWorld, nominalL1, nominalR1, activeInterval) {
-        const evaluator = AustralianRailCorpXYEvaluator_1.AustralianRailCorpXYEvaluator.create(nominalL1, nominalR1);
-        if (evaluator === undefined)
-            return undefined;
-        return new DirectSpiral3d(localToWorld.clone(), "AustralianRailCorp", undefined, nominalL1, nominalR1, activeInterval ? activeInterval.clone() : Segment1d_1.Segment1d.create(0, 1), evaluator);
-    }
-    static createDirectHalfCosine(localToWorld, nominalL1, nominalR1, activeInterval) {
-        return new this(localToWorld, "HalfCosine", undefined, nominalL1, nominalR1, activeInterval, new DirectHalfCosineSpiralEvaluator_1.DirectHalfCosineSpiralEvaluator(nominalL1, nominalR1));
-    }
-    /**
-     * Create an Arema spiral clothoid approximation
-     *   * X is 2 terms of the clothoid series as a function of nominal distance along
-     *   * Y is 2 terms f the clothoid series as a function of nominal distance along
-     *   * Remark: This is identical to the ChineseCubic
-     * @param localToWorld axes with inflection at origin, tangent along x axis
-     * @param nominalL1 nominal length as used in series LR terms.
-     * @param nominalR1 nominal final radius as used in series LR terms
-     * @param activeInterval fractional interval with (0, nominalL1) range for nominal distance along
-     */
-    static createArema(localToWorld, nominalL1, nominalR1, activeInterval) {
-        return this.createTruncatedClothoid("Arema", localToWorld, 2, 2, undefined, nominalL1, nominalR1, activeInterval);
-    }
-    /**
-     * Create a Chinese clothoid approximation
-     *   * X is 2 terms of the clothoid series as a function of nominal distance along
-     *   * Y is 2 terms f the clothoid series as a function of nominal distance along
-     *   * Remark: This is identical to the Arema spiral
-     * @param localToWorld axes with inflection at origin, tangent along x axis
-     * @param nominalL1 nominal length as used in series LR terms.
-     * @param nominalR1 nominal final radius as used in series LR terms
-     * @param activeInterval fractional interval with (0, nominalL1) range for nominal distance along
-     */
-    static createChineseCubic(localToWorld, nominalL1, nominalR1, activeInterval) {
-        return this.createTruncatedClothoid("ChineseCubic", localToWorld, 2, 2, undefined, nominalL1, nominalR1, activeInterval);
-    }
-    /**
-     * Create a Western Australian direct spiral.
-     *   * X is 2 terms of the clothoid series as a function of distance along
-     *   * Y is 1 term (cubic in nominal distance along)
-     * @param localToWorld axes with inflection at origin, tangent along x axis
-     * @param nominalL1 nominal length as used in series LR terms.
-     * @param nominalR1 nominal final radius as used in series LR terms
-     * @param activeInterval fractional interval with (0, nominalL1) range for nominal distance along
-     */
-    static createWesternAustralian(localToWorld, nominalL1, nominalR1, activeInterval) {
-        return this.createTruncatedClothoid("WesternAustralian", localToWorld, 2, 1, undefined, nominalL1, nominalR1, activeInterval);
-    }
-    /**
-     * Create (if possible) a DirectSpiral3d, applying various strict conditions appropriate to the spiral type.
-     * The parameter list includes extraneous values in order to directly match IntegratedSpiral3d.create, which has greater flexibility about
-     *    mixtures of values.
-     * * IMPORTANT RESTRICTIONS
-     *   * Direct spirals must have the inflection at the origin of their coordinate system, aligned with the x axis.
-     *      * hence bearing0 = 0
-     *      * hence radius0 = 0
-     *   * bearing1 is ignored
-     *   * radius1 must be given.
-     *   * arcLength must be given,
-     * @param spiralType one of the types in `DirectSpiralTypeNames`
-     * @param radius0 radius (or 0 for tangent to line) at start.   Must be ZERO or UNDEFINED
-     * @param radius1 radius (or 0 for tangent to line) at end.
-     * @param bearing0 bearing, measured CCW from x axis at start.   Must be ZERO or UNDEFINED
-     * @param bearing1 bearing, measured CCW from x axis at end.    IGNORED.
-     * @param fractionInterval optional fractional interval for an "active" portion of the curve.   if omitted, the full [0,1] is used.
-     * @param localToWorld placement transform
-     */
-    static createFromLengthAndRadius(spiralType, radius0, radius1, bearing0, _bearing1, arcLength, activeInterval, localToWorld) {
-        if (bearing0 !== undefined && !bearing0.isAlmostZero)
-            return undefined;
-        if (radius0 !== undefined && !Geometry_1.Geometry.isSmallMetricDistance(radius0))
-            return undefined;
-        if (radius1 === undefined || Geometry_1.Geometry.isSmallMetricDistance(radius1))
-            return undefined;
-        if (arcLength === undefined)
-            return undefined;
-        if (Geometry_1.Geometry.equalStringNoCase(spiralType, "Arema"))
-            return this.createArema(localToWorld, arcLength, radius1, activeInterval);
-        if (Geometry_1.Geometry.equalStringNoCase(spiralType, "ChineseCubic"))
-            return this.createChineseCubic(localToWorld, arcLength, radius1, activeInterval);
-        if (Geometry_1.Geometry.equalStringNoCase(spiralType, "JapaneseCubic"))
-            return this.createJapaneseCubic(localToWorld, arcLength, radius1, activeInterval);
-        if (Geometry_1.Geometry.equalStringNoCase(spiralType, "HalfCosine"))
-            return this.createDirectHalfCosine(localToWorld, arcLength, radius1, activeInterval);
-        if (Geometry_1.Geometry.equalStringNoCase(spiralType, "Czech"))
-            return this.createCzechCubic(localToWorld, arcLength, radius1, activeInterval);
-        if (Geometry_1.Geometry.equalStringNoCase(spiralType, "Italian"))
-            return this.createItalian(localToWorld, arcLength, radius1, activeInterval);
-        if (Geometry_1.Geometry.equalStringNoCase(spiralType, "AustralianRailCorp"))
-            return this.createAustralianRail(localToWorld, arcLength, radius1, activeInterval);
-        if (Geometry_1.Geometry.equalStringNoCase(spiralType, "MXCubicAlongArc"))
-            return this.createMXCubicAlongArc(localToWorld, arcLength, radius1, activeInterval);
-        if (Geometry_1.Geometry.equalStringNoCase(spiralType, "WesternAustralian"))
-            return this.createWesternAustralian(localToWorld, arcLength, radius1, activeInterval);
-        if (Geometry_1.Geometry.equalStringNoCase(spiralType, "PolishCubic"))
-            return this.createPolishCubic(localToWorld, arcLength, radius1, activeInterval);
-        return undefined;
-    }
-    /** Deep clone of this spiral */
-    clone() {
-        var _a, _b;
-        return new DirectSpiral3d(this.localToWorld.clone(), this._spiralType, (_a = this.designProperties) === null || _a === void 0 ? void 0 : _a.clone(), this._nominalL1, this._nominalR1, (_b = this._activeFractionInterval) === null || _b === void 0 ? void 0 : _b.clone(), this._evaluator.clone());
-    }
-    /** apply `transform` to this spiral's local to world transform. */
-    tryTransformInPlace(transformA) {
-        const rigidData = this.applyRigidPartOfTransform(transformA);
-        if (rigidData !== undefined) {
-            this._nominalL1 *= rigidData.scale;
-            this._nominalR1 *= rigidData.scale;
-            this.evaluator.scaleInPlace(rigidData.scale);
-        }
-        this.refreshComputedProperties();
-        return true;
-    }
-    /** Clone with a transform applied  */
-    cloneTransformed(transform) {
-        const result = this.clone();
-        result.tryTransformInPlace(transform); // ok, we're confident it will always work.
-        return result;
-    }
-    /** Return the spiral start point. */
-    startPoint() { return this.activeStrokes.startPoint(); }
-    /** return the spiral end point. */
-    endPoint() { return this.activeStrokes.endPoint(); }
-    /** test if the local to world transform places the spiral xy plane into `plane` */
-    isInPlane(plane) {
-        return plane.isPointInPlane(this.localToWorld.origin)
-            && Geometry_1.Geometry.isSameCoordinate(0.0, this.localToWorld.matrix.dotColumnX(plane.getNormalRef()))
-            && Geometry_1.Geometry.isSameCoordinate(0.0, this.localToWorld.matrix.dotColumnY(plane.getNormalRef()));
-    }
-    /** Return quick length of the spiral.
-     * The tangent vector of a true clothoid is length 1 everywhere, so simple proportion of nominalL1 is a good approximation.
-     */
-    quickLength() {
-        const distanceData = this._globalStrokes.packedUVParams;
-        const n = distanceData.length;
-        return distanceData.getYAtUncheckedPointIndex(n - 1);
-    }
-    /** Return length of the spiral.
-     * * True length is stored at back of uvParams . . .
-     */
-    curveLength() { return this.quickLength(); }
-    /** Test if `other` is an instance of `TransitionSpiral3d` */
-    isSameGeometryClass(other) { return other instanceof DirectSpiral3d; }
-    /** Add strokes from this spiral to `dest`.
-     * * Linestrings will usually stroke as just their points.
-     * * If maxEdgeLength is given, this will sub-stroke within the linestring -- not what we want.
-     */
-    emitStrokes(dest, options) { this.activeStrokes.emitStrokes(dest, options); }
-    /** emit stroke fragments to `dest` handler. */
-    emitStrokableParts(dest, options) {
-        const n = this.computeStrokeCountForOptions(options);
-        const activeStrokes = this.activeStrokes;
-        dest.startParentCurvePrimitive(this);
-        if (n <= activeStrokes.numPoints()) {
-            // this.activeStrokes.emitStrokableParts(dest, options);
-            dest.announceIntervalForUniformStepStrokes(this, 2 * activeStrokes.numPoints(), 0.0, 1.0);
-        }
-        else {
-            dest.announceIntervalForUniformStepStrokes(this, n, 0.0, 1.0);
-        }
-        dest.endParentCurvePrimitive(this);
-    }
-    /**
-     * return the stroke count required for given options.
-     * @param options StrokeOptions that determine count
-     */
-    computeStrokeCountForOptions(options) {
-        let numStroke;
-        const nominalRadians = this._nominalL1 / (2.0 * this._nominalR1);
-        if (options) {
-            const rMin = Math.abs(this._nominalR1);
-            numStroke = options.applyTolerancesToArc(rMin, nominalRadians);
-            numStroke = options.applyMaxEdgeLength(numStroke, this.quickLength());
-            numStroke = options.applyMinStrokesPerPrimitive(numStroke);
-        }
-        else {
-            numStroke = StrokeOptions_1.StrokeOptions.applyAngleTol(undefined, 4, nominalRadians, 0.02);
-        }
-        numStroke = Math.ceil(this._activeFractionInterval.absoluteDelta() * numStroke);
-        return numStroke;
-    }
-    /** Reverse the active interval and active strokes.
-     * * Primary defining data remains unchanged !!!
-     */
-    reverseInPlace() {
-        this._activeFractionInterval.reverseInPlace();
-        if (this._activeStrokes === undefined)
-            this._activeStrokes = this._globalStrokes.clone();
-        this._activeStrokes.reverseInPlace();
-    }
-    /** Evaluate curve point with respect to fraction. */
-    fractionToPoint(activeFraction, result) {
-        const globalFraction = this._activeFractionInterval.fractionToPoint(activeFraction);
-        result = this._evaluator.fractionToPoint(globalFraction, result);
-        this.localToWorld.multiplyPoint3d(result, result);
-        return result;
-    }
-    /** Evaluate curve point and derivative with respect to fraction. */
-    fractionToPointAndDerivative(activeFraction, result) {
-        const globalFraction = this._activeFractionInterval.fractionToPoint(activeFraction);
-        result = this._evaluator.fractionToPointAndDerivative(globalFraction, result);
-        result.direction.scaleInPlace(this._activeFractionInterval.signedDelta());
-        result.transformInPlace(this.localToWorld);
-        return result;
-    }
-    /** Return a plane with
-     *
-     * * origin at fractional position along the curve
-     * * vectorU is the first derivative, i.e. tangent vector with length equal to the rate of change with respect to the fraction.
-     * * vectorV is the second derivative, i.e.derivative of vectorU.
-     */
-    fractionToPointAnd2Derivatives(activeFraction, result) {
-        const globalFraction = this._activeFractionInterval.fractionToPoint(activeFraction);
-        result = this._evaluator.fractionToPointAnd2Derivatives(globalFraction, result);
-        const a = this._activeFractionInterval.signedDelta();
-        result.vectorU.scaleInPlace(a);
-        result.vectorV.scaleInPlace(a * a);
-        result.transformInPlace(this.localToWorld);
-        return result;
-    }
-    /** Second step of double dispatch:  call `handler.handleTransitionSpiral(this)` */
-    dispatchToGeometryHandler(handler) {
-        return handler.handleTransitionSpiral(this);
-    }
-    /** extend the range by the strokes of the spiral */
-    extendRange(rangeToExtend, transform) {
-        this.activeStrokes.extendRange(rangeToExtend, transform);
-    }
-    /** compare various coordinate quantities */
-    isAlmostEqual(other) {
-        if (other instanceof DirectSpiral3d) {
-            return Geometry_1.Geometry.isSameCoordinate(this._nominalL1, other._nominalL1)
-                && Geometry_1.Geometry.isSameCoordinate(this._nominalR1, other._nominalR1)
-                && this.localToWorld.isAlmostEqual(other.localToWorld)
-                && this._activeFractionInterval.isAlmostEqual(other._activeFractionInterval)
-                && this._evaluator.isAlmostEqual(other._evaluator);
-        }
-        return false;
-    }
-}
-exports.DirectSpiral3d = DirectSpiral3d;
-//# sourceMappingURL=DirectSpiral3d.js.map

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

@@ -1 +0,0 @@
-{"version":3,"file":"DirectSpiral3d.js","sourceRoot":"","sources":["../../../src/curve/spiral/DirectSpiral3d.ts"],"names":[],"mappings":";;;AAAA;;;+FAG+F;AAC/F;;GAEG;AACH,6CAA0C;AAO1C,0DAAuD;AAEvD,kDAA+C;AAC/C,oDAAiD;AAEjD,qDAA6D;AAC7D,iEAAsF;AACtF,uFAAoF;AACpF,mFAAgF;AAEhF,6DAA0D;AAE1D,qFAA4E;AAC5E,6EAAoE;AACpE;;;;;;;;;;;;;;;EAeE;AACF,MAAa,cAAe,SAAQ,uCAAkB;IA+BpD,sEAAsE;IACtE,iEAAiE;IACjE,YACE,YAAuB,EACvB,UAA8B,EAC9B,kBAA+D,EAC/D,SAAiB,EACjB,SAAiB,EACjB,sBAA6C,EAC7C,SAA2B;QAC3B,KAAK,CAAC,UAAU,EAAE,YAAY,EAAE,sBAAsB,EAAE,kBAAkB,CAAC,CAAC;QAvC9E,wCAAwC;QAExB,uBAAkB,GAAG,kBAAkB,CAAC;QAsCtD,IAAI,CAAC,UAAU,GAAG,SAAS,CAAC;QAC5B,IAAI,CAAC,UAAU,GAAG,SAAS,CAAC;QAC5B,IAAI,CAAC,UAAU,GAAG,SAAS,CAAC;QAC5B,IAAI,CAAC,cAAc,GAAG,2BAAY,CAAC,MAAM,EAAE,CAAC;QAC5C,IAAI,CAAC,cAAc,GAAG,2BAAY,CAAC,MAAM,EAAE,CAAC;QAC5C,yFAAyF;QACzF,IAAI,CAAC,yBAAyB,EAAE,CAAC;IACnC,CAAC;IApCD,0EAA0E;IAC1E,IAAW,aAAa,KAAmB,OAAO,IAAI,CAAC,cAAc,KAAK,SAAS,CAAC,CAAC,CAAC,IAAI,CAAC,cAAc,CAAC,CAAC,CAAC,IAAI,CAAC,cAAc,CAAC,CAAC,CAAC;IAMlI,qCAAqC;IACrC,IAAW,SAAS,KAAa,OAAO,IAAI,CAAC,UAAU,CAAC,CAAC,CAAC;IAC1D,+DAA+D;IAC/D,IAAW,SAAS,KAAa,OAAO,IAAI,CAAC,UAAU,CAAC,CAAC,CAAC;IAC1D,uCAAuC;IACvC,IAAW,iBAAiB,KAAa,OAAO,uCAAkB,CAAC,iBAAiB,CAAC,IAAI,CAAC,UAAU,CAAC,CAAC,CAAC,CAAC;IACxG;;OAEG;IACH,IAAW,SAAS,KAAuB,OAAO,IAAI,CAAC,UAAU,CAAC,CAAC,CAAC;IAqBpE;;;;;OAKG;IACK,cAAc,CAAC,OAAqB,EAAE,SAAiB,EAAE,SAAiB,EAAE,WAAmB;QACrG,IAAI,WAAW,GAAG,CAAC;YACjB,WAAW,GAAG,CAAC,CAAC;QAClB,OAAO,CAAC,KAAK,EAAE,CAAC;QAChB,OAAO,CAAC,mBAAmB,EAAE,CAAC;QAC9B,OAAO,CAAC,oBAAoB,EAAE,CAAC;QAC/B,MAAM,SAAS,GAAG,OAAO,CAAC,cAAe,CAAC;QAC1C,MAAM,qBAAqB,GAAG,IAAI,CAAC,GAAG,CAAC,SAAS,GAAG,SAAS,CAAC,GAAG,IAAI,CAAC,UAAU,CAAC;QAChF,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,WAAW,EAAE,CAAC,EAAE,EAAE;YACrC,MAAM,QAAQ,GAAG,mBAAQ,CAAC,WAAW,CAAC,SAAS,EAAE,CAAC,GAAG,WAAW,EAAE,SAAS,CAAC,CAAC;YAC7E,MAAM,oBAAoB,GAAG,QAAQ,GAAG,qBAAqB,CAAC;YAC9D,OAAO,CAAC,YAAY,CAAC,OAAO,CAAC,IAAI,CAAC,UAAU,CAAC,WAAW,CAAC,QAAQ,CAAC,EAChE,IAAI,CAAC,UAAU,CAAC,WAAW,CAAC,QAAQ,CAAC,EAAE,CAAC,CAAC,CAAC;YAC5C,SAAS,CAAC,MAAM,CAAC,QAAQ,EAAE,oBAAoB,CAAC,CAAC,CAAC,4CAA4C;SAC/F;QAED,IAAI,SAAS,GAAG,SAAS,CAAC,yBAAyB,CAAC,CAAC,CAAC,CAAC;QACvD,IAAI,aAAa,GAAG,SAAS,CAAC,yBAAyB,CAAC,CAAC,CAAC,CAAC,CAAC,kDAAkD;QAC9G,IAAI,aAAa,EAAE,SAAS,CAAC;QAC7B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,WAAW,EAAE,CAAC,EAAE,EAAE;YACrC,SAAS,GAAG,SAAS,CAAC,yBAAyB,CAAC,CAAC,CAAC,CAAC;YACnD,aAAa,GAAG,aAAa,GAAG,IAAI,CAAC,UAAU,CAAC,iCAAiC,CAAC,SAAS,EAAE,SAAS,CAAC,CAAC;YACxG,SAAS,CAAC,yBAAyB,CAAC,CAAC,EAAE,SAAS,EAAE,aAAa,CAAC,CAAC;YACjE,SAAS,GAAG,SAAS,CAAC;YACtB,aAAa,GAAG,aAAa,CAAC;SAC/B;IAEH,CAAC;IACD,wBAAwB;IACjB,yBAAyB;QAC9B,MAAM,YAAY,GAAG,IAAI,CAAC,SAAS,GAAG,CAAC,GAAG,GAAG,IAAI,CAAC,SAAS,CAAC,CAAC;QAC7D,MAAM,WAAW,GAAG,IAAI,CAAC;QACzB,MAAM,WAAW,GAAG,6BAAa,CAAC,aAAa,CAAC,SAAS,EAAE,CAAC,EAAE,YAAY,EAAE,WAAW,CAAC,CAAC;QACzF,IAAI,CAAC,cAAc,CAAC,IAAI,CAAC,cAAc,EAAE,CAAC,EAAE,CAAC,EAAE,WAAW,CAAC,CAAC;QAC5D,MAAM,iBAAiB,GAAG,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,uBAAuB,CAAC,aAAa,EAAE,GAAG,WAAW,CAAC,CAAC;QAChG,IAAI,CAAC,cAAc,GAAG,2BAAY,CAAC,MAAM,EAAE,CAAC;QAC5C,IAAI,CAAC,cAAc,CAAC,IAAI,CAAC,cAAc,EAAE,IAAI,CAAC,uBAAuB,CAAC,EAAE,EAAE,IAAI,CAAC,uBAAuB,CAAC,EAAE,EACvG,iBAAiB,CAAC,CAAC;IACvB,CAAC;IACD;;;;;;;;OAQG;IACI,MAAM,CAAC,uBAAuB,CACnC,UAAkB,EAClB,YAAuB,EACvB,QAAgB,EAAE,QAAgB,EAClC,kBAA+D,EAC/D,SAAiB,EAAE,SAAiB,EACpC,cAAqC;QACrC,IAAI,QAAQ,GAAG,CAAC;YACd,QAAQ,GAAG,CAAC,CAAC;QACf,IAAI,QAAQ,GAAG,CAAC;YACd,QAAQ,GAAG,CAAC,CAAC;QACf,MAAM,SAAS,GAAG,IAAI,0CAAyB,CAAC,SAAS,EAAE,GAAG,GAAG,CAAC,GAAG,GAAG,SAAS,GAAG,SAAS,CAAC,EAAE,QAAQ,EAAE,QAAQ,CAAC,CAAC;QACpH,OAAO,IAAI,cAAc,CACvB,YAAY,CAAC,KAAK,EAAE,EACpB,UAAU,EACV,kBAAkB,EAClB,SAAS,EAAE,SAAS,EACpB,cAAc,CAAC,CAAC,CAAC,cAAc,CAAC,KAAK,EAAE,CAAC,CAAC,CAAC,qBAAS,CAAC,MAAM,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,SAAS,CAAC,CAAC;IACjF,CAAC;IACD;;;;;;;;;OASG;IACI,MAAM,CAAC,mBAAmB,CAC/B,YAAuB,EACvB,SAAiB,EACjB,SAAiB,EACjB,cAA0B;QAC1B,OAAO,IAAI,CAAC,uBAAuB,CAAC,eAAe,EAAE,YAAY,EAAE,CAAC,EAAE,CAAC,EAAE,SAAS,EAAE,SAAS,EAAE,SAAS,EAAE,cAAc,CAAC,CAAC;IAC5H,CAAC;IACD;;;;;;;;;;;;OAYG;IACI,MAAM,CAAC,gBAAgB,CAC5B,YAAuB,EACvB,SAAiB,EACjB,SAAiB,EACjB,cAA0B;QAC1B,MAAM,SAAS,GAAG,2CAAoB,CAAC,MAAM,CAAC,SAAS,EAAE,SAAS,CAAC,CAAC;QACpE,IAAI,SAAS,KAAK,SAAS;YACzB,OAAO,SAAS,CAAC;QACnB,OAAO,IAAI,cAAc,CACvB,YAAY,CAAC,KAAK,EAAE,EACpB,OAAO,EACP,SAAS,EACT,SAAS,EAAE,SAAS,EACpB,cAAc,CAAC,CAAC,CAAC,cAAc,CAAC,KAAK,EAAE,CAAC,CAAC,CAAC,qBAAS,CAAC,MAAM,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,SAAS,CAAC,CAAC;IACjF,CAAC;IACD;;;;;;;;;;;;;OAaG;IACI,MAAM,CAAC,aAAa,CACzB,YAAuB,EACvB,SAAiB,EACjB,SAAiB,EACjB,cAA0B;QAC1B,MAAM,SAAS,GAAG,6CAAsB,CAAC,MAAM,CAAC,SAAS,EAAE,SAAS,CAAC,CAAC;QACtE,IAAI,SAAS,KAAK,SAAS;YACzB,OAAO,SAAS,CAAC;QACnB,OAAO,IAAI,cAAc,CACvB,YAAY,CAAC,KAAK,EAAE,EACpB,SAAS,EACT,SAAS,EACT,SAAS,EAAE,SAAS,EACpB,cAAc,CAAC,CAAC,CAAC,cAAc,CAAC,KAAK,EAAE,CAAC,CAAC,CAAC,qBAAS,CAAC,MAAM,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,SAAS,CAAC,CAAC;IACjF,CAAC;IAED;;;;;;;;;;;;OAYG;IACI,MAAM,CAAC,qBAAqB,CACjC,YAAuB,EACvB,SAAiB,EACjB,SAAiB,EACjB,cAA0B;QAC1B,MAAM,SAAS,GAAG,yDAAwB,CAAC,MAAM,CAAC,SAAS,EAAE,SAAS,CAAC,CAAC;QACxE,IAAI,SAAS,KAAK,SAAS;YACzB,OAAO,SAAS,CAAC;QACnB,OAAO,IAAI,cAAc,CACvB,YAAY,CAAC,KAAK,EAAE,EACpB,iBAAiB,EACjB,SAAS,EACT,SAAS,EAAE,SAAS,EACpB,cAAc,CAAC,CAAC,CAAC,cAAc,CAAC,KAAK,EAAE,CAAC,CAAC,CAAC,qBAAS,CAAC,MAAM,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,SAAS,CAAC,CAAC;IACjF,CAAC;IAED;;;;;;;;;OASG;IACI,MAAM,CAAC,iBAAiB,CAC7B,YAAuB,EACvB,SAAiB,EACjB,SAAiB,EACjB,cAA0B;QAC1B,MAAM,SAAS,GAAG,iDAAoB,CAAC,MAAM,CAAC,SAAS,EAAE,SAAS,CAAC,CAAC;QACpE,IAAI,SAAS,KAAK,SAAS;YACzB,OAAO,SAAS,CAAC;QACnB,OAAO,IAAI,cAAc,CACvB,YAAY,CAAC,KAAK,EAAE,EACpB,aAAa,EACb,SAAS,EACT,SAAS,EAAE,SAAS,EACpB,cAAc,CAAC,CAAC,CAAC,cAAc,CAAC,KAAK,EAAE,CAAC,CAAC,CAAC,qBAAS,CAAC,MAAM,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,SAAS,CAAC,CAAC;IACjF,CAAC;IAED;;;;;;;;;;;;OAYG;IACI,MAAM,CAAC,oBAAoB,CAChC,YAAuB,EACvB,SAAiB,EACjB,SAAiB,EACjB,cAA0B;QAC1B,MAAM,SAAS,GAAG,6DAA6B,CAAC,MAAM,CAAC,SAAS,EAAE,SAAS,CAAC,CAAC;QAC7E,IAAI,SAAS,KAAK,SAAS;YACzB,OAAO,SAAS,CAAC;QACnB,OAAO,IAAI,cAAc,CACvB,YAAY,CAAC,KAAK,EAAE,EACpB,oBAAoB,EACpB,SAAS,EACT,SAAS,EAAE,SAAS,EACpB,cAAc,CAAC,CAAC,CAAC,cAAc,CAAC,KAAK,EAAE,CAAC,CAAC,CAAC,qBAAS,CAAC,MAAM,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,SAAS,CAAC,CAAC;IACjF,CAAC;IAEM,MAAM,CAAC,sBAAsB,CAClC,YAAuB,EACvB,SAAiB,EACjB,SAAiB,EACjB,cAA0B;QAC1B,OAAO,IAAI,IAAI,CAAC,YAAY,EAAE,YAAY,EAAE,SAAS,EAAE,SAAS,EAAE,SAAS,EAAE,cAAc,EACzF,IAAI,iEAA+B,CAAC,SAAS,EAAE,SAAS,CAAC,CAAC,CAAC;IAC/D,CAAC;IACD;;;;;;;;;OASG;IACI,MAAM,CAAC,WAAW,CACvB,YAAuB,EACvB,SAAiB,EACjB,SAAiB,EACjB,cAA0B;QAC1B,OAAO,IAAI,CAAC,uBAAuB,CAAC,OAAO,EAAE,YAAY,EAAE,CAAC,EAAE,CAAC,EAAE,SAAS,EAAE,SAAS,EAAE,SAAS,EAAE,cAAc,CAAC,CAAC;IACpH,CAAC;IAED;;;;;;;;;OASG;IACI,MAAM,CAAC,kBAAkB,CAC9B,YAAuB,EACvB,SAAiB,EACjB,SAAiB,EACjB,cAA0B;QAC1B,OAAO,IAAI,CAAC,uBAAuB,CAAC,cAAc,EAAE,YAAY,EAAE,CAAC,EAAE,CAAC,EAAE,SAAS,EAAE,SAAS,EAAE,SAAS,EAAE,cAAc,CAAC,CAAC;IAC3H,CAAC;IACD;;;;;;;;OAQG;IACI,MAAM,CAAC,uBAAuB,CACnC,YAAuB,EACvB,SAAiB,EACjB,SAAiB,EACjB,cAA0B;QAC1B,OAAO,IAAI,CAAC,uBAAuB,CAAC,mBAAmB,EAAE,YAAY,EAAE,CAAC,EAAE,CAAC,EAAE,SAAS,EAAE,SAAS,EAAE,SAAS,EAAE,cAAc,CAAC,CAAC;IAChI,CAAC;IACD;;;;;;;;;;;;;;;;;;OAkBG;IACI,MAAM,CAAC,yBAAyB,CACrC,UAAkB,EAClB,OAA2B,EAC3B,OAA2B,EAC3B,QAA2B,EAC3B,SAA4B,EAC5B,SAA6B,EAC7B,cAAqC,EACrC,YAAuB;QACvB,IAAI,QAAQ,KAAK,SAAS,IAAI,CAAC,QAAQ,CAAC,YAAY;YAClD,OAAO,SAAS,CAAC;QACnB,IAAI,OAAO,KAAK,SAAS,IAAI,CAAC,mBAAQ,CAAC,qBAAqB,CAAC,OAAO,CAAC;YACnE,OAAO,SAAS,CAAC;QACnB,IAAI,OAAO,KAAK,SAAS,IAAI,mBAAQ,CAAC,qBAAqB,CAAC,OAAO,CAAC;YAClE,OAAO,SAAS,CAAC;QACnB,IAAI,SAAS,KAAK,SAAS;YACzB,OAAO,SAAS,CAAC;QACnB,IAAI,mBAAQ,CAAC,iBAAiB,CAAC,UAAU,EAAE,OAAO,CAAC;YACjD,OAAO,IAAI,CAAC,WAAW,CAAC,YAAY,EAAE,SAAS,EAAE,OAAO,EAAE,cAAc,CAAC,CAAC;QAC5E,IAAI,mBAAQ,CAAC,iBAAiB,CAAC,UAAU,EAAE,cAAc,CAAC;YACxD,OAAO,IAAI,CAAC,kBAAkB,CAAC,YAAY,EAAE,SAAS,EAAE,OAAO,EAAE,cAAc,CAAC,CAAC;QACnF,IAAI,mBAAQ,CAAC,iBAAiB,CAAC,UAAU,EAAE,eAAe,CAAC;YACzD,OAAO,IAAI,CAAC,mBAAmB,CAAC,YAAY,EAAE,SAAS,EAAE,OAAO,EAAE,cAAc,CAAC,CAAC;QACpF,IAAI,mBAAQ,CAAC,iBAAiB,CAAC,UAAU,EAAE,YAAY,CAAC;YACtD,OAAO,IAAI,CAAC,sBAAsB,CAAC,YAAY,EAAE,SAAS,EAAE,OAAO,EAAE,cAAc,CAAC,CAAC;QACvF,IAAI,mBAAQ,CAAC,iBAAiB,CAAC,UAAU,EAAE,OAAO,CAAC;YACjD,OAAO,IAAI,CAAC,gBAAgB,CAAC,YAAY,EAAE,SAAS,EAAE,OAAO,EAAE,cAAc,CAAC,CAAC;QACjF,IAAI,mBAAQ,CAAC,iBAAiB,CAAC,UAAU,EAAE,SAAS,CAAC;YACnD,OAAO,IAAI,CAAC,aAAa,CAAC,YAAY,EAAE,SAAS,EAAE,OAAO,EAAE,cAAc,CAAC,CAAC;QAC9E,IAAI,mBAAQ,CAAC,iBAAiB,CAAC,UAAU,EAAE,oBAAoB,CAAC;YAC9D,OAAO,IAAI,CAAC,oBAAoB,CAAC,YAAY,EAAE,SAAS,EAAE,OAAO,EAAE,cAAc,CAAC,CAAC;QACrF,IAAI,mBAAQ,CAAC,iBAAiB,CAAC,UAAU,EAAE,iBAAiB,CAAC;YAC3D,OAAO,IAAI,CAAC,qBAAqB,CAAC,YAAY,EAAE,SAAS,EAAE,OAAO,EAAE,cAAc,CAAC,CAAC;QACtF,IAAI,mBAAQ,CAAC,iBAAiB,CAAC,UAAU,EAAE,mBAAmB,CAAC;YAC7D,OAAO,IAAI,CAAC,uBAAuB,CAAC,YAAY,EAAE,SAAS,EAAE,OAAO,EAAE,cAAc,CAAC,CAAC;QACxF,IAAI,mBAAQ,CAAC,iBAAiB,CAAC,UAAU,EAAE,aAAa,CAAC;YACvD,OAAO,IAAI,CAAC,iBAAiB,CAAC,YAAY,EAAE,SAAS,EAAE,OAAO,EAAE,cAAc,CAAC,CAAC;QAClF,OAAO,SAAS,CAAC;IACnB,CAAC;IACD,gCAAgC;IACzB,KAAK;;QACV,OAAO,IAAI,cAAc,CACvB,IAAI,CAAC,YAAY,CAAC,KAAK,EAAE,EACzB,IAAI,CAAC,WAAW,EAChB,MAAA,IAAI,CAAC,gBAAgB,0CAAE,KAAK,EAAE,EAC9B,IAAI,CAAC,UAAU,EACf,IAAI,CAAC,UAAU,EACf,MAAA,IAAI,CAAC,uBAAuB,0CAAE,KAAK,EAAE,EACrC,IAAI,CAAC,UAAU,CAAC,KAAK,EAAE,CAAC,CAAC;IAC7B,CAAC;IACD,mEAAmE;IAC5D,mBAAmB,CAAC,UAAqB;QAC9C,MAAM,SAAS,GAAG,IAAI,CAAC,yBAAyB,CAAC,UAAU,CAAC,CAAC;QAC7D,IAAI,SAAS,KAAK,SAAS,EAAE;YAC3B,IAAI,CAAC,UAAU,IAAI,SAAS,CAAC,KAAK,CAAC;YACnC,IAAI,CAAC,UAAU,IAAI,SAAS,CAAC,KAAK,CAAC;YACnC,IAAI,CAAC,SAAS,CAAC,YAAY,CAAC,SAAS,CAAC,KAAK,CAAC,CAAC;SAC9C;QACD,IAAI,CAAC,yBAAyB,EAAE,CAAC;QACjC,OAAO,IAAI,CAAC;IACd,CAAC;IACD,sCAAsC;IAC/B,gBAAgB,CAAC,SAAoB;QAC1C,MAAM,MAAM,GAAG,IAAI,CAAC,KAAK,EAAE,CAAC;QAC5B,MAAM,CAAC,mBAAmB,CAAC,SAAS,CAAC,CAAC,CAAC,2CAA2C;QAClF,OAAO,MAAM,CAAC;IAChB,CAAC;IACD,qCAAqC;IACrB,UAAU,KAAc,OAAO,IAAI,CAAC,aAAa,CAAC,UAAU,EAAE,CAAC,CAAC,CAAC;IACjF,mCAAmC;IACnB,QAAQ,KAAc,OAAO,IAAI,CAAC,aAAa,CAAC,QAAQ,EAAE,CAAC,CAAC,CAAC;IAC7E,mFAAmF;IAC5E,SAAS,CAAC,KAAmC;QAClD,OAAO,KAAK,CAAC,cAAc,CAAC,IAAI,CAAC,YAAY,CAAC,MAAiB,CAAC;eAC3D,mBAAQ,CAAC,gBAAgB,CAAC,GAAG,EAAE,IAAI,CAAC,YAAY,CAAC,MAAM,CAAC,UAAU,CAAC,KAAK,CAAC,YAAY,EAAE,CAAC,CAAC;eACzF,mBAAQ,CAAC,gBAAgB,CAAC,GAAG,EAAE,IAAI,CAAC,YAAY,CAAC,MAAM,CAAC,UAAU,CAAC,KAAK,CAAC,YAAY,EAAE,CAAC,CAAC,CAAC;IACjG,CAAC;IACD;;OAEG;IACI,WAAW;QAChB,MAAM,YAAY,GAAG,IAAI,CAAC,cAAc,CAAC,cAAe,CAAC;QACzD,MAAM,CAAC,GAAG,YAAY,CAAC,MAAM,CAAC;QAC9B,OAAO,YAAY,CAAC,yBAAyB,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;IACvD,CAAC;IACD;;OAEG;IACa,WAAW,KAAK,OAAO,IAAI,CAAC,WAAW,EAAE,CAAC,CAAC,CAAC;IAC5D,6DAA6D;IACtD,mBAAmB,CAAC,KAAU,IAAa,OAAO,KAAK,YAAY,cAAc,CAAC,CAAC,CAAC;IAC3F;;;OAGG;IACI,WAAW,CAAC,IAAkB,EAAE,OAAuB,IAAU,IAAI,CAAC,aAAa,CAAC,WAAW,CAAC,IAAI,EAAE,OAAO,CAAC,CAAC,CAAC,CAAC;IACxH,+CAA+C;IACxC,kBAAkB,CAAC,IAAoB,EAAE,OAAuB;QACrE,MAAM,CAAC,GAAG,IAAI,CAAC,4BAA4B,CAAC,OAAO,CAAC,CAAC;QACrD,MAAM,aAAa,GAAG,IAAI,CAAC,aAAa,CAAC;QACzC,IAAI,CAAC,yBAAyB,CAAC,IAAI,CAAC,CAAC;QACrC,IAAI,CAAC,IAAI,aAAa,CAAC,SAAS,EAAE,EAAE;YAClC,wDAAwD;YACxD,IAAI,CAAC,qCAAqC,CAAC,IAAI,EAAE,CAAC,GAAG,aAAa,CAAC,SAAS,EAAE,EAAE,GAAG,EAAE,GAAG,CAAC,CAAC;SAC3F;aAAM;YACL,IAAI,CAAC,qCAAqC,CAAC,IAAI,EAAE,CAAC,EAAE,GAAG,EAAE,GAAG,CAAC,CAAC;SAC/D;QACD,IAAI,CAAC,uBAAuB,CAAC,IAAI,CAAC,CAAC;IACrC,CAAC;IAED;;;OAGG;IAEI,4BAA4B,CAAC,OAAuB;QACzD,IAAI,SAAS,CAAC;QACd,MAAM,cAAc,GAAG,IAAI,CAAC,UAAU,GAAG,CAAC,GAAG,GAAG,IAAI,CAAC,UAAU,CAAC,CAAC;QACjE,IAAI,OAAO,EAAE;YACX,MAAM,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,UAAU,CAAC,CAAC;YACvC,SAAS,GAAG,OAAO,CAAC,oBAAoB,CAAC,IAAI,EAAE,cAAc,CAAC,CAAC;YAC/D,SAAS,GAAG,OAAO,CAAC,kBAAkB,CAAC,SAAS,EAAE,IAAI,CAAC,WAAW,EAAE,CAAC,CAAC;YACtE,SAAS,GAAG,OAAO,CAAC,2BAA2B,CAAC,SAAS,CAAC,CAAC;SAC5D;aAAM;YACL,SAAS,GAAG,6BAAa,CAAC,aAAa,CAAC,SAAS,EAAE,CAAC,EAAE,cAAc,EAAE,IAAI,CAAC,CAAC;SAC7E;QACD,SAAS,GAAG,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,uBAAuB,CAAC,aAAa,EAAE,GAAG,SAAS,CAAC,CAAC;QAChF,OAAO,SAAS,CAAC;IACnB,CAAC;IAED;;OAEG;IAEI,cAAc;QACnB,IAAI,CAAC,uBAAuB,CAAC,cAAc,EAAE,CAAC;QAC9C,IAAI,IAAI,CAAC,cAAc,KAAK,SAAS;YACnC,IAAI,CAAC,cAAc,GAAG,IAAI,CAAC,cAAc,CAAC,KAAK,EAAE,CAAC;QACpD,IAAI,CAAC,cAAc,CAAC,cAAc,EAAE,CAAC;IACvC,CAAC;IACD,qDAAqD;IAC9C,eAAe,CAAC,cAAsB,EAAE,MAAgB;QAC7D,MAAM,cAAc,GAAG,IAAI,CAAC,uBAAuB,CAAC,eAAe,CAAC,cAAc,CAAC,CAAC;QACpF,MAAM,GAAG,IAAI,CAAC,UAAU,CAAC,eAAe,CAAC,cAAc,EAAE,MAAM,CAAC,CAAC;QACjE,IAAI,CAAC,YAAY,CAAC,eAAe,CAAC,MAAM,EAAE,MAAM,CAAC,CAAC;QAClD,OAAO,MAAM,CAAC;IAChB,CAAC;IACD,oEAAoE;IAC7D,4BAA4B,CAAC,cAAsB,EAAE,MAAc;QACxE,MAAM,cAAc,GAAG,IAAI,CAAC,uBAAuB,CAAC,eAAe,CAAC,cAAc,CAAC,CAAC;QACpF,MAAM,GAAG,IAAI,CAAC,UAAU,CAAC,4BAA4B,CAAC,cAAc,EAAE,MAAM,CAAC,CAAC;QAC9E,MAAM,CAAC,SAAS,CAAC,YAAY,CAAC,IAAI,CAAC,uBAAuB,CAAC,WAAW,EAAE,CAAC,CAAC;QAC1E,MAAM,CAAC,gBAAgB,CAAC,IAAI,CAAC,YAAY,CAAC,CAAC;QAC3C,OAAO,MAAM,CAAC;IAChB,CAAC;IAED;;;;;OAKG;IAEI,8BAA8B,CAAC,cAAsB,EAAE,MAAkC;QAC9F,MAAM,cAAc,GAAG,IAAI,CAAC,uBAAuB,CAAC,eAAe,CAAC,cAAc,CAAC,CAAC;QACpF,MAAM,GAAG,IAAI,CAAC,UAAU,CAAC,8BAA8B,CAAC,cAAc,EAAE,MAAM,CAAC,CAAC;QAChF,MAAM,CAAC,GAAG,IAAI,CAAC,uBAAuB,CAAC,WAAW,EAAE,CAAC;QACrD,MAAM,CAAC,OAAO,CAAC,YAAY,CAAC,CAAC,CAAC,CAAC;QAC/B,MAAM,CAAC,OAAO,CAAC,YAAY,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;QACnC,MAAM,CAAC,gBAAgB,CAAC,IAAI,CAAC,YAAY,CAAC,CAAC;QAC3C,OAAO,MAAM,CAAC;IAChB,CAAC;IACD,mFAAmF;IAC5E,yBAAyB,CAAC,OAAwB;QACvD,OAAO,OAAO,CAAC,sBAAsB,CAAC,IAAI,CAAC,CAAC;IAC9C,CAAC;IACD,oDAAoD;IAC7C,WAAW,CAAC,aAAsB,EAAE,SAAqB;QAC9D,IAAI,CAAC,aAAa,CAAC,WAAW,CAAC,aAAa,EAAE,SAAS,CAAC,CAAC;IAC3D,CAAC;IACD,4CAA4C;IAC5B,aAAa,CAAC,KAAU;QACtC,IAAI,KAAK,YAAY,cAAc,EAAE;YACnC,OAAO,mBAAQ,CAAC,gBAAgB,CAAC,IAAI,CAAC,UAAU,EAAE,KAAK,CAAC,UAAU,CAAC;mBAC9D,mBAAQ,CAAC,gBAAgB,CAAC,IAAI,CAAC,UAAU,EAAE,KAAK,CAAC,UAAU,CAAC;mBAC5D,IAAI,CAAC,YAAY,CAAC,aAAa,CAAC,KAAK,CAAC,YAAY,CAAC;mBACnD,IAAI,CAAC,uBAAuB,CAAC,aAAa,CAAC,KAAK,CAAC,uBAAuB,CAAC;mBACzE,IAAI,CAAC,UAAU,CAAC,aAAa,CAAC,KAAK,CAAC,UAAU,CAAC,CAAC;SACtD;QACD,OAAO,KAAK,CAAC;IACf,CAAC;CACF;AAviBD,wCAuiBC","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/** @packageDocumentation\r\n * @module Curve\r\n */\r\nimport { Geometry } from \"../../Geometry\";\r\nimport { GeometryHandler, IStrokeHandler } from \"../../geometry3d/GeometryHandler\";\r\nimport { Plane3dByOriginAndUnitNormal } from \"../../geometry3d/Plane3dByOriginAndUnitNormal\";\r\nimport { Plane3dByOriginAndVectors } from \"../../geometry3d/Plane3dByOriginAndVectors\";\r\nimport { Point3d } from \"../../geometry3d/Point3dVector3d\";\r\nimport { Range3d } from \"../../geometry3d/Range\";\r\nimport { Ray3d } from \"../../geometry3d/Ray3d\";\r\nimport { Segment1d } from \"../../geometry3d/Segment1d\";\r\nimport { Transform } from \"../../geometry3d/Transform\";\r\nimport { LineString3d } from \"../LineString3d\";\r\nimport { StrokeOptions } from \"../StrokeOptions\";\r\nimport { TransitionConditionalProperties } from \"./TransitionConditionalProperties\";\r\nimport { ClothoidSeriesRLEvaluator } from \"./ClothoidSeries\";\r\nimport { CzechSpiralEvaluator, ItalianSpiralEvaluator } from \"./CzechSpiralEvaluator\";\r\nimport { DirectHalfCosineSpiralEvaluator } from \"./DirectHalfCosineSpiralEvaluator\";\r\nimport { AustralianRailCorpXYEvaluator } from \"./AustralianRailCorpXYEvaluator\";\r\nimport { XYCurveEvaluator } from \"./XYCurveEvaluator\";\r\nimport { TransitionSpiral3d } from \"./TransitionSpiral3d\";\r\nimport { Angle } from \"../../geometry3d/Angle\";\r\nimport { MXCubicAlongArcEvaluator } from \"./MXCubicAlongArcSpiralEvaluator\";\r\nimport { PolishCubicEvaluator } from \"./PolishCubicSpiralEvaluator\";\r\n/**\r\n* DirectSpiral3d acts like a TransitionSpiral3d for serialization purposes, but implements spiral types that have \"direct\" xy calculations without the integrations required\r\n* for IntegratedSpiral3d.\r\n* * Each DirectSpiral3d carries an XYCurveEvaluator to give it specialized behavior.\r\n* * Direct spirals that flow through serialization to native imodel02 are create with these static methods:\r\n*   * createArema\r\n*   * createJapaneseCubic\r\n*   * createAustralianRail\r\n*   * createDirectHalfCosine\r\n*   * createChineseCubic\r\n*   * createCzechCubic\r\n*   * createPolishCubic\r\n*   * createItalian\r\n*   * createWesternAustralian\r\n* @public\r\n*/\r\nexport class DirectSpiral3d extends TransitionSpiral3d {\r\n\r\n  /** String name for schema properties */\r\n\r\n  public readonly curvePrimitiveType = \"transitionSpiral\";\r\n\r\n  /** stroked approximation of entire spiral. */\r\n  private _globalStrokes: LineString3d;\r\n  /** stroked approximation of active spiral.\r\n   * * Same count as global -- possibly overly fine, but it gives some consistency between same clothoid constructed as partial versus complete.\r\n   * * If no trimming, this points to the same place as the _globalStrokes !!!  Don't double transform!!!\r\n   */\r\n  private _activeStrokes?: LineString3d;\r\n  /** Return the internal stroked form of the (possibly partial) spiral   */\r\n  public get activeStrokes(): LineString3d { return this._activeStrokes !== undefined ? this._activeStrokes : this._globalStrokes; }\r\n\r\n  private _nominalL1: number;\r\n  private _nominalR1: number;\r\n  private _evaluator: XYCurveEvaluator;\r\n\r\n  /** Return the nominal end radius. */\r\n  public get nominalR1(): number { return this._nominalR1; }\r\n  /** Return the nominal distance from inflection to endpoint. */\r\n  public get nominalL1(): number { return this._nominalL1; }\r\n  /** Return the nominal end curvature */\r\n  public get nominalCurvature1(): number { return TransitionSpiral3d.radiusToCurvature(this._nominalR1); }\r\n  /** Return the low level evaluator\r\n   * @internal\r\n   */\r\n  public get evaluator(): XYCurveEvaluator { return this._evaluator; }\r\n\r\n  // constructor demands radius1 and distance1 for nominal construction.\r\n  // caller is responsible for managing intervals of partial spiral\r\n  constructor(\r\n    localToWorld: Transform,\r\n    spiralType: string | undefined,\r\n    originalProperties: TransitionConditionalProperties | undefined,\r\n    nominalL1: number,\r\n    nominalR1: number,\r\n    activeFractionInterval: Segment1d | undefined,\r\n    evaluator: XYCurveEvaluator) {\r\n    super(spiralType, localToWorld, activeFractionInterval, originalProperties);\r\n    this._nominalL1 = nominalL1;\r\n    this._nominalR1 = nominalR1;\r\n    this._evaluator = evaluator;\r\n    this._globalStrokes = LineString3d.create();\r\n    this._activeStrokes = LineString3d.create();\r\n    // initialize for compiler -- but this will be recomputed in refreshComputeProperties ...\r\n    this.refreshComputedProperties();\r\n  }\r\n  /**\r\n   * Compute stroke data in an interval.\r\n   * @param strokes strokes to clear and refill.\r\n   * @param fraction0 start fraction\r\n   * @param fraction1 end fraction\r\n   */\r\n  private computeStrokes(strokes: LineString3d, fractionA: number, fractionB: number, numInterval: number) {\r\n    if (numInterval < 1)\r\n      numInterval = 1;\r\n    strokes.clear();\r\n    strokes.ensureEmptyUVParams();\r\n    strokes.ensureEmptyFractions();\r\n    const distances = strokes.packedUVParams!;\r\n    const nominalIntervalLength = Math.abs(fractionB - fractionA) * this._nominalL1;\r\n    for (let i = 0; i <= numInterval; i++) {\r\n      const fraction = Geometry.interpolate(fractionA, i / numInterval, fractionB);\r\n      const nominalDistanceAlong = fraction * nominalIntervalLength;\r\n      strokes.packedPoints.pushXYZ(this._evaluator.fractionToX(fraction),\r\n        this._evaluator.fractionToY(fraction), 0);\r\n      distances.pushXY(fraction, nominalDistanceAlong); // the second distance will be updated below\r\n    }\r\n\r\n    let fraction0 = distances.getXAtUncheckedPointIndex(0);\r\n    let trueDistance0 = distances.getYAtUncheckedPointIndex(0); // whatever was assigned as start distance is fine\r\n    let trueDistance1, fraction1;\r\n    for (let i = 1; i <= numInterval; i++) {\r\n      fraction1 = distances.getXAtUncheckedPointIndex(i);\r\n      trueDistance1 = trueDistance0 + this._evaluator.integrateDistanceBetweenFractions(fraction0, fraction1);\r\n      distances.setXYZAtCheckedPointIndex(i, fraction1, trueDistance1);\r\n      fraction0 = fraction1;\r\n      trueDistance0 = trueDistance1;\r\n    }\r\n\r\n  }\r\n  /** Recompute strokes */\r\n  public refreshComputedProperties() {\r\n    const sweepRadians = this.nominalL1 / (2.0 * this.nominalR1);\r\n    const radiansStep = 0.02;\r\n    const numInterval = StrokeOptions.applyAngleTol(undefined, 4, sweepRadians, radiansStep);\r\n    this.computeStrokes(this._globalStrokes, 0, 1, numInterval);\r\n    const numActiveInterval = Math.ceil(this._activeFractionInterval.absoluteDelta() * numInterval);\r\n    this._activeStrokes = LineString3d.create();\r\n    this.computeStrokes(this._activeStrokes, this._activeFractionInterval.x0, this._activeFractionInterval.x1,\r\n      numActiveInterval);\r\n  }\r\n  /**\r\n   * Create a spiral object which uses numXTerm terms from the clothoid X series and numYTerm from the clothoid Y series.\r\n   * @param numXTerm  number of terms to use from X series\r\n   * @param numYTerm number of terms to use from Y series\r\n   * @param localToWorld placement frame.  Inflection point is at origin, initial direction is along x axis.\r\n   * @param nominalL1 design distance from inflection to end point.\r\n   * @param nominalR1 design radius at end point.\r\n   * @param activeInterval active interval (as fractions of nominalL1 !!!)\r\n   */\r\n  public static createTruncatedClothoid(\r\n    spiralType: string,\r\n    localToWorld: Transform,\r\n    numXTerm: number, numYTerm: number,\r\n    originalProperties: TransitionConditionalProperties | undefined,\r\n    nominalL1: number, nominalR1: number,\r\n    activeInterval: Segment1d | undefined): DirectSpiral3d | undefined {\r\n    if (numXTerm < 1)\r\n      numXTerm = 1;\r\n    if (numYTerm < 1)\r\n      numYTerm = 1;\r\n    const evaluator = new ClothoidSeriesRLEvaluator(nominalL1, 1.0 / (2.0 * nominalL1 * nominalR1), numXTerm, numYTerm);\r\n    return new DirectSpiral3d(\r\n      localToWorld.clone(),\r\n      spiralType,\r\n      originalProperties,\r\n      nominalL1, nominalR1,\r\n      activeInterval ? activeInterval.clone() : Segment1d.create(0, 1), evaluator);\r\n  }\r\n  /**\r\n   * Create an Japanese spiral clothoid approximation\r\n   *   * X is 1 terms of the clothoid series as a function of nominal distance along.\r\n   *   * Y is 1 terms f the clothoid series as a function of nominal distance along.\r\n   *   * Remark: This is identical to the ChineseCubic\r\n   * @param localToWorld axes with inflection at origin, tangent along x axis\r\n   * @param nominalL1 nominal length as used in series LR terms.\r\n   * @param nominalR1 nominal final radius as used in series LR terms\r\n   * @param activeInterval fractional interval with (0, nominalL1) range for nominal distance along\r\n   */\r\n  public static createJapaneseCubic(\r\n    localToWorld: Transform,\r\n    nominalL1: number,\r\n    nominalR1: number,\r\n    activeInterval?: Segment1d): DirectSpiral3d | undefined {\r\n    return this.createTruncatedClothoid(\"JapaneseCubic\", localToWorld, 1, 1, undefined, nominalL1, nominalR1, activeInterval);\r\n  }\r\n  /**\r\n   * Create a 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 up a little bit to simulate the lost series terms.\r\n   * @param localToWorld\r\n   * @param nominalLx nominal length along x axis\r\n   * @param nominalR1\r\n   * @param activeInterval\r\n   */\r\n  public static createCzechCubic(\r\n    localToWorld: Transform,\r\n    nominalLx: number,\r\n    nominalR1: number,\r\n    activeInterval?: Segment1d): DirectSpiral3d | undefined {\r\n    const evaluator = CzechSpiralEvaluator.create(nominalLx, nominalR1);\r\n    if (evaluator === undefined)\r\n      return undefined;\r\n    return new DirectSpiral3d(\r\n      localToWorld.clone(),\r\n      \"Czech\",\r\n      undefined,\r\n      nominalLx, nominalR1,\r\n      activeInterval ? activeInterval.clone() : Segment1d.create(0, 1), evaluator);\r\n  }\r\n  /**\r\n   * Create an italian spiral\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 up a little bit to simulate the lost series terms.\r\n   * * L in gamma and m is the\r\n   * @param localToWorld\r\n   * @param nominalL1 nominal length along the spiral\r\n   * @param nominalR1\r\n   * @param activeInterval\r\n   */\r\n  public static createItalian(\r\n    localToWorld: Transform,\r\n    nominalL1: number,\r\n    nominalR1: number,\r\n    activeInterval?: Segment1d): DirectSpiral3d | undefined {\r\n    const evaluator = ItalianSpiralEvaluator.create(nominalL1, nominalR1);\r\n    if (evaluator === undefined)\r\n      return undefined;\r\n    return new DirectSpiral3d(\r\n      localToWorld.clone(),\r\n      \"Italian\",\r\n      undefined,\r\n      nominalL1, nominalR1,\r\n      activeInterval ? activeInterval.clone() : Segment1d.create(0, 1), evaluator);\r\n  }\r\n\r\n  /**\r\n   * Create an MX Cubic whose nominal length is close to along the curve.\r\n   * This is y= m*x^3 with\r\n   * * m is 1/ (6RL1)\r\n   *    * 1/(6RL) is the leading term of the sine series.\r\n   * * L1 is an along-the-x-axis distance that is slightly LESS THAN the nominal length\r\n   * * x is axis position that is slightly LESS than nominal distance along\r\n   * * L1, x use the approximation   `x = s * ( 1 - s^4/ (40 R R L L))\r\n   * @param localToWorld\r\n   * @param nominalL1\r\n   * @param nominalR1\r\n   * @param activeInterval\r\n   */\r\n  public static createMXCubicAlongArc(\r\n    localToWorld: Transform,\r\n    nominalL1: number,\r\n    nominalR1: number,\r\n    activeInterval?: Segment1d): DirectSpiral3d | undefined {\r\n    const evaluator = MXCubicAlongArcEvaluator.create(nominalL1, nominalR1);\r\n    if (evaluator === undefined)\r\n      return undefined;\r\n    return new DirectSpiral3d(\r\n      localToWorld.clone(),\r\n      \"MXCubicAlongArc\",\r\n      undefined,\r\n      nominalL1, nominalR1,\r\n      activeInterval ? activeInterval.clone() : Segment1d.create(0, 1), evaluator);\r\n  }\r\n\r\n  /**\r\n   * Create a polish cubic\r\n   * This is y= m*x^3 with\r\n   * * m is 1/ (6RL)\r\n   *    * 1/(6RL) is the leading term of the sine series.\r\n   * * L is nominal length\r\n   * * R is nominal end radius.\r\n   * * x ranges up to the x axis distance for which the polish distance series produces f(x)=L\r\n   * * The support class PolishCubicEvaluator has static methods for the distance series and its inversion.\r\n   */\r\n  public static createPolishCubic(\r\n    localToWorld: Transform,\r\n    nominalL1: number,\r\n    nominalR1: number,\r\n    activeInterval?: Segment1d): DirectSpiral3d | undefined {\r\n    const evaluator = PolishCubicEvaluator.create(nominalL1, nominalR1);\r\n    if (evaluator === undefined)\r\n      return undefined;\r\n    return new DirectSpiral3d(\r\n      localToWorld.clone(),\r\n      \"PolishCubic\",\r\n      undefined,\r\n      nominalL1, nominalR1,\r\n      activeInterval ? activeInterval.clone() : Segment1d.create(0, 1), evaluator);\r\n  }\r\n\r\n  /**\r\n   * Create an AustralianRailCorp spiral\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 up a little bit to simulate the lost series terms.\r\n   * @param localToWorld\r\n   * @param nominalL1\r\n   * @param nominalR1\r\n   * @param activeInterval\r\n   */\r\n  public static createAustralianRail(\r\n    localToWorld: Transform,\r\n    nominalL1: number,\r\n    nominalR1: number,\r\n    activeInterval?: Segment1d): DirectSpiral3d | undefined {\r\n    const evaluator = AustralianRailCorpXYEvaluator.create(nominalL1, nominalR1);\r\n    if (evaluator === undefined)\r\n      return undefined;\r\n    return new DirectSpiral3d(\r\n      localToWorld.clone(),\r\n      \"AustralianRailCorp\",\r\n      undefined,\r\n      nominalL1, nominalR1,\r\n      activeInterval ? activeInterval.clone() : Segment1d.create(0, 1), evaluator);\r\n  }\r\n\r\n  public static createDirectHalfCosine(\r\n    localToWorld: Transform,\r\n    nominalL1: number,\r\n    nominalR1: number,\r\n    activeInterval?: Segment1d): DirectSpiral3d | undefined {\r\n    return new this(localToWorld, \"HalfCosine\", undefined, nominalL1, nominalR1, activeInterval,\r\n      new DirectHalfCosineSpiralEvaluator(nominalL1, nominalR1));\r\n  }\r\n  /**\r\n   * Create an Arema spiral clothoid approximation\r\n   *   * X is 2 terms of the clothoid series as a function of nominal distance along\r\n   *   * Y is 2 terms f the clothoid series as a function of nominal distance along\r\n   *   * Remark: This is identical to the ChineseCubic\r\n   * @param localToWorld axes with inflection at origin, tangent along x axis\r\n   * @param nominalL1 nominal length as used in series LR terms.\r\n   * @param nominalR1 nominal final radius as used in series LR terms\r\n   * @param activeInterval fractional interval with (0, nominalL1) range for nominal distance along\r\n   */\r\n  public static createArema(\r\n    localToWorld: Transform,\r\n    nominalL1: number,\r\n    nominalR1: number,\r\n    activeInterval?: Segment1d): DirectSpiral3d | undefined {\r\n    return this.createTruncatedClothoid(\"Arema\", localToWorld, 2, 2, undefined, nominalL1, nominalR1, activeInterval);\r\n  }\r\n\r\n  /**\r\n   * Create a Chinese clothoid approximation\r\n   *   * X is 2 terms of the clothoid series as a function of nominal distance along\r\n   *   * Y is 2 terms f the clothoid series as a function of nominal distance along\r\n   *   * Remark: This is identical to the Arema spiral\r\n   * @param localToWorld axes with inflection at origin, tangent along x axis\r\n   * @param nominalL1 nominal length as used in series LR terms.\r\n   * @param nominalR1 nominal final radius as used in series LR terms\r\n   * @param activeInterval fractional interval with (0, nominalL1) range for nominal distance along\r\n   */\r\n  public static createChineseCubic(\r\n    localToWorld: Transform,\r\n    nominalL1: number,\r\n    nominalR1: number,\r\n    activeInterval?: Segment1d): DirectSpiral3d | undefined {\r\n    return this.createTruncatedClothoid(\"ChineseCubic\", localToWorld, 2, 2, undefined, nominalL1, nominalR1, activeInterval);\r\n  }\r\n  /**\r\n   * Create a Western Australian direct spiral.\r\n   *   * X is 2 terms of the clothoid series as a function of distance along\r\n   *   * Y is 1 term (cubic in nominal distance along)\r\n   * @param localToWorld axes with inflection at origin, tangent along x axis\r\n   * @param nominalL1 nominal length as used in series LR terms.\r\n   * @param nominalR1 nominal final radius as used in series LR terms\r\n   * @param activeInterval fractional interval with (0, nominalL1) range for nominal distance along\r\n   */\r\n  public static createWesternAustralian(\r\n    localToWorld: Transform,\r\n    nominalL1: number,\r\n    nominalR1: number,\r\n    activeInterval?: Segment1d): DirectSpiral3d | undefined {\r\n    return this.createTruncatedClothoid(\"WesternAustralian\", localToWorld, 2, 1, undefined, nominalL1, nominalR1, activeInterval);\r\n  }\r\n  /**\r\n   * Create (if possible) a DirectSpiral3d, applying various strict conditions appropriate to the spiral type.\r\n   * The parameter list includes extraneous values in order to directly match IntegratedSpiral3d.create, which has greater flexibility about\r\n   *    mixtures of values.\r\n   * * IMPORTANT RESTRICTIONS\r\n   *   * Direct spirals must have the inflection at the origin of their coordinate system, aligned with the x axis.\r\n   *      * hence bearing0 = 0\r\n   *      * hence radius0 = 0\r\n   *   * bearing1 is ignored\r\n   *   * radius1 must be given.\r\n   *   * arcLength must be given,\r\n   * @param spiralType one of the types in `DirectSpiralTypeNames`\r\n   * @param radius0 radius (or 0 for tangent to line) at start.   Must be ZERO or UNDEFINED\r\n   * @param radius1 radius (or 0 for tangent to line) at end.\r\n   * @param bearing0 bearing, measured CCW from x axis at start.   Must be ZERO or UNDEFINED\r\n   * @param bearing1 bearing, measured CCW from x axis at end.    IGNORED.\r\n   * @param fractionInterval optional fractional interval for an \"active\" portion of the curve.   if omitted, the full [0,1] is used.\r\n   * @param localToWorld placement transform\r\n   */\r\n  public static createFromLengthAndRadius(\r\n    spiralType: string,\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    activeInterval: undefined | Segment1d,\r\n    localToWorld: Transform): TransitionSpiral3d | undefined {\r\n    if (bearing0 !== undefined && !bearing0.isAlmostZero)\r\n      return undefined;\r\n    if (radius0 !== undefined && !Geometry.isSmallMetricDistance(radius0))\r\n      return undefined;\r\n    if (radius1 === undefined || Geometry.isSmallMetricDistance(radius1))\r\n      return undefined;\r\n    if (arcLength === undefined)\r\n      return undefined;\r\n    if (Geometry.equalStringNoCase(spiralType, \"Arema\"))\r\n      return this.createArema(localToWorld, arcLength, radius1, activeInterval);\r\n    if (Geometry.equalStringNoCase(spiralType, \"ChineseCubic\"))\r\n      return this.createChineseCubic(localToWorld, arcLength, radius1, activeInterval);\r\n    if (Geometry.equalStringNoCase(spiralType, \"JapaneseCubic\"))\r\n      return this.createJapaneseCubic(localToWorld, arcLength, radius1, activeInterval);\r\n    if (Geometry.equalStringNoCase(spiralType, \"HalfCosine\"))\r\n      return this.createDirectHalfCosine(localToWorld, arcLength, radius1, activeInterval);\r\n    if (Geometry.equalStringNoCase(spiralType, \"Czech\"))\r\n      return this.createCzechCubic(localToWorld, arcLength, radius1, activeInterval);\r\n    if (Geometry.equalStringNoCase(spiralType, \"Italian\"))\r\n      return this.createItalian(localToWorld, arcLength, radius1, activeInterval);\r\n    if (Geometry.equalStringNoCase(spiralType, \"AustralianRailCorp\"))\r\n      return this.createAustralianRail(localToWorld, arcLength, radius1, activeInterval);\r\n    if (Geometry.equalStringNoCase(spiralType, \"MXCubicAlongArc\"))\r\n      return this.createMXCubicAlongArc(localToWorld, arcLength, radius1, activeInterval);\r\n    if (Geometry.equalStringNoCase(spiralType, \"WesternAustralian\"))\r\n      return this.createWesternAustralian(localToWorld, arcLength, radius1, activeInterval);\r\n    if (Geometry.equalStringNoCase(spiralType, \"PolishCubic\"))\r\n      return this.createPolishCubic(localToWorld, arcLength, radius1, activeInterval);\r\n    return undefined;\r\n  }\r\n  /** Deep clone of this spiral */\r\n  public clone(): DirectSpiral3d {\r\n    return new DirectSpiral3d(\r\n      this.localToWorld.clone(),\r\n      this._spiralType,\r\n      this.designProperties?.clone(),\r\n      this._nominalL1,\r\n      this._nominalR1,\r\n      this._activeFractionInterval?.clone(),\r\n      this._evaluator.clone());\r\n  }\r\n  /** apply `transform` to this spiral's local to world transform. */\r\n  public tryTransformInPlace(transformA: Transform): boolean {\r\n    const rigidData = this.applyRigidPartOfTransform(transformA);\r\n    if (rigidData !== undefined) {\r\n      this._nominalL1 *= rigidData.scale;\r\n      this._nominalR1 *= rigidData.scale;\r\n      this.evaluator.scaleInPlace(rigidData.scale);\r\n    }\r\n    this.refreshComputedProperties();\r\n    return true;\r\n  }\r\n  /** Clone with a transform applied  */\r\n  public cloneTransformed(transform: Transform): DirectSpiral3d {\r\n    const result = this.clone();\r\n    result.tryTransformInPlace(transform); // ok, we're confident it will always work.\r\n    return result;\r\n  }\r\n  /** Return the spiral start point. */\r\n  public override startPoint(): Point3d { return this.activeStrokes.startPoint(); }\r\n  /** return the spiral end point. */\r\n  public override endPoint(): Point3d { return this.activeStrokes.endPoint(); }\r\n  /** test if the local to world transform places the spiral xy plane into `plane` */\r\n  public isInPlane(plane: Plane3dByOriginAndUnitNormal): boolean {\r\n    return plane.isPointInPlane(this.localToWorld.origin as Point3d)\r\n      && Geometry.isSameCoordinate(0.0, this.localToWorld.matrix.dotColumnX(plane.getNormalRef()))\r\n      && Geometry.isSameCoordinate(0.0, this.localToWorld.matrix.dotColumnY(plane.getNormalRef()));\r\n  }\r\n  /** Return quick length of the spiral.\r\n   * The tangent vector of a true clothoid is length 1 everywhere, so simple proportion of nominalL1 is a good approximation.\r\n   */\r\n  public quickLength() {\r\n    const distanceData = this._globalStrokes.packedUVParams!;\r\n    const n = distanceData.length;\r\n    return distanceData.getYAtUncheckedPointIndex(n - 1);\r\n  }\r\n  /** Return length of the spiral.\r\n   * * True length is stored at back of uvParams . . .\r\n   */\r\n  public override curveLength() { return this.quickLength(); }\r\n  /** Test if `other` is an instance of `TransitionSpiral3d` */\r\n  public isSameGeometryClass(other: any): boolean { return other instanceof DirectSpiral3d; }\r\n  /** Add strokes from this spiral to `dest`.\r\n   * * Linestrings will usually stroke as just their points.\r\n   * * If maxEdgeLength is given, this will sub-stroke within the linestring -- not what we want.\r\n   */\r\n  public emitStrokes(dest: LineString3d, options?: StrokeOptions): void { this.activeStrokes.emitStrokes(dest, options); }\r\n  /** emit stroke fragments to `dest` handler. */\r\n  public emitStrokableParts(dest: IStrokeHandler, options?: StrokeOptions): void {\r\n    const n = this.computeStrokeCountForOptions(options);\r\n    const activeStrokes = this.activeStrokes;\r\n    dest.startParentCurvePrimitive(this);\r\n    if (n <= activeStrokes.numPoints()) {\r\n      // this.activeStrokes.emitStrokableParts(dest, options);\r\n      dest.announceIntervalForUniformStepStrokes(this, 2 * activeStrokes.numPoints(), 0.0, 1.0);\r\n    } else {\r\n      dest.announceIntervalForUniformStepStrokes(this, n, 0.0, 1.0);\r\n    }\r\n    dest.endParentCurvePrimitive(this);\r\n  }\r\n\r\n  /**\r\n   * return the stroke count required for given options.\r\n   * @param options StrokeOptions that determine count\r\n   */\r\n\r\n  public computeStrokeCountForOptions(options?: StrokeOptions): number {\r\n    let numStroke;\r\n    const nominalRadians = this._nominalL1 / (2.0 * this._nominalR1);\r\n    if (options) {\r\n      const rMin = Math.abs(this._nominalR1);\r\n      numStroke = options.applyTolerancesToArc(rMin, nominalRadians);\r\n      numStroke = options.applyMaxEdgeLength(numStroke, this.quickLength());\r\n      numStroke = options.applyMinStrokesPerPrimitive(numStroke);\r\n    } else {\r\n      numStroke = StrokeOptions.applyAngleTol(undefined, 4, nominalRadians, 0.02);\r\n    }\r\n    numStroke = Math.ceil(this._activeFractionInterval.absoluteDelta() * numStroke);\r\n    return numStroke;\r\n  }\r\n\r\n  /** Reverse the active interval and active strokes.\r\n   * * Primary defining data remains unchanged !!!\r\n   */\r\n\r\n  public reverseInPlace(): void {\r\n    this._activeFractionInterval.reverseInPlace();\r\n    if (this._activeStrokes === undefined)\r\n      this._activeStrokes = this._globalStrokes.clone();\r\n    this._activeStrokes.reverseInPlace();\r\n  }\r\n  /** Evaluate curve point with respect to fraction. */\r\n  public fractionToPoint(activeFraction: number, result?: Point3d): Point3d {\r\n    const globalFraction = this._activeFractionInterval.fractionToPoint(activeFraction);\r\n    result = this._evaluator.fractionToPoint(globalFraction, result);\r\n    this.localToWorld.multiplyPoint3d(result, result);\r\n    return result;\r\n  }\r\n  /** Evaluate curve point and derivative with respect to fraction. */\r\n  public fractionToPointAndDerivative(activeFraction: number, result?: Ray3d): Ray3d {\r\n    const globalFraction = this._activeFractionInterval.fractionToPoint(activeFraction);\r\n    result = this._evaluator.fractionToPointAndDerivative(globalFraction, result);\r\n    result.direction.scaleInPlace(this._activeFractionInterval.signedDelta());\r\n    result.transformInPlace(this.localToWorld);\r\n    return result;\r\n  }\r\n\r\n  /** Return a plane with\r\n   *\r\n   * * origin at fractional position along the curve\r\n   * * vectorU is the first derivative, i.e. tangent vector with length equal to the rate of change with respect to the fraction.\r\n   * * vectorV is the second derivative, i.e.derivative of vectorU.\r\n   */\r\n\r\n  public fractionToPointAnd2Derivatives(activeFraction: number, result?: Plane3dByOriginAndVectors): Plane3dByOriginAndVectors | undefined {\r\n    const globalFraction = this._activeFractionInterval.fractionToPoint(activeFraction);\r\n    result = this._evaluator.fractionToPointAnd2Derivatives(globalFraction, result);\r\n    const a = this._activeFractionInterval.signedDelta();\r\n    result.vectorU.scaleInPlace(a);\r\n    result.vectorV.scaleInPlace(a * a);\r\n    result.transformInPlace(this.localToWorld);\r\n    return result;\r\n  }\r\n  /** Second step of double dispatch:  call `handler.handleTransitionSpiral(this)` */\r\n  public dispatchToGeometryHandler(handler: GeometryHandler): any {\r\n    return handler.handleTransitionSpiral(this);\r\n  }\r\n  /** extend the range by the strokes of the spiral */\r\n  public extendRange(rangeToExtend: Range3d, transform?: Transform): void {\r\n    this.activeStrokes.extendRange(rangeToExtend, transform);\r\n  }\r\n  /** compare various coordinate quantities */\r\n  public override isAlmostEqual(other: any): boolean {\r\n    if (other instanceof DirectSpiral3d) {\r\n      return Geometry.isSameCoordinate(this._nominalL1, other._nominalL1)\r\n        && Geometry.isSameCoordinate(this._nominalR1, other._nominalR1)\r\n        && this.localToWorld.isAlmostEqual(other.localToWorld)\r\n        && this._activeFractionInterval.isAlmostEqual(other._activeFractionInterval)\r\n        && this._evaluator.isAlmostEqual(other._evaluator);\r\n    }\r\n    return false;\r\n  }\r\n}\r\n"]}