@itwin/core-geometry 5.0.0-dev.9 → 5.0.0-dev.90

package/lib/cjs/curve/CurveFactory.js

@@ -8,7 +8,6 @@ exports.CurveFactory = exports.MiteredSweepOutputSelect = void 0;
 /** @packageDocumentation
  * @module Curve
  */
-// import { Geometry, Angle, AngleSweep } from "../Geometry";
 const Geometry_1 = require("../Geometry");
 const Angle_1 = require("../geometry3d/Angle");
 const AngleSweep_1 = require("../geometry3d/AngleSweep");
@@ -33,6 +32,7 @@ const LineSegment3d_1 = require("./LineSegment3d");
 const LineString3d_1 = require("./LineString3d");
 const Loop_1 = require("./Loop");
 const Path_1 = require("./Path");
+const RegionOps_1 = require("./RegionOps");
 const IntegratedSpiral3d_1 = require("./spiral/IntegratedSpiral3d");
 /**
  * Enumeration of geometric output for [CurveFactory.createMiteredSweepSections].
@@ -47,12 +47,13 @@ var MiteredSweepOutputSelect;
     /** Output planes and sections, as well as the assembled ruled sweep and its stroked mesh. */
     MiteredSweepOutputSelect[MiteredSweepOutputSelect["AlsoMesh"] = 2] = "AlsoMesh";
 })(MiteredSweepOutputSelect || (exports.MiteredSweepOutputSelect = MiteredSweepOutputSelect = {}));
+;
 /**
  * The `CurveFactory` class contains methods for specialized curve constructions.
  * @public
  */
 class CurveFactory {
-    /** (cautiously) construct and save a line segment between fractional positions. */
+    /** (Cautiously) construct and save a line segment between fractional positions. */
     static addPartialSegment(path, allowBackup, pointA, pointB, fraction0, fraction1) {
         if (allowBackup || (fraction1 > fraction0)) {
             if (pointA !== undefined && pointB !== undefined && !Geometry_1.Geometry.isAlmostEqualNumber(fraction0, fraction1))
@@ -149,10 +150,11 @@ class CurveFactory {
         }
         return path;
     }
-    /** Create a `Loop` with given xy corners and fixed z.
+    /**
+     * Create a `Loop` with given xy corners and fixed z.
      * * The corners always proceed counter clockwise from lower left.
      * * If the radius is too large for the outer rectangle size, it is reduced to half of the the smaller x or y size.
-    */
+     */
     static createRectangleXY(x0, y0, x1, y1, z = 0, filletRadius) {
         let radius = Geometry_1.Geometry.correctSmallMetricDistance(filletRadius);
         const xMin = Math.min(x0, x1);
@@ -161,7 +163,13 @@ class CurveFactory {
         const yMax = Math.max(y0, y1);
         radius = Math.min(Math.abs(radius), 0.5 * (xMax - xMin), 0.5 * (yMax - yMin));
         if (radius === 0.0)
-            return Loop_1.Loop.createPolygon([Point3dVector3d_1.Point3d.create(xMin, yMin, z), Point3dVector3d_1.Point3d.create(xMax, yMin, z), Point3dVector3d_1.Point3d.create(xMax, yMax, z), Point3dVector3d_1.Point3d.create(xMin, yMax, z), Point3dVector3d_1.Point3d.create(xMin, yMin, z)]);
+            return Loop_1.Loop.createPolygon([
+                Point3dVector3d_1.Point3d.create(xMin, yMin, z),
+                Point3dVector3d_1.Point3d.create(xMax, yMin, z),
+                Point3dVector3d_1.Point3d.create(xMax, yMax, z),
+                Point3dVector3d_1.Point3d.create(xMin, yMax, z),
+                Point3dVector3d_1.Point3d.create(xMin, yMin, z),
+            ]);
         else {
             const vectorU = Point3dVector3d_1.Vector3d.create(radius, 0, 0);
             const vectorV = Point3dVector3d_1.Vector3d.create(0, radius, 0);
@@ -169,7 +177,12 @@ class CurveFactory {
             const y0A = yMin + radius;
             const x1A = xMax - radius;
             const y1A = yMax - radius;
-            const centers = [Point3dVector3d_1.Point3d.create(x1A, y1A, z), Point3dVector3d_1.Point3d.create(x0A, y1A, z), Point3dVector3d_1.Point3d.create(x0A, y0A, z), Point3dVector3d_1.Point3d.create(x1A, y0A, z)];
+            const centers = [
+                Point3dVector3d_1.Point3d.create(x1A, y1A, z),
+                Point3dVector3d_1.Point3d.create(x0A, y1A, z),
+                Point3dVector3d_1.Point3d.create(x0A, y0A, z),
+                Point3dVector3d_1.Point3d.create(x1A, y0A, z),
+            ];
             const loop = Loop_1.Loop.create();
             for (let i = 0; i < 4; i++) {
                 const center = centers[i];
@@ -189,10 +202,10 @@ class CurveFactory {
     /**
      * If `arcB` is a continuation of `arcA`, extend `arcA` (in place) to include the range of `arcB`
      * * This only succeeds if the two arcs are part of identical complete arcs and end of `arcA` matches the beginning of `arcB`.
-     * @param arcA first arc, modified in place
-     * @param arcB second arc, unmodified
-     * @param allowReversed whether to consolidate even when second arc is reversed
-     * @returns whether `arcA` was modified
+     * @param arcA first arc, modified in place.
+     * @param arcB second arc, unmodified.
+     * @param allowReversed whether to consolidate even when second arc is reversed.
+     * @returns whether `arcA` was modified.
      */
     static appendToArcInPlace(arcA, arcB, allowReverse = false) {
         if (arcA.center.isAlmostEqual(arcB.center)) {
@@ -208,7 +221,7 @@ class CurveFactory {
                 arcA.sweep.setStartEndRadians(arcA.sweep.startRadians, arcA.sweep.startRadians + arcA.sweep.sweepRadians + sweepSign * arcB.sweep.sweepRadians);
                 return true;
             }
-            // Also ok if negated tangent . ..
+            // Also ok if negated tangent
             if (allowReverse) {
                 startB.direction.scaleInPlace(-1.0);
                 if (endA.isAlmostEqual(startB)) {
@@ -221,7 +234,8 @@ class CurveFactory {
     }
     /**
      * Return a `Path` containing arcs are on the surface of an ellipsoid and pass through a sequence of points.
-     * * Each arc passes through the two given endpoints and in the plane containing the true surface normal at given `fractionForIntermediateNormal`
+     * * Each arc passes through the two given endpoints and in the plane containing the true surface normal at
+     * given `fractionForIntermediateNormal`
      * @param ellipsoid
      * @param pathPoints
      * @param fractionForIntermediateNormal fractional position for surface normal used to create the section plane.
@@ -244,7 +258,7 @@ class CurveFactory {
     }
     /**
      * Create solid primitives for pipe segments (e.g. Cone or TorusPipe) around line and arc primitives.
-     * @param centerline centerline geometry/
+     * @param centerline centerline geometry.
      * @param pipeRadius radius of pipe.
      */
     static createPipeSegments(centerline, pipeRadius) {
@@ -269,122 +283,201 @@ class CurveFactory {
         }
         return undefined;
     }
-    /**
-     * * Create section arcs for mitered pipe.
-     * * At each end of each pipe, the pipe is cut by the plane that bisects the angle between successive pipe centerlines.
-     * * The arc definitions are constructed so that lines between corresponding fractional positions on the arcs are
-     *     axial lines on the pipes.
-     * * This means that each arc definition axes (aka vector0 and vector90) are _not_ perpendicular to each other.
-     * * Circular or elliptical pipe cross sections can be specified by supplying either a radius, a pair of semi-axis lengths, or a full Arc3d.
-     *    * For semi-axis length input, x corresponds to an ellipse local axis nominally situated parallel to the xy-plane.
-     *    * The center of Arc3d input is translated to the centerline start point to act as initial cross section.
-     * @param centerline centerline of pipe
-     * @param sectionData circle radius, ellipse semi-axis lengths, or full Arc3d
-     */
-    static createMiteredPipeSections(centerline, sectionData) {
-        const arcs = [];
-        if (centerline.length < 2)
-            return [];
-        const vector0 = Point3dVector3d_1.Vector3d.create();
-        const vector90 = Point3dVector3d_1.Vector3d.create();
-        const vectorBC = Point3dVector3d_1.Vector3d.create();
-        const currentCenter = Point3dVector3d_1.Point3d.create();
-        centerline.vectorIndexIndex(0, 1, vectorBC);
-        centerline.getPoint3dAtUncheckedPointIndex(0, currentCenter);
-        let initialSection;
+    /** Get start point and tangent for variant curve data. */
+    static startPointAndTangent(curve) {
+        if (curve instanceof CurvePrimitive_1.CurvePrimitive)
+            return curve.fractionToPointAndDerivative(0.0);
+        if (curve.length < 2)
+            return undefined;
+        if (Array.isArray(curve))
+            curve = new Point3dArrayCarrier_1.Point3dArrayCarrier(curve);
+        const ray = Ray3d_1.Ray3d.createZero();
+        curve.getPoint3dAtUncheckedPointIndex(0, ray.origin);
+        curve.vectorIndexIndex(0, 1, ray.direction);
+        return ray;
+    }
+    ;
+    /** Create an [[Arc3d]] from `sectionData` that has its center at the start point of the centerline. */
+    static createArcFromSectionData(centerline, sectionData) {
+        const ray = CurveFactory.startPointAndTangent(centerline);
+        if (!ray)
+            return undefined;
+        let arc;
         if (sectionData instanceof Arc3d_1.Arc3d) {
-            initialSection = sectionData.clone();
-            initialSection.center.setFrom(currentCenter);
-            vector0.setFrom(sectionData.vector0);
-            vector90.setFrom(sectionData.vector90);
+            arc = sectionData.clone();
+            arc.center = ray.origin;
         }
-        else if (typeof sectionData === "number" || Point3dVector3d_1.Point3d.isXAndY(sectionData)) {
+        else {
+            const vector0 = Point3dVector3d_1.Vector3d.create();
+            const vector90 = Point3dVector3d_1.Vector3d.create();
             const length0 = (typeof sectionData === "number") ? sectionData : sectionData.x;
             const length90 = (typeof sectionData === "number") ? sectionData : sectionData.y;
-            const baseFrame = Matrix3d_1.Matrix3d.createRigidHeadsUp(vectorBC, Geometry_1.AxisOrder.ZXY);
+            const baseFrame = Matrix3d_1.Matrix3d.createRigidHeadsUp(ray.direction, Geometry_1.AxisOrder.ZXY);
             baseFrame.columnX(vector0).scaleInPlace(length0);
             baseFrame.columnY(vector90).scaleInPlace(length90);
-            initialSection = Arc3d_1.Arc3d.create(currentCenter, vector0, vector90, AngleSweep_1.AngleSweep.create360());
+            arc = Arc3d_1.Arc3d.create(ray.origin, vector0, vector90);
         }
-        else {
+        if (arc.binormalVector().dotProduct(ray.direction) < 0.0)
+            arc.reverseInPlace();
+        return arc;
+    }
+    /**
+     * Create section arcs for mitered pipe.
+     * * At the end of each pipe segment, the pipe is mitered by the plane that bisects the angle between successive
+     * centerline segments.
+     * * The section arcs are constructed so that lines between corresponding fractional positions on the arcs are
+     *   axial lines on the pipes.
+     * * This means that the initial arc's vector0 and vector90 lengths and angular separation are _not_ preserved in
+     * the section arcs.
+     * * Circular or elliptical pipe cross sections can be specified by supplying either a radius, a pair of semi-axis
+     * lengths, or an Arc3d:
+     *    * For semi-axis length input, x and y correspond to ellipse local axes perpendicular to each other and to the
+     * start tangent.
+     *    * For Arc3d input, the center is translated to the centerline start point. For best results, ensure this arc
+     * is perpendicular to the centerline start tangent.
+     * * This function internally calls [[createMiteredSweepSections]] with default options.
+     * @param centerline centerline of pipe. For best results, ensure no successive duplicate points with e.g.,
+     * [[GrowableXYZArray.createCompressed]].
+     * @param sectionData circle radius, ellipse semi-axis lengths, or full Arc3d (if not full, function makes it full).
+     * @returns array of sections or empty array if section creation failed.
+     */
+    static createMiteredPipeSections(centerline, sectionData) {
+        const arc = CurveFactory.createArcFromSectionData(centerline, sectionData);
+        if (!arc)
             return [];
+        const miteredSweeps = CurveFactory.createMiteredSweepSections(centerline, arc);
+        if (miteredSweeps)
+            return miteredSweeps.sections;
+        return [];
+    }
+    /** For a smooth curve, stroke and return unnormalized start/end tangents. */
+    static strokeSmoothCurve(curve, options) {
+        let startTangent, endTangent;
+        if (curve instanceof CurvePrimitive_1.CurvePrimitive) {
+            startTangent = curve.fractionToPointAndDerivative(0.0).direction;
+            endTangent = curve.fractionToPointAndDerivative(1.0).direction;
+            const strokes = LineString3d_1.LineString3d.create();
+            curve.emitStrokes(strokes, options);
+            curve = strokes.packedPoints;
         }
-        arcs.push(initialSection);
-        const vectorAB = Point3dVector3d_1.Vector3d.create();
-        const bisector = Point3dVector3d_1.Vector3d.create();
-        for (let i = 1; i < centerline.length; i++) {
-            vectorAB.setFromVector3d(vectorBC);
-            centerline.getPoint3dAtUncheckedPointIndex(i, currentCenter);
-            if (i + 1 < centerline.length) {
-                centerline.vectorIndexIndex(i, i + 1, vectorBC);
-            }
-            else {
-                vectorBC.setFromVector3d(vectorAB);
+        else if (curve instanceof CurveCollection_1.CurveChain) {
+            startTangent = curve.startPointAndDerivative()?.direction;
+            endTangent = curve.endPointAndDerivative()?.direction;
+            const strokes = curve.getPackedStrokes(options);
+            if (!strokes)
+                return undefined;
+            curve = strokes;
+        }
+        else if (Array.isArray(curve))
+            curve = new Point3dArrayCarrier_1.Point3dArrayCarrier(curve);
+        return { strokes: curve, startTangent, endTangent };
+    }
+    /**
+     * Align two bisector plane normals to given smooth rail tangents or optional overrides.
+     * * Optionally average the normals for physically closed rail.
+     */
+    static alignFirstAndLastBisectorPlanes(firstPlane, lastPlane, smoothRailData, options) {
+        const normal0 = options?.startTangent ?? (smoothRailData?.startTangent ?? firstPlane.getNormalRef());
+        const normal1 = options?.endTangent ?? (smoothRailData?.endTangent ?? lastPlane.getNormalRef());
+        if (options?.wrapIfPhysicallyClosed && firstPlane.getOriginRef().isAlmostEqual(lastPlane.getOriginRef())) {
+            const avgNormal = normal0.plus(normal1);
+            if (avgNormal.tryNormalizeInPlace()) { // ignore cusp at seam
+                firstPlane.getNormalRef().setFrom(avgNormal);
+                lastPlane.getNormalRef().setFrom(avgNormal);
+                return;
             }
-            if (vectorAB.normalizeInPlace() && vectorBC.normalizeInPlace()) {
-                vectorAB.interpolate(0.5, vectorBC, bisector);
-                // On the end ellipse for this pipe section. ..
-                // center comes directly from centerline[i]
-                // vector0 and vector90 are obtained by sweeping the corresponding vectors of the start ellipse to the split plane.
-                moveVectorToPlane(vector0, vectorAB, bisector, vector0);
-                moveVectorToPlane(vector90, vectorAB, bisector, vector90);
-                arcs.push(Arc3d_1.Arc3d.create(currentCenter, vector0, vector90, AngleSweep_1.AngleSweep.create360()));
+        }
+        if (normal0.tryNormalizeInPlace())
+            firstPlane.getNormalRef().setFrom(normal0);
+        if (normal1.tryNormalizeInPlace())
+            lastPlane.getNormalRef().setFrom(normal1);
+    }
+    /** Reverse a closed curve or region if necessary so that its orientation is CCW with respect to the plane normal. */
+    static alignClosedCurveToPlane(curve, planeNormal) {
+        let closedCurve;
+        if (curve instanceof CurvePrimitive_1.CurvePrimitive) {
+            if (curve.startPoint().isAlmostEqual(curve.endPoint()))
+                closedCurve = Loop_1.Loop.create(curve);
+        }
+        else if (curve.isAnyRegion())
+            closedCurve = curve;
+        if (closedCurve) {
+            // The alignment condition is equivalent to positive projected curve area computed wrt to the plane normal.
+            const toLocal = Matrix3d_1.Matrix3d.createRigidHeadsUp(planeNormal).transpose();
+            const projection = closedCurve.cloneTransformed(Transform_1.Transform.createOriginAndMatrix(undefined, toLocal));
+            if (projection) { // now we can ignore z-coords
+                const areaXY = RegionOps_1.RegionOps.computeXYArea(projection);
+                if (areaXY && areaXY < 0)
+                    curve.reverseInPlace();
             }
         }
-        return arcs;
     }
     /**
-     * Sweep the initialSection along each segment of the centerLine until it hits the bisector plane at the next vertex.
-     * * The caller should place the initialSection on a plane perpendicular to the first edge.
-     *   * This plane is commonly (but not necessarily) through the start point itself.
-     *   * If the geometry is not "on a perpendicular plane", the output geometry will still be flattened onto the various planes.
-     * * In the "open path" case (i.e when wrapIfPhysicallyClosed is false or the path does not have matched first and last points)
-     *       the first/last output plane will be at the start/end of the first/last edge and on a perpendicular plane.
-     * * In the "closed path" case, the output plane for the first and last point is the bisector of the start and end planes from the "open path" case,
-     *    and the first/last section geometry may be different from `initialSection`.
-     * * The centerline path does NOT have to be planar, however twisting effects effects will appear in the various bisector planes.
-     * @param centerline sweep path, e.g., as stroked from a smooth centerline curve
-     * @param initialSection profile curve to be swept. As noted above, this should be on a plane perpendicular to the first segment of the centerline.
-     * @param options options for computation and output
-     * @return array of sections, starting with `initialSection` projected along the first edge to the first plane.
+     * Projection to target plane, constructing sweep direction from two given planes.
+     * * If successful, push the target plane and swept section to the output arrays and return the swept section.
+     * * If unsuccessful, leave the output arrays alone and return the input section.
+     */
+    static doSweepToPlane(output, edgePlane0, edgePlane1, targetPlane, section) {
+        const sweepVector = Point3dVector3d_1.Vector3d.createStartEnd(edgePlane0.getOriginRef(), edgePlane1.getOriginRef());
+        const transform = Transform_1.Transform.createFlattenAlongVectorToPlane(sweepVector, targetPlane.getOriginRef(), targetPlane.getNormalRef());
+        if (transform === undefined)
+            return section;
+        const transformedSection = section.cloneTransformed(transform);
+        if (transformedSection === undefined)
+            return section;
+        output.planes.push(targetPlane);
+        output.sections.push(transformedSection);
+        return transformedSection;
+    }
+    /**
+     * Sweep the `initialSection` along each segment of the (stroked) `centerline` until it hits the bisector plane at
+     * the next vertex.
+     * * For best results, the caller should place `initialSection` in a plane perpendicular to the `centerline`
+     * start tangent.
+     *   * This plane is commonly (but not necessarily) through the centerline start point itself.
+     *   * To compute the sections, `initialSection` is projected in the direction of the centerline start tangent onto
+     * the first bisector plane at the centerline start. The result of this projection will be likewise projected onto
+     * the second plane, and so on in sequence.
+     * * By default, the first/last bisector plane normals are set to the centerline start/end tangents. The caller can
+     * override these with tangents supplied in `options`. If the centerline is physically closed and
+     * `options.wrapIfPhysicallyClosed` is true, the first and last plane normals are averaged and equated.
+     * * The centerline path does NOT have to be planar, however non-planarity will result in twisting of the sections
+     * in the bisector planes.
+     * @param centerline sweep path. Will be stroked if smooth.
+     * @param initialSection profile curve to be swept. As noted above, this should be on a plane perpendicular to the
+     * centerline start tangent.
+     * @param options options for computation and output.
+     * @return array of sections, formed from projecting `initialSection` successively onto the bisector planes.
      */
     static createMiteredSweepSections(centerline, initialSection, options) {
+        if (!options)
+            options = {};
+        const rail = this.strokeSmoothCurve(centerline, options.strokeOptions);
+        if (!rail)
+            return undefined;
+        const planes = PolylineOps_1.PolylineOps.createBisectorPlanesForDistinctPoints(rail.strokes);
+        if (!planes || planes.length < 2)
+            return undefined;
+        this.alignFirstAndLastBisectorPlanes(planes[0], planes[planes.length - 1], rail, options);
+        // RuledSweep facet construction assumes the contours are oriented CCW with respect to the sweep direction so that
+        // facet normals point outward. We only have to align the first contour; the rest will inherit its orientation.
+        this.alignClosedCurveToPlane(initialSection, planes[0].getNormalRef());
         const sectionData = { sections: [], planes: [] };
-        const planes = PolylineOps_1.PolylineOps.createBisectorPlanesForDistinctPoints(centerline, options.wrapIfPhysicallyClosed);
-        if (planes !== undefined && planes.length > 1) {
-            // Projection to target plane, constructing sweep direction from two given planes.
-            // If successful, push the target plane and swept section to the output arrays and return the swept section.
-            // If unsuccessful, leave the output arrays alone and return the input section.
-            const doSweepToPlane = function (edgePlane0, edgePlane1, targetPlane, section) {
-                const sweepVector = Point3dVector3d_1.Vector3d.createStartEnd(edgePlane0.getOriginRef(), edgePlane1.getOriginRef());
-                const transform = Transform_1.Transform.createFlattenAlongVectorToPlane(sweepVector, targetPlane.getOriginRef(), targetPlane.getNormalRef());
-                if (transform === undefined)
-                    return section;
-                const section1 = section.cloneTransformed(transform);
-                if (section1 === undefined)
-                    return section;
-                sectionData.planes.push(targetPlane);
-                sectionData.sections.push(section1);
-                return section1;
-            };
-            let currentSection = doSweepToPlane(planes[0], planes[1], planes[0], initialSection);
-            for (let i = 1; i < planes.length; i++) {
-                currentSection = doSweepToPlane(planes[i - 1], planes[i], planes[i], currentSection);
-            }
-            if (options.outputSelect) {
-                const ruledSweep = RuledSweep_1.RuledSweep.create(sectionData.sections, options.capped ?? false);
-                if (ruledSweep) {
-                    sectionData.ruledSweep = ruledSweep;
-                    if (MiteredSweepOutputSelect.AlsoMesh === options.outputSelect) {
-                        const builder = PolyfaceBuilder_1.PolyfaceBuilder.create(options.strokeOptions);
-                        builder.addRuledSweep(ruledSweep);
-                        sectionData.mesh = builder.claimPolyface();
-                    }
+        let currentSection = this.doSweepToPlane(sectionData, planes[0], planes[1], planes[0], initialSection);
+        for (let i = 1; i < planes.length; i++)
+            currentSection = this.doSweepToPlane(sectionData, planes[i - 1], planes[i], planes[i], currentSection);
+        if (options.outputSelect) {
+            const ruledSweep = RuledSweep_1.RuledSweep.create(sectionData.sections, options.capped ?? false);
+            if (ruledSweep) {
+                sectionData.ruledSweep = ruledSweep;
+                if (MiteredSweepOutputSelect.AlsoMesh === options.outputSelect) {
+                    const builder = PolyfaceBuilder_1.PolyfaceBuilder.create(options.strokeOptions);
+                    builder.addRuledSweep(ruledSweep);
+                    sectionData.mesh = builder.claimPolyface();
                 }
             }
-            return sectionData;
         }
-        return undefined;
+        return sectionData;
     }
     /**
      * Create a circular arc from start point, tangent at start, radius, optional plane normal, arc sweep.
@@ -401,11 +494,11 @@ class CurveFactory {
     /**
      * Compute 2 spirals (all in XY) for a symmetric line-to-line transition.
      * * First spiral begins at given start point.
-     * * first tangent aims at shoulder
+     * * first tangent aims at shoulder.
      * * outbound spiral joins line from shoulder to target.
-     * @param spiralType name of spiral type.  THIS MUST BE AN "Integrated" SPIRAL TYPE
+     * @param spiralType name of spiral type. THIS MUST BE AN "Integrated" SPIRAL TYPE.
      * @param startPoint inbound start point.
-     * @param shoulder point target point for (both) spiral-to-line tangencies
+     * @param shoulder point target point for (both) spiral-to-line tangencies.
      * @return array with the computed spirals, or undefined if failure.
      */
     static createLineSpiralSpiralLine(spiralType, startPoint, shoulderPoint, targetPoint) {
@@ -443,9 +536,9 @@ class CurveFactory {
      * Compute 2 spirals (all in XY) for a symmetric line-to-line transition.
      * * Spiral length is given.
      * * tangency points float on both lines.
-     * @param spiralType name of spiral type.  THIS MUST BE AN "Integrated" SPIRAL TYPE
+     * @param spiralType name of spiral type.  THIS MUST BE AN "Integrated" SPIRAL TYPE.
      * @param pointA inbound start point.
-     * @param shoulder point target point for (both) spiral-to-line tangencies
+     * @param shoulder point target point for (both) spiral-to-line tangencies.
      * @param spiralLength for each part of the spiral pair.
      * @return array with the computed spirals, or undefined if failure.
      */
@@ -483,12 +576,12 @@ class CurveFactory {
     }
     /**
      * Compute 2 spirals and an arc (all in XY) for a symmetric line-to-line transition.
-     * Spiral lengths and arc radius are given.   (e.g. from design speed standards.)
-     * @param spiralType name of spiral type.  THIS MUST BE AN "Integrated" SPIRAL TYPE
+     * Spiral lengths and arc radius are given (e.g., from design speed standards).
+     * @param spiralType name of spiral type. THIS MUST BE AN "Integrated" SPIRAL TYPE.
      * @param pointA inbound start point.
-     * @param pointB shoulder (target)  point for (both) spiral-to-line tangencies
-     * @param lengthA inbound spiral length
-     * @param lengthB outbound spiral length
+     * @param pointB shoulder (target)  point for (both) spiral-to-line tangencies.
+     * @param lengthA inbound spiral length.
+     * @param lengthB outbound spiral length.
      * @return array with the computed spirals, or undefined if failure.
      */
     static createLineSpiralArcSpiralLine(spiralType, pointA, pointB, pointC, lengthA, lengthB, arcRadius) {
@@ -538,12 +631,7 @@ class CurveFactory {
         }
         return undefined;
     }
-    /**
-     * Return the intersection point of 3 planes.
-     * @param planeA
-     * @param planeB
-     * @param planeC
-     */
+    /** Return the intersection point of 3 planes. */
     static planePlaneIntersectionRay(planeA, planeB) {
         const altitudeA = planeA.altitudeXYZ(0, 0, 0);
         const altitudeB = planeB.altitudeXYZ(0, 0, 0);
@@ -564,14 +652,4 @@ class CurveFactory {
     }
 }
 exports.CurveFactory = CurveFactory;
-/**
- * Starting at vectorR, move parallel to vectorV until perpendicular to planeNormal
- */
-function moveVectorToPlane(vectorR, vectorV, planeNormal, result) {
-    // find s such that (vectorR + s * vectorV) DOT planeNormal = 0.
-    const dotRN = vectorR.dotProduct(planeNormal);
-    const dotVN = vectorV.dotProduct(planeNormal);
-    const s = Geometry_1.Geometry.safeDivideFraction(dotRN, dotVN, 0.0);
-    return vectorR.plusScaled(vectorV, -s, result);
-}
 //# sourceMappingURL=CurveFactory.js.map