@itwin/core-geometry 5.0.0-dev.10 → 5.0.0-dev.100

package/lib/cjs/geometry3d/Ellipsoid.js

@@ -44,6 +44,20 @@ const Transform_1 = require("./Transform");
  * @internal
  */
 class EllipsoidComponentExtrema {
+    c;
+    u;
+    v;
+    w;
+    axisIndex;
+    theta0Radians;
+    phi0Radians;
+    cosTheta0;
+    sinTheta0;
+    cosPhi0;
+    sinPhi0;
+    // temp vars used in arc range evaluation:
+    _axisRange;
+    _trigForm;
     constructor(c, u, v, w, axisIndex) {
         this.c = c;
         this.u = u;
@@ -117,19 +131,25 @@ class EllipsoidComponentExtrema {
     }
 }
 /**
- * * An Ellipsoid is a (complete) unit sphere with an arbitrary (possibly skewed) `Transform` to 3d.
+ * A complete unit sphere mapped by an arbitrary [[Transform]].
  * * The (unit) sphere parameterization with respect to longitude `theta` and latitude `phi` is
  *    * `u = cos(theta) * cos (phi)`
  *    * `v = sin(theta) * cos(phi)`
  *    * `w = sin(phi)`
- *  * The sphere (u,v,w) multiply the x,y,z columns of the Ellipsoid transform.
+ * * The sphere (u,v,w) multiply the x,y,z columns of the Ellipsoid transform.
+ * * Compare to [[Sphere]], which has the same parameterization, but is a [[SolidPrimitive]] with latitude sweep.
  * @public
  */
 class Ellipsoid {
+    _transform;
+    _workUnitVectorA;
+    _workUnitVectorB;
+    _workPointA;
+    _workPointB;
     constructor(transform) {
         this._transform = transform;
-        this._unitVectorA = Point3dVector3d_1.Vector3d.create();
-        this._unitVectorB = Point3dVector3d_1.Vector3d.create();
+        this._workUnitVectorA = Point3dVector3d_1.Vector3d.create();
+        this._workUnitVectorB = Point3dVector3d_1.Vector3d.create();
         this._workPointA = Point3dVector3d_1.Point3d.create();
         this._workPointB = Point3dVector3d_1.Point3d.create();
     }
@@ -341,11 +361,11 @@ class Ellipsoid {
      * @param result optional preallocated result
      */
     radiansPairToGreatArc(thetaARadians, phiARadians, thetaBRadians, phiBRadians, result) {
-        Polynomials_1.SphereImplicit.radiansToUnitSphereXYZ(thetaARadians, phiARadians, this._unitVectorA);
-        Polynomials_1.SphereImplicit.radiansToUnitSphereXYZ(thetaBRadians, phiBRadians, this._unitVectorB);
-        const sweepAngle = this._unitVectorA.angleTo(this._unitVectorB);
+        Polynomials_1.SphereImplicit.radiansToUnitSphereXYZ(thetaARadians, phiARadians, this._workUnitVectorA);
+        Polynomials_1.SphereImplicit.radiansToUnitSphereXYZ(thetaBRadians, phiBRadians, this._workUnitVectorB);
+        const sweepAngle = this._workUnitVectorA.angleTo(this._workUnitVectorB);
         // the unit vectors (on unit sphere) are never 0, so this cannot fail.
-        const matrix = Matrix3d_1.Matrix3d.createRigidFromColumns(this._unitVectorA, this._unitVectorB, Geometry_1.AxisOrder.XYZ);
+        const matrix = Matrix3d_1.Matrix3d.createRigidFromColumns(this._workUnitVectorA, this._workUnitVectorB, Geometry_1.AxisOrder.XYZ);
         if (matrix !== undefined) {
             const matrix1 = this._transform.matrix.multiplyMatrixMatrix(matrix);
             return Arc3d_1.Arc3d.create(this._transform.getOrigin(), matrix1.columnX(), matrix1.columnY(), AngleSweep_1.AngleSweep.createStartEndRadians(0.0, sweepAngle.radians), result);
@@ -430,9 +450,9 @@ class Ellipsoid {
      * @param result optional preallocated result
      */
     radiansPairToEquatorialEllipsoid(thetaARadians, phiARadians, thetaBRadians, phiBRadians, result) {
-        Polynomials_1.SphereImplicit.radiansToUnitSphereXYZ(thetaARadians, phiARadians, this._unitVectorA);
-        Polynomials_1.SphereImplicit.radiansToUnitSphereXYZ(thetaBRadians, phiBRadians, this._unitVectorB);
-        const matrix = Matrix3d_1.Matrix3d.createRigidFromColumns(this._unitVectorA, this._unitVectorB, Geometry_1.AxisOrder.XYZ);
+        Polynomials_1.SphereImplicit.radiansToUnitSphereXYZ(thetaARadians, phiARadians, this._workUnitVectorA);
+        Polynomials_1.SphereImplicit.radiansToUnitSphereXYZ(thetaBRadians, phiBRadians, this._workUnitVectorB);
+        const matrix = Matrix3d_1.Matrix3d.createRigidFromColumns(this._workUnitVectorA, this._workUnitVectorB, Geometry_1.AxisOrder.XYZ);
         if (matrix) {
             if (result) {
                 this._transform.multiplyTransformMatrix3d(matrix, result._transform);
@@ -609,14 +629,14 @@ class Ellipsoid {
         result.direction.setFromPoint3d(result.origin);
         return result;
     }
-    /** Implement the `isPointInOnOrOutside` test fom the `interface` */
+    /** Implementation of [[Clipper.isPointOnOrInside]]. */
     isPointOnOrInside(point) {
         const localPoint = this._transform.multiplyInversePoint3d(point, this._workPointA);
         if (localPoint !== undefined)
             return localPoint.magnitude() <= 1.0;
         return false;
     }
-    /** Announce "in" portions of a line segment.  See `Clipper.announceClippedSegmentIntervals` */
+    /** Announce "in" portions of a line segment. Implementation of [[Clipper.announceClippedSegmentIntervals]]. */
     announceClippedSegmentIntervals(f0, f1, pointA, pointB, announce) {
         const localA = this._transform.multiplyInversePoint3d(pointA, this._workPointA);
         const localB = this._transform.multiplyInversePoint3d(pointB, this._workPointB);
@@ -655,7 +675,7 @@ class Ellipsoid {
         }
         return false;
     }
-    /** Announce "in" portions of a line segment.  See `Clipper.announceClippedSegmentIntervals` */
+    /** Announce "in" portions of a line segment. Implementation of [[Clipper.announceClippedArcIntervals]] */
     announceClippedArcIntervals(arc, announce) {
         const arcData = arc.toVectors();
         let numAnnounce = 0;
@@ -708,6 +728,9 @@ exports.Ellipsoid = Ellipsoid;
  * @public
  */
 class EllipsoidPatch {
+    ellipsoid;
+    longitudeSweep;
+    latitudeSweep;
     /**
      * CAPTURE ellipsoid and sweeps as an EllipsoidPatch.
      * @param ellipsoid
@@ -818,6 +841,15 @@ exports.EllipsoidPatch = EllipsoidPatch;
  * Internal class for searching for the closest point (projection of spacePoint) on an ellipsoid.
  */
 class EllipsoidClosestPoint extends Newton_1.NewtonEvaluatorRRtoRRD {
+    _ellipsoid;
+    _spacePoint;
+    _surfacePoint;
+    _d1Theta;
+    _d2Theta;
+    _d1Phi;
+    _d2Phi;
+    _d2ThetaPhi;
+    _delta;
     constructor(ellipsoid) {
         super();
         this._ellipsoid = ellipsoid;
@@ -863,6 +895,17 @@ class EllipsoidClosestPoint extends Newton_1.NewtonEvaluatorRRtoRRD {
  * @public
  */
 class GeodesicPathPoint {
+    /** First angle, in radians */
+    thetaRadians;
+    /** Second angle, in radians */
+    phiRadians;
+    point;
+    dTheta;
+    dPhi;
+    d2Theta;
+    d2Phi;
+    d2ThetaPhi;
+    d1Cross;
     constructor() {
         this.thetaRadians = 0;
         this.phiRadians = 0;
@@ -879,6 +922,9 @@ class GeodesicPathPoint {
         ellipsoid.radiansToPointAnd2Derivatives(this.thetaRadians, this.phiRadians, this.point, this.dTheta, this.dPhi, this.d2Theta, this.d2Phi, this.d2ThetaPhi);
         this.dTheta.crossProduct(this.dPhi, this.d1Cross);
     }
+    static _vectorAB;
+    static _vectorCB;
+    static _vectorCross;
     /** Evaluate the newton function and derivatives:
      *          `(UAB cross UCB) dot d1cross`
      * with as the central data, UAB = vector from pointA to pointB, UCB = vector from pointC to pointB.
@@ -920,6 +966,9 @@ exports.GeodesicPathPoint = GeodesicPathPoint;
  * @public
  */
 class GeodesicPathSolver {
+    _defaultArc;
+    _pathPoints;
+    _tridiagonalSolver;
     constructor(defaultArc) {
         this._pathPoints = [];
         this._defaultArc = defaultArc;