@itwin/core-geometry 5.0.0-dev.11 → 5.0.0-dev.110

package/lib/cjs/curve/Arc3d.js

@@ -57,6 +57,13 @@ var EllipticalArcSampleMethod;
  * @public
  */
 class EllipticalArcApproximationOptions {
+    _sampleMethod;
+    _numSamplesInQuadrant;
+    _maxError;
+    _remapFunction;
+    _forcePath;
+    /** Default error tolerance. */
+    static defaultMaxError = Constant_1.Constant.oneCentimeter;
     constructor(method, numSamplesInQuadrant, maxError, remapFunction, forcePath) {
         this._sampleMethod = method;
         this._numSamplesInQuadrant = numSamplesInQuadrant;
@@ -134,8 +141,6 @@ class EllipticalArcApproximationOptions {
     }
 }
 exports.EllipticalArcApproximationOptions = EllipticalArcApproximationOptions;
-/** Default error tolerance. */
-EllipticalArcApproximationOptions.defaultMaxError = Constant_1.Constant.oneCentimeter;
 /**
  * Circular or elliptic arc.
  * * The angle to point equation is:
@@ -152,14 +157,32 @@ EllipticalArcApproximationOptions.defaultMaxError = Constant_1.Constant.oneCenti
  * @public
  */
 class Arc3d extends CurvePrimitive_1.CurvePrimitive {
+    /** String name for schema properties. */
+    curvePrimitiveType = "arc";
     /** Test if this and other are both instances of Arc3d. */
     isSameGeometryClass(other) {
         return other instanceof Arc3d;
     }
-    /** Read property for (clone of) center. */
+    _center;
+    _matrix; // columns are [vector0, vector90, unitNormal]
+    _sweep; // sweep limits
+    static _workPointA = Point3dVector3d_1.Point3d.create();
+    static _workPointB = Point3dVector3d_1.Point3d.create();
+    static _workPointC = Point3dVector3d_1.Point3d.create();
+    static _workVectorU = Point3dVector3d_1.Vector3d.create();
+    static _workVectorV = Point3dVector3d_1.Vector3d.create();
+    static _workVectorW = Point3dVector3d_1.Vector3d.create();
+    /** Read/write the center. Getter returns clone. */
     get center() {
         return this._center.clone();
     }
+    set center(center) {
+        this._center.setFrom(center);
+    }
+    /** Read property for (reference to) center. */
+    get centerRef() {
+        return this._center;
+    }
     /** Read property for (clone of) vector0. */
     get vector0() {
         return this._matrix.columnX();
@@ -168,7 +191,19 @@ class Arc3d extends CurvePrimitive_1.CurvePrimitive {
     get vector90() {
         return this._matrix.columnY();
     }
-    /** Read property for (clone of) plane normal, with arbitrary length. */
+    /**
+     * Compute an arc binormal vector with arbitrary length.
+     * * The arc parameterization is counter-clockwise with respect to this vector.
+     * * This vector is parallel to [[perpendicularVector]] and possibly opposite.
+     */
+    binormalVector(result) {
+        const plane = this.fractionToPointAnd2Derivatives(0.0);
+        return plane.vectorU.crossProduct(plane.vectorV, result);
+    }
+    /**
+     * Read property for (clone of) plane normal, with arbitrary length.
+     * * Does not take arc sweep direction into account. See also [[binormalVector]].
+     */
     get perpendicularVector() {
         return this._matrix.columnZ();
     }
@@ -180,7 +215,7 @@ class Arc3d extends CurvePrimitive_1.CurvePrimitive {
     get matrixRef() {
         return this._matrix;
     }
-    /** Sweep of the angle. */
+    /** Read/write the sweep. Getter returns reference. */
     get sweep() {
         return this._sweep;
     }
@@ -194,11 +229,9 @@ class Arc3d extends CurvePrimitive_1.CurvePrimitive {
     // constructor copies the pointers
     constructor(center, matrix, sweep) {
         super();
-        /** String name for schema properties. */
-        this.curvePrimitiveType = "arc";
         this._center = center;
         this._matrix = matrix;
-        this._sweep = sweep;
+        this._sweep = sweep.clampToFullCircle(sweep);
     }
     /**
      * Return a clone of the arc, with transform applied.
@@ -315,7 +348,7 @@ class Arc3d extends CurvePrimitive_1.CurvePrimitive {
         const v1DotV90 = vector1.dotProduct(vector90);
         // 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))
+        if (Geometry_1.Geometry.isSmallMetricDistanceSquared(v90Len)) // tighter than smallMetricDistance to allow flatter long elliptical arcs
             return undefined;
         vector90.scaleInPlace(v90Len);
         return Arc3d.create(center, vector0, vector90, sweep, result);
@@ -596,6 +629,12 @@ class Arc3d extends CurvePrimitive_1.CurvePrimitive {
     curveLength() {
         return this.curveLengthBetweenFractions(0, 1);
     }
+    // !! 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). */
+    static quadratureGuassCount = 5;
+    /** In quadrature for arc length, use this interval (divided by quickEccentricity). */
+    static quadratureIntervalAngleDegrees = 10.0;
     /**
      * * If this is a circular arc, return the simple length derived from radius and sweep.
      * * Otherwise (i.e. if this elliptical) fall through CurvePrimitive integrator.
@@ -734,6 +773,42 @@ class Arc3d extends CurvePrimitive_1.CurvePrimitive {
         }
         return result;
     }
+    /** Override of [[CurvePrimitive.emitTangents]] for Arc3d. */
+    emitTangents(spacePoint, announceTangent, options) {
+        const centerToPoint = Point3dVector3d_1.Vector3d.createStartEnd(this.centerRef, spacePoint);
+        let centerToLocalPoint;
+        if (options?.vectorToEye) {
+            const arcToView = Matrix3d_1.Matrix3d.createColumns(this.matrixRef.getColumn(0), this.matrixRef.getColumn(1), options.vectorToEye);
+            centerToLocalPoint = arcToView.multiplyInverse(centerToPoint);
+        }
+        else {
+            centerToLocalPoint = this.matrixRef.multiplyInverse(centerToPoint);
+        }
+        if (centerToLocalPoint === undefined)
+            return;
+        // centerToLocalPoint is a vector in the local coordinate system of the as-viewed arc.
+        // In other words, the local arc is the unit circle.
+        // alpha is the angle from the local x-axis to centerToLocalPoint.
+        // beta is the nonnegative angle from centerToLocalPoint to a tangency radial.
+        // Tangency angles are preserved by local <-> world transformation.
+        if (centerToLocalPoint !== undefined) {
+            const hypotenuseSquared = centerToLocalPoint.magnitudeSquaredXY();
+            if (hypotenuseSquared >= 1.0) { // localPoint lies outside or on the unit circle...
+                // ...and forms a right triangle with unit radial leg to tangent point
+                const distanceToTangency = Math.sqrt(hypotenuseSquared - 1.0);
+                const alpha = Math.atan2(centerToLocalPoint.y, centerToLocalPoint.x);
+                const beta = Math.atan2(distanceToTangency, 1);
+                const angles = Geometry_1.Geometry.isSmallAngleRadians(beta) ? [alpha] : [alpha + beta, alpha - beta];
+                for (const theta of angles) {
+                    const f = CurveExtendMode_1.CurveExtendOptions.resolveRadiansToValidSweepFraction(options?.extend ?? false, theta, this.sweep);
+                    if (f.isValid) {
+                        const tangent = CurveLocationDetail_1.CurveLocationDetail.createCurveFractionPoint(this, f.fraction, this.fractionToPoint(f.fraction));
+                        announceTangent(tangent);
+                    }
+                }
+            }
+        }
+    }
     /** Reverse the sweep  of the arc. */
     reverseInPlace() {
         this._sweep.reverseInPlace();
@@ -1236,6 +1311,8 @@ class Arc3d extends CurvePrimitive_1.CurvePrimitive {
     }
     /**
      * Construct a circular arc chain approximation to the instance elliptical arc.
+     * * Visualization can be found at https://www.itwinjs.org/sandbox/SaeedTorabi/ArcApproximationGeneral and
+     *  https://www.itwinjs.org/sandbox/SaeedTorabi/ArcApproximation
      * @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.
      */
@@ -1250,16 +1327,4 @@ class Arc3d extends CurvePrimitive_1.CurvePrimitive {
     }
 }
 exports.Arc3d = Arc3d;
-Arc3d._workPointA = Point3dVector3d_1.Point3d.create();
-Arc3d._workPointB = Point3dVector3d_1.Point3d.create();
-Arc3d._workPointC = Point3dVector3d_1.Point3d.create();
-Arc3d._workVectorU = Point3dVector3d_1.Vector3d.create();
-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). */
-Arc3d.quadratureGuassCount = 5;
-/** In quadrature for arc length, use this interval (divided by quickEccentricity). */
-Arc3d.quadratureIntervalAngleDegrees = 10.0;
 //# sourceMappingURL=Arc3d.js.map