@itwin/core-geometry 4.10.0-dev.9 → 5.0.0-dev.0

package/lib/cjs/curve/Arc3d.js

@@ -15,6 +15,7 @@ const Point3dVector3d_1 = require("../geometry3d/Point3dVector3d");
 const Range_1 = require("../geometry3d/Range");
 const Ray3d_1 = require("../geometry3d/Ray3d");
 const Polynomials_1 = require("../numerics/Polynomials");
+const SmallSystem_1 = require("../numerics/SmallSystem");
 const CurveExtendMode_1 = require("./CurveExtendMode");
 const CurveLocationDetail_1 = require("./CurveLocationDetail");
 const CurvePrimitive_1 = require("./CurvePrimitive");
@@ -154,15 +155,15 @@ class Arc3d extends CurvePrimitive_1.CurvePrimitive {
     isSameGeometryClass(other) {
         return other instanceof Arc3d;
     }
-    /** Read property for (clone of) center */
+    /** Read property for (clone of) center. */
     get center() {
         return this._center.clone();
     }
-    /** Read property for (clone of) vector0 */
+    /** Read property for (clone of) vector0. */
     get vector0() {
         return this._matrix.columnX();
     }
-    /** Read property for (clone of) vector90 */
+    /** Read property for (clone of) vector90. */
     get vector90() {
         return this._matrix.columnY();
     }
@@ -170,11 +171,11 @@ class Arc3d extends CurvePrimitive_1.CurvePrimitive {
     get perpendicularVector() {
         return this._matrix.columnZ();
     }
-    /** Read property for (clone of) matrix of vector0, vector90, unit normal */
+    /** Read property for (clone of) matrix of vector0, vector90, unit normal. */
     matrixClone() {
         return this._matrix.clone();
     }
-    /** Read property for (reference to) matrix of vector0, vector90, unit normal */
+    /** Read property for (reference to) matrix of vector0, vector90, unit normal. */
     get matrixRef() {
         return this._matrix;
     }
@@ -185,21 +186,21 @@ class Arc3d extends CurvePrimitive_1.CurvePrimitive {
     set sweep(value) {
         this._sweep.setFrom(value);
     }
-    /** An Arc3d extends along its complete elliptic arc */
+    /** An Arc3d extends along its complete elliptic arc. */
     get isExtensibleFractionSpace() {
         return true;
     }
     // constructor copies the pointers
     constructor(center, matrix, sweep) {
         super();
-        /** String name for schema properties */
+        /** String name for schema properties. */
         this.curvePrimitiveType = "arc";
         this._center = center;
         this._matrix = matrix;
         this._sweep = sweep;
     }
     /**
-     * Return a clone of the arc, with transform applied
+     * Return a clone of the arc, with transform applied.
      * @param transform
      */
     cloneTransformed(transform) {
@@ -209,9 +210,9 @@ class Arc3d extends CurvePrimitive_1.CurvePrimitive {
     }
     /**
      * Redefine the arc with (captured references to) given data.
-     * @param center arc center
-     * @param matrix matrix with columns vector0, vector 90, and their unit cross product
-     * @param sweep angle sweep
+     * @param center arc center.
+     * @param matrix matrix with columns vector0, vector 90, and their unit cross product.
+     * @param sweep angle sweep.
      */
     setRefs(center, matrix, sweep) {
         this._center = center;
@@ -220,14 +221,14 @@ class Arc3d extends CurvePrimitive_1.CurvePrimitive {
     }
     /**
      * Redefine the arc with (clones of) given data.
-     * @param center arc center
-     * @param matrix matrix with columns vector0, vector 90, and their unit cross product
-     * @param sweep angle sweep
+     * @param center arc center.
+     * @param matrix matrix with columns vector0, vector 90, and their unit cross product.
+     * @param sweep angle sweep.
      */
     set(center, matrix, sweep) {
         this.setRefs(center.clone(), matrix.clone(), sweep ? sweep.clone() : AngleSweep_1.AngleSweep.create360());
     }
-    /** Copy center, matrix, and sweep from other Arc3d */
+    /** Copy center, matrix, and sweep from other Arc3d. */
     setFrom(other) {
         this._center.setFrom(other._center);
         this._matrix.setFrom(other._matrix);
@@ -239,9 +240,9 @@ class Arc3d extends CurvePrimitive_1.CurvePrimitive {
     }
     /**
      * Create an arc, capturing references to center, matrix and sweep.
-     * @param center center point
-     * @param matrix matrix with columns vector0, vector90, and unit cross product
-     * @param sweep sweep limits
+     * @param center center point.
+     * @param matrix matrix with columns vector0, vector90, and unit cross product.
+     * @param sweep sweep limits.
      * @param result optional preallocated result.
      */
     static createRefs(center, matrix, sweep, result) {
@@ -253,11 +254,11 @@ class Arc3d extends CurvePrimitive_1.CurvePrimitive {
     }
     /**
      * Create an arc from center, x column to be scaled, and y column to be scaled.
-     * @param center center of ellipse
-     * @param matrix matrix whose x and y columns are unit vectors to be scaled by radius0 and radius90
+     * @param center center of ellipse.
+     * @param matrix matrix whose x and y columns are scaled by radius0 and radius90.
      * @param radius0 radius in x direction.
      * @param radius90 radius in y direction.
-     * @param sweep sweep limits
+     * @param sweep sweep limits.
      * @param result optional preallocated result.
      */
     static createScaledXYColumns(center, matrix, radius0, radius90, sweep, result) {
@@ -266,9 +267,9 @@ class Arc3d extends CurvePrimitive_1.CurvePrimitive {
         return Arc3d.create(center, vector0.scale(radius0, vector0), vector90.scale(radius90, vector90), sweep, result);
     }
     /**
-     * Create a (full circular) arc from center, normal and radius
-     * @param center center of ellipse.  If undefined, center at 000
-     * @param normal normal vector
+     * Create a (full circular) arc from center, normal and radius.
+     * @param center center of ellipse. If undefined, center at 000.
+     * @param normal normal vector.
      * @param radius radius in x direction.
      * @param result optional preallocated result.
      */
@@ -278,11 +279,11 @@ class Arc3d extends CurvePrimitive_1.CurvePrimitive {
     }
     /**
      * Create an arc by center with vectors to points at 0 and 90 degrees in parameter space.
-     * @param center arc center
-     * @param vector0 vector to 0 degrees (commonly major axis)
-     * @param vector90 vector to 90 degree point (commonly minor axis)
-     * @param sweep sweep limits (default full sweep)
-     * @param result optional preallocated result
+     * @param center arc center.
+     * @param vector0 vector to 0 degrees (commonly major axis).
+     * @param vector90 vector to 90 degree point (commonly minor axis).
+     * @param sweep sweep limits; defaults to full sweep.
+     * @param result optional preallocated result.
      */
     static create(center, vector0, vector90, sweep, result) {
         const normal = vector0.unitCrossProductWithDefault(vector90, 0, 0, 0); // normal will be 000 for degenerate case
@@ -291,19 +292,19 @@ class Arc3d extends CurvePrimitive_1.CurvePrimitive {
     }
     /**
      * Create an elliptical arc from three points on the ellipse: two points on an axis and one in between.
-     * @param point0 start of arc, on an axis
-     * @param point1 point on arc somewhere between `point0` and `point2`
-     * @param point2 point on arc directly opposite `point0`
-     * @param sweep angular sweep, measured from `point0` in the direction of `point1`.
-     *  For a half-ellipse from `point0` to `point2` passing through `point1`, pass `AngleSweep.createStartEndDegrees(0,180)`.
+     * @param start start of arc, on an axis.
+     * @param middle point on arc somewhere between `start` and `end`.
+     * @param end point on arc directly opposite `start`.
+     * @param sweep angular sweep, measured from `start` in the direction of `middle`.
+     *  For a half-ellipse from `start` to `end` passing through `middle`, pass `AngleSweep.createStartEndDegrees(0,180)`.
      *  Default value is full sweep to create the entire ellipse.
-     * @param result optional preallocated result
+     * @param result optional preallocated result.
      * @returns elliptical arc, or undefined if construction impossible.
      */
-    static createStartMiddleEnd(point0, point1, point2, sweep, result) {
-        const center = Point3dVector3d_1.Point3d.createAdd2Scaled(point0, 0.5, point2, 0.5);
-        const vector0 = Point3dVector3d_1.Vector3d.createStartEnd(center, point0);
-        const vector1 = Point3dVector3d_1.Vector3d.createStartEnd(center, point1);
+    static createStartMiddleEnd(start, middle, end, sweep, result) {
+        const center = Point3dVector3d_1.Point3d.createAdd2Scaled(start, 0.5, end, 0.5);
+        const vector0 = Point3dVector3d_1.Vector3d.createStartEnd(center, start);
+        const vector1 = Point3dVector3d_1.Vector3d.createStartEnd(center, middle);
         const v0DotV1 = vector0.dotProduct(vector1);
         const v0Len2 = vector0.magnitudeSquared();
         if (Math.abs(v0DotV1) >= v0Len2)
@@ -311,7 +312,7 @@ class Arc3d extends CurvePrimitive_1.CurvePrimitive {
         const normal = vector0.crossProduct(vector1);
         const vector90 = normal.unitCrossProductWithDefault(vector0, 0, 0, 0);
         const v1DotV90 = vector1.dotProduct(vector90);
-        // Solve the standard ellipse equation for the unknown axis length, given local coords of point1 (v0.v1/||v0||, v90.v1)
+        // solve the standard ellipse equation for the unknown axis length, given local coords of middle (v0.v1/||v0||, v90.v1)
         const v90Len = Geometry_1.Geometry.safeDivideFraction(v0Len2 * v1DotV90, Math.sqrt(v0Len2 * v0Len2 - v0DotV1 * v0DotV1), 0);
         if (Geometry_1.Geometry.isSmallMetricDistanceSquared(v90Len))
             return undefined;
@@ -320,39 +321,61 @@ class Arc3d extends CurvePrimitive_1.CurvePrimitive {
     }
     /**
      * Create a circular arc defined by start point, tangent at start point, and end point.
-     * If tangent is parallel to line segment from start to end, return the line segment.
+     * * The circular arc is swept from `start` to `end` in the direction of `tangentAtStart`.
+     * * If `tangentAtStart` is parallel to the line segment from `start` to `end`, return the line segment.
      */
     static createCircularStartTangentEnd(start, tangentAtStart, end, result) {
-        // To find the circle passing through start and end with tangentAtStart at start:
-        //      - find line 1: the perpendicular bisector of the line from start to end.
-        //      - find line 2: the perpendicular to the tangentAtStart.
-        //      - intersection of the two lines would be the circle center.
-        const vector = Point3dVector3d_1.Vector3d.createStartEnd(start, end);
-        const normal = tangentAtStart.crossProduct(vector).normalize();
-        if (normal) {
-            const vectorPerp = normal.crossProduct(vector);
-            const tangentPerp = normal.crossProduct(tangentAtStart);
-            const midPoint = start.plusScaled(vector, 0.5);
-            const lineSeg1 = LineSegment3d_1.LineSegment3d.create(start, start.plusScaled(tangentPerp, 1));
-            const lineSeg2 = LineSegment3d_1.LineSegment3d.create(midPoint, midPoint.plusScaled(vectorPerp, 1));
-            const intersection = LineSegment3d_1.LineSegment3d.closestApproach(lineSeg1, true, lineSeg2, true);
-            if (intersection) {
-                const center = intersection.detailA.point;
-                const vector0 = Point3dVector3d_1.Vector3d.createStartEnd(center, start);
-                const vector90 = normal.crossProduct(vector0);
-                const endVector = Point3dVector3d_1.Vector3d.createStartEnd(center, end);
-                const sweep = AngleSweep_1.AngleSweep.create(vector0.signedAngleTo(endVector, normal));
+        // see itwinjs-core\core\geometry\internaldocs\Arc3d.md to clarify below algorithm
+        const startToEnd = Point3dVector3d_1.Vector3d.createStartEnd(start, end);
+        const frame = Matrix3d_1.Matrix3d.createRigidFromColumns(tangentAtStart, startToEnd, Geometry_1.AxisOrder.XYZ);
+        if (frame !== undefined) {
+            const vv = startToEnd.dotProduct(startToEnd);
+            const vw = frame.dotColumnY(startToEnd);
+            const radius = Geometry_1.Geometry.conditionalDivideCoordinate(vv, 2 * vw);
+            if (radius !== undefined) {
+                const vector0 = frame.columnY();
+                vector0.scaleInPlace(-radius); // center to start
+                const vector90 = frame.columnX();
+                vector90.scaleInPlace(radius);
+                const centerToEnd = vector0.plus(startToEnd);
+                const sweepAngle = vector0.angleTo(centerToEnd);
+                let sweepRadians = sweepAngle.radians; // always positive and less than PI
+                if (tangentAtStart.dotProduct(centerToEnd) < 0.0) // sweepRadians is the wrong way
+                    sweepRadians = 2.0 * Math.PI - sweepRadians;
+                const center = start.plusScaled(vector0, -1.0);
+                const sweep = AngleSweep_1.AngleSweep.createStartEndRadians(0.0, sweepRadians);
                 return Arc3d.create(center, vector0, vector90, sweep, result);
             }
         }
         return LineSegment3d_1.LineSegment3d.create(start, end);
     }
+    /**
+     * Create a circular arc from start point, tangent at start, radius, optional plane normal, arc sweep.
+     * * The vector from start point to center is in the direction of upVector crossed with tangentA.
+     * @param start start point.
+     * @param tangentAtStart vector in tangent direction at the start.
+     * @param radius signed radius.
+     * @param upVector optional out-of-plane vector. Defaults to positive Z.
+     * @param sweep angular range. If single `Angle` is given, start angle is at 0 degrees (the start point).
+     */
+    static createCircularStartTangentRadius(start, tangentAtStart, radius, upVector, sweep) {
+        if (upVector === undefined)
+            upVector = Point3dVector3d_1.Vector3d.unitZ();
+        const vector0 = upVector.unitCrossProduct(tangentAtStart);
+        if (vector0 === undefined)
+            return undefined;
+        const center = start.plusScaled(vector0, radius);
+        // reverse the A-to-center vector and bring it up to scale
+        vector0.scaleInPlace(-radius);
+        const vector90 = tangentAtStart.scaleToLength(Math.abs(radius)); // cannot fail; prior unitCrossProduct would have failed first
+        return Arc3d.create(center, vector0, vector90, AngleSweep_1.AngleSweep.create(sweep));
+    }
     /**
      * Create a circular arc defined by start and end points and radius.
-     * @param start start point of the arc
-     * @param end end point of the arc
-     * @param helper a third point near the arc in its plane, or a vector in the direction of the arc normal
-     * @returns the constructed arc, or undefined if desired arc cannot be constructed
+     * @param start start point of the arc.
+     * @param end end point of the arc.
+     * @param helper a third point near the arc in its plane, or a vector in the direction of the arc normal.
+     * @returns the constructed arc, or undefined if desired arc cannot be constructed.
      */
     static createCircularStartEndRadius(start, end, radius, helper) {
         // Construct a line segment from start to end. It is a chord of the circle,
@@ -384,7 +407,7 @@ class Arc3d extends CurvePrimitive_1.CurvePrimitive {
      * Return a clone of this arc, projected to given z value.
      * * If `z` is omitted, the clone is at the z of the center.
      * * This function projects the arc into a plane parallel to xy-plane.
-     * * Note that projection to fixed z can change circle into ellipse (and (rarely) ellipse to circle)
+     * * Note that projection to fixed z can change circle into ellipse (and (rarely) ellipse to circle).
      */
     cloneAtZ(z) {
         if (z === undefined)
@@ -394,7 +417,7 @@ class Arc3d extends CurvePrimitive_1.CurvePrimitive {
     /**
      * Create an arc by center (cx,cy,xz) with vectors (ux,uy,uz) and (vx,vy,vz) to points at 0 and 90 degrees in
      * parameter space.
-     * @param result optional preallocated result
+     * @param result optional preallocated result.
      */
     static createXYZXYZXYZ(cx, cy, cz, ux, uy, uz, vx, vy, vz, sweep, result) {
         return Arc3d.create(Point3dVector3d_1.Point3d.create(cx, cy, cz), Point3dVector3d_1.Vector3d.create(ux, uy, uz), Point3dVector3d_1.Vector3d.create(vx, vy, vz), sweep, result);
@@ -423,7 +446,7 @@ class Arc3d extends CurvePrimitive_1.CurvePrimitive {
         const ac2 = vectorAC.magnitudeSquared();
         const normal = vectorAB.sizedCrossProduct(vectorAC, Math.sqrt(Math.sqrt(ab2 * ac2)));
         if (normal) {
-            const vectorToCenter = Polynomials_1.SmallSystem.linearSystem3d(normal.x, normal.y, normal.z, vectorAB.x, vectorAB.y, vectorAB.z, vectorAC.x, vectorAC.y, vectorAC.z, 0, // vectorToCenter DOT normal = 0 (ensure normal is perp to the plane of the 3 points)
+            const vectorToCenter = SmallSystem_1.SmallSystem.linearSystem3d(normal.x, normal.y, normal.z, vectorAB.x, vectorAB.y, vectorAB.z, vectorAC.x, vectorAC.y, vectorAC.z, 0, // vectorToCenter DOT normal = 0 (ensure normal is perp to the plane of the 3 points)
             0.5 * ab2, // vectorToCenter DOT vectorAB = ab2 / 2 (ensure the projection of vectorToCenter on AB bisects AB)
             0.5 * ac2);
             if (vectorToCenter) { // i.e., the negative of vectorX
@@ -449,18 +472,18 @@ class Arc3d extends CurvePrimitive_1.CurvePrimitive {
         return undefined;
     }
     /**
-     * Convert a fractional position to xyz coordinates
-     * @param fraction fractional position on arc
-     * @param result optional preallocated result
+     * Convert a fractional position to xyz coordinates.
+     * @param fraction fractional position on arc.
+     * @param result optional preallocated result.
      */
     fractionToPoint(fraction, result) {
         const radians = this._sweep.fractionToRadians(fraction);
         return this._matrix.originPlusMatrixTimesXY(this._center, Math.cos(radians), Math.sin(radians), result);
     }
     /**
-     * Convert fractional arc and radial positions to xyz coordinates
-     * @param fraction fractional position on arc
-     * @param result optional preallocated result
+     * Convert fractional arc and radial positions to xyz coordinates.
+     * @param fraction fractional position on arc.
+     * @param result optional preallocated result.
      */
     fractionAndRadialFractionToPoint(arcFraction, radialFraction, result) {
         const radians = this._sweep.fractionToRadians(arcFraction);
@@ -468,8 +491,8 @@ class Arc3d extends CurvePrimitive_1.CurvePrimitive {
     }
     /**
      * Convert a fractional position to xyz coordinates and derivative with respect to fraction.
-     * @param fraction fractional position on arc
-     * @param result optional preallocated result
+     * @param fraction fractional position on arc.
+     * @param result optional preallocated result.
      */
     fractionToPointAndDerivative(fraction, result) {
         result = this.radiansToPointAndDerivative(this._sweep.fractionToRadians(fraction), result);
@@ -478,10 +501,10 @@ class Arc3d extends CurvePrimitive_1.CurvePrimitive {
     }
     /**
      * Construct a plane with
-     * * origin at the fractional position along the arc
-     * * x axis is the first derivative, i.e. tangent along the arc
+     * * origin at the fractional position along the arc.
+     * * x axis is the first derivative, i.e. tangent along the arc.
      * * y axis is the second derivative, i.e. in the plane and on the center side of the tangent.
-     * If the arc is circular, the second derivative is directly towards the center
+     * If the arc is circular, the second derivative is directly towards the center.
      */
     fractionToPointAnd2Derivatives(fraction, result) {
         const radians = this._sweep.fractionToRadians(fraction);
@@ -497,8 +520,8 @@ class Arc3d extends CurvePrimitive_1.CurvePrimitive {
         return result;
     }
     /**
-     * Evaluate the point and derivative with respect to the angle (in radians)
-     * @param radians angular position
+     * Evaluate the point and derivative with respect to the angle (in radians).
+     * @param radians angular position.
      * @param result optional preallocated ray.
      */
     radiansToPointAndDerivative(radians, result) {
@@ -510,8 +533,8 @@ class Arc3d extends CurvePrimitive_1.CurvePrimitive {
         return result;
     }
     /**
-     * Evaluate the point with respect to the angle (in radians)
-     * @param radians angular position
+     * Evaluate the point with respect to the angle (in radians).
+     * @param radians angular position.
      * @param result optional preallocated ray.
      */
     radiansToPoint(radians, result) {
@@ -523,11 +546,11 @@ class Arc3d extends CurvePrimitive_1.CurvePrimitive {
     }
     /**
      * Return a parametric plane with
-     * * origin at arc center
-     * * vectorU from center to arc at angle (in radians)
+     * * origin at arc center.
+     * * vectorU from center to arc at angle (in radians).
      * * vectorV from center to arc at 90 degrees past the angle.
-     * @param radians angular position
-     * @param result optional preallocated plane
+     * @param radians angular position.
+     * @param result optional preallocated plane.
      */
     radiansToRotatedBasis(radians, result) {
         result = result ? result : Plane3dByOriginAndVectors_1.Plane3dByOriginAndVectors.createXYPlane();
@@ -539,8 +562,8 @@ class Arc3d extends CurvePrimitive_1.CurvePrimitive {
         return result;
     }
     /**
-     * Evaluate the point and derivative with respect to the angle (in radians)
-     * @param theta angular position
+     * Evaluate the point and derivative with respect to the angle (in radians).
+     * @param theta angular position.
      * @param result optional preallocated ray.
      */
     angleToPointAndDerivative(theta, result) {
@@ -553,14 +576,14 @@ class Arc3d extends CurvePrimitive_1.CurvePrimitive {
     }
     /**
      * Return the start point of the arc.
-     * @param result optional preallocated result
+     * @param result optional preallocated result.
      */
     startPoint(result) {
         return this.fractionToPoint(0.0, result);
     }
     /**
      * Return the end point of the arc.
-     * @param result optional preallocated result
+     * @param result optional preallocated result.
      */
     endPoint(result) {
         return this.fractionToPoint(1.0, result);
@@ -632,9 +655,9 @@ class Arc3d extends CurvePrimitive_1.CurvePrimitive {
         return chordSum * factor;
     }
     /**
-     * * See extended comments on `CurvePrimitive.moveSignedDistanceFromFraction`
-     * * A zero length line generates `CurveSearchStatus.error`
-     * * Nonzero length line generates `CurveSearchStatus.success` or `CurveSearchStatus.stoppedAtBoundary`
+     * * See extended comments on `CurvePrimitive.moveSignedDistanceFromFraction`.
+     * * A zero length line generates `CurveSearchStatus.error`.
+     * * Nonzero length line generates `CurveSearchStatus.success` or `CurveSearchStatus.stoppedAtBoundary`.
      */
     moveSignedDistanceFromFraction(startFraction, signedDistance, allowExtension, result) {
         if (!this.isCircular) // suppress extension !!!
@@ -648,7 +671,7 @@ class Arc3d extends CurvePrimitive_1.CurvePrimitive {
     }
     /**
      * Return all radian angles where the ellipse tangent is perpendicular to the vector to a spacePoint.
-     * @param spacePoint point of origin of vectors to the ellipse
+     * @param spacePoint point of origin of vectors to the ellipse.
      * @param _extend always true. Sweep is ignored: perpendiculars for the full ellipse are returned.
      * @param endpoints if true, force the end radians into the result.
      */
@@ -717,7 +740,7 @@ class Arc3d extends CurvePrimitive_1.CurvePrimitive {
     /**
      * Apply a transform to the arc basis vectors.
      * * nonuniform (i.e. skewing) transforms are allowed.
-     * * The transformed vector0 and vector90 are NOT squared up as major minor axes.  (This is a good feature!!)
+     * * The transformed vector0 and vector90 are NOT squared up as major minor axes (this is a good feature).
      */
     tryTransformInPlace(transform) {
         this._center = transform.multiplyPoint3d(this._center, this._center);
@@ -726,10 +749,7 @@ class Arc3d extends CurvePrimitive_1.CurvePrimitive {
         this.setVector0Vector90(this._matrix.columnX(), this._matrix.columnY());
         return true;
     }
-    /**
-     * Return true if the ellipse center and basis vectors are in the plane
-     * @param plane
-     */
+    /** Return true if the ellipse center and basis vectors are in the plane. */
     isInPlane(plane) {
         const normal = plane.getNormalRef();
         // The ellipse vectors are full-length  -- true distance comparisons say things.
@@ -757,17 +777,17 @@ class Arc3d extends CurvePrimitive_1.CurvePrimitive {
             return Geometry_1.Geometry.hypotenuseXY(ux, uy);
         return undefined;
     }
-    /** If the arc is circular, return its radius. Otherwise return undefined */
+    /** If the arc is circular, return its radius. Otherwise return undefined. */
     circularRadius() {
         return this.isCircular ? this._matrix.columnXMagnitude() : undefined;
     }
-    /** Return the larger of the two defining vectors. */
+    /** Return the larger length of the two defining vectors. */
     maxVectorLength() {
         return Math.max(this._matrix.columnXMagnitude(), this._matrix.columnYMagnitude());
     }
     /**
      * Compute intersections with a plane.
-     * @param plane plane to intersect
+     * @param plane plane to intersect.
      * @param result array of locations on the curve.
      */
     appendPlaneIntersectionPoints(plane, result) {
@@ -830,8 +850,8 @@ class Arc3d extends CurvePrimitive_1.CurvePrimitive {
         range.extend(highPoint);
     }
     /**
-     * Returns a (high accuracy) range of the curve between fractional positions
-     * * Default implementation returns teh range of the curve from clonePartialCurve
+     * Returns a (high accuracy) range of the curve between fractional positions.
+     * * Default implementation returns teh range of the curve from clonePartialCurve.
      */
     rangeBetweenFractions(fraction0, fraction1, transform) {
         const sweep = AngleSweep_1.AngleSweep.createStartEndRadians(this.sweep.fractionToRadians(fraction0), this.sweep.fractionToRadians(fraction1));
@@ -841,7 +861,7 @@ class Arc3d extends CurvePrimitive_1.CurvePrimitive {
     }
     /**
      * Set up a SineCosinePolynomial as the function c+u*cos(theta)+v*sin(theta) where
-     *  c,u,v are coefficients obtained by evaluating altitude and velocity relative to the plane.
+     * c,u,v are coefficients obtained by evaluating altitude and velocity relative to the plane.
      * @param plane plane for altitude calculation.
      * @param result optional result.
      * @internal
@@ -861,32 +881,32 @@ class Arc3d extends CurvePrimitive_1.CurvePrimitive {
     }
     /**
      * Create a new arc which is parallel to the xy plane, with given center and radius and optional angle sweep.
-     * @param center center of arc
-     * @param radius radius of arc
-     * @param sweep sweep limits.  defaults to full circle.
+     * @param center center of arc.
+     * @param radius radius of arc.
+     * @param sweep sweep limits; defaults to full circle.
      */
     static createXY(center, radius, sweep = AngleSweep_1.AngleSweep.create360()) {
         return new Arc3d(center.clone(), Matrix3d_1.Matrix3d.createScale(radius, radius, 1.0), sweep.clone());
     }
     /**
      * Create a new arc which is parallel to the xy plane, with given center and x,y radii, and optional angle sweep
-     * @param center center of ellipse
-     * @param radiusA x axis radius
-     * @param radiusB y axis radius
-     * @param sweep angle sweep
+     * @param center center of ellipse.
+     * @param radiusA x axis radius.
+     * @param radiusB y axis radius.
+     * @param sweep angle sweep.
      */
     static createXYEllipse(center, radiusA, radiusB, sweep = AngleSweep_1.AngleSweep.create360()) {
         return new Arc3d(center.clone(), Matrix3d_1.Matrix3d.createScale(radiusA, radiusB, 1.0), sweep.clone());
     }
     /**
      * Replace the arc's 0 and 90 degree vectors.
-     * @param vector0 vector from center to ellipse point at 0 degrees in parameter space
-     * @param vector90 vector from center to ellipse point at 90 degrees in parameter space
+     * @param vector0 vector from center to ellipse point at 0 degrees in parameter space.
+     * @param vector90 vector from center to ellipse point at 90 degrees in parameter space.
      */
     setVector0Vector90(vector0, vector90) {
         this._matrix.setColumns(vector0, vector90, vector0.unitCrossProductWithDefault(vector90, 0, 0, 0));
     }
-    /** Return the arc definition with rigid matrix form with axis radii */
+    /** Return the arc definition with rigid matrix form with axis radii. */
     toScaledMatrix3d() {
         const angleData = Angle_1.Angle.dotProductsToHalfAngleTrigValues(this._matrix.columnXMagnitudeSquared(), this._matrix.columnYMagnitudeSquared(), this._matrix.columnXDotColumnY(), true);
         const vector0A = this._matrix.multiplyXY(angleData.c, angleData.s);
@@ -900,7 +920,7 @@ class Arc3d extends CurvePrimitive_1.CurvePrimitive {
             sweep: this.sweep.cloneMinusRadians(angleData.radians),
         };
     }
-    /** Return the arc definition with center, two vectors, and angle sweep; */
+    /** Return the arc definition with center, two vectors, and angle sweep. */
     toVectors() {
         return {
             center: this.center.clone(),
@@ -935,9 +955,9 @@ class Arc3d extends CurvePrimitive_1.CurvePrimitive {
     }
     /**
      * Set this arc from a json object with these values:
-     * * center center point
-     * * vector0 vector from center to 0 degree point in parameter space (commonly but not always the major axis vector)
-     * * vector90 vector from center to 90 degree point in parameter space (commonly but not always the minor axis vector)
+     * * center center point.
+     * * vector0 vector from center to 0 degree point in parameter space (commonly but not always the major axis vector).
+     * * vector90 vector from center to 90 degree point in parameter space (commonly but not always the minor axis vector).
      * @param json
      */
     setFromJSON(json) {
@@ -968,22 +988,22 @@ class Arc3d extends CurvePrimitive_1.CurvePrimitive {
             vector90: this._matrix.columnY().toJSON(),
         };
     }
-    /** Test if this arc is almost equal to another GeometryQuery object */
-    isAlmostEqual(otherGeometry) {
+    /** Test if this arc is almost equal to another GeometryQuery object. */
+    isAlmostEqual(otherGeometry, distanceTol = Geometry_1.Geometry.smallMetricDistance, radianTol = Geometry_1.Geometry.smallAngleRadians) {
         if (otherGeometry instanceof Arc3d) {
             const other = otherGeometry;
-            return this._center.isAlmostEqual(other._center)
-                && this._matrix.isAlmostEqual(other._matrix)
-                && this._sweep.isAlmostEqualAllowPeriodShift(other._sweep);
+            return this._center.isAlmostEqual(other._center, distanceTol)
+                && this._matrix.isAlmostEqual(other._matrix, distanceTol)
+                && this._sweep.isAlmostEqualAllowPeriodShift(other._sweep, radianTol);
         }
         return false;
     }
-    /** Emit strokes to caller-supplied linestring */
+    /** Emit strokes to caller-supplied linestring. */
     emitStrokes(dest, options) {
         const numStrokes = this.computeStrokeCountForOptions(options);
         dest.appendFractionalStrokePoints(this, numStrokes, 0.0, 1.0, true);
     }
-    /** Emit strokes to caller-supplied handler */
+    /** Emit strokes to caller-supplied handler. */
     emitStrokableParts(handler, options) {
         const numStrokes = this.computeStrokeCountForOptions(options);
         handler.startCurvePrimitive(this);
@@ -992,7 +1012,7 @@ class Arc3d extends CurvePrimitive_1.CurvePrimitive {
     }
     /**
      * Return the stroke count required for given options.
-     * @param options StrokeOptions that determine count
+     * @param options StrokeOptions that determine count.
      */
     computeStrokeCountForOptions(options) {
         let numStroke;
@@ -1005,14 +1025,14 @@ class Arc3d extends CurvePrimitive_1.CurvePrimitive {
         }
         return numStroke;
     }
-    /** Second step of double dispatch:  call `handler.handleArc3d(this)` */
+    /** Second step of double dispatch: call `handler.handleArc3d(this)` */
     dispatchToGeometryHandler(handler) {
         return handler.handleArc3d(this);
     }
     /**
      * Return (if possible) an arc which is a portion of this curve.
-     * @param fractionA [in] start fraction
-     * @param fractionB [in] end fraction
+     * @param fractionA start fraction.
+     * @param fractionB end fraction.
      */
     clonePartialCurve(fractionA, fractionB) {
         if (fractionB < fractionA) {
@@ -1026,8 +1046,8 @@ class Arc3d extends CurvePrimitive_1.CurvePrimitive {
     }
     /**
      * Return an arc whose basis vectors are rotated by given angle within the current basis space.
-     * * the result arc will have its zero-degree point (new `vector0`) at the current
-     * `vector0 * cos(theta) + vector90 * sin(theta)`
+     * * the result arc will have its zero-degree point (new `vector0`) at the current.
+     * `vector0 * cos(theta) + vector90 * sin(theta)`.
      * * the result sweep is adjusted so all fractional coordinates (e.g. start and end) evaluate to the same xyz.
      *   * Specifically, theta is subtracted from the original start and end angles.
      * @param theta the angle (in the input arc space) which is to become the 0-degree point in the new arc.
@@ -1042,10 +1062,10 @@ class Arc3d extends CurvePrimitive_1.CurvePrimitive {
         return arcB;
     }
     /**
-     * Find intervals of this CurvePrimitive that are interior to a clipper
-     * @param clipper clip structure (e.g.clip planes)
-     * @param announce(optional) function to be called announcing fractional intervals"
-     * ` announce(fraction0, fraction1, curvePrimitive)`
+     * Find intervals of this CurvePrimitive that are interior to a clipper.
+     * @param clipper clip structure (e.g.clip planes).
+     * @param announce (optional) function to be called announcing fractional intervals
+     * `announce(fraction0, fraction1, curvePrimitive)`.
      * @returns true if any "in" segments are announced.
      */
     announceClipIntervals(clipper, announce) {
@@ -1057,7 +1077,12 @@ class Arc3d extends CurvePrimitive_1.CurvePrimitive {
         const otherVector0 = this._matrix.multiplyInverse(other.vector0);
         const otherVector90 = this._matrix.multiplyInverse(other.vector90);
         if (otherOrigin && otherVector0 && otherVector90) {
-            return { center: otherOrigin, vector0: otherVector0, vector90: otherVector90, sweep: this.sweep.clone() };
+            return {
+                center: otherOrigin,
+                vector0: otherVector0,
+                vector90: otherVector90,
+                sweep: this.sweep.clone(),
+            };
         }
         return undefined;
     }
@@ -1065,20 +1090,20 @@ class Arc3d extends CurvePrimitive_1.CurvePrimitive {
      * Determine an arc "at a point of inflection" of a point sequence.
      * * Return the arc along with the fractional positions of the tangency points.
      * * In the returned object:
-     *   * `arc` is the (bounded) arc
+     *   * `arc` is the (bounded) arc.
      *   * `fraction10` is the tangency point's position as an interpolating fraction of the line segment from
-     * `point1` (backwards) to `point0`
+     * `point1` (backwards) to `point0`.
      *   * `fraction12` is the tangency point's position as an interpolating fraction of the line segment from
-     * `point1` (forward) to `point2`
+     * `point1` (forward) to `point2`.
      *   * `point` is the `point1` input.
      * * If unable to construct the arc:
      *   * `point` is the `point1` input.
-     *   * both fractions are zero
+     *   * both fractions are zero.
      *   * `arc` is undefined.
-     * @param point0 first point of path. (the point before the point of inflection)
-     * @param point1 second point of path (the point of inflection)
-     * @param point2 third point of path (the point after the point of inflection)
-     * @param radius arc radius
+     * @param point0 first point of path (the point before the point of inflection).
+     * @param point1 second point of path (the point of inflection).
+     * @param point2 third point of path (the point after the point of inflection).
+     * @param radius arc radius.
      *
      */
     static createFilletArc(point0, point1, point2, radius) {
@@ -1113,11 +1138,11 @@ class Arc3d extends CurvePrimitive_1.CurvePrimitive {
         }
         return { fraction10: 0.0, fraction12: 0.0, point: point1.clone() };
     }
-    /** Scale the vector0 and vector90 vectors by `scaleFactor` */
+    /** Scale the vector0 and vector90 vectors by `scaleFactor`. */
     scaleAboutCenterInPlace(scaleFactor) {
         this._matrix.scaleColumnsInPlace(scaleFactor, scaleFactor, 1.0);
     }
-    /** Return the (signed!) area between (a fractional portion of) the arc and the chord between those points */
+    /** Return the (signed!) area between (a fractional portion of) the arc and the chord between those points. */
     areaToChordXY(fraction0, fraction1) {
         let detJ = Geometry_1.Geometry.crossProductXYXY(this._matrix.coffs[0], this._matrix.coffs[3], this._matrix.coffs[1], this._matrix.coffs[4]);
         // areas in arc of unit circle with radians limits
@@ -1132,7 +1157,7 @@ class Arc3d extends CurvePrimitive_1.CurvePrimitive {
     }
     /**
      * Construct an offset of the instance curve as viewed in the xy-plane (ignoring z).
-     * @param offsetDistanceOrOptions offset distance (positive to left of the instance curve), or options object
+     * @param offsetDistanceOrOptions offset distance (positive to left of the instance curve), or options object.
      */
     constructOffsetXY(offsetDistanceOrOptions) {
         const options = OffsetOptions_1.OffsetOptions.create(offsetDistanceOrOptions);
@@ -1164,7 +1189,7 @@ class Arc3d extends CurvePrimitive_1.CurvePrimitive {
     /**
      * Project instance geometry (via dispatch) onto the given ray, and return the extreme fractional parameters of projection.
      * @param ray ray onto which the instance is projected. A `Vector3d` is treated as a `Ray3d` with zero origin.
-     * @param lowHigh optional receiver for output
+     * @param lowHigh optional receiver for output.
      * @returns range of fractional projection parameters onto the ray, where 0.0 is start of the ray and 1.0 is the end of the ray.
      */
     projectedParameterRange(ray, lowHigh) {
@@ -1172,7 +1197,7 @@ class Arc3d extends CurvePrimitive_1.CurvePrimitive {
     }
     /**
      * Construct a circular arc chain approximation to the instance elliptical arc.
-     * @param options bundle of options for sampling an elliptical arc (use default options if undefined)
+     * @param options bundle of options for sampling an elliptical arc (use default options if undefined).
      * @returns the approximating curve chain, the circular instance, or undefined if construction fails.
      */
     constructCircularArcChainApproximation(options) {
@@ -1194,8 +1219,8 @@ Arc3d._workVectorV = Point3dVector3d_1.Vector3d.create();
 Arc3d._workVectorW = Point3dVector3d_1.Vector3d.create();
 // !! misspelled Gauss in the published static !! Declare it ok.
 // cspell::word Guass
-/** Gauss point quadrature count for evaluating curve length. (The number of intervals is adjusted to the arc sweep) */
+/** Gauss point quadrature count for evaluating curve length. (The number of intervals is adjusted to the arc sweep). */
 Arc3d.quadratureGuassCount = 5;
-/** In quadrature for arc length, use this interval (divided by quickEccentricity) */
+/** In quadrature for arc length, use this interval (divided by quickEccentricity). */
 Arc3d.quadratureIntervalAngleDegrees = 10.0;
 //# sourceMappingURL=Arc3d.js.map