@itwin/core-geometry 5.0.0-dev.99 → 5.0.0

package/lib/cjs/geometry4d/MomentData.js

@@ -28,8 +28,8 @@ const Point4d_1 = require("./Point4d");
  *         * e.g. entry 03 is summed x.
  *    * In this level:
  *        * the `absoluteQuantity` member is undefined.
- *        * the `localToWorldMap` and `radiiOfGyration` are created by have undefined contents.
- *  * Second level: after a call to inertiaProductsToPrincipalAxes, the `localToWorldMap`, `absoluteQuantity` and
+ *        * the `localToWorldMap` and `radiiOfGyration` are created but have undefined contents.
+ *  * Second level: after a call to `inertiaProductsToPrincipalAxes`, the `localToWorldMap`, `absoluteQuantity` and
  *    `radiiOfGyration` are filled in.
  * @public
  */
@@ -49,18 +49,35 @@ class MomentData {
      * * This set up with its inverse already constructed.
      */
     localToWorldMap;
+    /** Radii of gyration (square roots of principal second moments). */
+    radiusOfGyration;
+    /**
+     * Principal quantity (e.g. length, area, or volume). This is undefined in raw moments, and becomes defined by
+     * inertiaProductsToPrincipalAxes.
+     */
+    absoluteQuantity;
     // private variables
     static _vectorA;
     static _vectorB;
     static _vectorC;
     _point0 = Point3dVector3d_1.Point3d.create();
     _point1 = Point3dVector3d_1.Point3d.create();
+    /** Constructor. */
+    constructor() {
+        this.origin = Point3dVector3d_1.Point3d.createZero();
+        this.needOrigin = false;
+        this.sums = Matrix4d_1.Matrix4d.createZero();
+        this.localToWorldMap = Transform_1.Transform.createIdentity();
+        this.radiusOfGyration = Point3dVector3d_1.Vector3d.create();
+        this.absoluteQuantity = 0.1; // so optimizer sees its type
+        this.absoluteQuantity = undefined;
+    }
     /**
      * Return the lower-right (3,3) entry in the sums.
      * * This is the quantity (i.e. length, area, or volume) summed.
      */
     get quantitySum() {
-        return this.sums.atIJ(3, 3);
+        return this.sums.weight();
     }
     /**
      * Return a scale factor to make these sums match the target orientation sign.
@@ -91,23 +108,6 @@ class MomentData {
             this.needOrigin = false;
         }
     }
-    /** Radii of gyration (square roots of principal second moments). */
-    radiusOfGyration;
-    /**
-     * Principal quantity (e.g. length, area, or volume). This is undefined in raw moments, and becomes defined by
-     * inertiaProductsToPrincipalAxes.
-     */
-    absoluteQuantity;
-    /** Constructor. */
-    constructor() {
-        this.origin = Point3dVector3d_1.Point3d.createZero();
-        this.sums = Matrix4d_1.Matrix4d.createZero();
-        this.localToWorldMap = Transform_1.Transform.createIdentity();
-        this.radiusOfGyration = Point3dVector3d_1.Vector3d.create();
-        this.needOrigin = false;
-        this.absoluteQuantity = 0.1; // so optimizer sees its type
-        this.absoluteQuantity = undefined;
-    }
     /**
      * Create moments with optional origin.
      * * Origin and needOrigin are quirky.
@@ -151,7 +151,7 @@ class MomentData {
         axes.setColumnsPoint4dXYZ(points[0], points[1], points[2]);
         if (axes.determinant() < 0)
             axes.scaleColumnsInPlace(-1.0, -1.0, -1.0);
-        // prefer x and z positive -- y falls wherever . ..
+        // prefer x and z positive; y falls wherever
         if (axes.at(0, 0) < 0.0)
             axes.scaleColumnsInPlace(-1.0, -1.0, 1.0);
         if (axes.at(2, 2) < 0.0)
@@ -176,7 +176,8 @@ class MomentData {
      * * Hence x axis is long direction.
      * * Hence planar data generates large moment as Z.
      * @param origin The origin used for the inertia products.
-     * @param inertiaProducts The inertia products: sums or integrals of [xx,xy,xz,xw; yx,yy, yz,yw; zx,zy,zz,zw; wx,wy,wz,w].
+     * @param inertiaProducts The inertia products: sums or integrals of
+     * [xx,xy,xz,xw; yx,yy,yz,yw; zx,zy,zz,zw; wx,wy,wz,w].
      */
     static inertiaProductsToPrincipalAxes(origin, inertiaProducts) {
         const moments = new MomentData();
@@ -220,23 +221,21 @@ class MomentData {
      */
     static areEquivalentPrincipalAxes(dataA, dataB) {
         if (dataA && dataB
-            && Geometry_1.Geometry.isSameCoordinate(dataA.quantitySum, dataB.quantitySum)) { // um.. need different tolerance for area, volume?)
+            && Geometry_1.Geometry.isSameCoordinate(dataA.quantitySum, dataB.quantitySum)) { // TODO: need different tolerance for area, volume?
             if (dataA.localToWorldMap.getOrigin().isAlmostEqual(dataB.localToWorldMap.getOrigin())
                 && dataA.radiusOfGyration.isAlmostEqual(dataB.radiusOfGyration)) {
                 if (Geometry_1.Geometry.isSameCoordinate(dataA.radiusOfGyration.x, dataA.radiusOfGyration.y)) {
-                    // We have at least xy symmetry ....
+                    // we have at least xy symmetry
                     if (Geometry_1.Geometry.isSameCoordinate(dataA.radiusOfGyration.x, dataA.radiusOfGyration.z))
                         return true;
-                    // just xy.
-                    // allow opposite z directions.
-                    // If the z's are aligned, x an dy can spin freely.
+                    // just xy; allow opposite z directions; if the z's are aligned, x and y can spin freely
                     const zA = dataA.localToWorldMap.matrix.columnZ();
                     const zB = dataB.localToWorldMap.matrix.columnZ();
                     if (zA.isParallelTo(zB, true))
                         return true;
                     return false;
                 }
-                // no symmetry.  Test all three axes.
+                // no symmetry; test all three axes
                 const vectorA = Point3dVector3d_1.Vector3d.create();
                 const vectorB = Point3dVector3d_1.Vector3d.create();
                 for (let i = 0; i < 3; i++) {
@@ -266,7 +265,7 @@ class MomentData {
     }
     /** Revise the accumulated sums to be "around the centroid". */
     shiftOriginAndSumsToCentroidOfSums() {
-        const xyz = this.sums.columnW().realPoint();
+        const xyz = this.sums.columnW().realPoint(); // centroid of the geometry
         if (xyz) {
             this.shiftOriginAndSumsByXYZ(xyz.x, xyz.y, xyz.z);
             return true;
@@ -275,9 +274,9 @@ class MomentData {
     }
     /**
      * Revise the accumulated sums.
-     * * add ax,ay,ax to the origin coordinates.
-     * * apply the negative translation to the sums.
-    */
+     * * Add (ax,ay,az) to the origin coordinates.
+     * * Apply the negative translation to the sums.
+     */
     shiftOriginAndSumsByXYZ(ax, ay, az) {
         this.origin.addXYZInPlace(ax, ay, az);
         this.sums.multiplyTranslationSandwichInPlace(-ax, -ay, -az);
@@ -287,23 +286,24 @@ class MomentData {
         this.shiftOriginAndSumsByXYZ(newOrigin.x - this.origin.x, newOrigin.y - this.origin.y, newOrigin.z - this.origin.z);
     }
     /**
-     * Compute moments of a triangle from the origin to the given line.
-     * Accumulate them to this.sums.
-     * * If `pointA` is undefined, use `this.origin` as pointA.
-     * * If `this.needOrigin` is set, pointB is used
-    */
+     * Compute moments of a triangle from the origin. Accumulate them to `this.sums`.
+     * * If `this.needOrigin` is set, `this.origin` is set to `pointB`.
+     * * If `pointA` is undefined, use `this.origin` as `pointA`.
+     */
     accumulateTriangleMomentsXY(pointA, pointB, pointC) {
         this.setOriginXYZIfNeeded(pointB.x, pointB.y, 0.0);
         const x0 = this.origin.x;
         const y0 = this.origin.y;
-        const vectorA = MomentData._vectorA =
-            pointA !== undefined ? Point4d_1.Point4d.create(pointA.x - x0, pointA.y - y0, 0.0, 1.0, MomentData._vectorA)
-                : Point4d_1.Point4d.create(this.origin.x, this.origin.y, 0.0, 1.0, MomentData._vectorA);
+        const vectorA = MomentData._vectorA = (pointA !== undefined) ?
+            Point4d_1.Point4d.create(pointA.x - x0, pointA.y - y0, 0.0, 1.0, MomentData._vectorA) :
+            Point4d_1.Point4d.create(0.0, 0.0, 0.0, 1.0, MomentData._vectorA);
         const vectorB = MomentData._vectorB = Point4d_1.Point4d.create(pointB.x - x0, pointB.y - y0, 0.0, 1.0, MomentData._vectorB);
         const vectorC = MomentData._vectorC = Point4d_1.Point4d.create(pointC.x - x0, pointC.y - y0, 0.0, 1.0, MomentData._vectorC);
-        // accumulate Return product integrals I(0<=u<=1) I (0<=v<= u)  (w*W + u *U + v * V)(w*W + u *U + v * V)^  du dv
-        //  where w = 1-u-v
-        //  W = column vector (point00.x, point00.y, point00.z, 1.0) etc.
+        // Below we calculate 16 double integrals: \iint_T [x y 0 1]^ [x y 0 1] dT over triangle T=(A,B,C).
+        // Each accumulates contributions from 9 scaled outer products. Integration computations use the barycentric
+        // change of variables [B-A C-A][u,v]^ = [x,y]^ with Jacobian detJ = B-A x C-A = twice the area of T.
+        // This converts the integration domain from T to the triangle bounded by u=0, v=0 and v=1-u, yielding e.g.,
+        // \iint_T x^2 dT = detJ \int_0^1 \int_0^{1-u} u^2 dv du = detJ / 12, and similarly \iint_T xy dT = detJ / 24.
         const detJ = Geometry_1.Geometry.crossProductXYXY(vectorB.x - vectorA.x, vectorB.y - vectorA.y, vectorC.x - vectorA.x, vectorC.y - vectorA.y);
         if (detJ !== 0.0) {
             const r1_12 = detJ / 12.0;
@@ -319,7 +319,7 @@ class MomentData {
             this.sums.addScaledOuterProductInPlace(vectorC, vectorC, r1_12);
         }
     }
-    /** Add scaled outer product of (4d, unit weight) point to this.sums. */
+    /** Add scaled outer product of (4d, unit weight) point to `this.sums`. */
     accumulateScaledOuterProduct(point, scaleFactor) {
         this.setOriginXYZIfNeeded(point.x, point.y, 0.0);
         const vectorA = MomentData._vectorA = Point4d_1.Point4d.create(point.x - this.origin.x, point.y - this.origin.y, point.z - this.origin.z, 1.0, MomentData._vectorA);
@@ -342,16 +342,15 @@ class MomentData {
         this.sums.addScaledOuterProductInPlace(vectorB, vectorB, r1_3);
     }
     /**
-     * Compute moments of triangles from a base point to the given linestring.
-     * Accumulate them to this.sums.
-     * * If `pointA` is undefined, use `this.origin` as pointA.
-     * * If `this.needOrigin` is set, the first point of the array is captured as local origin for subsequent sums.
-     *
+     * Compute moments of triangles from a base point to the given linestring. Accumulate them to `this.sums`.
+     * * If `this.needOrigin` is set, `this.origin` is set to the first point of the array.
+     * * If `sweepBase` is undefined, use `this.origin` as `sweepBase`.
      */
     accumulateTriangleToLineStringMomentsXY(sweepBase, points) {
         const n = points.length;
         if (n > 1) {
             points.getPoint3dAtUncheckedPointIndex(0, this._point0);
+            // The linestring forms a polygon with sweepBase. Integrate over this polygon using Shoelace algorithm.
             for (let i = 1; i < n; i++) {
                 points.getPoint3dAtUncheckedPointIndex(i, this._point1);
                 this.accumulateTriangleMomentsXY(sweepBase, this._point0, this._point1);
@@ -360,17 +359,17 @@ class MomentData {
         }
     }
     /**
-     * * Assemble XX, YY, XY products into a full matrix form [xx,xy,0,0; xy,yy,0,0;0,0,0,0;0,0,0,1].
-     * * Sandwich this between transforms with columns [vectorU, vectorV, 0000, origin].  (Column weights 0001) (only xy
-     *   parts of vectors).
-     * * scale by detJ for the xy-only determinant of the vectors.
+     * Assemble XX, YY, XY products into a full matrix form [xx,xy,0,0; xy,yy,0,0; 0,0,0,0; 0,0,0,1].
+     * * Sandwich this between transforms with columns [vectorU, vectorV, 0000, origin].
+     *   (column weights 0001; only xy parts of vectors).
+     * * Scale by detJ for the xy-only determinant of the vectors.
      * @param productXX
      * @param productXY
      * @param productYY
-     * @param area Area in caller's system.
-     * @param origin Caller's origin.
-     * @param vectorU Caller's U axis (not necessarily unit).
-     * @param vectorV Caller's V axis (not necessarily unit).
+     * @param area area in caller's system.
+     * @param origin caller's origin.
+     * @param vectorU caller's U axis (not necessarily unit).
+     * @param vectorV caller's V axis (not necessarily unit).
      */
     accumulateXYProductsInCentroidalFrame(productXX, productXY, productYY, area, origin, vectorU, vectorV) {
         const centroidalProducts = Matrix4d_1.Matrix4d.createRowValues(productXX, productXY, 0, 0, productXY, productYY, 0, 0, 0, 0, 0, 0, 0, 0, 0, area);

package/lib/cjs/geometry4d/MomentData.js.map

@@ -1 +1 @@
-{"version":3,"file":"MomentData.js","sourceRoot":"","sources":["../../../src/geometry4d/MomentData.ts"],"names":[],"mappings":";AAAA;;;+FAG+F;;;AAE/F;;GAEG;AAEH,yDAAyD;AACzD,mBAAmB;AAEnB,0CAAuC;AAEvC,qDAAkD;AAClD,mEAAuE;AACvE,uDAAoD;AAEpD,yCAAsC;AACtC,uCAAoC;AAEpC;;;;;;;;;;;;;;;;GAgBG;AACH,MAAa,UAAU;IACrB,4BAA4B;IACrB,MAAM,CAAU;IACvB,6CAA6C;IACtC,UAAU,CAAU;IAC3B;;;;OAIG;IACI,IAAI,CAAW;IACtB;;;OAGG;IACI,eAAe,CAAY;IAClC,oBAAoB;IACZ,MAAM,CAAC,QAAQ,CAAW;IAC1B,MAAM,CAAC,QAAQ,CAAW;IAC1B,MAAM,CAAC,QAAQ,CAAW;IAC1B,OAAO,GAAG,yBAAO,CAAC,MAAM,EAAE,CAAC;IAC3B,OAAO,GAAG,yBAAO,CAAC,MAAM,EAAE,CAAC;IACnC;;;OAGG;IACH,IAAW,WAAW;QACpB,OAAO,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;IAC9B,CAAC;IACD;;;;OAIG;IACI,UAAU,CAAC,UAAkB;QAClC,OAAO,UAAU,GAAG,IAAI,CAAC,WAAW,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC;IACxD,CAAC;IACD,2FAA2F;IACpF,iBAAiB,CAAC,MAAe;QACtC,IAAI,IAAI,CAAC,UAAU,EAAE,CAAC;YACpB,IAAI,CAAC,MAAM,CAAC,cAAc,CAAC,MAAM,CAAC,CAAC;YACnC,IAAI,CAAC,UAAU,GAAG,KAAK,CAAC;QAC1B,CAAC;IACH,CAAC;IACD,2FAA2F;IACpF,qCAAqC,CAAC,MAAwB;QACnE,IAAI,IAAI,CAAC,UAAU,IAAI,MAAM,CAAC,MAAM,GAAG,CAAC,EAAE,CAAC;YACzC,MAAM,CAAC,6BAA6B,CAAC,CAAC,EAAE,IAAI,CAAC,MAAM,CAAC,CAAC;YACrD,IAAI,CAAC,UAAU,GAAG,KAAK,CAAC;QAC1B,CAAC;IACH,CAAC;IACD,2FAA2F;IACpF,oBAAoB,CAAC,CAAS,EAAE,CAAS,EAAE,CAAS;QACzD,IAAI,IAAI,CAAC,UAAU,EAAE,CAAC;YACpB,IAAI,CAAC,MAAM,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;YACzB,IAAI,CAAC,UAAU,GAAG,KAAK,CAAC;QAC1B,CAAC;IACH,CAAC;IACD,oEAAoE;IAC7D,gBAAgB,CAAW;IAClC;;;OAGG;IACI,gBAAgB,CAAU;IACjC,mBAAmB;IACnB;QACE,IAAI,CAAC,MAAM,GAAG,yBAAO,CAAC,UAAU,EAAE,CAAC;QACnC,IAAI,CAAC,IAAI,GAAG,mBAAQ,CAAC,UAAU,EAAE,CAAC;QAClC,IAAI,CAAC,eAAe,GAAG,qBAAS,CAAC,cAAc,EAAE,CAAC;QAClD,IAAI,CAAC,gBAAgB,GAAG,0BAAQ,CAAC,MAAM,EAAE,CAAC;QAC1C,IAAI,CAAC,UAAU,GAAG,KAAK,CAAC;QACxB,IAAI,CAAC,gBAAgB,GAAG,GAAG,CAAC,CAAC,6BAA6B;QAC1D,IAAI,CAAC,gBAAgB,GAAG,SAAS,CAAC;IACpC,CAAC;IACD;;;;;;;MAOE;IACK,MAAM,CAAC,MAAM,CAAC,MAA4B,EAAE,aAAsB,KAAK;QAC5E,MAAM,IAAI,GAAG,IAAI,UAAU,EAAE,CAAC;QAC9B,IAAI,CAAC,UAAU,GAAG,UAAU,CAAC;QAC7B,IAAI,MAAM,EAAE,CAAC;YACX,IAAI,CAAC,MAAM,CAAC,cAAc,CAAC,MAAM,CAAC,CAAC;YACnC,IAAI,CAAC,UAAU,GAAG,KAAK,CAAC;QAC1B,CAAC;QACD,OAAO,IAAI,CAAC;IACd,CAAC;IACD;;;OAGG;IACI,MAAM,CAAC,+BAA+B,CAAC,QAAkB;QAC9D,MAAM,EAAE,GAAG,QAAQ,CAAC,WAAW,EAAE,CAAC;QAClC,MAAM,MAAM,GAAG,mBAAQ,CAAC,WAAW,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC;QAChD,MAAM,CAAC,gBAAgB,CAAC,QAAQ,EAAE,CAAC,GAAG,CAAC,CAAC;QACxC,OAAO,MAAM,CAAC;IAChB,CAAC;IACD,6DAA6D;IACtD,MAAM,CAAC,+BAA+B,CAAC,IAAc,EAAE,OAAiB;QAC7E,MAAM,MAAM,GAAG;YACb,IAAI,CAAC,uBAAuB,CAAC,CAAC,EAAE,OAAO,CAAC,CAAC,CAAC;YAC1C,IAAI,CAAC,uBAAuB,CAAC,CAAC,EAAE,OAAO,CAAC,CAAC,CAAC;YAC1C,IAAI,CAAC,uBAAuB,CAAC,CAAC,EAAE,OAAO,CAAC,CAAC,CAAC;SAC3C,CAAC,IAAI,CAAC,CAAC,KAAc,EAAE,KAAc,EAAU,EAAE;YAChD,IAAI,KAAK,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC;gBAAE,OAAO,CAAC,CAAC,CAAC;YACjC,IAAI,KAAK,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC;gBAAE,OAAO,CAAC,CAAC;YAChC,OAAO,CAAC,CAAC;QACX,CAAC,CAAC,CAAC;QACH,IAAI,CAAC,oBAAoB,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC;QAC3D,IAAI,IAAI,CAAC,WAAW,EAAE,GAAG,CAAC;YACxB,IAAI,CAAC,mBAAmB,CAAC,CAAC,GAAG,EAAE,CAAC,GAAG,EAAE,CAAC,GAAG,CAAC,CAAC;QAC7C,mDAAmD;QACnD,IAAI,IAAI,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,GAAG;YACrB,IAAI,CAAC,mBAAmB,CAAC,CAAC,GAAG,EAAE,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC;QAC5C,IAAI,IAAI,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,GAAG;YACrB,IAAI,CAAC,mBAAmB,CAAC,GAAG,EAAE,CAAC,GAAG,EAAE,CAAC,GAAG,CAAC,CAAC;QAC5C,OAAO,CAAC,GAAG,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IACrD,CAAC;IACD;;;OAGG;IACI,MAAM,CAAC,qBAAqB,CAAC,MAAiB;QACnD,MAAM,OAAO,GAAG,IAAI,UAAU,EAAE,CAAC;QACjC,IAAI,MAAM,CAAC,MAAM,KAAK,CAAC;YACrB,OAAO,OAAO,CAAC;QACjB,OAAO,CAAC,SAAS,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC;QAC7B,OAAO,CAAC,gCAAgC,CAAC,MAAM,CAAC,CAAC;QACjD,OAAO,IAAI,CAAC,8BAA8B,CAAC,OAAO,CAAC,MAAM,EAAE,OAAO,CAAC,IAAI,CAAC,CAAC;IAC3E,CAAC;IACD;;;;;;;OAOG;IACI,MAAM,CAAC,8BAA8B,CAAC,MAAW,EAAE,eAAyB;QACjF,MAAM,OAAO,GAAG,IAAI,UAAU,EAAE,CAAC;QACjC,OAAO,CAAC,IAAI,CAAC,OAAO,CAAC,eAAe,CAAC,CAAC;QACtC,OAAO,CAAC,MAAM,CAAC,OAAO,CAAC,MAAM,CAAC,CAAC;QAC/B,IAAI,CAAC,OAAO,CAAC,kCAAkC,EAAE;YAC/C,OAAO,SAAS,CAAC;QACnB,MAAM,QAAQ,GAAG,OAAO,CAAC,IAAI,CAAC,UAAU,EAAE,CAAC;QAC3C,MAAM,CAAC,GAAG,OAAO,CAAC,IAAI,CAAC,MAAM,EAAE,CAAC;QAChC,IAAI,CAAC,GAAG,GAAG;YACT,QAAQ,CAAC,mBAAmB,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;QAC3C,MAAM,MAAM,GAAG,UAAU,CAAC,+BAA+B,CAAC,QAAQ,CAAC,CAAC;QACpE,MAAM,OAAO,GAAG,0BAAQ,CAAC,MAAM,EAAE,CAAC;QAClC,MAAM,WAAW,GAAG,mBAAQ,CAAC,UAAU,EAAE,CAAC;QAC1C,MAAM,CAAC,wBAAwB,CAAC,WAAW,EAAE,OAAO,CAAC,CAAC;QACtD,IAAI,OAAO,CAAC,CAAC,GAAG,GAAG;YACjB,OAAO,SAAS,CAAC;QACnB,UAAU,CAAC,+BAA+B,CAAC,WAAW,EAAE,OAAO,CAAC,CAAC;QACjE,IAAI,CAAC,GAAG,GAAG;YACT,WAAW,CAAC,mBAAmB,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;QAC7C,OAAO,CAAC,eAAe,GAAG,qBAAS,CAAC,qBAAqB,CAAC,OAAO,CAAC,MAAM,EAAE,WAAW,CAAC,CAAC;QACvF,OAAO,CAAC,gBAAgB,CAAC,GAAG,CAC1B,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,OAAO,CAAC,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,OAAO,CAAC,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QAClG,OAAO,CAAC,gBAAgB,CAAC,YAAY,CAAC,GAAG,GAAG,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QACpE,OAAO,CAAC,gBAAgB,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;QACvC,OAAO,OAAO,CAAC;IACjB,CAAC;IACD;;;;;;;;;;;;;;OAcG;IACI,MAAM,CAAC,0BAA0B,CAAC,KAA6B,EAAE,KAA6B;QACnG,IAAI,KAAK,IAAI,KAAK;eACb,mBAAQ,CAAC,gBAAgB,CAAC,KAAK,CAAC,WAAW,EAAE,KAAK,CAAC,WAAW,CAAC,EAAE,CAAC,CAAE,mDAAmD;YAC1H,IAAI,KAAK,CAAC,eAAe,CAAC,SAAS,EAAE,CAAC,aAAa,CAAC,KAAK,CAAC,eAAe,CAAC,SAAS,EAAE,CAAC;mBACjF,KAAK,CAAC,gBAAgB,CAAC,aAAa,CAAC,KAAK,CAAC,gBAAgB,CAAC,EAAE,CAAC;gBAClE,IAAI,mBAAQ,CAAC,gBAAgB,CAAC,KAAK,CAAC,gBAAgB,CAAC,CAAC,EAAE,KAAK,CAAC,gBAAgB,CAAC,CAAC,CAAC,EAAE,CAAC;oBAClF,oCAAoC;oBACpC,IAAI,mBAAQ,CAAC,gBAAgB,CAAC,KAAK,CAAC,gBAAgB,CAAC,CAAC,EAAE,KAAK,CAAC,gBAAgB,CAAC,CAAC,CAAC;wBAC/E,OAAO,IAAI,CAAC;oBACd,WAAW;oBACX,+BAA+B;oBAC/B,mDAAmD;oBACnD,MAAM,EAAE,GAAG,KAAK,CAAC,eAAe,CAAC,MAAM,CAAC,OAAO,EAAE,CAAC;oBAClD,MAAM,EAAE,GAAG,KAAK,CAAC,eAAe,CAAC,MAAM,CAAC,OAAO,EAAE,CAAC;oBAClD,IAAI,EAAE,CAAC,YAAY,CAAC,EAAE,EAAE,IAAI,CAAC;wBAC3B,OAAO,IAAI,CAAC;oBACd,OAAO,KAAK,CAAC;gBACf,CAAC;gBACD,qCAAqC;gBACrC,MAAM,OAAO,GAAG,0BAAQ,CAAC,MAAM,EAAE,CAAC;gBAClC,MAAM,OAAO,GAAG,0BAAQ,CAAC,MAAM,EAAE,CAAC;gBAClC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;oBAC3B,KAAK,CAAC,eAAe,CAAC,MAAM,CAAC,SAAS,CAAC,CAAC,EAAE,OAAO,CAAC,CAAC;oBACnD,KAAK,CAAC,eAAe,CAAC,MAAM,CAAC,SAAS,CAAC,CAAC,EAAE,OAAO,CAAC,CAAC;oBACnD,IAAI,CAAC,OAAO,CAAC,YAAY,CAAC,OAAO,EAAE,IAAI,CAAC;wBACtC,OAAO,KAAK,CAAC;gBACjB,CAAC;gBACD,OAAO,IAAI,CAAC;YACd,CAAC;QACH,CAAC;QACD,OAAO,KAAK,CAAC;IACf,CAAC;IACD,qEAAqE;IAC9D,SAAS,CAAC,MAAgB;QAC/B,IAAI,CAAC,IAAI,CAAC,OAAO,EAAE,CAAC;QACpB,IAAI,MAAM;YACR,IAAI,CAAC,MAAM,CAAC,OAAO,CAAC,MAAM,CAAC,CAAC;;YAE5B,IAAI,CAAC,MAAM,CAAC,OAAO,EAAE,CAAC;IAC1B,CAAC;IACD,0DAA0D;IACnD,gCAAgC,CAAC,MAAiB;QACvD,KAAK,MAAM,CAAC,IAAI,MAAM,EAAE,CAAC;YACvB,IAAI,CAAC,IAAI,CAAC,iBAAiB,CACzB,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,MAAM,CAAC,CAAC,EACnB,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,MAAM,CAAC,CAAC,EACnB,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,MAAM,CAAC,CAAC,EACnB,GAAG,CACJ,CAAC;QACJ,CAAC;IACH,CAAC;IACD,+DAA+D;IACxD,kCAAkC;QACvC,MAAM,GAAG,GAAG,IAAI,CAAC,IAAI,CAAC,OAAO,EAAE,CAAC,SAAS,EAAE,CAAC;QAC5C,IAAI,GAAG,EAAE,CAAC;YACR,IAAI,CAAC,uBAAuB,CAAC,GAAG,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC;YAClD,OAAO,IAAI,CAAC;QACd,CAAC;QACD,OAAO,KAAK,CAAC;IACf,CAAC;IAED;;;;MAIE;IACK,uBAAuB,CAAC,EAAU,EAAE,EAAU,EAAE,EAAU;QAC/D,IAAI,CAAC,MAAM,CAAC,aAAa,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC;QACtC,IAAI,CAAC,IAAI,CAAC,kCAAkC,CAAC,CAAC,EAAE,EAAE,CAAC,EAAE,EAAE,CAAC,EAAE,CAAC,CAAC;IAC9D,CAAC;IACD,2EAA2E;IACpE,6BAA6B,CAAC,SAAiB;QACpD,IAAI,CAAC,uBAAuB,CAAC,SAAS,CAAC,CAAC,GAAG,IAAI,CAAC,MAAM,CAAC,CAAC,EAAE,SAAS,CAAC,CAAC,GAAG,IAAI,CAAC,MAAM,CAAC,CAAC,EAAE,SAAS,CAAC,CAAC,GAAG,IAAI,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC;IACtH,CAAC;IACD;;;;;MAKE;IACK,2BAA2B,CAAC,MAAyB,EAAE,MAAa,EAAE,MAAa;QACxF,IAAI,CAAC,oBAAoB,CAAC,MAAM,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC;QACnD,MAAM,EAAE,GAAG,IAAI,CAAC,MAAM,CAAC,CAAC,CAAC;QACzB,MAAM,EAAE,GAAG,IAAI,CAAC,MAAM,CAAC,CAAC,CAAC;QACzB,MAAM,OAAO,GAAG,UAAU,CAAC,QAAQ;YACjC,MAAM,KAAK,SAAS,CAAC,CAAC,CAAC,iBAAO,CAAC,MAAM,CAAC,MAAM,CAAC,CAAC,GAAG,EAAE,EAAE,MAAM,CAAC,CAAC,GAAG,EAAE,EAAE,GAAG,EAAE,GAAG,EAAE,UAAU,CAAC,QAAQ,CAAC;gBAChG,CAAC,CAAC,iBAAO,CAAC,MAAM,CAAC,IAAI,CAAC,MAAM,CAAC,CAAC,EAAE,IAAI,CAAC,MAAM,CAAC,CAAC,EAAE,GAAG,EAAE,GAAG,EAAE,UAAU,CAAC,QAAQ,CAAC,CAAC;QAClF,MAAM,OAAO,GAAG,UAAU,CAAC,QAAQ,GAAG,iBAAO,CAAC,MAAM,CAAC,MAAM,CAAC,CAAC,GAAG,EAAE,EAAE,MAAM,CAAC,CAAC,GAAG,EAAE,EAAE,GAAG,EAAE,GAAG,EAAE,UAAU,CAAC,QAAQ,CAAC,CAAC;QAClH,MAAM,OAAO,GAAG,UAAU,CAAC,QAAQ,GAAG,iBAAO,CAAC,MAAM,CAAC,MAAM,CAAC,CAAC,GAAG,EAAE,EAAE,MAAM,CAAC,CAAC,GAAG,EAAE,EAAE,GAAG,EAAE,GAAG,EAAE,UAAU,CAAC,QAAQ,CAAC,CAAC;QAClH,gHAAgH;QAChH,mBAAmB;QACnB,iEAAiE;QACjE,MAAM,IAAI,GAAG,mBAAQ,CAAC,gBAAgB,CAAC,OAAO,CAAC,CAAC,GAAG,OAAO,CAAC,CAAC,EAAE,OAAO,CAAC,CAAC,GAAG,OAAO,CAAC,CAAC,EAAE,OAAO,CAAC,CAAC,GAAG,OAAO,CAAC,CAAC,EAAE,OAAO,CAAC,CAAC,GAAG,OAAO,CAAC,CAAC,CAAC,CAAC;QACnI,IAAI,IAAI,KAAK,GAAG,EAAE,CAAC;YACjB,MAAM,KAAK,GAAG,IAAI,GAAG,IAAI,CAAC;YAC1B,MAAM,KAAK,GAAG,IAAI,GAAG,IAAI,CAAC;YAE1B,IAAI,CAAC,IAAI,CAAC,4BAA4B,CAAC,OAAO,EAAE,OAAO,EAAE,KAAK,CAAC,CAAC;YAChE,IAAI,CAAC,IAAI,CAAC,4BAA4B,CAAC,OAAO,EAAE,OAAO,EAAE,KAAK,CAAC,CAAC;YAChE,IAAI,CAAC,IAAI,CAAC,4BAA4B,CAAC,OAAO,EAAE,OAAO,EAAE,KAAK,CAAC,CAAC;YAEhE,IAAI,CAAC,IAAI,CAAC,4BAA4B,CAAC,OAAO,EAAE,OAAO,EAAE,KAAK,CAAC,CAAC;YAChE,IAAI,CAAC,IAAI,CAAC,4BAA4B,CAAC,OAAO,EAAE,OAAO,EAAE,KAAK,CAAC,CAAC;YAChE,IAAI,CAAC,IAAI,CAAC,4BAA4B,CAAC,OAAO,EAAE,OAAO,EAAE,KAAK,CAAC,CAAC;YAEhE,IAAI,CAAC,IAAI,CAAC,4BAA4B,CAAC,OAAO,EAAE,OAAO,EAAE,KAAK,CAAC,CAAC;YAChE,IAAI,CAAC,IAAI,CAAC,4BAA4B,CAAC,OAAO,EAAE,OAAO,EAAE,KAAK,CAAC,CAAC;YAChE,IAAI,CAAC,IAAI,CAAC,4BAA4B,CAAC,OAAO,EAAE,OAAO,EAAE,KAAK,CAAC,CAAC;QAClE,CAAC;IACH,CAAC;IACD,wEAAwE;IACjE,4BAA4B,CAAC,KAAa,EAAE,WAAmB;QACpE,IAAI,CAAC,oBAAoB,CAAC,KAAK,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC;QACjD,MAAM,OAAO,GAAG,UAAU,CAAC,QAAQ,GAAG,iBAAO,CAAC,MAAM,CAClD,KAAK,CAAC,CAAC,GAAG,IAAI,CAAC,MAAM,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,GAAG,IAAI,CAAC,MAAM,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,GAAG,IAAI,CAAC,MAAM,CAAC,CAAC,EAAE,GAAG,EAAE,UAAU,CAAC,QAAQ,CACpG,CAAC;QACF,IAAI,CAAC,IAAI,CAAC,4BAA4B,CAAC,OAAO,EAAE,OAAO,EAAE,WAAW,CAAC,CAAC;IACxE,CAAC;IACD,4DAA4D;IACrD,wBAAwB,CAAC,MAAe,EAAE,MAAe;QAC9D,IAAI,CAAC,oBAAoB,CAAC,MAAM,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC,CAAC,CAAC;QACxD,MAAM,EAAE,GAAG,IAAI,CAAC,MAAM,CAAC,CAAC,CAAC;QACzB,MAAM,EAAE,GAAG,IAAI,CAAC,MAAM,CAAC,CAAC,CAAC;QACzB,MAAM,EAAE,GAAG,IAAI,CAAC,MAAM,CAAC,CAAC,CAAC;QACzB,MAAM,OAAO,GAAG,UAAU,CAAC,QAAQ,GAAG,iBAAO,CAAC,MAAM,CAClD,MAAM,CAAC,CAAC,GAAG,EAAE,EAAE,MAAM,CAAC,CAAC,GAAG,EAAE,EAAE,MAAM,CAAC,CAAC,GAAG,EAAE,EAAE,GAAG,EAAE,UAAU,CAAC,QAAQ,CACtE,CAAC;QACF,MAAM,OAAO,GAAG,UAAU,CAAC,QAAQ,GAAG,iBAAO,CAAC,MAAM,CAClD,MAAM,CAAC,CAAC,GAAG,EAAE,EAAE,MAAM,CAAC,CAAC,GAAG,EAAE,EAAE,MAAM,CAAC,CAAC,GAAG,EAAE,EAAE,GAAG,EAAE,UAAU,CAAC,QAAQ,CACtE,CAAC;QACF,MAAM,IAAI,GAAG,MAAM,CAAC,QAAQ,CAAC,MAAM,CAAC,CAAC;QACrC,MAAM,IAAI,GAAG,IAAI,GAAG,GAAG,CAAC;QACxB,MAAM,IAAI,GAAG,IAAI,GAAG,GAAG,CAAC;QACxB,IAAI,CAAC,IAAI,CAAC,4BAA4B,CAAC,OAAO,EAAE,OAAO,EAAE,IAAI,CAAC,CAAC;QAC/D,IAAI,CAAC,IAAI,CAAC,4BAA4B,CAAC,OAAO,EAAE,OAAO,EAAE,IAAI,CAAC,CAAC;QAC/D,IAAI,CAAC,IAAI,CAAC,4BAA4B,CAAC,OAAO,EAAE,OAAO,EAAE,IAAI,CAAC,CAAC;QAC/D,IAAI,CAAC,IAAI,CAAC,4BAA4B,CAAC,OAAO,EAAE,OAAO,EAAE,IAAI,CAAC,CAAC;IACjE,CAAC;IACD;;;;;;OAMG;IACI,uCAAuC,CAAC,SAA4B,EAAE,MAAwB;QACnG,MAAM,CAAC,GAAG,MAAM,CAAC,MAAM,CAAC;QACxB,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC;YACV,MAAM,CAAC,+BAA+B,CAAC,CAAC,EAAE,IAAI,CAAC,OAAO,CAAC,CAAC;YACxD,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;gBAC3B,MAAM,CAAC,+BAA+B,CAAC,CAAC,EAAE,IAAI,CAAC,OAAO,CAAC,CAAC;gBACxD,IAAI,CAAC,2BAA2B,CAAC,SAAS,EAAE,IAAI,CAAC,OAAO,EAAE,IAAI,CAAC,OAAO,CAAC,CAAC;gBACxE,IAAI,CAAC,OAAO,CAAC,cAAc,CAAC,IAAI,CAAC,OAAO,CAAC,CAAC;YAC5C,CAAC;QACH,CAAC;IACH,CAAC;IACD;;;;;;;;;;;;OAYG;IACI,qCAAqC,CAC1C,SAAiB,EAAE,SAAiB,EAAE,SAAiB,EACvD,IAAY,EAAE,MAAa,EAAE,OAAc,EAAE,OAAc;QAE3D,MAAM,kBAAkB,GAAG,mBAAQ,CAAC,eAAe,CACjD,SAAS,EAAE,SAAS,EAAE,CAAC,EAAE,CAAC,EAC1B,SAAS,EAAE,SAAS,EAAE,CAAC,EAAE,CAAC,EAC1B,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EACV,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,IAAI,CACd,CAAC;QACF,MAAM,IAAI,GAAG,mBAAQ,CAAC,gBAAgB,CAAC,OAAO,CAAC,CAAC,EAAE,OAAO,CAAC,CAAC,EAAE,OAAO,CAAC,CAAC,EAAE,OAAO,CAAC,CAAC,CAAC,CAAC;QACnF,MAAM,SAAS,GAAG,mBAAQ,CAAC,eAAe,CACxC,OAAO,CAAC,CAAC,EAAE,OAAO,CAAC,CAAC,EAAE,CAAC,EAAE,MAAM,CAAC,CAAC,GAAG,IAAI,CAAC,MAAM,CAAC,CAAC,EACjD,OAAO,CAAC,CAAC,EAAE,OAAO,CAAC,CAAC,EAAE,CAAC,EAAE,MAAM,CAAC,CAAC,GAAG,IAAI,CAAC,MAAM,CAAC,CAAC,EACjD,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EACV,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CACX,CAAC;QACF,MAAM,EAAE,GAAG,SAAS,CAAC,oBAAoB,CAAC,kBAAkB,CAAC,CAAC;QAC9D,MAAM,IAAI,GAAG,EAAE,CAAC,6BAA6B,CAAC,SAAS,CAAC,CAAC;QACzD,IAAI,CAAC,IAAI,CAAC,gBAAgB,CAAC,IAAI,EAAE,IAAI,CAAC,CAAC;IACzC,CAAC;IACD;;;;OAIG;IACI,kBAAkB,CAAC,KAAiB,EAAE,KAAa;QACxD,IAAI,CAAC,iBAAiB,CAAC,KAAK,CAAC,MAAM,CAAC,CAAC;QACrC,IAAI,CAAC,IAAI,CAAC,6BAA6B,CACrC,KAAK,CAAC,IAAI,EAAE,IAAI,CAAC,MAAM,CAAC,CAAC,GAAG,KAAK,CAAC,MAAM,CAAC,CAAC,EAAE,IAAI,CAAC,MAAM,CAAC,CAAC,GAAG,KAAK,CAAC,MAAM,CAAC,CAAC,EAAE,IAAI,CAAC,MAAM,CAAC,CAAC,GAAG,KAAK,CAAC,MAAM,CAAC,CAAC,EAAE,KAAK,CAClH,CAAC;IACJ,CAAC;IAED;;;;OAIG;IACI,4BAA4B,CAAC,MAAe,EAAE,QAAkB,EAAE,KAAa;QACpF,IAAI,CAAC,iBAAiB,CAAC,MAAM,CAAC,CAAC;QAC/B,IAAI,CAAC,IAAI,CAAC,6BAA6B,CACrC,QAAQ,EAAE,IAAI,CAAC,MAAM,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,EAAE,IAAI,CAAC,MAAM,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,EAAE,IAAI,CAAC,MAAM,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,EAAE,KAAK,CAC9F,CAAC;IACJ,CAAC;IACD,qCAAqC;IAC9B,MAAM;QACX,OAAO;YACL,MAAM,EAAE,IAAI,CAAC,MAAM;YACnB,IAAI,EAAE,IAAI,CAAC,IAAI,CAAC,MAAM,EAAE;YACxB,gBAAgB,EAAE,IAAI,CAAC,gBAAgB,CAAC,MAAM,EAAE;YAChD,YAAY,EAAE,IAAI,CAAC,eAAe,CAAC,MAAM,EAAE;SAC5C,CAAC;IACJ,CAAC;CACF;AAtZD,gCAsZC","sourcesContent":["/*---------------------------------------------------------------------------------------------\n* Copyright (c) Bentley Systems, Incorporated. All rights reserved.\n* See LICENSE.md in the project root for license terms and full copyright notice.\n*--------------------------------------------------------------------------------------------*/\n\n/** @packageDocumentation\n * @module Numerics\n */\n\n/* eslint-disable @typescript-eslint/naming-convention */\n// cspell:word ABAT\n\nimport { Geometry } from \"../Geometry\";\nimport { GrowableXYZArray } from \"../geometry3d/GrowableXYZArray\";\nimport { Matrix3d } from \"../geometry3d/Matrix3d\";\nimport { Point3d, Vector3d, XYZ } from \"../geometry3d/Point3dVector3d\";\nimport { Transform } from \"../geometry3d/Transform\";\nimport { XAndY, XYAndZ } from \"../geometry3d/XYZProps\";\nimport { Matrix4d } from \"./Matrix4d\";\nimport { Point4d } from \"./Point4d\";\n\n/**\n * A MomentData structure exists in several levels:\n * * First level: as a carrier of sums of inertial products that determine moments.\n *   * origin = local origin used as moments are summed.\n *   * sums = array of summed moments.\n *     * The [i,j] entry of the sums is a summed or integrated moment for product of axis i and j.\n *       * axes 0,1,2 are x,y,z.\n *         * e.g. entry [0,1] is summed product xy\n *       * axis 3 is \"w\", which is 1 in sums.\n *         * e.g. entry 03 is summed x.\n *    * In this level:\n *        * the `absoluteQuantity` member is undefined.\n *        * the `localToWorldMap` and `radiiOfGyration` are created by have undefined contents.\n *  * Second level: after a call to inertiaProductsToPrincipalAxes, the `localToWorldMap`, `absoluteQuantity` and\n *    `radiiOfGyration` are filled in.\n * @public\n */\nexport class MomentData {\n  /** Origin used for sums. */\n  public origin: Point3d;\n  /** Flag to request deferred origin setup. */\n  public needOrigin: boolean;\n  /**\n   * Moment sums.\n   * * Set to zero at initialization and if requested later.\n   * * Accumulated during data entry phase.\n   */\n  public sums: Matrix4d;\n  /**\n   * The mapping between principal and world system.\n   * * This set up with its inverse already constructed.\n   */\n  public localToWorldMap: Transform;\n  // private variables\n  private static _vectorA?: Point4d;\n  private static _vectorB?: Point4d;\n  private static _vectorC?: Point4d;\n  private _point0 = Point3d.create();\n  private _point1 = Point3d.create();\n  /**\n   * Return the lower-right (3,3) entry in the sums.\n   * * This is the quantity (i.e. length, area, or volume) summed.\n   */\n  public get quantitySum(): number {\n    return this.sums.atIJ(3, 3);\n  }\n  /**\n   * Return a scale factor to make these sums match the target orientation sign.\n   * * 1.0 if `this.quantitySum` has the same sign as `targetSign`.\n   * * -1.0 if `this.quantitySum` has the opposite sign from `targetSign`\n   */\n  public signFactor(targetSign: number): number {\n    return targetSign * this.quantitySum > 0 ? 1.0 : -1.0;\n  }\n  /** If `this.needOrigin` flag is set, copy `origin` to `this.origin` and clear the flag. */\n  public setOriginIfNeeded(origin: Point3d) {\n    if (this.needOrigin) {\n      this.origin.setFromPoint3d(origin);\n      this.needOrigin = false;\n    }\n  }\n  /** If `this.needOrigin` flag is set, copy `origin` to `this.origin` and clear the flag. */\n  public setOriginFromGrowableXYZArrayIfNeeded(points: GrowableXYZArray) {\n    if (this.needOrigin && points.length > 0) {\n      points.getPoint3dAtCheckedPointIndex(0, this.origin);\n      this.needOrigin = false;\n    }\n  }\n  /** If `this.needOrigin` flag is set, copy `origin` to `this.origin` and clear the flag. */\n  public setOriginXYZIfNeeded(x: number, y: number, z: number) {\n    if (this.needOrigin) {\n      this.origin.set(x, y, z);\n      this.needOrigin = false;\n    }\n  }\n  /** Radii of gyration (square roots of principal second moments). */\n  public radiusOfGyration: Vector3d;\n  /**\n   * Principal quantity (e.g. length, area, or volume). This is undefined in raw moments, and becomes defined by\n   * inertiaProductsToPrincipalAxes.\n   */\n  public absoluteQuantity?: number;\n  /** Constructor. */\n  private constructor() {\n    this.origin = Point3d.createZero();\n    this.sums = Matrix4d.createZero();\n    this.localToWorldMap = Transform.createIdentity();\n    this.radiusOfGyration = Vector3d.create();\n    this.needOrigin = false;\n    this.absoluteQuantity = 0.1; // so optimizer sees its type\n    this.absoluteQuantity = undefined;\n  }\n  /**\n   * Create moments with optional origin.\n   * * Origin and needOrigin are quirky.\n   *   * (undefined, true) sets up to use first incoming point as origin.\n   *   * (origin) sets up to use that durable origin, set needsOrigin flag false.\n   *   * (origin, true) the \"true\" is meaningless.\n   *   * (undefined, false) makes 000 the durable origin.\n  */\n  public static create(origin?: Point3d | undefined, needOrigin: boolean = false): MomentData {\n    const data = new MomentData();\n    data.needOrigin = needOrigin;\n    if (origin) {\n      data.origin.setFromPoint3d(origin);\n      data.needOrigin = false;\n    }\n    return data;\n  }\n  /**\n   * Return the formal tensor of integrated values `[yy+zz,xy,xz][yx,xx+zz,yz][zx,xy,xx+yy]`.\n   * @param products matrix of (integrated) `[xx,xy,xz][yx,yy,yz][zx,xy,zz]`.\n   */\n  public static momentTensorFromInertiaProducts(products: Matrix3d): Matrix3d {\n    const rr = products.sumDiagonal();\n    const result = Matrix3d.createScale(rr, rr, rr);\n    result.addScaledInPlace(products, -1.0);\n    return result;\n  }\n  /** Sort the columns of the matrix for increasing moments. */\n  public static sortColumnsForIncreasingMoments(axes: Matrix3d, moments: Vector3d) {\n    const points = [\n      axes.indexedColumnWithWeight(0, moments.x),\n      axes.indexedColumnWithWeight(1, moments.y),\n      axes.indexedColumnWithWeight(2, moments.z),\n    ].sort((dataA: Point4d, dataB: Point4d): number => {\n      if (dataA.w < dataB.w) return -1;\n      if (dataA.w > dataB.w) return 1;\n      return 0;\n    });\n    axes.setColumnsPoint4dXYZ(points[0], points[1], points[2]);\n    if (axes.determinant() < 0)\n      axes.scaleColumnsInPlace(-1.0, -1.0, -1.0);\n    // prefer x and z positive -- y falls wherever . ..\n    if (axes.at(0, 0) < 0.0)\n      axes.scaleColumnsInPlace(-1.0, -1.0, 1.0);\n    if (axes.at(2, 2) < 0.0)\n      axes.scaleColumnsInPlace(1.0, -1.0, -1.0);\n    moments.set(points[0].w, points[1].w, points[2].w);\n  }\n  /**\n   * Return the principal moment data for an array of points.\n   * @param points array of points.\n   */\n  public static pointsToPrincipalAxes(points: Point3d[]): MomentData | undefined {\n    const moments = new MomentData();\n    if (points.length === 0)\n      return moments;\n    moments.clearSums(points[0]);\n    moments.accumulatePointMomentsFromOrigin(points);\n    return this.inertiaProductsToPrincipalAxes(moments.origin, moments.sums);\n  }\n  /**\n   * Compute principal axes from inertial products.\n   * * The radii of gyration are sorted smallest to largest.\n   * * Hence x axis is long direction.\n   * * Hence planar data generates large moment as Z.\n   * @param origin The origin used for the inertia products.\n   * @param inertiaProducts The inertia products: sums or integrals of [xx,xy,xz,xw; yx,yy, yz,yw; zx,zy,zz,zw; wx,wy,wz,w].\n   */\n  public static inertiaProductsToPrincipalAxes(origin: XYZ, inertiaProducts: Matrix4d): MomentData | undefined {\n    const moments = new MomentData();\n    moments.sums.setFrom(inertiaProducts);\n    moments.origin.setFrom(origin);\n    if (!moments.shiftOriginAndSumsToCentroidOfSums())\n      return undefined;\n    const products = moments.sums.matrixPart();\n    const w = moments.sums.weight();\n    if (w < 0.0)\n      products.scaleColumnsInPlace(-1, -1, -1);\n    const tensor = MomentData.momentTensorFromInertiaProducts(products);\n    const moment2 = Vector3d.create();\n    const axisVectors = Matrix3d.createZero();\n    tensor.fastSymmetricEigenvalues(axisVectors, moment2);\n    if (moment2.x < 0.0)\n      return undefined;\n    MomentData.sortColumnsForIncreasingMoments(axisVectors, moment2);\n    if (w < 0.0)\n      axisVectors.scaleColumnsInPlace(1, -1, -1);\n    moments.localToWorldMap = Transform.createOriginAndMatrix(moments.origin, axisVectors);\n    moments.radiusOfGyration.set(\n      Math.sqrt(Math.abs(moment2.x)), Math.sqrt(Math.abs(moment2.y)), Math.sqrt(Math.abs(moment2.z)));\n    moments.radiusOfGyration.scaleInPlace(1.0 / Math.sqrt(Math.abs(w)));\n    moments.absoluteQuantity = Math.abs(w);\n    return moments;\n  }\n  /**\n   * Test for match among selected members as they exist after `inertiaProductsToPrincipalAxes`.\n   * * The members considered are:\n   *   * origin of local to world map (i.e. centroid),\n   *   * radius of gyration,\n   *   * axes of localToWorldMap.\n   * * Axis direction tests allow these quirks:\n   *   * opposite orientation is considered matched.\n   * * Full xyz symmetry: If x,y,z radii are matched, axes are not tested.\n   * * Symmetry in xy plane: If x and y radii are matched, the x and y axes area allowed to spin freely. Only Z direction\n   *   is tested.\n   * * If either or both are undefined, returns false.\n   * @param dataA first set of moments.\n   * @param dataB second set of moments.\n   */\n  public static areEquivalentPrincipalAxes(dataA: MomentData | undefined, dataB: MomentData | undefined): boolean {\n    if (dataA && dataB\n      && Geometry.isSameCoordinate(dataA.quantitySum, dataB.quantitySum)) {  // um.. need different tolerance for area, volume?)\n      if (dataA.localToWorldMap.getOrigin().isAlmostEqual(dataB.localToWorldMap.getOrigin())\n        && dataA.radiusOfGyration.isAlmostEqual(dataB.radiusOfGyration)) {\n        if (Geometry.isSameCoordinate(dataA.radiusOfGyration.x, dataA.radiusOfGyration.y)) {\n          // We have at least xy symmetry ....\n          if (Geometry.isSameCoordinate(dataA.radiusOfGyration.x, dataA.radiusOfGyration.z))\n            return true;\n          // just xy.\n          // allow opposite z directions.\n          // If the z's are aligned, x an dy can spin freely.\n          const zA = dataA.localToWorldMap.matrix.columnZ();\n          const zB = dataB.localToWorldMap.matrix.columnZ();\n          if (zA.isParallelTo(zB, true))\n            return true;\n          return false;\n        }\n        // no symmetry.  Test all three axes.\n        const vectorA = Vector3d.create();\n        const vectorB = Vector3d.create();\n        for (let i = 0; i < 3; i++) {\n          dataA.localToWorldMap.matrix.getColumn(i, vectorA);\n          dataB.localToWorldMap.matrix.getColumn(i, vectorB);\n          if (!vectorA.isParallelTo(vectorB, true))\n            return false;\n        }\n        return true;\n      }\n    }\n    return false;\n  }\n  /** Clear the MomentData sums to zero, and establish a new origin. */\n  public clearSums(origin?: Point3d) {\n    this.sums.setZero();\n    if (origin)\n      this.origin.setFrom(origin);\n    else\n      this.origin.setZero();\n  }\n  /** Accumulate products-of-components for given points. */\n  public accumulatePointMomentsFromOrigin(points: Point3d[]) {\n    for (const p of points) {\n      this.sums.addMomentsInPlace(\n        p.x - this.origin.x,\n        p.y - this.origin.y,\n        p.z - this.origin.z,\n        1.0,\n      );\n    }\n  }\n  /** Revise the accumulated sums to be \"around the centroid\". */\n  public shiftOriginAndSumsToCentroidOfSums(): boolean {\n    const xyz = this.sums.columnW().realPoint();\n    if (xyz) {\n      this.shiftOriginAndSumsByXYZ(xyz.x, xyz.y, xyz.z);\n      return true;\n    }\n    return false;\n  }\n\n  /**\n   * Revise the accumulated sums.\n   * * add ax,ay,ax to the origin coordinates.\n   * * apply the negative translation to the sums.\n  */\n  public shiftOriginAndSumsByXYZ(ax: number, ay: number, az: number) {\n    this.origin.addXYZInPlace(ax, ay, az);\n    this.sums.multiplyTranslationSandwichInPlace(-ax, -ay, -az);\n  }\n  /** Revise the accumulated sums so they are based at a specified origin. */\n  public shiftOriginAndSumsToNewOrigin(newOrigin: XYAndZ) {\n    this.shiftOriginAndSumsByXYZ(newOrigin.x - this.origin.x, newOrigin.y - this.origin.y, newOrigin.z - this.origin.z);\n  }\n  /**\n   * Compute moments of a triangle from the origin to the given line.\n   * Accumulate them to this.sums.\n   * * If `pointA` is undefined, use `this.origin` as pointA.\n   * * If `this.needOrigin` is set, pointB is used\n  */\n  public accumulateTriangleMomentsXY(pointA: XAndY | undefined, pointB: XAndY, pointC: XAndY) {\n    this.setOriginXYZIfNeeded(pointB.x, pointB.y, 0.0);\n    const x0 = this.origin.x;\n    const y0 = this.origin.y;\n    const vectorA = MomentData._vectorA =\n      pointA !== undefined ? Point4d.create(pointA.x - x0, pointA.y - y0, 0.0, 1.0, MomentData._vectorA)\n        : Point4d.create(this.origin.x, this.origin.y, 0.0, 1.0, MomentData._vectorA);\n    const vectorB = MomentData._vectorB = Point4d.create(pointB.x - x0, pointB.y - y0, 0.0, 1.0, MomentData._vectorB);\n    const vectorC = MomentData._vectorC = Point4d.create(pointC.x - x0, pointC.y - y0, 0.0, 1.0, MomentData._vectorC);\n    // accumulate Return product integrals I(0<=u<=1) I (0<=v<= u)  (w*W + u *U + v * V)(w*W + u *U + v * V)^  du dv\n    //  where w = 1-u-v\n    //  W = column vector (point00.x, point00.y, point00.z, 1.0) etc.\n    const detJ = Geometry.crossProductXYXY(vectorB.x - vectorA.x, vectorB.y - vectorA.y, vectorC.x - vectorA.x, vectorC.y - vectorA.y);\n    if (detJ !== 0.0) {\n      const r1_12 = detJ / 12.0;\n      const r1_24 = detJ / 24.0;\n\n      this.sums.addScaledOuterProductInPlace(vectorA, vectorA, r1_12);\n      this.sums.addScaledOuterProductInPlace(vectorA, vectorB, r1_24);\n      this.sums.addScaledOuterProductInPlace(vectorA, vectorC, r1_24);\n\n      this.sums.addScaledOuterProductInPlace(vectorB, vectorA, r1_24);\n      this.sums.addScaledOuterProductInPlace(vectorB, vectorB, r1_12);\n      this.sums.addScaledOuterProductInPlace(vectorB, vectorC, r1_24);\n\n      this.sums.addScaledOuterProductInPlace(vectorC, vectorA, r1_24);\n      this.sums.addScaledOuterProductInPlace(vectorC, vectorB, r1_24);\n      this.sums.addScaledOuterProductInPlace(vectorC, vectorC, r1_12);\n    }\n  }\n  /** Add scaled outer product of (4d, unit weight) point to this.sums. */\n  public accumulateScaledOuterProduct(point: XYAndZ, scaleFactor: number) {\n    this.setOriginXYZIfNeeded(point.x, point.y, 0.0);\n    const vectorA = MomentData._vectorA = Point4d.create(\n      point.x - this.origin.x, point.y - this.origin.y, point.z - this.origin.z, 1.0, MomentData._vectorA,\n    );\n    this.sums.addScaledOuterProductInPlace(vectorA, vectorA, scaleFactor);\n  }\n  /** Accumulate wire moment integral from pointA to pointB */\n  public accumulateLineMomentsXYZ(pointA: Point3d, pointB: Point3d) {\n    this.setOriginXYZIfNeeded(pointA.x, pointA.y, pointA.z);\n    const x0 = this.origin.x;\n    const y0 = this.origin.y;\n    const z0 = this.origin.z;\n    const vectorA = MomentData._vectorA = Point4d.create(\n      pointA.x - x0, pointA.y - y0, pointA.z - z0, 1.0, MomentData._vectorA,\n    );\n    const vectorB = MomentData._vectorB = Point4d.create(\n      pointB.x - x0, pointB.y - y0, pointB.z - z0, 1.0, MomentData._vectorB,\n    );\n    const detJ = pointA.distance(pointB);\n    const r1_3 = detJ / 3.0;\n    const r1_6 = detJ / 6.0;\n    this.sums.addScaledOuterProductInPlace(vectorA, vectorA, r1_3);\n    this.sums.addScaledOuterProductInPlace(vectorA, vectorB, r1_6);\n    this.sums.addScaledOuterProductInPlace(vectorB, vectorA, r1_6);\n    this.sums.addScaledOuterProductInPlace(vectorB, vectorB, r1_3);\n  }\n  /**\n   * Compute moments of triangles from a base point to the given linestring.\n   * Accumulate them to this.sums.\n   * * If `pointA` is undefined, use `this.origin` as pointA.\n   * * If `this.needOrigin` is set, the first point of the array is captured as local origin for subsequent sums.\n   *\n   */\n  public accumulateTriangleToLineStringMomentsXY(sweepBase: XAndY | undefined, points: GrowableXYZArray) {\n    const n = points.length;\n    if (n > 1) {\n      points.getPoint3dAtUncheckedPointIndex(0, this._point0);\n      for (let i = 1; i < n; i++) {\n        points.getPoint3dAtUncheckedPointIndex(i, this._point1);\n        this.accumulateTriangleMomentsXY(sweepBase, this._point0, this._point1);\n        this._point0.setFromPoint3d(this._point1);\n      }\n    }\n  }\n  /**\n   * * Assemble XX, YY, XY products into a full matrix form [xx,xy,0,0; xy,yy,0,0;0,0,0,0;0,0,0,1].\n   * * Sandwich this between transforms with columns [vectorU, vectorV, 0000, origin].  (Column weights 0001) (only xy\n   *   parts of vectors).\n   * * scale by detJ for the xy-only determinant of the vectors.\n   * @param productXX\n   * @param productXY\n   * @param productYY\n   * @param area Area in caller's system.\n   * @param origin Caller's origin.\n   * @param vectorU Caller's U axis (not necessarily unit).\n   * @param vectorV Caller's V axis (not necessarily unit).\n   */\n  public accumulateXYProductsInCentroidalFrame(\n    productXX: number, productXY: number, productYY: number,\n    area: number, origin: XAndY, vectorU: XAndY, vectorV: XAndY,\n  ) {\n    const centroidalProducts = Matrix4d.createRowValues(\n      productXX, productXY, 0, 0,\n      productXY, productYY, 0, 0,\n      0, 0, 0, 0,\n      0, 0, 0, area,\n    );\n    const detJ = Geometry.crossProductXYXY(vectorU.x, vectorV.x, vectorU.y, vectorV.y);\n    const placement = Matrix4d.createRowValues(\n      vectorU.x, vectorV.x, 0, origin.x - this.origin.x,\n      vectorU.y, vectorV.y, 0, origin.y - this.origin.y,\n      0, 0, 0, 0,\n      0, 0, 0, 1,\n    );\n    const AB = placement.multiplyMatrixMatrix(centroidalProducts);\n    const ABAT = AB.multiplyMatrixMatrixTranspose(placement);\n    this.sums.addScaledInPlace(ABAT, detJ);\n  }\n  /**\n   * Accumulate sums from other moments.\n   * * Scale by given scaleFactor (e.g. sign to correct orientation).\n   * * Pull the origin from `other` if `this` needs an origin.\n   */\n  public accumulateProducts(other: MomentData, scale: number) {\n    this.setOriginIfNeeded(other.origin);\n    this.sums.addTranslationSandwichInPlace(\n      other.sums, this.origin.x - other.origin.x, this.origin.y - other.origin.y, this.origin.z - other.origin.z, scale,\n    );\n  }\n\n  /**\n   * Accumulate sums from Matrix4d and origin.\n   * * Scale by given scaleFactor (e.g. sign to correct orientation).\n   * * Trap the origin if `this` needs an origin.\n   */\n  public accumulateProductsFromOrigin(origin: Point3d, products: Matrix4d, scale: number) {\n    this.setOriginIfNeeded(origin);\n    this.sums.addTranslationSandwichInPlace(\n      products, this.origin.x - origin.x, this.origin.y - origin.y, this.origin.z - origin.z, scale,\n    );\n  }\n  /** Convert to a json data object. */\n  public toJSON(): any {\n    return {\n      origin: this.origin,\n      sums: this.sums.toJSON(),\n      radiusOfGyration: this.radiusOfGyration.toJSON(),\n      localToWorld: this.localToWorldMap.toJSON(),\n    };\n  }\n}\n"]}
+{"version":3,"file":"MomentData.js","sourceRoot":"","sources":["../../../src/geometry4d/MomentData.ts"],"names":[],"mappings":";AAAA;;;+FAG+F;;;AAE/F;;GAEG;AAEH,yDAAyD;AACzD,mBAAmB;AAEnB,0CAAuC;AAEvC,qDAAkD;AAClD,mEAAuE;AACvE,uDAAoD;AAEpD,yCAAsC;AACtC,uCAAoC;AAEpC;;;;;;;;;;;;;;;;GAgBG;AACH,MAAa,UAAU;IACrB,4BAA4B;IACrB,MAAM,CAAU;IACvB,6CAA6C;IACtC,UAAU,CAAU;IAC3B;;;;OAIG;IACI,IAAI,CAAW;IACtB;;;OAGG;IACI,eAAe,CAAY;IAClC,oEAAoE;IAC7D,gBAAgB,CAAW;IAClC;;;OAGG;IACI,gBAAgB,CAAU;IACjC,oBAAoB;IACZ,MAAM,CAAC,QAAQ,CAAW;IAC1B,MAAM,CAAC,QAAQ,CAAW;IAC1B,MAAM,CAAC,QAAQ,CAAW;IAC1B,OAAO,GAAG,yBAAO,CAAC,MAAM,EAAE,CAAC;IAC3B,OAAO,GAAG,yBAAO,CAAC,MAAM,EAAE,CAAC;IACnC,mBAAmB;IACnB;QACE,IAAI,CAAC,MAAM,GAAG,yBAAO,CAAC,UAAU,EAAE,CAAC;QACnC,IAAI,CAAC,UAAU,GAAG,KAAK,CAAC;QACxB,IAAI,CAAC,IAAI,GAAG,mBAAQ,CAAC,UAAU,EAAE,CAAC;QAClC,IAAI,CAAC,eAAe,GAAG,qBAAS,CAAC,cAAc,EAAE,CAAC;QAClD,IAAI,CAAC,gBAAgB,GAAG,0BAAQ,CAAC,MAAM,EAAE,CAAC;QAC1C,IAAI,CAAC,gBAAgB,GAAG,GAAG,CAAC,CAAC,6BAA6B;QAC1D,IAAI,CAAC,gBAAgB,GAAG,SAAS,CAAC;IACpC,CAAC;IACD;;;OAGG;IACH,IAAW,WAAW;QACpB,OAAO,IAAI,CAAC,IAAI,CAAC,MAAM,EAAE,CAAC;IAC5B,CAAC;IACD;;;;OAIG;IACI,UAAU,CAAC,UAAkB;QAClC,OAAO,UAAU,GAAG,IAAI,CAAC,WAAW,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC;IACxD,CAAC;IACD,2FAA2F;IACpF,iBAAiB,CAAC,MAAe;QACtC,IAAI,IAAI,CAAC,UAAU,EAAE,CAAC;YACpB,IAAI,CAAC,MAAM,CAAC,cAAc,CAAC,MAAM,CAAC,CAAC;YACnC,IAAI,CAAC,UAAU,GAAG,KAAK,CAAC;QAC1B,CAAC;IACH,CAAC;IACD,2FAA2F;IACpF,qCAAqC,CAAC,MAAwB;QACnE,IAAI,IAAI,CAAC,UAAU,IAAI,MAAM,CAAC,MAAM,GAAG,CAAC,EAAE,CAAC;YACzC,MAAM,CAAC,6BAA6B,CAAC,CAAC,EAAE,IAAI,CAAC,MAAM,CAAC,CAAC;YACrD,IAAI,CAAC,UAAU,GAAG,KAAK,CAAC;QAC1B,CAAC;IACH,CAAC;IACD,2FAA2F;IACpF,oBAAoB,CAAC,CAAS,EAAE,CAAS,EAAE,CAAS;QACzD,IAAI,IAAI,CAAC,UAAU,EAAE,CAAC;YACpB,IAAI,CAAC,MAAM,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;YACzB,IAAI,CAAC,UAAU,GAAG,KAAK,CAAC;QAC1B,CAAC;IACH,CAAC;IACD;;;;;;;MAOE;IACK,MAAM,CAAC,MAAM,CAAC,MAA4B,EAAE,aAAsB,KAAK;QAC5E,MAAM,IAAI,GAAG,IAAI,UAAU,EAAE,CAAC;QAC9B,IAAI,CAAC,UAAU,GAAG,UAAU,CAAC;QAC7B,IAAI,MAAM,EAAE,CAAC;YACX,IAAI,CAAC,MAAM,CAAC,cAAc,CAAC,MAAM,CAAC,CAAC;YACnC,IAAI,CAAC,UAAU,GAAG,KAAK,CAAC;QAC1B,CAAC;QACD,OAAO,IAAI,CAAC;IACd,CAAC;IACD;;;OAGG;IACI,MAAM,CAAC,+BAA+B,CAAC,QAAkB;QAC9D,MAAM,EAAE,GAAG,QAAQ,CAAC,WAAW,EAAE,CAAC;QAClC,MAAM,MAAM,GAAG,mBAAQ,CAAC,WAAW,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC;QAChD,MAAM,CAAC,gBAAgB,CAAC,QAAQ,EAAE,CAAC,GAAG,CAAC,CAAC;QACxC,OAAO,MAAM,CAAC;IAChB,CAAC;IACD,6DAA6D;IACtD,MAAM,CAAC,+BAA+B,CAAC,IAAc,EAAE,OAAiB;QAC7E,MAAM,MAAM,GAAG;YACb,IAAI,CAAC,uBAAuB,CAAC,CAAC,EAAE,OAAO,CAAC,CAAC,CAAC;YAC1C,IAAI,CAAC,uBAAuB,CAAC,CAAC,EAAE,OAAO,CAAC,CAAC,CAAC;YAC1C,IAAI,CAAC,uBAAuB,CAAC,CAAC,EAAE,OAAO,CAAC,CAAC,CAAC;SAC3C,CAAC,IAAI,CAAC,CAAC,KAAc,EAAE,KAAc,EAAU,EAAE;YAChD,IAAI,KAAK,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC;gBAAE,OAAO,CAAC,CAAC,CAAC;YACjC,IAAI,KAAK,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC;gBAAE,OAAO,CAAC,CAAC;YAChC,OAAO,CAAC,CAAC;QACX,CAAC,CAAC,CAAC;QACH,IAAI,CAAC,oBAAoB,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC;QAC3D,IAAI,IAAI,CAAC,WAAW,EAAE,GAAG,CAAC;YACxB,IAAI,CAAC,mBAAmB,CAAC,CAAC,GAAG,EAAE,CAAC,GAAG,EAAE,CAAC,GAAG,CAAC,CAAC;QAC7C,4CAA4C;QAC5C,IAAI,IAAI,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,GAAG;YACrB,IAAI,CAAC,mBAAmB,CAAC,CAAC,GAAG,EAAE,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC;QAC5C,IAAI,IAAI,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,GAAG;YACrB,IAAI,CAAC,mBAAmB,CAAC,GAAG,EAAE,CAAC,GAAG,EAAE,CAAC,GAAG,CAAC,CAAC;QAC5C,OAAO,CAAC,GAAG,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IACrD,CAAC;IACD;;;OAGG;IACI,MAAM,CAAC,qBAAqB,CAAC,MAAiB;QACnD,MAAM,OAAO,GAAG,IAAI,UAAU,EAAE,CAAC;QACjC,IAAI,MAAM,CAAC,MAAM,KAAK,CAAC;YACrB,OAAO,OAAO,CAAC;QACjB,OAAO,CAAC,SAAS,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC;QAC7B,OAAO,CAAC,gCAAgC,CAAC,MAAM,CAAC,CAAC;QACjD,OAAO,IAAI,CAAC,8BAA8B,CAAC,OAAO,CAAC,MAAM,EAAE,OAAO,CAAC,IAAI,CAAC,CAAC;IAC3E,CAAC;IACD;;;;;;;;OAQG;IACI,MAAM,CAAC,8BAA8B,CAAC,MAAW,EAAE,eAAyB;QACjF,MAAM,OAAO,GAAG,IAAI,UAAU,EAAE,CAAC;QACjC,OAAO,CAAC,IAAI,CAAC,OAAO,CAAC,eAAe,CAAC,CAAC;QACtC,OAAO,CAAC,MAAM,CAAC,OAAO,CAAC,MAAM,CAAC,CAAC;QAC/B,IAAI,CAAC,OAAO,CAAC,kCAAkC,EAAE;YAC/C,OAAO,SAAS,CAAC;QACnB,MAAM,QAAQ,GAAG,OAAO,CAAC,IAAI,CAAC,UAAU,EAAE,CAAC;QAC3C,MAAM,CAAC,GAAG,OAAO,CAAC,IAAI,CAAC,MAAM,EAAE,CAAC;QAChC,IAAI,CAAC,GAAG,GAAG;YACT,QAAQ,CAAC,mBAAmB,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;QAC3C,MAAM,MAAM,GAAG,UAAU,CAAC,+BAA+B,CAAC,QAAQ,CAAC,CAAC;QACpE,MAAM,OAAO,GAAG,0BAAQ,CAAC,MAAM,EAAE,CAAC;QAClC,MAAM,WAAW,GAAG,mBAAQ,CAAC,UAAU,EAAE,CAAC;QAC1C,MAAM,CAAC,wBAAwB,CAAC,WAAW,EAAE,OAAO,CAAC,CAAC;QACtD,IAAI,OAAO,CAAC,CAAC,GAAG,GAAG;YACjB,OAAO,SAAS,CAAC;QACnB,UAAU,CAAC,+BAA+B,CAAC,WAAW,EAAE,OAAO,CAAC,CAAC;QACjE,IAAI,CAAC,GAAG,GAAG;YACT,WAAW,CAAC,mBAAmB,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;QAC7C,OAAO,CAAC,eAAe,GAAG,qBAAS,CAAC,qBAAqB,CAAC,OAAO,CAAC,MAAM,EAAE,WAAW,CAAC,CAAC;QACvF,OAAO,CAAC,gBAAgB,CAAC,GAAG,CAC1B,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,OAAO,CAAC,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,OAAO,CAAC,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,OAAO,CAAC,CAAC,CAAC,CAAC,CAC/F,CAAC;QACF,OAAO,CAAC,gBAAgB,CAAC,YAAY,CAAC,GAAG,GAAG,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QACpE,OAAO,CAAC,gBAAgB,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;QACvC,OAAO,OAAO,CAAC;IACjB,CAAC;IACD;;;;;;;;;;;;;;OAcG;IACI,MAAM,CAAC,0BAA0B,CAAC,KAA6B,EAAE,KAA6B;QACnG,IAAI,KAAK,IAAI,KAAK;eACb,mBAAQ,CAAC,gBAAgB,CAAC,KAAK,CAAC,WAAW,EAAE,KAAK,CAAC,WAAW,CAAC,EAAE,CAAC,CAAC,mDAAmD;YACzH,IAAI,KAAK,CAAC,eAAe,CAAC,SAAS,EAAE,CAAC,aAAa,CAAC,KAAK,CAAC,eAAe,CAAC,SAAS,EAAE,CAAC;mBACjF,KAAK,CAAC,gBAAgB,CAAC,aAAa,CAAC,KAAK,CAAC,gBAAgB,CAAC,EAAE,CAAC;gBAClE,IAAI,mBAAQ,CAAC,gBAAgB,CAAC,KAAK,CAAC,gBAAgB,CAAC,CAAC,EAAE,KAAK,CAAC,gBAAgB,CAAC,CAAC,CAAC,EAAE,CAAC;oBAClF,+BAA+B;oBAC/B,IAAI,mBAAQ,CAAC,gBAAgB,CAAC,KAAK,CAAC,gBAAgB,CAAC,CAAC,EAAE,KAAK,CAAC,gBAAgB,CAAC,CAAC,CAAC;wBAC/E,OAAO,IAAI,CAAC;oBACd,wFAAwF;oBACxF,MAAM,EAAE,GAAG,KAAK,CAAC,eAAe,CAAC,MAAM,CAAC,OAAO,EAAE,CAAC;oBAClD,MAAM,EAAE,GAAG,KAAK,CAAC,eAAe,CAAC,MAAM,CAAC,OAAO,EAAE,CAAC;oBAClD,IAAI,EAAE,CAAC,YAAY,CAAC,EAAE,EAAE,IAAI,CAAC;wBAC3B,OAAO,IAAI,CAAC;oBACd,OAAO,KAAK,CAAC;gBACf,CAAC;gBACD,mCAAmC;gBACnC,MAAM,OAAO,GAAG,0BAAQ,CAAC,MAAM,EAAE,CAAC;gBAClC,MAAM,OAAO,GAAG,0BAAQ,CAAC,MAAM,EAAE,CAAC;gBAClC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;oBAC3B,KAAK,CAAC,eAAe,CAAC,MAAM,CAAC,SAAS,CAAC,CAAC,EAAE,OAAO,CAAC,CAAC;oBACnD,KAAK,CAAC,eAAe,CAAC,MAAM,CAAC,SAAS,CAAC,CAAC,EAAE,OAAO,CAAC,CAAC;oBACnD,IAAI,CAAC,OAAO,CAAC,YAAY,CAAC,OAAO,EAAE,IAAI,CAAC;wBACtC,OAAO,KAAK,CAAC;gBACjB,CAAC;gBACD,OAAO,IAAI,CAAC;YACd,CAAC;QACH,CAAC;QACD,OAAO,KAAK,CAAC;IACf,CAAC;IACD,qEAAqE;IAC9D,SAAS,CAAC,MAAgB;QAC/B,IAAI,CAAC,IAAI,CAAC,OAAO,EAAE,CAAC;QACpB,IAAI,MAAM;YACR,IAAI,CAAC,MAAM,CAAC,OAAO,CAAC,MAAM,CAAC,CAAC;;YAE5B,IAAI,CAAC,MAAM,CAAC,OAAO,EAAE,CAAC;IAC1B,CAAC;IACD,0DAA0D;IACnD,gCAAgC,CAAC,MAAiB;QACvD,KAAK,MAAM,CAAC,IAAI,MAAM,EAAE,CAAC;YACvB,IAAI,CAAC,IAAI,CAAC,iBAAiB,CACzB,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,MAAM,CAAC,CAAC,EACnB,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,MAAM,CAAC,CAAC,EACnB,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,MAAM,CAAC,CAAC,EACnB,GAAG,CACJ,CAAC;QACJ,CAAC;IACH,CAAC;IACD,+DAA+D;IACxD,kCAAkC;QACvC,MAAM,GAAG,GAAG,IAAI,CAAC,IAAI,CAAC,OAAO,EAAE,CAAC,SAAS,EAAE,CAAC,CAAC,2BAA2B;QACxE,IAAI,GAAG,EAAE,CAAC;YACR,IAAI,CAAC,uBAAuB,CAAC,GAAG,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC;YAClD,OAAO,IAAI,CAAC;QACd,CAAC;QACD,OAAO,KAAK,CAAC;IACf,CAAC;IACD;;;;OAIG;IACI,uBAAuB,CAAC,EAAU,EAAE,EAAU,EAAE,EAAU;QAC/D,IAAI,CAAC,MAAM,CAAC,aAAa,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC;QACtC,IAAI,CAAC,IAAI,CAAC,kCAAkC,CAAC,CAAC,EAAE,EAAE,CAAC,EAAE,EAAE,CAAC,EAAE,CAAC,CAAC;IAC9D,CAAC;IACD,2EAA2E;IACpE,6BAA6B,CAAC,SAAiB;QACpD,IAAI,CAAC,uBAAuB,CAAC,SAAS,CAAC,CAAC,GAAG,IAAI,CAAC,MAAM,CAAC,CAAC,EAAE,SAAS,CAAC,CAAC,GAAG,IAAI,CAAC,MAAM,CAAC,CAAC,EAAE,SAAS,CAAC,CAAC,GAAG,IAAI,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC;IACtH,CAAC;IACD;;;;OAIG;IACI,2BAA2B,CAAC,MAAyB,EAAE,MAAa,EAAE,MAAa;QACxF,IAAI,CAAC,oBAAoB,CAAC,MAAM,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC;QACnD,MAAM,EAAE,GAAG,IAAI,CAAC,MAAM,CAAC,CAAC,CAAC;QACzB,MAAM,EAAE,GAAG,IAAI,CAAC,MAAM,CAAC,CAAC,CAAC;QACzB,MAAM,OAAO,GAAG,UAAU,CAAC,QAAQ,GAAG,CAAC,MAAM,KAAK,SAAS,CAAC,CAAC,CAAC;YAC5D,iBAAO,CAAC,MAAM,CAAC,MAAM,CAAC,CAAC,GAAG,EAAE,EAAE,MAAM,CAAC,CAAC,GAAG,EAAE,EAAE,GAAG,EAAE,GAAG,EAAE,UAAU,CAAC,QAAQ,CAAC,CAAC,CAAC;YAC7E,iBAAO,CAAC,MAAM,CAAC,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,UAAU,CAAC,QAAQ,CAAC,CAAC;QAC1D,MAAM,OAAO,GAAG,UAAU,CAAC,QAAQ,GAAG,iBAAO,CAAC,MAAM,CAAC,MAAM,CAAC,CAAC,GAAG,EAAE,EAAE,MAAM,CAAC,CAAC,GAAG,EAAE,EAAE,GAAG,EAAE,GAAG,EAAE,UAAU,CAAC,QAAQ,CAAC,CAAC;QAClH,MAAM,OAAO,GAAG,UAAU,CAAC,QAAQ,GAAG,iBAAO,CAAC,MAAM,CAAC,MAAM,CAAC,CAAC,GAAG,EAAE,EAAE,MAAM,CAAC,CAAC,GAAG,EAAE,EAAE,GAAG,EAAE,GAAG,EAAE,UAAU,CAAC,QAAQ,CAAC,CAAC;QAClH,mGAAmG;QACnG,4GAA4G;QAC5G,qGAAqG;QACrG,4GAA4G;QAC5G,8GAA8G;QAC9G,MAAM,IAAI,GAAG,mBAAQ,CAAC,gBAAgB,CACpC,OAAO,CAAC,CAAC,GAAG,OAAO,CAAC,CAAC,EAAE,OAAO,CAAC,CAAC,GAAG,OAAO,CAAC,CAAC,EAAE,OAAO,CAAC,CAAC,GAAG,OAAO,CAAC,CAAC,EAAE,OAAO,CAAC,CAAC,GAAG,OAAO,CAAC,CAAC,CAC3F,CAAC;QACF,IAAI,IAAI,KAAK,GAAG,EAAE,CAAC;YACjB,MAAM,KAAK,GAAG,IAAI,GAAG,IAAI,CAAC;YAC1B,MAAM,KAAK,GAAG,IAAI,GAAG,IAAI,CAAC;YAE1B,IAAI,CAAC,IAAI,CAAC,4BAA4B,CAAC,OAAO,EAAE,OAAO,EAAE,KAAK,CAAC,CAAC;YAChE,IAAI,CAAC,IAAI,CAAC,4BAA4B,CAAC,OAAO,EAAE,OAAO,EAAE,KAAK,CAAC,CAAC;YAChE,IAAI,CAAC,IAAI,CAAC,4BAA4B,CAAC,OAAO,EAAE,OAAO,EAAE,KAAK,CAAC,CAAC;YAEhE,IAAI,CAAC,IAAI,CAAC,4BAA4B,CAAC,OAAO,EAAE,OAAO,EAAE,KAAK,CAAC,CAAC;YAChE,IAAI,CAAC,IAAI,CAAC,4BAA4B,CAAC,OAAO,EAAE,OAAO,EAAE,KAAK,CAAC,CAAC;YAChE,IAAI,CAAC,IAAI,CAAC,4BAA4B,CAAC,OAAO,EAAE,OAAO,EAAE,KAAK,CAAC,CAAC;YAEhE,IAAI,CAAC,IAAI,CAAC,4BAA4B,CAAC,OAAO,EAAE,OAAO,EAAE,KAAK,CAAC,CAAC;YAChE,IAAI,CAAC,IAAI,CAAC,4BAA4B,CAAC,OAAO,EAAE,OAAO,EAAE,KAAK,CAAC,CAAC;YAChE,IAAI,CAAC,IAAI,CAAC,4BAA4B,CAAC,OAAO,EAAE,OAAO,EAAE,KAAK,CAAC,CAAC;QAClE,CAAC;IACH,CAAC;IACD,0EAA0E;IACnE,4BAA4B,CAAC,KAAa,EAAE,WAAmB;QACpE,IAAI,CAAC,oBAAoB,CAAC,KAAK,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC;QACjD,MAAM,OAAO,GAAG,UAAU,CAAC,QAAQ,GAAG,iBAAO,CAAC,MAAM,CAClD,KAAK,CAAC,CAAC,GAAG,IAAI,CAAC,MAAM,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,GAAG,IAAI,CAAC,MAAM,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,GAAG,IAAI,CAAC,MAAM,CAAC,CAAC,EAAE,GAAG,EAAE,UAAU,CAAC,QAAQ,CACpG,CAAC;QACF,IAAI,CAAC,IAAI,CAAC,4BAA4B,CAAC,OAAO,EAAE,OAAO,EAAE,WAAW,CAAC,CAAC;IACxE,CAAC;IACD,4DAA4D;IACrD,wBAAwB,CAAC,MAAe,EAAE,MAAe;QAC9D,IAAI,CAAC,oBAAoB,CAAC,MAAM,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC,CAAC,CAAC;QACxD,MAAM,EAAE,GAAG,IAAI,CAAC,MAAM,CAAC,CAAC,CAAC;QACzB,MAAM,EAAE,GAAG,IAAI,CAAC,MAAM,CAAC,CAAC,CAAC;QACzB,MAAM,EAAE,GAAG,IAAI,CAAC,MAAM,CAAC,CAAC,CAAC;QACzB,MAAM,OAAO,GAAG,UAAU,CAAC,QAAQ,GAAG,iBAAO,CAAC,MAAM,CAClD,MAAM,CAAC,CAAC,GAAG,EAAE,EAAE,MAAM,CAAC,CAAC,GAAG,EAAE,EAAE,MAAM,CAAC,CAAC,GAAG,EAAE,EAAE,GAAG,EAAE,UAAU,CAAC,QAAQ,CACtE,CAAC;QACF,MAAM,OAAO,GAAG,UAAU,CAAC,QAAQ,GAAG,iBAAO,CAAC,MAAM,CAClD,MAAM,CAAC,CAAC,GAAG,EAAE,EAAE,MAAM,CAAC,CAAC,GAAG,EAAE,EAAE,MAAM,CAAC,CAAC,GAAG,EAAE,EAAE,GAAG,EAAE,UAAU,CAAC,QAAQ,CACtE,CAAC;QACF,MAAM,IAAI,GAAG,MAAM,CAAC,QAAQ,CAAC,MAAM,CAAC,CAAC;QACrC,MAAM,IAAI,GAAG,IAAI,GAAG,GAAG,CAAC;QACxB,MAAM,IAAI,GAAG,IAAI,GAAG,GAAG,CAAC;QACxB,IAAI,CAAC,IAAI,CAAC,4BAA4B,CAAC,OAAO,EAAE,OAAO,EAAE,IAAI,CAAC,CAAC;QAC/D,IAAI,CAAC,IAAI,CAAC,4BAA4B,CAAC,OAAO,EAAE,OAAO,EAAE,IAAI,CAAC,CAAC;QAC/D,IAAI,CAAC,IAAI,CAAC,4BAA4B,CAAC,OAAO,EAAE,OAAO,EAAE,IAAI,CAAC,CAAC;QAC/D,IAAI,CAAC,IAAI,CAAC,4BAA4B,CAAC,OAAO,EAAE,OAAO,EAAE,IAAI,CAAC,CAAC;IACjE,CAAC;IACD;;;;OAIG;IACI,uCAAuC,CAAC,SAA4B,EAAE,MAAwB;QACnG,MAAM,CAAC,GAAG,MAAM,CAAC,MAAM,CAAC;QACxB,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC;YACV,MAAM,CAAC,+BAA+B,CAAC,CAAC,EAAE,IAAI,CAAC,OAAO,CAAC,CAAC;YACxD,uGAAuG;YACvG,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;gBAC3B,MAAM,CAAC,+BAA+B,CAAC,CAAC,EAAE,IAAI,CAAC,OAAO,CAAC,CAAC;gBACxD,IAAI,CAAC,2BAA2B,CAAC,SAAS,EAAE,IAAI,CAAC,OAAO,EAAE,IAAI,CAAC,OAAO,CAAC,CAAC;gBACxE,IAAI,CAAC,OAAO,CAAC,cAAc,CAAC,IAAI,CAAC,OAAO,CAAC,CAAC;YAC5C,CAAC;QACH,CAAC;IACH,CAAC;IACD;;;;;;;;;;;;OAYG;IACI,qCAAqC,CAC1C,SAAiB,EAAE,SAAiB,EAAE,SAAiB,EACvD,IAAY,EAAE,MAAa,EAAE,OAAc,EAAE,OAAc;QAE3D,MAAM,kBAAkB,GAAG,mBAAQ,CAAC,eAAe,CACjD,SAAS,EAAE,SAAS,EAAE,CAAC,EAAE,CAAC,EAC1B,SAAS,EAAE,SAAS,EAAE,CAAC,EAAE,CAAC,EAC1B,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EACV,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,IAAI,CACd,CAAC;QACF,MAAM,IAAI,GAAG,mBAAQ,CAAC,gBAAgB,CAAC,OAAO,CAAC,CAAC,EAAE,OAAO,CAAC,CAAC,EAAE,OAAO,CAAC,CAAC,EAAE,OAAO,CAAC,CAAC,CAAC,CAAC;QACnF,MAAM,SAAS,GAAG,mBAAQ,CAAC,eAAe,CACxC,OAAO,CAAC,CAAC,EAAE,OAAO,CAAC,CAAC,EAAE,CAAC,EAAE,MAAM,CAAC,CAAC,GAAG,IAAI,CAAC,MAAM,CAAC,CAAC,EACjD,OAAO,CAAC,CAAC,EAAE,OAAO,CAAC,CAAC,EAAE,CAAC,EAAE,MAAM,CAAC,CAAC,GAAG,IAAI,CAAC,MAAM,CAAC,CAAC,EACjD,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EACV,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CACX,CAAC;QACF,MAAM,EAAE,GAAG,SAAS,CAAC,oBAAoB,CAAC,kBAAkB,CAAC,CAAC;QAC9D,MAAM,IAAI,GAAG,EAAE,CAAC,6BAA6B,CAAC,SAAS,CAAC,CAAC;QACzD,IAAI,CAAC,IAAI,CAAC,gBAAgB,CAAC,IAAI,EAAE,IAAI,CAAC,CAAC;IACzC,CAAC;IACD;;;;OAIG;IACI,kBAAkB,CAAC,KAAiB,EAAE,KAAa;QACxD,IAAI,CAAC,iBAAiB,CAAC,KAAK,CAAC,MAAM,CAAC,CAAC;QACrC,IAAI,CAAC,IAAI,CAAC,6BAA6B,CACrC,KAAK,CAAC,IAAI,EAAE,IAAI,CAAC,MAAM,CAAC,CAAC,GAAG,KAAK,CAAC,MAAM,CAAC,CAAC,EAAE,IAAI,CAAC,MAAM,CAAC,CAAC,GAAG,KAAK,CAAC,MAAM,CAAC,CAAC,EAAE,IAAI,CAAC,MAAM,CAAC,CAAC,GAAG,KAAK,CAAC,MAAM,CAAC,CAAC,EAAE,KAAK,CAClH,CAAC;IACJ,CAAC;IACD;;;;OAIG;IACI,4BAA4B,CAAC,MAAe,EAAE,QAAkB,EAAE,KAAa;QACpF,IAAI,CAAC,iBAAiB,CAAC,MAAM,CAAC,CAAC;QAC/B,IAAI,CAAC,IAAI,CAAC,6BAA6B,CACrC,QAAQ,EAAE,IAAI,CAAC,MAAM,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,EAAE,IAAI,CAAC,MAAM,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,EAAE,IAAI,CAAC,MAAM,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,EAAE,KAAK,CAC9F,CAAC;IACJ,CAAC;IACD,qCAAqC;IAC9B,MAAM;QACX,OAAO;YACL,MAAM,EAAE,IAAI,CAAC,MAAM;YACnB,IAAI,EAAE,IAAI,CAAC,IAAI,CAAC,MAAM,EAAE;YACxB,gBAAgB,EAAE,IAAI,CAAC,gBAAgB,CAAC,MAAM,EAAE;YAChD,YAAY,EAAE,IAAI,CAAC,eAAe,CAAC,MAAM,EAAE;SAC5C,CAAC;IACJ,CAAC;CACF;AAtZD,gCAsZC","sourcesContent":["/*---------------------------------------------------------------------------------------------\n* Copyright (c) Bentley Systems, Incorporated. All rights reserved.\n* See LICENSE.md in the project root for license terms and full copyright notice.\n*--------------------------------------------------------------------------------------------*/\n\n/** @packageDocumentation\n * @module Numerics\n */\n\n/* eslint-disable @typescript-eslint/naming-convention */\n// cspell:word ABAT\n\nimport { Geometry } from \"../Geometry\";\nimport { GrowableXYZArray } from \"../geometry3d/GrowableXYZArray\";\nimport { Matrix3d } from \"../geometry3d/Matrix3d\";\nimport { Point3d, Vector3d, XYZ } from \"../geometry3d/Point3dVector3d\";\nimport { Transform } from \"../geometry3d/Transform\";\nimport { XAndY, XYAndZ } from \"../geometry3d/XYZProps\";\nimport { Matrix4d } from \"./Matrix4d\";\nimport { Point4d } from \"./Point4d\";\n\n/**\n * A MomentData structure exists in several levels:\n * * First level: as a carrier of sums of inertial products that determine moments.\n *   * origin = local origin used as moments are summed.\n *   * sums = array of summed moments.\n *     * The [i,j] entry of the sums is a summed or integrated moment for product of axis i and j.\n *       * axes 0,1,2 are x,y,z.\n *         * e.g. entry [0,1] is summed product xy\n *       * axis 3 is \"w\", which is 1 in sums.\n *         * e.g. entry 03 is summed x.\n *    * In this level:\n *        * the `absoluteQuantity` member is undefined.\n *        * the `localToWorldMap` and `radiiOfGyration` are created but have undefined contents.\n *  * Second level: after a call to `inertiaProductsToPrincipalAxes`, the `localToWorldMap`, `absoluteQuantity` and\n *    `radiiOfGyration` are filled in.\n * @public\n */\nexport class MomentData {\n  /** Origin used for sums. */\n  public origin: Point3d;\n  /** Flag to request deferred origin setup. */\n  public needOrigin: boolean;\n  /**\n   * Moment sums.\n   * * Set to zero at initialization and if requested later.\n   * * Accumulated during data entry phase.\n   */\n  public sums: Matrix4d;\n  /**\n   * The mapping between principal and world system.\n   * * This set up with its inverse already constructed.\n   */\n  public localToWorldMap: Transform;\n  /** Radii of gyration (square roots of principal second moments). */\n  public radiusOfGyration: Vector3d;\n  /**\n   * Principal quantity (e.g. length, area, or volume). This is undefined in raw moments, and becomes defined by\n   * inertiaProductsToPrincipalAxes.\n   */\n  public absoluteQuantity?: number;\n  // private variables\n  private static _vectorA?: Point4d;\n  private static _vectorB?: Point4d;\n  private static _vectorC?: Point4d;\n  private _point0 = Point3d.create();\n  private _point1 = Point3d.create();\n  /** Constructor. */\n  private constructor() {\n    this.origin = Point3d.createZero();\n    this.needOrigin = false;\n    this.sums = Matrix4d.createZero();\n    this.localToWorldMap = Transform.createIdentity();\n    this.radiusOfGyration = Vector3d.create();\n    this.absoluteQuantity = 0.1; // so optimizer sees its type\n    this.absoluteQuantity = undefined;\n  }\n  /**\n   * Return the lower-right (3,3) entry in the sums.\n   * * This is the quantity (i.e. length, area, or volume) summed.\n   */\n  public get quantitySum(): number {\n    return this.sums.weight();\n  }\n  /**\n   * Return a scale factor to make these sums match the target orientation sign.\n   * * 1.0 if `this.quantitySum` has the same sign as `targetSign`.\n   * * -1.0 if `this.quantitySum` has the opposite sign from `targetSign`\n   */\n  public signFactor(targetSign: number): number {\n    return targetSign * this.quantitySum > 0 ? 1.0 : -1.0;\n  }\n  /** If `this.needOrigin` flag is set, copy `origin` to `this.origin` and clear the flag. */\n  public setOriginIfNeeded(origin: Point3d) {\n    if (this.needOrigin) {\n      this.origin.setFromPoint3d(origin);\n      this.needOrigin = false;\n    }\n  }\n  /** If `this.needOrigin` flag is set, copy `origin` to `this.origin` and clear the flag. */\n  public setOriginFromGrowableXYZArrayIfNeeded(points: GrowableXYZArray) {\n    if (this.needOrigin && points.length > 0) {\n      points.getPoint3dAtCheckedPointIndex(0, this.origin);\n      this.needOrigin = false;\n    }\n  }\n  /** If `this.needOrigin` flag is set, copy `origin` to `this.origin` and clear the flag. */\n  public setOriginXYZIfNeeded(x: number, y: number, z: number) {\n    if (this.needOrigin) {\n      this.origin.set(x, y, z);\n      this.needOrigin = false;\n    }\n  }\n  /**\n   * Create moments with optional origin.\n   * * Origin and needOrigin are quirky.\n   *   * (undefined, true) sets up to use first incoming point as origin.\n   *   * (origin) sets up to use that durable origin, set needsOrigin flag false.\n   *   * (origin, true) the \"true\" is meaningless.\n   *   * (undefined, false) makes 000 the durable origin.\n  */\n  public static create(origin?: Point3d | undefined, needOrigin: boolean = false): MomentData {\n    const data = new MomentData();\n    data.needOrigin = needOrigin;\n    if (origin) {\n      data.origin.setFromPoint3d(origin);\n      data.needOrigin = false;\n    }\n    return data;\n  }\n  /**\n   * Return the formal tensor of integrated values `[yy+zz,xy,xz][yx,xx+zz,yz][zx,xy,xx+yy]`.\n   * @param products matrix of (integrated) `[xx,xy,xz][yx,yy,yz][zx,xy,zz]`.\n   */\n  public static momentTensorFromInertiaProducts(products: Matrix3d): Matrix3d {\n    const rr = products.sumDiagonal();\n    const result = Matrix3d.createScale(rr, rr, rr);\n    result.addScaledInPlace(products, -1.0);\n    return result;\n  }\n  /** Sort the columns of the matrix for increasing moments. */\n  public static sortColumnsForIncreasingMoments(axes: Matrix3d, moments: Vector3d) {\n    const points = [\n      axes.indexedColumnWithWeight(0, moments.x),\n      axes.indexedColumnWithWeight(1, moments.y),\n      axes.indexedColumnWithWeight(2, moments.z),\n    ].sort((dataA: Point4d, dataB: Point4d): number => {\n      if (dataA.w < dataB.w) return -1;\n      if (dataA.w > dataB.w) return 1;\n      return 0;\n    });\n    axes.setColumnsPoint4dXYZ(points[0], points[1], points[2]);\n    if (axes.determinant() < 0)\n      axes.scaleColumnsInPlace(-1.0, -1.0, -1.0);\n    // prefer x and z positive; y falls wherever\n    if (axes.at(0, 0) < 0.0)\n      axes.scaleColumnsInPlace(-1.0, -1.0, 1.0);\n    if (axes.at(2, 2) < 0.0)\n      axes.scaleColumnsInPlace(1.0, -1.0, -1.0);\n    moments.set(points[0].w, points[1].w, points[2].w);\n  }\n  /**\n   * Return the principal moment data for an array of points.\n   * @param points array of points.\n   */\n  public static pointsToPrincipalAxes(points: Point3d[]): MomentData | undefined {\n    const moments = new MomentData();\n    if (points.length === 0)\n      return moments;\n    moments.clearSums(points[0]);\n    moments.accumulatePointMomentsFromOrigin(points);\n    return this.inertiaProductsToPrincipalAxes(moments.origin, moments.sums);\n  }\n  /**\n   * Compute principal axes from inertial products.\n   * * The radii of gyration are sorted smallest to largest.\n   * * Hence x axis is long direction.\n   * * Hence planar data generates large moment as Z.\n   * @param origin The origin used for the inertia products.\n   * @param inertiaProducts The inertia products: sums or integrals of\n   * [xx,xy,xz,xw; yx,yy,yz,yw; zx,zy,zz,zw; wx,wy,wz,w].\n   */\n  public static inertiaProductsToPrincipalAxes(origin: XYZ, inertiaProducts: Matrix4d): MomentData | undefined {\n    const moments = new MomentData();\n    moments.sums.setFrom(inertiaProducts);\n    moments.origin.setFrom(origin);\n    if (!moments.shiftOriginAndSumsToCentroidOfSums())\n      return undefined;\n    const products = moments.sums.matrixPart();\n    const w = moments.sums.weight();\n    if (w < 0.0)\n      products.scaleColumnsInPlace(-1, -1, -1);\n    const tensor = MomentData.momentTensorFromInertiaProducts(products);\n    const moment2 = Vector3d.create();\n    const axisVectors = Matrix3d.createZero();\n    tensor.fastSymmetricEigenvalues(axisVectors, moment2);\n    if (moment2.x < 0.0)\n      return undefined;\n    MomentData.sortColumnsForIncreasingMoments(axisVectors, moment2);\n    if (w < 0.0)\n      axisVectors.scaleColumnsInPlace(1, -1, -1);\n    moments.localToWorldMap = Transform.createOriginAndMatrix(moments.origin, axisVectors);\n    moments.radiusOfGyration.set(\n      Math.sqrt(Math.abs(moment2.x)), Math.sqrt(Math.abs(moment2.y)), Math.sqrt(Math.abs(moment2.z)),\n    );\n    moments.radiusOfGyration.scaleInPlace(1.0 / Math.sqrt(Math.abs(w)));\n    moments.absoluteQuantity = Math.abs(w);\n    return moments;\n  }\n  /**\n   * Test for match among selected members as they exist after `inertiaProductsToPrincipalAxes`.\n   * * The members considered are:\n   *   * origin of local to world map (i.e. centroid),\n   *   * radius of gyration,\n   *   * axes of localToWorldMap.\n   * * Axis direction tests allow these quirks:\n   *   * opposite orientation is considered matched.\n   * * Full xyz symmetry: If x,y,z radii are matched, axes are not tested.\n   * * Symmetry in xy plane: If x and y radii are matched, the x and y axes area allowed to spin freely. Only Z direction\n   *   is tested.\n   * * If either or both are undefined, returns false.\n   * @param dataA first set of moments.\n   * @param dataB second set of moments.\n   */\n  public static areEquivalentPrincipalAxes(dataA: MomentData | undefined, dataB: MomentData | undefined): boolean {\n    if (dataA && dataB\n      && Geometry.isSameCoordinate(dataA.quantitySum, dataB.quantitySum)) { // TODO: need different tolerance for area, volume?\n      if (dataA.localToWorldMap.getOrigin().isAlmostEqual(dataB.localToWorldMap.getOrigin())\n        && dataA.radiusOfGyration.isAlmostEqual(dataB.radiusOfGyration)) {\n        if (Geometry.isSameCoordinate(dataA.radiusOfGyration.x, dataA.radiusOfGyration.y)) {\n          // we have at least xy symmetry\n          if (Geometry.isSameCoordinate(dataA.radiusOfGyration.x, dataA.radiusOfGyration.z))\n            return true;\n          // just xy; allow opposite z directions; if the z's are aligned, x and y can spin freely\n          const zA = dataA.localToWorldMap.matrix.columnZ();\n          const zB = dataB.localToWorldMap.matrix.columnZ();\n          if (zA.isParallelTo(zB, true))\n            return true;\n          return false;\n        }\n        // no symmetry; test all three axes\n        const vectorA = Vector3d.create();\n        const vectorB = Vector3d.create();\n        for (let i = 0; i < 3; i++) {\n          dataA.localToWorldMap.matrix.getColumn(i, vectorA);\n          dataB.localToWorldMap.matrix.getColumn(i, vectorB);\n          if (!vectorA.isParallelTo(vectorB, true))\n            return false;\n        }\n        return true;\n      }\n    }\n    return false;\n  }\n  /** Clear the MomentData sums to zero, and establish a new origin. */\n  public clearSums(origin?: Point3d) {\n    this.sums.setZero();\n    if (origin)\n      this.origin.setFrom(origin);\n    else\n      this.origin.setZero();\n  }\n  /** Accumulate products-of-components for given points. */\n  public accumulatePointMomentsFromOrigin(points: Point3d[]) {\n    for (const p of points) {\n      this.sums.addMomentsInPlace(\n        p.x - this.origin.x,\n        p.y - this.origin.y,\n        p.z - this.origin.z,\n        1.0,\n      );\n    }\n  }\n  /** Revise the accumulated sums to be \"around the centroid\". */\n  public shiftOriginAndSumsToCentroidOfSums(): boolean {\n    const xyz = this.sums.columnW().realPoint(); // centroid of the geometry\n    if (xyz) {\n      this.shiftOriginAndSumsByXYZ(xyz.x, xyz.y, xyz.z);\n      return true;\n    }\n    return false;\n  }\n  /**\n   * Revise the accumulated sums.\n   * * Add (ax,ay,az) to the origin coordinates.\n   * * Apply the negative translation to the sums.\n   */\n  public shiftOriginAndSumsByXYZ(ax: number, ay: number, az: number) {\n    this.origin.addXYZInPlace(ax, ay, az);\n    this.sums.multiplyTranslationSandwichInPlace(-ax, -ay, -az);\n  }\n  /** Revise the accumulated sums so they are based at a specified origin. */\n  public shiftOriginAndSumsToNewOrigin(newOrigin: XYAndZ) {\n    this.shiftOriginAndSumsByXYZ(newOrigin.x - this.origin.x, newOrigin.y - this.origin.y, newOrigin.z - this.origin.z);\n  }\n  /**\n   * Compute moments of a triangle from the origin. Accumulate them to `this.sums`.\n   * * If `this.needOrigin` is set, `this.origin` is set to `pointB`.\n   * * If `pointA` is undefined, use `this.origin` as `pointA`.\n   */\n  public accumulateTriangleMomentsXY(pointA: XAndY | undefined, pointB: XAndY, pointC: XAndY) {\n    this.setOriginXYZIfNeeded(pointB.x, pointB.y, 0.0);\n    const x0 = this.origin.x;\n    const y0 = this.origin.y;\n    const vectorA = MomentData._vectorA = (pointA !== undefined) ?\n      Point4d.create(pointA.x - x0, pointA.y - y0, 0.0, 1.0, MomentData._vectorA) :\n      Point4d.create(0.0, 0.0, 0.0, 1.0, MomentData._vectorA);\n    const vectorB = MomentData._vectorB = Point4d.create(pointB.x - x0, pointB.y - y0, 0.0, 1.0, MomentData._vectorB);\n    const vectorC = MomentData._vectorC = Point4d.create(pointC.x - x0, pointC.y - y0, 0.0, 1.0, MomentData._vectorC);\n    // Below we calculate 16 double integrals: \\iint_T [x y 0 1]^ [x y 0 1] dT over triangle T=(A,B,C).\n    // Each accumulates contributions from 9 scaled outer products. Integration computations use the barycentric\n    // change of variables [B-A C-A][u,v]^ = [x,y]^ with Jacobian detJ = B-A x C-A = twice the area of T.\n    // This converts the integration domain from T to the triangle bounded by u=0, v=0 and v=1-u, yielding e.g.,\n    // \\iint_T x^2 dT = detJ \\int_0^1 \\int_0^{1-u} u^2 dv du = detJ / 12, and similarly \\iint_T xy dT = detJ / 24.\n    const detJ = Geometry.crossProductXYXY(\n      vectorB.x - vectorA.x, vectorB.y - vectorA.y, vectorC.x - vectorA.x, vectorC.y - vectorA.y,\n    );\n    if (detJ !== 0.0) {\n      const r1_12 = detJ / 12.0;\n      const r1_24 = detJ / 24.0;\n\n      this.sums.addScaledOuterProductInPlace(vectorA, vectorA, r1_12);\n      this.sums.addScaledOuterProductInPlace(vectorA, vectorB, r1_24);\n      this.sums.addScaledOuterProductInPlace(vectorA, vectorC, r1_24);\n\n      this.sums.addScaledOuterProductInPlace(vectorB, vectorA, r1_24);\n      this.sums.addScaledOuterProductInPlace(vectorB, vectorB, r1_12);\n      this.sums.addScaledOuterProductInPlace(vectorB, vectorC, r1_24);\n\n      this.sums.addScaledOuterProductInPlace(vectorC, vectorA, r1_24);\n      this.sums.addScaledOuterProductInPlace(vectorC, vectorB, r1_24);\n      this.sums.addScaledOuterProductInPlace(vectorC, vectorC, r1_12);\n    }\n  }\n  /** Add scaled outer product of (4d, unit weight) point to `this.sums`. */\n  public accumulateScaledOuterProduct(point: XYAndZ, scaleFactor: number) {\n    this.setOriginXYZIfNeeded(point.x, point.y, 0.0);\n    const vectorA = MomentData._vectorA = Point4d.create(\n      point.x - this.origin.x, point.y - this.origin.y, point.z - this.origin.z, 1.0, MomentData._vectorA,\n    );\n    this.sums.addScaledOuterProductInPlace(vectorA, vectorA, scaleFactor);\n  }\n  /** Accumulate wire moment integral from pointA to pointB */\n  public accumulateLineMomentsXYZ(pointA: Point3d, pointB: Point3d) {\n    this.setOriginXYZIfNeeded(pointA.x, pointA.y, pointA.z);\n    const x0 = this.origin.x;\n    const y0 = this.origin.y;\n    const z0 = this.origin.z;\n    const vectorA = MomentData._vectorA = Point4d.create(\n      pointA.x - x0, pointA.y - y0, pointA.z - z0, 1.0, MomentData._vectorA,\n    );\n    const vectorB = MomentData._vectorB = Point4d.create(\n      pointB.x - x0, pointB.y - y0, pointB.z - z0, 1.0, MomentData._vectorB,\n    );\n    const detJ = pointA.distance(pointB);\n    const r1_3 = detJ / 3.0;\n    const r1_6 = detJ / 6.0;\n    this.sums.addScaledOuterProductInPlace(vectorA, vectorA, r1_3);\n    this.sums.addScaledOuterProductInPlace(vectorA, vectorB, r1_6);\n    this.sums.addScaledOuterProductInPlace(vectorB, vectorA, r1_6);\n    this.sums.addScaledOuterProductInPlace(vectorB, vectorB, r1_3);\n  }\n  /**\n   * Compute moments of triangles from a base point to the given linestring. Accumulate them to `this.sums`.\n   * * If `this.needOrigin` is set, `this.origin` is set to the first point of the array.\n   * * If `sweepBase` is undefined, use `this.origin` as `sweepBase`.\n   */\n  public accumulateTriangleToLineStringMomentsXY(sweepBase: XAndY | undefined, points: GrowableXYZArray) {\n    const n = points.length;\n    if (n > 1) {\n      points.getPoint3dAtUncheckedPointIndex(0, this._point0);\n      // The linestring forms a polygon with sweepBase. Integrate over this polygon using Shoelace algorithm.\n      for (let i = 1; i < n; i++) {\n        points.getPoint3dAtUncheckedPointIndex(i, this._point1);\n        this.accumulateTriangleMomentsXY(sweepBase, this._point0, this._point1);\n        this._point0.setFromPoint3d(this._point1);\n      }\n    }\n  }\n  /**\n   * Assemble XX, YY, XY products into a full matrix form [xx,xy,0,0; xy,yy,0,0; 0,0,0,0; 0,0,0,1].\n   * * Sandwich this between transforms with columns [vectorU, vectorV, 0000, origin].\n   *   (column weights 0001; only xy parts of vectors).\n   * * Scale by detJ for the xy-only determinant of the vectors.\n   * @param productXX\n   * @param productXY\n   * @param productYY\n   * @param area area in caller's system.\n   * @param origin caller's origin.\n   * @param vectorU caller's U axis (not necessarily unit).\n   * @param vectorV caller's V axis (not necessarily unit).\n   */\n  public accumulateXYProductsInCentroidalFrame(\n    productXX: number, productXY: number, productYY: number,\n    area: number, origin: XAndY, vectorU: XAndY, vectorV: XAndY,\n  ) {\n    const centroidalProducts = Matrix4d.createRowValues(\n      productXX, productXY, 0, 0,\n      productXY, productYY, 0, 0,\n      0, 0, 0, 0,\n      0, 0, 0, area,\n    );\n    const detJ = Geometry.crossProductXYXY(vectorU.x, vectorV.x, vectorU.y, vectorV.y);\n    const placement = Matrix4d.createRowValues(\n      vectorU.x, vectorV.x, 0, origin.x - this.origin.x,\n      vectorU.y, vectorV.y, 0, origin.y - this.origin.y,\n      0, 0, 0, 0,\n      0, 0, 0, 1,\n    );\n    const AB = placement.multiplyMatrixMatrix(centroidalProducts);\n    const ABAT = AB.multiplyMatrixMatrixTranspose(placement);\n    this.sums.addScaledInPlace(ABAT, detJ);\n  }\n  /**\n   * Accumulate sums from other moments.\n   * * Scale by given scaleFactor (e.g. sign to correct orientation).\n   * * Pull the origin from `other` if `this` needs an origin.\n   */\n  public accumulateProducts(other: MomentData, scale: number) {\n    this.setOriginIfNeeded(other.origin);\n    this.sums.addTranslationSandwichInPlace(\n      other.sums, this.origin.x - other.origin.x, this.origin.y - other.origin.y, this.origin.z - other.origin.z, scale,\n    );\n  }\n  /**\n   * Accumulate sums from Matrix4d and origin.\n   * * Scale by given scaleFactor (e.g. sign to correct orientation).\n   * * Trap the origin if `this` needs an origin.\n   */\n  public accumulateProductsFromOrigin(origin: Point3d, products: Matrix4d, scale: number) {\n    this.setOriginIfNeeded(origin);\n    this.sums.addTranslationSandwichInPlace(\n      products, this.origin.x - origin.x, this.origin.y - origin.y, this.origin.z - origin.z, scale,\n    );\n  }\n  /** Convert to a json data object. */\n  public toJSON(): any {\n    return {\n      origin: this.origin,\n      sums: this.sums.toJSON(),\n      radiusOfGyration: this.radiusOfGyration.toJSON(),\n      localToWorld: this.localToWorldMap.toJSON(),\n    };\n  }\n}\n"]}

package/lib/cjs/geometry4d/Point4d.d.ts

@@ -61,7 +61,7 @@ export declare class Point4d extends Plane3d implements BeJSONFunctions {
     /** Create a new point with coordinates from a json array `[x,y,z,w]` */
     static fromJSON(json?: Point4dProps): Point4d;
     /** Near-equality test, using `Geometry.isSameCoordinate` on all 4 x,y,z,w */
-    isAlmostEqual(other: Point4d): boolean;
+    isAlmostEqual(other: Point4d, tolerance?: number): boolean;
     /**
      * Test for same coordinate by direct x,y,z,w args
      * @param x x to test
@@ -69,7 +69,7 @@ export declare class Point4d extends Plane3d implements BeJSONFunctions {
      * @param z z to test
      * @param w w to test
      */
-    isAlmostEqualXYZW(x: number, y: number, z: number, w: number): boolean;
+    isAlmostEqualXYZW(x: number, y: number, z: number, w: number, tolerance?: number): boolean;
     /**
      * Convert an Angle to a JSON object.
      * @return {*} [x,y,z,w]
@@ -83,9 +83,10 @@ export declare class Point4d extends Plane3d implements BeJSONFunctions {
      * * x,y,z,w all participate without normalization.
      */
     distanceSquaredXYZW(other: Point4d): number;
-    /** Return the distance between the instance and other after normalizing by weights
-     */
+    /** Return the xy distance between the instance and `other` after normalizing by weights */
     realDistanceXY(other: Point4d): number | undefined;
+    /** Return the xy squared distance between the instance and `other` after normalizing by weights */
+    realDistanceSquaredXY(other: Point4d): number | undefined;
     /** Return the largest absolute distance between corresponding components
      * * x,y,z,w all participate without normalization.
      */
@@ -130,7 +131,7 @@ export declare class Point4d extends Plane3d implements BeJSONFunctions {
      * extract 4 consecutive numbers from a Float64Array into a Point4d.
      * @param data buffer of numbers
      * @param xIndex first index for x,y,z,w sequence. Assumed to be a valid index!
-     * @deprecated in 4.x. Use createFromPacked instead.
+     * @deprecated in 4.x - will not be removed until after 2026-06-13. Use createFromPacked instead.
      */
     static createFromPackedXYZW(data: Float64Array, xIndex?: number, result?: Point4d): Point4d;
     /** Create a `Point4d` with x,y,z from an `XYAndZ` input, and w from a separate number. */
@@ -199,6 +200,8 @@ export declare class Point4d extends Plane3d implements BeJSONFunctions {
      * * If the xyz part of `this` are all zero, (a clone of) `spacePoint` is returned.
      */
     projectPointToPlane(spacePoint: Point3d, result?: Point3d): Point3d;
+    /** Return the projection of (x,y,z) onto the plane. */
+    projectXYZToPlane(x: number, y: number, z: number, result?: Point3d): Point3d;
     /** scale all components (including w!!) */
     scale(scale: number, result?: Point4d): Point4d;
     /** Negate components (including w!!) */

package/lib/cjs/geometry4d/Point4d.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"Point4d.d.ts","sourceRoot":"","sources":["../../../src/geometry4d/Point4d.ts"],"names":[],"mappings":"AAKA;;GAEG;AACH,OAAO,EAAE,eAAe,EAAY,sBAAsB,EAAE,MAAM,aAAa,CAAC;AAChF,OAAO,EAAE,OAAO,EAAE,MAAM,uBAAuB,CAAC;AAChD,OAAO,EAAE,4BAA4B,EAAE,MAAM,4CAA4C,CAAC;AAC1F,OAAO,EAAE,yBAAyB,EAAE,MAAM,yCAAyC,CAAC;AAEpF,OAAO,EAAE,OAAO,EAAE,QAAQ,EAAE,MAAM,+BAA+B,CAAC;AAClE,OAAO,EAAE,KAAK,EAAE,MAAM,qBAAqB,CAAC;AAC5C,OAAO,EAAE,KAAK,EAAE,MAAM,EAAE,MAAM,wBAAwB,CAAC;AAEvD;;;GAGG;AACH,MAAM,MAAM,YAAY,GAAG,MAAM,EAAE,CAAC;AAgBpC;;;;;;;GAOG;AACH,qBAAa,OAAQ,SAAQ,OAAQ,YAAW,eAAe;IAC7D,6CAA6C;IACtC,IAAI,EAAE,YAAY,CAAC;IAC1B,kCAAkC;IAC3B,GAAG,CAAC,CAAC,GAAE,MAAU,EAAE,CAAC,GAAE,MAAU,EAAE,CAAC,GAAE,MAAU,EAAE,CAAC,GAAE,MAAU,GAAG,OAAO;IAO/E;;OAEG;IACI,YAAY,CAAC,KAAK,EAAE,MAAM,EAAE,KAAK,EAAE,MAAM;IAKhD,uBAAuB;IACvB,IAAW,CAAC,IACM,MAAM,CADe;IACvC,IAAW,CAAC,CAAC,GAAG,EAAE,MAAM,EAAyB;IACjD,uBAAuB;IACvB,IAAW,CAAC,IACM,MAAM,CADe;IACvC,IAAW,CAAC,CAAC,GAAG,EAAE,MAAM,EAAyB;IACjD,uBAAuB;IACvB,IAAW,CAAC,IACM,MAAM,CADe;IACvC,IAAW,CAAC,CAAC,GAAG,EAAE,MAAM,EAAyB;IACjD,qCAAqC;IACrC,IAAW,CAAC,IACM,MAAM,CADe;IACvC,IAAW,CAAC,CAAC,GAAG,EAAE,MAAM,EAAyB;IACjD,kCAAkC;IAClC,SAAS,aAAa,CAAC,GAAE,MAAU,EAAE,CAAC,GAAE,MAAU,EAAE,CAAC,GAAE,MAAU,EAAE,CAAC,GAAE,MAAU;IAQhF,8CAA8C;WAChC,MAAM,CAAC,CAAC,GAAE,MAAU,EAAE,CAAC,GAAE,MAAU,EAAE,CAAC,GAAE,MAAU,EAAE,CAAC,GAAE,MAAU,EAAE,MAAM,CAAC,EAAE,OAAO,GAAG,OAAO;IAG3G;;;;OAIG;WACW,eAAe,CAAC,MAAM,EAAE,sBAAsB,GAAG,OAAO,GAAG,SAAS;IAGlF,qCAAqC;IAC9B,OAAO,CAAC,KAAK,EAAE,OAAO,GAAG,OAAO;IAOvC,uBAAuB;IAChB,KAAK,CAAC,MAAM,CAAC,EAAE,OAAO,GAAG,OAAO;IAGvC,0DAA0D;IACnD,WAAW,CAAC,IAAI,CAAC,EAAE,YAAY;IAOtC,wEAAwE;WAC1D,QAAQ,CAAC,IAAI,CAAC,EAAE,YAAY,GAAG,OAAO;IAKpD,6EAA6E;IACtE,aAAa,CAAC,KAAK,EAAE,OAAO,GAAG,OAAO;IAM7C;;;;;;OAMG;IACI,iBAAiB,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,GAAG,OAAO;IAO7E;;;OAGG;IACI,MAAM,IAAI,YAAY;IAG7B;;OAEG;IACI,YAAY,CAAC,KAAK,EAAE,OAAO,GAAG,MAAM;IAG3C;;OAEG;IACI,mBAAmB,CAAC,KAAK,EAAE,OAAO,GAAG,MAAM;IAGlD;OACG;IACI,cAAc,CAAC,KAAK,EAAE,OAAO,GAAG,MAAM,GAAG,SAAS;IAOzD;;OAEG;IACI,OAAO,CAAC,KAAK,EAAE,OAAO,GAAG,MAAM;IAGtC,oEAAoE;IAC7D,MAAM,IAAI,MAAM;IAGvB,+DAA+D;IACxD,aAAa,IAAI,MAAM;IAG9B,wHAAwH;IACjH,mBAAmB,IAAI,MAAM;IAIpC,wEAAwE;IACjE,KAAK,CAAC,KAAK,EAAE,OAAO,EAAE,MAAM,CAAC,EAAE,OAAO,GAAG,OAAO;IAGvD,sDAAsD;IAC/C,kBAAkB,CAAC,KAAK,EAAE,OAAO,EAAE,MAAM,CAAC,EAAE,QAAQ,GAAG,QAAQ;IAKtE,kFAAkF;IAC3E,yBAAyB,CAAC,KAAK,EAAE,OAAO,EAAE,MAAM,CAAC,EAAE,QAAQ,GAAG,QAAQ;IAI7E,uEAAuE;IAChE,IAAI,CAAC,KAAK,EAAE,OAAO,EAAE,MAAM,CAAC,EAAE,OAAO,GAAG,OAAO;IAGtD,8CAA8C;IAC9C,IAAW,YAAY,IAAI,OAAO,CAEjC;IACD,mDAAmD;WACrC,UAAU,IAAI,OAAO;IACnC;;;;OAIG;WACW,sBAAsB,CAAC,MAAM,EAAE,OAAO,EAAE,MAAM,EAAE,OAAO,EAAE,MAAM,CAAC,EAAE,OAAO;IAGvF;;;;OAIG;WACW,mBAAmB,CAAC,IAAI,EAAE,YAAY,EAAE,MAAM,GAAE,MAAU,EAAE,MAAM,CAAC,EAAE,OAAO,GAAG,OAAO,GAAG,SAAS;IAKhH;;;;OAIG;WACW,gBAAgB,CAAC,IAAI,EAAE,YAAY,EAAE,MAAM,GAAE,MAAU,EAAE,MAAM,CAAC,EAAE,OAAO,GAAG,OAAO,GAAG,SAAS;IAK7G;;;;;OAKG;WACW,oBAAoB,CAAC,IAAI,EAAE,YAAY,EAAE,MAAM,GAAE,MAAU,EAAE,MAAM,CAAC,EAAE,OAAO,GAAG,OAAO;IAGrG,0FAA0F;WAC5E,wBAAwB,CAAC,GAAG,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,GAAG,OAAO;IAGvE;;;;;;OAMG;WACW,eAAe,CAAC,KAAK,EAAE,KAAK,GAAG,MAAM,GAAG,OAAO,GAAG,MAAM,EAAE,GAAG,OAAO;IAsBlF,uCAAuC;IAChC,UAAU,CAAC,MAAM,EAAE,OAAO,EAAE,WAAW,EAAE,MAAM,EAAE,MAAM,CAAC,EAAE,OAAO,GAAG,OAAO;IAGlF;OACG;IACI,WAAW,CAAC,QAAQ,EAAE,MAAM,EAAE,MAAM,EAAE,OAAO,EAAE,MAAM,CAAC,EAAE,OAAO,GAAG,OAAO;IAIhF,6DAA6D;IACtD,WAAW,CAAC,OAAO,EAAE,OAAO,EAAE,OAAO,EAAE,MAAM,EAAE,OAAO,EAAE,OAAO,EAAE,OAAO,EAAE,MAAM,EAAE,MAAM,CAAC,EAAE,OAAO,GAAG,OAAO;IAGnH,iFAAiF;IAC1E,WAAW,CAAC,OAAO,EAAE,OAAO,EAAE,OAAO,EAAE,MAAM,EAAE,OAAO,EAAE,OAAO,EAAE,OAAO,EAAE,MAAM,EAAE,OAAO,EAAE,OAAO,EAAE,OAAO,EAAE,MAAM,EAAE,MAAM,CAAC,EAAE,OAAO,GAAG,OAAO;IAGtJ,6DAA6D;WAC/C,gBAAgB,CAAC,OAAO,EAAE,OAAO,EAAE,OAAO,EAAE,MAAM,EAAE,OAAO,EAAE,OAAO,EAAE,OAAO,EAAE,MAAM,EAAE,MAAM,CAAC,EAAE,OAAO,GAAG,OAAO;IAG/H,iFAAiF;WACnE,gBAAgB,CAAC,OAAO,EAAE,OAAO,EAAE,OAAO,EAAE,MAAM,EAAE,OAAO,EAAE,OAAO,EAAE,OAAO,EAAE,MAAM,EAAE,OAAO,EAAE,OAAO,EAAE,OAAO,EAAE,MAAM,EAAE,MAAM,CAAC,EAAE,OAAO,GAAG,OAAO;IAGlK,kFAAkF;IAC3E,mBAAmB,CAAC,OAAO,EAAE,OAAO,EAAE,OAAO,EAAE,OAAO,GAAG,MAAM;IAMtE,+DAA+D;IACxD,UAAU,CAAC,KAAK,EAAE,OAAO,GAAG,MAAM;IAGzC,wDAAwD;IACjD,cAAc,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,GAAG,MAAM;IAGzE,8FAA8F;IACvF,QAAQ,CAAC,KAAK,EAAE,OAAO,GAAG,MAAM;IAGvC,4EAA4E;IACrE,WAAW,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,GAAG,MAAM;IAG3D,oGAAoG;IAC7F,gBAAgB,CAAC,KAAK,EAAE,OAAO,GAAG,MAAM;IAG/C,mGAAmG;IAC5F,QAAQ,CAAC,MAAM,EAAE,QAAQ,GAAG,MAAM;IAGzC,4EAA4E;IACrE,WAAW,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,GAAG,MAAM;IAG3D;;OAEG;IACI,OAAO,IAAI,MAAM;IACxB;;MAEE;IACK,OAAO,IAAI,MAAM;IACxB;;MAEE;IACK,OAAO,IAAI,MAAM;IAExB,oBAAoB;WACN,KAAK,IAAI,OAAO;IAC9B,oBAAoB;WACN,KAAK,IAAI,OAAO;IAC9B,oBAAoB;WACN,KAAK,IAAI,OAAO;IAC9B,oBAAoB;WACN,KAAK,IAAI,OAAO;IAC9B,0EAA0E;IACnE,gBAAgB,CAAC,WAAW,EAAE,MAAM,EAAE,MAAM,CAAC,EAAE,OAAO,GAAG,OAAO,GAAG,SAAS;IAMnF;;;OAGG;IACI,mBAAmB,CAAC,UAAU,EAAE,OAAO,EAAE,MAAM,CAAC,EAAE,OAAO,GAAG,OAAO;IAS1E,2CAA2C;IACpC,KAAK,CAAC,KAAK,EAAE,MAAM,EAAE,MAAM,CAAC,EAAE,OAAO,GAAG,OAAO;IAQtD,wCAAwC;IACjC,MAAM,CAAC,MAAM,CAAC,EAAE,OAAO,GAAG,OAAO;IAQxC;;;;OAIG;IACI,eAAe,CAAC,MAAM,CAAC,EAAE,OAAO,GAAG,OAAO,GAAG,SAAS;IAK7D;;;;OAIG;IACI,SAAS,CAAC,MAAM,CAAC,EAAE,OAAO,GAAG,OAAO,GAAG,SAAS;IAQvD;;;OAGG;IACI,iBAAiB,IAAI,OAAO,GAAG,QAAQ;IAO9C;;;;;;;;OAQG;WACW,2BAA2B,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,MAAM,CAAC,EAAE,OAAO,GAAG,OAAO;IAKhH;;;;;;;;;;;;OAYG;WACW,mCAAmC,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,EAAE,EAAE,MAAM,EAAE,EAAE,EAAE,MAAM,EAAE,EAAE,EAAE,MAAM,EAAE,EAAE,EAAE,MAAM,EAAE,MAAM,CAAC,EAAE,KAAK,GAAG,KAAK;IAQpK;;;;;;;;;;;;OAYG;WACW,uDAAuD,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,EAAE,EAAE,MAAM,EAAE,EAAE,EAAE,MAAM,EAAE,EAAE,EAAE,MAAM,EAAE,EAAE,EAAE,MAAM,EAAE,GAAG,EAAE,MAAM,EAAE,GAAG,EAAE,MAAM,EAAE,GAAG,EAAE,MAAM,EAAE,GAAG,EAAE,MAAM,EAAE,MAAM,CAAC,EAAE,yBAAyB,GAAG,yBAAyB;IAcpR;;;OAGG;IACI,mBAAmB,CAAC,MAAM,CAAC,EAAE,OAAO,GAAG,OAAO;IAQrD;;;;OAIG;IACI,aAAa,CAAC,MAAM,CAAC,EAAE,OAAO,GAAG,OAAO,GAAG,SAAS;IAM3D;;OAEG;WACW,qBAAqB,CAAC,MAAM,EAAE,OAAO,EAAE,MAAM,EAAE,OAAO,EAAE,MAAM,EAAE,OAAO,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM;IAMtH;;;;;OAKG;WACW,yBAAyB,CAAC,MAAM,EAAE,OAAO,EAAE,MAAM,EAAE,OAAO,EAAE,MAAM,EAAE,OAAO,GAAG,OAAO;IAOnG,sFAAsF;IAC/E,8BAA8B,CAAC,MAAM,CAAC,EAAE,4BAA4B,GAAG,4BAA4B,GAAG,SAAS;IAGtH,6DAA6D;IACtD,mBAAmB;IAY1B,6EAA6E;WAC/D,sBAAsB,CAAC,WAAW,EAAE,OAAO,EAAE,iBAAiB,EAAE,MAAM,EAAE,WAAW,EAAE,OAAO,EAAE,MAAM,CAAC,EAAE,OAAO,GAAG,OAAO;IAqDtI;;OAEG;IACI,oBAAoB,CAAC,KAAK,EAAE,OAAO,GAAG,MAAM,GAAG,SAAS;CAShE"}
+{"version":3,"file":"Point4d.d.ts","sourceRoot":"","sources":["../../../src/geometry4d/Point4d.ts"],"names":[],"mappings":"AAKA;;GAEG;AACH,OAAO,EAAE,eAAe,EAAY,sBAAsB,EAAE,MAAM,aAAa,CAAC;AAChF,OAAO,EAAE,OAAO,EAAE,MAAM,uBAAuB,CAAC;AAChD,OAAO,EAAE,4BAA4B,EAAE,MAAM,4CAA4C,CAAC;AAC1F,OAAO,EAAE,yBAAyB,EAAE,MAAM,yCAAyC,CAAC;AAEpF,OAAO,EAAE,OAAO,EAAE,QAAQ,EAAE,MAAM,+BAA+B,CAAC;AAClE,OAAO,EAAE,KAAK,EAAE,MAAM,qBAAqB,CAAC;AAC5C,OAAO,EAAE,KAAK,EAAE,MAAM,EAAE,MAAM,wBAAwB,CAAC;AAEvD;;;GAGG;AACH,MAAM,MAAM,YAAY,GAAG,MAAM,EAAE,CAAC;AAgBpC;;;;;;;GAOG;AACH,qBAAa,OAAQ,SAAQ,OAAQ,YAAW,eAAe;IAC7D,6CAA6C;IACtC,IAAI,EAAE,YAAY,CAAC;IAC1B,kCAAkC;IAC3B,GAAG,CAAC,CAAC,GAAE,MAAU,EAAE,CAAC,GAAE,MAAU,EAAE,CAAC,GAAE,MAAU,EAAE,CAAC,GAAE,MAAU,GAAG,OAAO;IAO/E;;OAEG;IACI,YAAY,CAAC,KAAK,EAAE,MAAM,EAAE,KAAK,EAAE,MAAM;IAKhD,uBAAuB;IACvB,IAAW,CAAC,IACM,MAAM,CADe;IACvC,IAAW,CAAC,CAAC,GAAG,EAAE,MAAM,EAAyB;IACjD,uBAAuB;IACvB,IAAW,CAAC,IACM,MAAM,CADe;IACvC,IAAW,CAAC,CAAC,GAAG,EAAE,MAAM,EAAyB;IACjD,uBAAuB;IACvB,IAAW,CAAC,IACM,MAAM,CADe;IACvC,IAAW,CAAC,CAAC,GAAG,EAAE,MAAM,EAAyB;IACjD,qCAAqC;IACrC,IAAW,CAAC,IACM,MAAM,CADe;IACvC,IAAW,CAAC,CAAC,GAAG,EAAE,MAAM,EAAyB;IACjD,kCAAkC;IAClC,SAAS,aAAa,CAAC,GAAE,MAAU,EAAE,CAAC,GAAE,MAAU,EAAE,CAAC,GAAE,MAAU,EAAE,CAAC,GAAE,MAAU;IAQhF,8CAA8C;WAChC,MAAM,CAAC,CAAC,GAAE,MAAU,EAAE,CAAC,GAAE,MAAU,EAAE,CAAC,GAAE,MAAU,EAAE,CAAC,GAAE,MAAU,EAAE,MAAM,CAAC,EAAE,OAAO,GAAG,OAAO;IAG3G;;;;OAIG;WACW,eAAe,CAAC,MAAM,EAAE,sBAAsB,GAAG,OAAO,GAAG,SAAS;IAGlF,qCAAqC;IAC9B,OAAO,CAAC,KAAK,EAAE,OAAO,GAAG,OAAO;IAOvC,uBAAuB;IAChB,KAAK,CAAC,MAAM,CAAC,EAAE,OAAO,GAAG,OAAO;IAGvC,0DAA0D;IACnD,WAAW,CAAC,IAAI,CAAC,EAAE,YAAY;IAOtC,wEAAwE;WAC1D,QAAQ,CAAC,IAAI,CAAC,EAAE,YAAY,GAAG,OAAO;IAKpD,6EAA6E;IACtE,aAAa,CAAC,KAAK,EAAE,OAAO,EAAE,SAAS,GAAE,MAAqC,GAAG,OAAO;IAM/F;;;;;;OAMG;IACI,iBAAiB,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,SAAS,GAAE,MAAqC,GAAG,OAAO;IAO/H;;;OAGG;IACI,MAAM,IAAI,YAAY;IAG7B;;OAEG;IACI,YAAY,CAAC,KAAK,EAAE,OAAO,GAAG,MAAM;IAG3C;;OAEG;IACI,mBAAmB,CAAC,KAAK,EAAE,OAAO,GAAG,MAAM;IAGlD,2FAA2F;IACpF,cAAc,CAAC,KAAK,EAAE,OAAO,GAAG,MAAM,GAAG,SAAS;IAOzD,mGAAmG;IAC5F,qBAAqB,CAAC,KAAK,EAAE,OAAO,GAAG,MAAM,GAAG,SAAS;IAOhE;;OAEG;IACI,OAAO,CAAC,KAAK,EAAE,OAAO,GAAG,MAAM;IAGtC,oEAAoE;IAC7D,MAAM,IAAI,MAAM;IAGvB,+DAA+D;IACxD,aAAa,IAAI,MAAM;IAG9B,wHAAwH;IACjH,mBAAmB,IAAI,MAAM;IAIpC,wEAAwE;IACjE,KAAK,CAAC,KAAK,EAAE,OAAO,EAAE,MAAM,CAAC,EAAE,OAAO,GAAG,OAAO;IAGvD,sDAAsD;IAC/C,kBAAkB,CAAC,KAAK,EAAE,OAAO,EAAE,MAAM,CAAC,EAAE,QAAQ,GAAG,QAAQ;IAKtE,kFAAkF;IAC3E,yBAAyB,CAAC,KAAK,EAAE,OAAO,EAAE,MAAM,CAAC,EAAE,QAAQ,GAAG,QAAQ;IAI7E,uEAAuE;IAChE,IAAI,CAAC,KAAK,EAAE,OAAO,EAAE,MAAM,CAAC,EAAE,OAAO,GAAG,OAAO;IAGtD,8CAA8C;IAC9C,IAAW,YAAY,IAAI,OAAO,CAEjC;IACD,mDAAmD;WACrC,UAAU,IAAI,OAAO;IACnC;;;;OAIG;WACW,sBAAsB,CAAC,MAAM,EAAE,OAAO,EAAE,MAAM,EAAE,OAAO,EAAE,MAAM,CAAC,EAAE,OAAO;IAGvF;;;;OAIG;WACW,mBAAmB,CAAC,IAAI,EAAE,YAAY,EAAE,MAAM,GAAE,MAAU,EAAE,MAAM,CAAC,EAAE,OAAO,GAAG,OAAO,GAAG,SAAS;IAKhH;;;;OAIG;WACW,gBAAgB,CAAC,IAAI,EAAE,YAAY,EAAE,MAAM,GAAE,MAAU,EAAE,MAAM,CAAC,EAAE,OAAO,GAAG,OAAO,GAAG,SAAS;IAK7G;;;;;OAKG;WACW,oBAAoB,CAAC,IAAI,EAAE,YAAY,EAAE,MAAM,GAAE,MAAU,EAAE,MAAM,CAAC,EAAE,OAAO,GAAG,OAAO;IAGrG,0FAA0F;WAC5E,wBAAwB,CAAC,GAAG,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,GAAG,OAAO;IAGvE;;;;;;OAMG;WACW,eAAe,CAAC,KAAK,EAAE,KAAK,GAAG,MAAM,GAAG,OAAO,GAAG,MAAM,EAAE,GAAG,OAAO;IAsBlF,uCAAuC;IAChC,UAAU,CAAC,MAAM,EAAE,OAAO,EAAE,WAAW,EAAE,MAAM,EAAE,MAAM,CAAC,EAAE,OAAO,GAAG,OAAO;IAGlF;OACG;IACI,WAAW,CAAC,QAAQ,EAAE,MAAM,EAAE,MAAM,EAAE,OAAO,EAAE,MAAM,CAAC,EAAE,OAAO,GAAG,OAAO;IAIhF,6DAA6D;IACtD,WAAW,CAAC,OAAO,EAAE,OAAO,EAAE,OAAO,EAAE,MAAM,EAAE,OAAO,EAAE,OAAO,EAAE,OAAO,EAAE,MAAM,EAAE,MAAM,CAAC,EAAE,OAAO,GAAG,OAAO;IAGnH,iFAAiF;IAC1E,WAAW,CAAC,OAAO,EAAE,OAAO,EAAE,OAAO,EAAE,MAAM,EAAE,OAAO,EAAE,OAAO,EAAE,OAAO,EAAE,MAAM,EAAE,OAAO,EAAE,OAAO,EAAE,OAAO,EAAE,MAAM,EAAE,MAAM,CAAC,EAAE,OAAO,GAAG,OAAO;IAGtJ,6DAA6D;WAC/C,gBAAgB,CAAC,OAAO,EAAE,OAAO,EAAE,OAAO,EAAE,MAAM,EAAE,OAAO,EAAE,OAAO,EAAE,OAAO,EAAE,MAAM,EAAE,MAAM,CAAC,EAAE,OAAO,GAAG,OAAO;IAG/H,iFAAiF;WACnE,gBAAgB,CAAC,OAAO,EAAE,OAAO,EAAE,OAAO,EAAE,MAAM,EAAE,OAAO,EAAE,OAAO,EAAE,OAAO,EAAE,MAAM,EAAE,OAAO,EAAE,OAAO,EAAE,OAAO,EAAE,MAAM,EAAE,MAAM,CAAC,EAAE,OAAO,GAAG,OAAO;IAGlK,kFAAkF;IAC3E,mBAAmB,CAAC,OAAO,EAAE,OAAO,EAAE,OAAO,EAAE,OAAO,GAAG,MAAM;IAMtE,+DAA+D;IACxD,UAAU,CAAC,KAAK,EAAE,OAAO,GAAG,MAAM;IAGzC,wDAAwD;IACjD,cAAc,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,GAAG,MAAM;IAGzE,8FAA8F;IACvF,QAAQ,CAAC,KAAK,EAAE,OAAO,GAAG,MAAM;IAGvC,4EAA4E;IACrE,WAAW,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,GAAG,MAAM;IAG3D,oGAAoG;IAC7F,gBAAgB,CAAC,KAAK,EAAE,OAAO,GAAG,MAAM;IAG/C,mGAAmG;IAC5F,QAAQ,CAAC,MAAM,EAAE,QAAQ,GAAG,MAAM;IAGzC,4EAA4E;IACrE,WAAW,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,GAAG,MAAM;IAG3D;;OAEG;IACI,OAAO,IAAI,MAAM;IACxB;;MAEE;IACK,OAAO,IAAI,MAAM;IACxB;;MAEE;IACK,OAAO,IAAI,MAAM;IAExB,oBAAoB;WACN,KAAK,IAAI,OAAO;IAC9B,oBAAoB;WACN,KAAK,IAAI,OAAO;IAC9B,oBAAoB;WACN,KAAK,IAAI,OAAO;IAC9B,oBAAoB;WACN,KAAK,IAAI,OAAO;IAC9B,0EAA0E;IACnE,gBAAgB,CAAC,WAAW,EAAE,MAAM,EAAE,MAAM,CAAC,EAAE,OAAO,GAAG,OAAO,GAAG,SAAS;IAMnF;;;OAGG;IACI,mBAAmB,CAAC,UAAU,EAAE,OAAO,EAAE,MAAM,CAAC,EAAE,OAAO,GAAG,OAAO;IAG1E,uDAAuD;IACvC,iBAAiB,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,MAAM,CAAC,EAAE,OAAO,GAAG,OAAO;IAS7F,2CAA2C;IACpC,KAAK,CAAC,KAAK,EAAE,MAAM,EAAE,MAAM,CAAC,EAAE,OAAO,GAAG,OAAO;IAQtD,wCAAwC;IACjC,MAAM,CAAC,MAAM,CAAC,EAAE,OAAO,GAAG,OAAO;IAQxC;;;;OAIG;IACI,eAAe,CAAC,MAAM,CAAC,EAAE,OAAO,GAAG,OAAO,GAAG,SAAS;IAK7D;;;;OAIG;IACI,SAAS,CAAC,MAAM,CAAC,EAAE,OAAO,GAAG,OAAO,GAAG,SAAS;IAQvD;;;OAGG;IACI,iBAAiB,IAAI,OAAO,GAAG,QAAQ;IAO9C;;;;;;;;OAQG;WACW,2BAA2B,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,MAAM,CAAC,EAAE,OAAO,GAAG,OAAO;IAKhH;;;;;;;;;;;;OAYG;WACW,mCAAmC,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,EAAE,EAAE,MAAM,EAAE,EAAE,EAAE,MAAM,EAAE,EAAE,EAAE,MAAM,EAAE,EAAE,EAAE,MAAM,EAAE,MAAM,CAAC,EAAE,KAAK,GAAG,KAAK;IAQpK;;;;;;;;;;;;OAYG;WACW,uDAAuD,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,EAAE,EAAE,MAAM,EAAE,EAAE,EAAE,MAAM,EAAE,EAAE,EAAE,MAAM,EAAE,EAAE,EAAE,MAAM,EAAE,GAAG,EAAE,MAAM,EAAE,GAAG,EAAE,MAAM,EAAE,GAAG,EAAE,MAAM,EAAE,GAAG,EAAE,MAAM,EAAE,MAAM,CAAC,EAAE,yBAAyB,GAAG,yBAAyB;IAcpR;;;OAGG;IACI,mBAAmB,CAAC,MAAM,CAAC,EAAE,OAAO,GAAG,OAAO;IAQrD;;;;OAIG;IACI,aAAa,CAAC,MAAM,CAAC,EAAE,OAAO,GAAG,OAAO,GAAG,SAAS;IAM3D;;OAEG;WACW,qBAAqB,CAAC,MAAM,EAAE,OAAO,EAAE,MAAM,EAAE,OAAO,EAAE,MAAM,EAAE,OAAO,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM;IAMtH;;;;;OAKG;WACW,yBAAyB,CAAC,MAAM,EAAE,OAAO,EAAE,MAAM,EAAE,OAAO,EAAE,MAAM,EAAE,OAAO,GAAG,OAAO;IAOnG,sFAAsF;IAC/E,8BAA8B,CAAC,MAAM,CAAC,EAAE,4BAA4B,GAAG,4BAA4B,GAAG,SAAS;IAGtH,6DAA6D;IACtD,mBAAmB;IAY1B,6EAA6E;WAC/D,sBAAsB,CAAC,WAAW,EAAE,OAAO,EAAE,iBAAiB,EAAE,MAAM,EAAE,WAAW,EAAE,OAAO,EAAE,MAAM,CAAC,EAAE,OAAO,GAAG,OAAO;IAqDtI;;OAEG;IACI,oBAAoB,CAAC,KAAK,EAAE,OAAO,GAAG,MAAM,GAAG,SAAS;CAShE"}

package/lib/cjs/geometry4d/Point4d.js

@@ -114,11 +114,11 @@ class Point4d extends Plane3d_1.Plane3d {
         return result;
     }
     /** Near-equality test, using `Geometry.isSameCoordinate` on all 4 x,y,z,w */
-    isAlmostEqual(other) {
-        return Geometry_1.Geometry.isSameCoordinate(this.x, other.x)
-            && Geometry_1.Geometry.isSameCoordinate(this.y, other.y)
-            && Geometry_1.Geometry.isSameCoordinate(this.z, other.z)
-            && Geometry_1.Geometry.isSameCoordinate(this.w, other.w);
+    isAlmostEqual(other, tolerance = Geometry_1.Geometry.smallMetricDistance) {
+        return Geometry_1.Geometry.isSameCoordinate(this.x, other.x, tolerance)
+            && Geometry_1.Geometry.isSameCoordinate(this.y, other.y, tolerance)
+            && Geometry_1.Geometry.isSameCoordinate(this.z, other.z, tolerance)
+            && Geometry_1.Geometry.isSameCoordinate(this.w, other.w, tolerance);
     }
     /**
      * Test for same coordinate by direct x,y,z,w args
@@ -127,11 +127,11 @@ class Point4d extends Plane3d_1.Plane3d {
      * @param z z to test
      * @param w w to test
      */
-    isAlmostEqualXYZW(x, y, z, w) {
-        return Geometry_1.Geometry.isSameCoordinate(this.x, x)
-            && Geometry_1.Geometry.isSameCoordinate(this.y, y)
-            && Geometry_1.Geometry.isSameCoordinate(this.z, z)
-            && Geometry_1.Geometry.isSameCoordinate(this.w, w);
+    isAlmostEqualXYZW(x, y, z, w, tolerance = Geometry_1.Geometry.smallMetricDistance) {
+        return Geometry_1.Geometry.isSameCoordinate(this.x, x, tolerance)
+            && Geometry_1.Geometry.isSameCoordinate(this.y, y, tolerance)
+            && Geometry_1.Geometry.isSameCoordinate(this.z, z, tolerance)
+            && Geometry_1.Geometry.isSameCoordinate(this.w, w, tolerance);
     }
     /**
      * Convert an Angle to a JSON object.
@@ -152,8 +152,7 @@ class Point4d extends Plane3d_1.Plane3d {
     distanceSquaredXYZW(other) {
         return Geometry_1.Geometry.hypotenuseSquaredXYZW(other.xyzw[0] - this.xyzw[0], other.xyzw[1] - this.xyzw[1], other.xyzw[2] - this.xyzw[2], other.xyzw[3] - this.xyzw[3]);
     }
-    /** Return the distance between the instance and other after normalizing by weights
-     */
+    /** Return the xy distance between the instance and `other` after normalizing by weights */
     realDistanceXY(other) {
         const wA = this.w;
         const wB = other.w;
@@ -161,6 +160,14 @@ class Point4d extends Plane3d_1.Plane3d {
             return undefined;
         return Geometry_1.Geometry.hypotenuseXY(other.xyzw[0] / wB - this.xyzw[0] / wA, other.xyzw[1] / wB - this.xyzw[1] / wA);
     }
+    /** Return the xy squared distance between the instance and `other` after normalizing by weights */
+    realDistanceSquaredXY(other) {
+        const wA = this.w;
+        const wB = other.w;
+        if (Geometry_1.Geometry.isSmallMetricDistance(wA) || Geometry_1.Geometry.isSmallMetricDistance(wB))
+            return undefined;
+        return Geometry_1.Geometry.hypotenuseSquaredXY(other.xyzw[0] / wB - this.xyzw[0] / wA, other.xyzw[1] / wB - this.xyzw[1] / wA);
+    }
     /** Return the largest absolute distance between corresponding components
      * * x,y,z,w all participate without normalization.
      */
@@ -236,7 +243,7 @@ class Point4d extends Plane3d_1.Plane3d {
      * extract 4 consecutive numbers from a Float64Array into a Point4d.
      * @param data buffer of numbers
      * @param xIndex first index for x,y,z,w sequence. Assumed to be a valid index!
-     * @deprecated in 4.x. Use createFromPacked instead.
+     * @deprecated in 4.x - will not be removed until after 2026-06-13. Use createFromPacked instead.
      */
     static createFromPackedXYZW(data, xIndex = 0, result) {
         return Point4d.create(data[xIndex], data[xIndex + 1], data[xIndex + 2], data[xIndex + 3], result);
@@ -366,13 +373,17 @@ class Point4d extends Plane3d_1.Plane3d {
      * * If the xyz part of `this` are all zero, (a clone of) `spacePoint` is returned.
      */
     projectPointToPlane(spacePoint, result) {
-        const h = this.altitude(spacePoint);
+        return this.projectXYZToPlane(spacePoint.x, spacePoint.y, spacePoint.z, result);
+    }
+    /** Return the projection of (x,y,z) onto the plane. */
+    projectXYZToPlane(x, y, z, result) {
+        const h = this.altitudeXYZ(x, y, z);
         const nn = this.magnitudeSquaredXYZ();
         // this unusual tol is needed so that toPlane3dByOriginAndUnitNormal agrees with its original implementation
         const alpha = Geometry_1.Geometry.conditionalDivideCoordinate(-h, nn, Geometry_1.Geometry.largeFractionResult * Geometry_1.Geometry.largeFractionResult);
         if (alpha === undefined)
-            return spacePoint.clone(result);
-        return spacePoint.plusXYZ(alpha * this.x, alpha * this.y, alpha * this.z, result);
+            return Point3dVector3d_1.Point3d.create(x, y, z, result);
+        return Point3dVector3d_1.Point3d.create(x + alpha * this.x, y + alpha * this.y, z + alpha * this.z, result);
     }
     /** scale all components (including w!!) */
     scale(scale, result) {