@itwin/core-geometry 5.1.0-dev.1 → 5.1.0-dev.4

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/esm/curve/RegionMomentsXY.d.ts

@@ -21,41 +21,34 @@ export declare class RegionMomentsXY extends NullGeometryHandler {
     private _activeMomentData?;
     private _point0;
     private _point1;
-    /** Accumulate (independent) integrations over
-     * * origin to chord of the arc.
-     * * origin to the "cap" between the chord and arc.
+    /**
+     * Accumulate (independent) integrations over:
+     * * The area between the arc and the chord connecting its endpoints.
+     * * The triangle with vertices: origin, arc start, arc end.
      */
     handleArc3d(arc: Arc3d): void;
-    /** Accumulate integrals over the (triangular) areas from the origin to each line segment */
+    /** Accumulate integrals over the (triangular) areas from the origin to each line segment. */
     handleLineString3d(ls: LineString3d): void;
-    /** Accumulate integrals over the (triangular) area from the origin to this line segment */
+    /** Accumulate integrals over the (triangular) area from the origin to this line segment. */
     handleLineSegment3d(segment: LineSegment3d): void;
     /** Accumulate integrals from origin to all primitives in the chain. */
     handleLoop(loop: Loop): MomentData | undefined;
-    /**
-     * ASSUMPTIONS FOR ORIENTATION AND CONTAINMENT ISSUES
-     * * Largest area is outer
-     * * All others are interior (and not overlapping)
-     * Hence
-     * * Outer area sign must be positive -- negate all integrations as needed
-     * * Outer area signs must be positive -- negate all integrations as needed
-     * @param region
-     */
+    /** Accumulate integrals from origin to the components of the parity region. */
     handleParityRegion(region: ParityRegion): MomentData | undefined;
-    /** Accumulate (as simple addition) products over each component of the union region. */
+    /** Accumulate integrals from origin to the components of the union region. */
     handleUnionRegion(region: UnionRegion): MomentData | undefined;
     private _strokeOptions?;
     private getStrokeOptions;
-    /** Single curve primitive (not loop . . .).
-     * * stroke the curve
-     * * accumulate stroke array.
+    /**
+     * Handle a single curve primitive (not loop).
+     * * Stroke the curve and accumulate stroke array.
      */
     handleCurvePrimitive(cp: CurvePrimitive): void;
-    /** handle strongly typed  BSplineCurve3d  as generic curve primitive */
+    /** Handle strongly typed  BSplineCurve3d  as generic curve primitive. */
     handleBSplineCurve3d(g: BSplineCurve3d): void;
-    /** handle strongly typed  BSplineCurve3dH  as generic curve primitive */
+    /** Handle strongly typed  BSplineCurve3dH  as generic curve primitive. */
     handleBSplineCurve3dH(g: BSplineCurve3dH): void;
-    /** handle strongly typed  TransitionSpiral as generic curve primitive  */
+    /** Handle strongly typed  TransitionSpiral as generic curve primitive. */
     handleTransitionSpiral(g: TransitionSpiral3d): void;
 }
 //# sourceMappingURL=RegionMomentsXY.d.ts.map

package/lib/esm/curve/RegionMomentsXY.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"RegionMomentsXY.d.ts","sourceRoot":"","sources":["../../../src/curve/RegionMomentsXY.ts"],"names":[],"mappings":"AAKA;;GAEG;AAIH,OAAO,EAAE,cAAc,EAAE,MAAM,yBAAyB,CAAC;AACzD,OAAO,EAAE,eAAe,EAAE,MAAM,4BAA4B,CAAC;AAE7D,OAAO,EAAE,mBAAmB,EAAE,MAAM,+BAA+B,CAAC;AAEpE,OAAO,EAAE,UAAU,EAAE,MAAM,0BAA0B,CAAC;AACtD,OAAO,EAAE,KAAK,EAAE,MAAM,SAAS,CAAC;AAChC,OAAO,EAAE,cAAc,EAAE,MAAM,kBAAkB,CAAC;AAClD,OAAO,EAAE,aAAa,EAAE,MAAM,iBAAiB,CAAC;AAChD,OAAO,EAAE,YAAY,EAAE,MAAM,gBAAgB,CAAC;AAC9C,OAAO,EAAE,IAAI,EAAE,MAAM,QAAQ,CAAC;AAC9B,OAAO,EAAE,YAAY,EAAE,MAAM,gBAAgB,CAAC;AAE9C,OAAO,EAAE,kBAAkB,EAAE,MAAM,6BAA6B,CAAC;AACjE,OAAO,EAAE,WAAW,EAAE,MAAM,eAAe,CAAC;AAE5C;;;GAGG;AACH,qBAAa,eAAgB,SAAQ,mBAAmB;IACtD,OAAO,CAAC,iBAAiB,CAAC,CAAa;IACvC,OAAO,CAAC,OAAO,CAAoB;IACnC,OAAO,CAAC,OAAO,CAAoB;IAEnC;;;OAGG;IACa,WAAW,CAAC,GAAG,EAAE,KAAK,GAAG,IAAI;IAgC7C,4FAA4F;IAC5E,kBAAkB,CAAC,EAAE,EAAE,YAAY,GAAG,IAAI;IAI1D,2FAA2F;IAC3E,mBAAmB,CAAC,OAAO,EAAE,aAAa,GAAG,IAAI;IAMjE,uEAAuE;IACvD,UAAU,CAAC,IAAI,EAAE,IAAI,GAAG,UAAU,GAAG,SAAS;IAQ9D;;;;;;;;OAQG;IACa,kBAAkB,CAAC,MAAM,EAAE,YAAY,GAAG,UAAU,GAAG,SAAS;IAgChF,wFAAwF;IACxE,iBAAiB,CAAC,MAAM,EAAE,WAAW,GAAG,UAAU,GAAG,SAAS;IAY9E,OAAO,CAAC,cAAc,CAAC,CAAgB;IACvC,OAAO,CAAC,gBAAgB;IASxB;;;OAGG;IACI,oBAAoB,CAAC,EAAE,EAAE,cAAc;IAM9C,wEAAwE;IACxD,oBAAoB,CAAC,CAAC,EAAE,cAAc;IACtD,yEAAyE;IACzD,qBAAqB,CAAC,CAAC,EAAE,eAAe;IACxD,0EAA0E;IAC1D,sBAAsB,CAAC,CAAC,EAAE,kBAAkB;CAE7D"}
+{"version":3,"file":"RegionMomentsXY.d.ts","sourceRoot":"","sources":["../../../src/curve/RegionMomentsXY.ts"],"names":[],"mappings":"AAKA;;GAEG;AAIH,OAAO,EAAE,cAAc,EAAE,MAAM,yBAAyB,CAAC;AACzD,OAAO,EAAE,eAAe,EAAE,MAAM,4BAA4B,CAAC;AAE7D,OAAO,EAAE,mBAAmB,EAAE,MAAM,+BAA+B,CAAC;AAEpE,OAAO,EAAE,UAAU,EAAE,MAAM,0BAA0B,CAAC;AACtD,OAAO,EAAE,KAAK,EAAE,MAAM,SAAS,CAAC;AAChC,OAAO,EAAE,cAAc,EAAE,MAAM,kBAAkB,CAAC;AAClD,OAAO,EAAE,aAAa,EAAE,MAAM,iBAAiB,CAAC;AAChD,OAAO,EAAE,YAAY,EAAE,MAAM,gBAAgB,CAAC;AAC9C,OAAO,EAAE,IAAI,EAAE,MAAM,QAAQ,CAAC;AAC9B,OAAO,EAAE,YAAY,EAAE,MAAM,gBAAgB,CAAC;AAE9C,OAAO,EAAE,kBAAkB,EAAE,MAAM,6BAA6B,CAAC;AACjE,OAAO,EAAE,WAAW,EAAE,MAAM,eAAe,CAAC;AAE5C;;;GAGG;AACH,qBAAa,eAAgB,SAAQ,mBAAmB;IACtD,OAAO,CAAC,iBAAiB,CAAC,CAAa;IACvC,OAAO,CAAC,OAAO,CAAoB;IACnC,OAAO,CAAC,OAAO,CAAoB;IACnC;;;;OAIG;IACa,WAAW,CAAC,GAAG,EAAE,KAAK,GAAG,IAAI;IAgC7C,6FAA6F;IAC7E,kBAAkB,CAAC,EAAE,EAAE,YAAY,GAAG,IAAI;IAI1D,4FAA4F;IAC5E,mBAAmB,CAAC,OAAO,EAAE,aAAa,GAAG,IAAI;IAMjE,uEAAuE;IACvD,UAAU,CAAC,IAAI,EAAE,IAAI,GAAG,UAAU,GAAG,SAAS;IAQ9D,+EAA+E;IAC/D,kBAAkB,CAAC,MAAM,EAAE,YAAY,GAAG,UAAU,GAAG,SAAS;IA+BhF,8EAA8E;IAC9D,iBAAiB,CAAC,MAAM,EAAE,WAAW,GAAG,UAAU,GAAG,SAAS;IAW9E,OAAO,CAAC,cAAc,CAAC,CAAgB;IACvC,OAAO,CAAC,gBAAgB;IASxB;;;OAGG;IACI,oBAAoB,CAAC,EAAE,EAAE,cAAc;IAM9C,yEAAyE;IACzD,oBAAoB,CAAC,CAAC,EAAE,cAAc;IAGtD,0EAA0E;IAC1D,qBAAqB,CAAC,CAAC,EAAE,eAAe;IAGxD,0EAA0E;IAC1D,sBAAsB,CAAC,CAAC,EAAE,kBAAkB;CAG7D"}

package/lib/esm/curve/RegionMomentsXY.js

@@ -17,15 +17,15 @@ export class RegionMomentsXY extends NullGeometryHandler {
     _activeMomentData;
     _point0 = Point3d.create();
     _point1 = Point3d.create();
-    /** Accumulate (independent) integrations over
-     * * origin to chord of the arc.
-     * * origin to the "cap" between the chord and arc.
+    /**
+     * Accumulate (independent) integrations over:
+     * * The area between the arc and the chord connecting its endpoints.
+     * * The triangle with vertices: origin, arc start, arc end.
      */
     handleArc3d(arc) {
         const momentData = this._activeMomentData;
         const sweepRadians = arc.sweep.sweepRadians;
         const alphaRadians = sweepRadians * 0.5;
-        // from https://apps.dtic.mil/dtic/tr/fulltext/u2/274936.pdf page 71  for radius = 1
         let s = Math.sin(alphaRadians);
         let c = Math.cos(alphaRadians);
         let s1 = Math.sin(sweepRadians);
@@ -52,12 +52,12 @@ export class RegionMomentsXY extends NullGeometryHandler {
         const pointC = arc.fractionToPoint(1.0);
         momentData.accumulateTriangleMomentsXY(undefined, pointB, pointC);
     }
-    /** Accumulate integrals over the (triangular) areas from the origin to each line segment */
+    /** Accumulate integrals over the (triangular) areas from the origin to each line segment. */
     handleLineString3d(ls) {
         const momentData = this._activeMomentData;
         momentData.accumulateTriangleToLineStringMomentsXY(undefined, ls.packedPoints);
     }
-    /** Accumulate integrals over the (triangular) area from the origin to this line segment */
+    /** Accumulate integrals over the (triangular) area from the origin to this line segment. */
     handleLineSegment3d(segment) {
         const momentData = this._activeMomentData;
         segment.startPoint(this._point0);
@@ -73,15 +73,7 @@ export class RegionMomentsXY extends NullGeometryHandler {
         this._activeMomentData = undefined;
         return momentData;
     }
-    /**
-     * ASSUMPTIONS FOR ORIENTATION AND CONTAINMENT ISSUES
-     * * Largest area is outer
-     * * All others are interior (and not overlapping)
-     * Hence
-     * * Outer area sign must be positive -- negate all integrations as needed
-     * * Outer area signs must be positive -- negate all integrations as needed
-     * @param region
-     */
+    /** Accumulate integrals from origin to the components of the parity region. */
     handleParityRegion(region) {
         const allChildMoments = [];
         let maxAbsArea = 0.0;
@@ -113,7 +105,7 @@ export class RegionMomentsXY extends NullGeometryHandler {
         }
         return undefined;
     }
-    /** Accumulate (as simple addition) products over each component of the union region. */
+    /** Accumulate integrals from origin to the components of the union region. */
     handleUnionRegion(region) {
         const summedMoments = MomentData.create();
         for (const child of region.children) {
@@ -135,9 +127,9 @@ export class RegionMomentsXY extends NullGeometryHandler {
         this._strokeOptions = options;
         return options;
     }
-    /** Single curve primitive (not loop . . .).
-     * * stroke the curve
-     * * accumulate stroke array.
+    /**
+     * Handle a single curve primitive (not loop).
+     * * Stroke the curve and accumulate stroke array.
      */
     handleCurvePrimitive(cp) {
         const strokes = LineString3d.create();
@@ -145,11 +137,17 @@ export class RegionMomentsXY extends NullGeometryHandler {
         cp.emitStrokes(strokes, options);
         this.handleLineString3d(strokes);
     }
-    /** handle strongly typed  BSplineCurve3d  as generic curve primitive */
-    handleBSplineCurve3d(g) { return this.handleCurvePrimitive(g); }
-    /** handle strongly typed  BSplineCurve3dH  as generic curve primitive */
-    handleBSplineCurve3dH(g) { return this.handleCurvePrimitive(g); }
-    /** handle strongly typed  TransitionSpiral as generic curve primitive  */
-    handleTransitionSpiral(g) { return this.handleCurvePrimitive(g); }
+    /** Handle strongly typed  BSplineCurve3d  as generic curve primitive. */
+    handleBSplineCurve3d(g) {
+        return this.handleCurvePrimitive(g);
+    }
+    /** Handle strongly typed  BSplineCurve3dH  as generic curve primitive. */
+    handleBSplineCurve3dH(g) {
+        return this.handleCurvePrimitive(g);
+    }
+    /** Handle strongly typed  TransitionSpiral as generic curve primitive. */
+    handleTransitionSpiral(g) {
+        return this.handleCurvePrimitive(g);
+    }
 }
 //# sourceMappingURL=RegionMomentsXY.js.map

package/lib/esm/curve/RegionMomentsXY.js.map

@@ -1 +1 @@
-{"version":3,"file":"RegionMomentsXY.js","sourceRoot":"","sources":["../../../src/curve/RegionMomentsXY.ts"],"names":[],"mappings":"AAAA;;;+FAG+F;AAU/F,OAAO,EAAE,KAAK,EAAE,MAAM,qBAAqB,CAAC;AAC5C,OAAO,EAAE,mBAAmB,EAAE,MAAM,+BAA+B,CAAC;AACpE,OAAO,EAAE,OAAO,EAAE,MAAM,+BAA+B,CAAC;AACxD,OAAO,EAAE,UAAU,EAAE,MAAM,0BAA0B,CAAC;AAItD,OAAO,EAAE,YAAY,EAAE,MAAM,gBAAgB,CAAC;AAC9C,OAAO,EAAE,IAAI,EAAE,MAAM,QAAQ,CAAC;AAE9B,OAAO,EAAE,aAAa,EAAE,MAAM,iBAAiB,CAAC;AAIhD;;;GAGG;AACH,MAAM,OAAO,eAAgB,SAAQ,mBAAmB;IAC9C,iBAAiB,CAAc;IAC/B,OAAO,GAAG,OAAO,CAAC,MAAM,EAAE,CAAC;IAC3B,OAAO,GAAG,OAAO,CAAC,MAAM,EAAE,CAAC;IAEnC;;;OAGG;IACa,WAAW,CAAC,GAAU;QACpC,MAAM,UAAU,GAAG,IAAI,CAAC,iBAAkB,CAAC;QAC3C,MAAM,YAAY,GAAG,GAAG,CAAC,KAAK,CAAC,YAAY,CAAC;QAC5C,MAAM,YAAY,GAAG,YAAY,GAAG,GAAG,CAAC;QACxC,oFAAoF;QACpF,IAAI,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,YAAY,CAAC,CAAC;QAC/B,IAAI,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,YAAY,CAAC,CAAC;QAC/B,IAAI,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,YAAY,CAAC,CAAC;QAChC,IAAI,KAAK,CAAC,mBAAmB,CAAC,YAAY,CAAC,EAAE,CAAC;YAC5C,CAAC,GAAG,GAAG,CAAC;YACR,CAAC,GAAG,CAAC,GAAG,CAAC;YACT,EAAE,GAAG,GAAG,CAAC;QACX,CAAC;QACD,MAAM,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,YAAY,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC;QACrD,MAAM,EAAE,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;QACrB,MAAM,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC;QACnB,MAAM,IAAI,GAAG,GAAG,GAAG,CAAC,YAAY,GAAG,EAAE,CAAC,CAAC;QACvC,MAAM,SAAS,GAAG,IAAI,GAAG,IAAI,GAAG,CAAC,GAAG,GAAG,CAAC,GAAG,GAAG,CAAC,CAAC;QAChD,MAAM,UAAU,GAAG,IAAI,GAAG,IAAI,GAAG,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC;QAC3C,MAAM,SAAS,GAAG,UAAU,GAAG,GAAG,GAAG,EAAE,GAAG,CAAC,GAAG,GAAG,IAAI,CAAC,CAAC;QACvD,MAAM,SAAS,GAAG,SAAS,CAAC;QAC5B,MAAM,SAAS,GAAG,SAAS,CAAC;QAC5B,MAAM,gBAAgB,GAAG,GAAG,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC,YAAY,GAAG,EAAE,CAAC,CAAC,CAAC;QACvE,MAAM,UAAU,GAAG,GAAG,CAAC,KAAK,CAAC,iBAAiB,CAAC,GAAG,CAAC,CAAC;QACpD,MAAM,YAAY,GAAG,GAAG,CAAC,qBAAqB,CAAC,UAAU,CAAC,CAAC;QAC3D,MAAM,QAAQ,GAAG,YAAY,CAAC,MAAM,CAAC,UAAU,CAAC,YAAY,CAAC,OAAO,EAAE,gBAAgB,CAAC,CAAC;QACxF,UAAU,CAAC,qCAAqC,CAAC,SAAS,EAAE,GAAG,EAAE,SAAS,EAAE,IAAI,EAAE,QAAQ,EAAE,YAAY,CAAC,OAAO,EAAE,YAAY,CAAC,OAAO,CAAC,CAAC;QACxI,MAAM,MAAM,GAAG,GAAG,CAAC,eAAe,CAAC,GAAG,CAAC,CAAC;QACxC,MAAM,MAAM,GAAG,GAAG,CAAC,eAAe,CAAC,GAAG,CAAC,CAAC;QACxC,UAAU,CAAC,2BAA2B,CAAC,SAAS,EAAE,MAAM,EAAE,MAAM,CAAC,CAAC;IACpE,CAAC;IAED,4FAA4F;IAC5E,kBAAkB,CAAC,EAAgB;QACjD,MAAM,UAAU,GAAG,IAAI,CAAC,iBAAkB,CAAC;QAC3C,UAAU,CAAC,uCAAuC,CAAC,SAAS,EAAE,EAAE,CAAC,YAAY,CAAC,CAAC;IACjF,CAAC;IACD,2FAA2F;IAC3E,mBAAmB,CAAC,OAAsB;QACxD,MAAM,UAAU,GAAG,IAAI,CAAC,iBAAkB,CAAC;QAC3C,OAAO,CAAC,UAAU,CAAC,IAAI,CAAC,OAAO,CAAC,CAAC;QACjC,OAAO,CAAC,QAAQ,CAAC,IAAI,CAAC,OAAO,CAAC,CAAC;QAC/B,UAAU,CAAC,2BAA2B,CAAC,SAAS,EAAE,IAAI,CAAC,OAAO,EAAE,IAAI,CAAC,OAAO,CAAC,CAAC;IAChF,CAAC;IACD,uEAAuE;IACvD,UAAU,CAAC,IAAU;QACnC,MAAM,UAAU,GAAG,IAAI,CAAC,iBAAiB,GAAG,UAAU,CAAC,MAAM,EAAE,CAAC;QAChE,UAAU,CAAC,UAAU,GAAG,KAAK,CAAC;QAC9B,KAAK,MAAM,KAAK,IAAI,IAAI,CAAC,QAAQ;YAC/B,KAAK,CAAC,yBAAyB,CAAC,IAAI,CAAC,CAAC;QACxC,IAAI,CAAC,iBAAiB,GAAG,SAAS,CAAC;QACnC,OAAO,UAAU,CAAC;IACpB,CAAC;IACD;;;;;;;;OAQG;IACa,kBAAkB,CAAC,MAAoB;QACrD,MAAM,eAAe,GAAiB,EAAE,CAAC;QACzC,IAAI,UAAU,GAAG,GAAG,CAAC;QACrB,IAAI,mBAA2C,CAAC;QAChD,KAAK,MAAM,KAAK,IAAI,MAAM,CAAC,QAAQ,EAAE,CAAC;YACpC,IAAI,KAAK,YAAY,IAAI,EAAE,CAAC;gBAC1B,MAAM,YAAY,GAAG,IAAI,CAAC,UAAU,CAAC,KAAK,CAAC,CAAC;gBAC5C,IAAI,YAAY,EAAE,CAAC;oBACjB,eAAe,CAAC,IAAI,CAAC,YAAY,CAAC,CAAC;oBACnC,MAAM,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,YAAY,CAAC,WAAW,CAAC,CAAC;oBAC7C,IAAI,CAAC,GAAG,UAAU,EAAE,CAAC;wBACnB,UAAU,GAAG,CAAC,CAAC;wBACf,mBAAmB,GAAG,YAAY,CAAC;oBACrC,CAAC;gBAEH,CAAC;YACH,CAAC;QACH,CAAC;QACD,IAAI,mBAAmB,EAAE,CAAC;YACxB,MAAM,aAAa,GAAG,UAAU,CAAC,MAAM,EAAE,CAAC;YAC1C,MAAM,KAAK,GAAG,mBAAmB,CAAC,UAAU,CAAC,GAAG,CAAC,CAAC;YAClD,aAAa,CAAC,kBAAkB,CAAC,mBAAmB,EAAE,KAAK,CAAC,CAAC;YAC7D,KAAK,MAAM,YAAY,IAAI,eAAe,EAAE,CAAC;gBAC3C,IAAI,YAAY,KAAK,mBAAmB,EAAE,CAAC;oBACzC,MAAM,KAAK,GAAG,YAAY,CAAC,UAAU,CAAC,CAAC,GAAG,CAAC,CAAC;oBAC5C,aAAa,CAAC,kBAAkB,CAAC,YAAY,EAAE,KAAK,CAAC,CAAC;gBACxD,CAAC;YACH,CAAC;YACD,OAAO,aAAa,CAAC;QACvB,CAAC;QACD,OAAO,SAAS,CAAC;IACnB,CAAC;IACD,wFAAwF;IACxE,iBAAiB,CAAC,MAAmB;QACnD,MAAM,aAAa,GAAG,UAAU,CAAC,MAAM,EAAE,CAAC;QAC1C,KAAK,MAAM,KAAK,IAAI,MAAM,CAAC,QAAQ,EAAE,CAAC;YACpC,MAAM,YAAY,GAAG,KAAK,CAAC,yBAAyB,CAAC,IAAI,CAAC,CAAC;YAC3D,IAAI,YAAY,EAAE,CAAC;gBACjB,MAAM,KAAK,GAAG,YAAY,CAAC,UAAU,CAAC,GAAG,CAAC,CAAC;gBAC3C,aAAa,CAAC,kBAAkB,CAAC,YAAY,EAAE,KAAK,CAAC,CAAC;YACxD,CAAC;QACH,CAAC;QACD,OAAO,aAAa,CAAC;IACvB,CAAC;IAEO,cAAc,CAAiB;IAC/B,gBAAgB;QACtB,IAAI,IAAI,CAAC,cAAc;YACrB,OAAO,IAAI,CAAC,cAAc,CAAC;QAC7B,MAAM,OAAO,GAAG,aAAa,CAAC,eAAe,EAAE,CAAC;QAChD,gEAAgE;QAChE,OAAO,CAAC,QAAQ,GAAG,KAAK,CAAC,aAAa,CAAC,GAAG,CAAC,CAAC;QAC5C,IAAI,CAAC,cAAc,GAAG,OAAO,CAAC;QAC9B,OAAO,OAAO,CAAC;IACjB,CAAC;IACD;;;OAGG;IACI,oBAAoB,CAAC,EAAkB;QAC5C,MAAM,OAAO,GAAG,YAAY,CAAC,MAAM,EAAE,CAAC;QACtC,MAAM,OAAO,GAAG,IAAI,CAAC,gBAAgB,EAAE,CAAC;QACxC,EAAE,CAAC,WAAW,CAAC,OAAO,EAAE,OAAO,CAAC,CAAC;QACjC,IAAI,CAAC,kBAAkB,CAAC,OAAO,CAAC,CAAC;IACnC,CAAC;IACD,wEAAwE;IACxD,oBAAoB,CAAC,CAAiB,IAAI,OAAO,IAAI,CAAC,oBAAoB,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IAChG,yEAAyE;IACzD,qBAAqB,CAAC,CAAkB,IAAI,OAAO,IAAI,CAAC,oBAAoB,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IAClG,0EAA0E;IAC1D,sBAAsB,CAAC,CAAqB,IAAI,OAAO,IAAI,CAAC,oBAAoB,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;CAEvG","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 Curve\n */\n\n// import { Geometry, Angle, AngleSweep } from \"../Geometry\";\n\nimport { BSplineCurve3d } from \"../bspline/BSplineCurve\";\nimport { BSplineCurve3dH } from \"../bspline/BSplineCurve3dH\";\nimport { Angle } from \"../geometry3d/Angle\";\nimport { NullGeometryHandler } from \"../geometry3d/GeometryHandler\";\nimport { Point3d } from \"../geometry3d/Point3dVector3d\";\nimport { MomentData } from \"../geometry4d/MomentData\";\nimport { Arc3d } from \"./Arc3d\";\nimport { CurvePrimitive } from \"./CurvePrimitive\";\nimport { LineSegment3d } from \"./LineSegment3d\";\nimport { LineString3d } from \"./LineString3d\";\nimport { Loop } from \"./Loop\";\nimport { ParityRegion } from \"./ParityRegion\";\nimport { StrokeOptions } from \"./StrokeOptions\";\nimport { TransitionSpiral3d } from \"./spiral/TransitionSpiral3d\";\nimport { UnionRegion } from \"./UnionRegion\";\n\n/**\n * Implementation class for computing XY area moments.\n * @internal\n */\nexport class RegionMomentsXY extends NullGeometryHandler {\n  private _activeMomentData?: MomentData;\n  private _point0 = Point3d.create();\n  private _point1 = Point3d.create();\n\n  /** Accumulate (independent) integrations over\n   * * origin to chord of the arc.\n   * * origin to the \"cap\" between the chord and arc.\n   */\n  public override handleArc3d(arc: Arc3d): void {\n    const momentData = this._activeMomentData!;\n    const sweepRadians = arc.sweep.sweepRadians;\n    const alphaRadians = sweepRadians * 0.5;\n    // from https://apps.dtic.mil/dtic/tr/fulltext/u2/274936.pdf page 71  for radius = 1\n    let s = Math.sin(alphaRadians);\n    let c = Math.cos(alphaRadians);\n    let s1 = Math.sin(sweepRadians);\n    if (Angle.isFullCircleRadians(sweepRadians)) {\n      s = 0.0;\n      c = -1.0;\n      s1 = 0.0;\n    }\n    const q = 2 * s * s * s * c / (alphaRadians - s * c);\n    const s3 = s * s * s;\n    const s6 = s3 * s3;\n    const area = 0.5 * (sweepRadians - s1);\n    const inertiaXX = 0.25 * area * (1.0 - q / 3.0);\n    const inertiaYY1 = 0.25 * area * (1.0 + q);\n    const inertiaYY = inertiaYY1 - 4.0 * s6 / (9.0 * area);\n    const productXX = inertiaYY;\n    const productYY = inertiaXX;\n    const centerToCentroid = 4.0 * s * s * s / (3.0 * (sweepRadians - s1));\n    const midRadians = arc.sweep.fractionToRadians(0.5);\n    const centralPlane = arc.radiansToRotatedBasis(midRadians);\n    const centroid = centralPlane.origin.plusScaled(centralPlane.vectorU, centerToCentroid);\n    momentData.accumulateXYProductsInCentroidalFrame(productXX, 0.0, productYY, area, centroid, centralPlane.vectorU, centralPlane.vectorV);\n    const pointB = arc.fractionToPoint(0.0);\n    const pointC = arc.fractionToPoint(1.0);\n    momentData.accumulateTriangleMomentsXY(undefined, pointB, pointC);\n  }\n\n  /** Accumulate integrals over the (triangular) areas from the origin to each line segment */\n  public override handleLineString3d(ls: LineString3d): void {\n    const momentData = this._activeMomentData!;\n    momentData.accumulateTriangleToLineStringMomentsXY(undefined, ls.packedPoints);\n  }\n  /** Accumulate integrals over the (triangular) area from the origin to this line segment */\n  public override handleLineSegment3d(segment: LineSegment3d): void {\n    const momentData = this._activeMomentData!;\n    segment.startPoint(this._point0);\n    segment.endPoint(this._point1);\n    momentData.accumulateTriangleMomentsXY(undefined, this._point0, this._point1);\n  }\n  /** Accumulate integrals from origin to all primitives in the chain. */\n  public override handleLoop(loop: Loop): MomentData | undefined {\n    const momentData = this._activeMomentData = MomentData.create();\n    momentData.needOrigin = false;\n    for (const child of loop.children)\n      child.dispatchToGeometryHandler(this);\n    this._activeMomentData = undefined;\n    return momentData;\n  }\n  /**\n   * ASSUMPTIONS FOR ORIENTATION AND CONTAINMENT ISSUES\n   * * Largest area is outer\n   * * All others are interior (and not overlapping)\n   * Hence\n   * * Outer area sign must be positive -- negate all integrations as needed\n   * * Outer area signs must be positive -- negate all integrations as needed\n   * @param region\n   */\n  public override handleParityRegion(region: ParityRegion): MomentData | undefined {\n    const allChildMoments: MomentData[] = [];\n    let maxAbsArea = 0.0;\n    let largestChildMoments: MomentData | undefined;\n    for (const child of region.children) {\n      if (child instanceof Loop) {\n        const childMoments = this.handleLoop(child);\n        if (childMoments) {\n          allChildMoments.push(childMoments);\n          const q = Math.abs(childMoments.quantitySum);\n          if (q > maxAbsArea) {\n            maxAbsArea = q;\n            largestChildMoments = childMoments;\n          }\n\n        }\n      }\n    }\n    if (largestChildMoments) {\n      const summedMoments = MomentData.create();\n      const sign0 = largestChildMoments.signFactor(1.0);\n      summedMoments.accumulateProducts(largestChildMoments, sign0);\n      for (const childMoments of allChildMoments) {\n        if (childMoments !== largestChildMoments) {\n          const sign1 = childMoments.signFactor(-1.0);\n          summedMoments.accumulateProducts(childMoments, sign1);\n        }\n      }\n      return summedMoments;\n    }\n    return undefined;\n  }\n  /** Accumulate (as simple addition) products over each component of the union region. */\n  public override handleUnionRegion(region: UnionRegion): MomentData | undefined {\n    const summedMoments = MomentData.create();\n    for (const child of region.children) {\n      const childMoments = child.dispatchToGeometryHandler(this);\n      if (childMoments) {\n        const sign0 = childMoments.signFactor(1.0);\n        summedMoments.accumulateProducts(childMoments, sign0);\n      }\n    }\n    return summedMoments;\n  }\n\n  private _strokeOptions?: StrokeOptions;\n  private getStrokeOptions(): StrokeOptions {\n    if (this._strokeOptions)\n      return this._strokeOptions;\n    const options = StrokeOptions.createForCurves();\n    // this is unusually fine for stroking, but appropriate for sum.\n    options.angleTol = Angle.createDegrees(5.0);\n    this._strokeOptions = options;\n    return options;\n  }\n  /** Single curve primitive (not loop . . .).\n   * * stroke the curve\n   * * accumulate stroke array.\n   */\n  public handleCurvePrimitive(cp: CurvePrimitive) {\n    const strokes = LineString3d.create();\n    const options = this.getStrokeOptions();\n    cp.emitStrokes(strokes, options);\n    this.handleLineString3d(strokes);\n  }\n  /** handle strongly typed  BSplineCurve3d  as generic curve primitive */\n  public override handleBSplineCurve3d(g: BSplineCurve3d) { return this.handleCurvePrimitive(g); }\n  /** handle strongly typed  BSplineCurve3dH  as generic curve primitive */\n  public override handleBSplineCurve3dH(g: BSplineCurve3dH) { return this.handleCurvePrimitive(g); }\n  /** handle strongly typed  TransitionSpiral as generic curve primitive  */\n  public override handleTransitionSpiral(g: TransitionSpiral3d) { return this.handleCurvePrimitive(g); }\n\n}\n"]}
+{"version":3,"file":"RegionMomentsXY.js","sourceRoot":"","sources":["../../../src/curve/RegionMomentsXY.ts"],"names":[],"mappings":"AAAA;;;+FAG+F;AAU/F,OAAO,EAAE,KAAK,EAAE,MAAM,qBAAqB,CAAC;AAC5C,OAAO,EAAE,mBAAmB,EAAE,MAAM,+BAA+B,CAAC;AACpE,OAAO,EAAE,OAAO,EAAE,MAAM,+BAA+B,CAAC;AACxD,OAAO,EAAE,UAAU,EAAE,MAAM,0BAA0B,CAAC;AAItD,OAAO,EAAE,YAAY,EAAE,MAAM,gBAAgB,CAAC;AAC9C,OAAO,EAAE,IAAI,EAAE,MAAM,QAAQ,CAAC;AAE9B,OAAO,EAAE,aAAa,EAAE,MAAM,iBAAiB,CAAC;AAIhD;;;GAGG;AACH,MAAM,OAAO,eAAgB,SAAQ,mBAAmB;IAC9C,iBAAiB,CAAc;IAC/B,OAAO,GAAG,OAAO,CAAC,MAAM,EAAE,CAAC;IAC3B,OAAO,GAAG,OAAO,CAAC,MAAM,EAAE,CAAC;IACnC;;;;OAIG;IACa,WAAW,CAAC,GAAU;QACpC,MAAM,UAAU,GAAG,IAAI,CAAC,iBAAkB,CAAC;QAC3C,MAAM,YAAY,GAAG,GAAG,CAAC,KAAK,CAAC,YAAY,CAAC;QAC5C,MAAM,YAAY,GAAG,YAAY,GAAG,GAAG,CAAC;QACxC,IAAI,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,YAAY,CAAC,CAAC;QAC/B,IAAI,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,YAAY,CAAC,CAAC;QAC/B,IAAI,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,YAAY,CAAC,CAAC;QAChC,IAAI,KAAK,CAAC,mBAAmB,CAAC,YAAY,CAAC,EAAE,CAAC;YAC5C,CAAC,GAAG,GAAG,CAAC;YACR,CAAC,GAAG,CAAC,GAAG,CAAC;YACT,EAAE,GAAG,GAAG,CAAC;QACX,CAAC;QACD,MAAM,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,YAAY,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC;QACrD,MAAM,EAAE,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;QACrB,MAAM,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC;QACnB,MAAM,IAAI,GAAG,GAAG,GAAG,CAAC,YAAY,GAAG,EAAE,CAAC,CAAC;QACvC,MAAM,SAAS,GAAG,IAAI,GAAG,IAAI,GAAG,CAAC,GAAG,GAAG,CAAC,GAAG,GAAG,CAAC,CAAC;QAChD,MAAM,UAAU,GAAG,IAAI,GAAG,IAAI,GAAG,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC;QAC3C,MAAM,SAAS,GAAG,UAAU,GAAG,GAAG,GAAG,EAAE,GAAG,CAAC,GAAG,GAAG,IAAI,CAAC,CAAC;QACvD,MAAM,SAAS,GAAG,SAAS,CAAC;QAC5B,MAAM,SAAS,GAAG,SAAS,CAAC;QAC5B,MAAM,gBAAgB,GAAG,GAAG,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC,YAAY,GAAG,EAAE,CAAC,CAAC,CAAC;QACvE,MAAM,UAAU,GAAG,GAAG,CAAC,KAAK,CAAC,iBAAiB,CAAC,GAAG,CAAC,CAAC;QACpD,MAAM,YAAY,GAAG,GAAG,CAAC,qBAAqB,CAAC,UAAU,CAAC,CAAC;QAC3D,MAAM,QAAQ,GAAG,YAAY,CAAC,MAAM,CAAC,UAAU,CAAC,YAAY,CAAC,OAAO,EAAE,gBAAgB,CAAC,CAAC;QACxF,UAAU,CAAC,qCAAqC,CAC9C,SAAS,EAAE,GAAG,EAAE,SAAS,EAAE,IAAI,EAAE,QAAQ,EAAE,YAAY,CAAC,OAAO,EAAE,YAAY,CAAC,OAAO,CACtF,CAAC;QACF,MAAM,MAAM,GAAG,GAAG,CAAC,eAAe,CAAC,GAAG,CAAC,CAAC;QACxC,MAAM,MAAM,GAAG,GAAG,CAAC,eAAe,CAAC,GAAG,CAAC,CAAC;QACxC,UAAU,CAAC,2BAA2B,CAAC,SAAS,EAAE,MAAM,EAAE,MAAM,CAAC,CAAC;IACpE,CAAC;IACD,6FAA6F;IAC7E,kBAAkB,CAAC,EAAgB;QACjD,MAAM,UAAU,GAAG,IAAI,CAAC,iBAAkB,CAAC;QAC3C,UAAU,CAAC,uCAAuC,CAAC,SAAS,EAAE,EAAE,CAAC,YAAY,CAAC,CAAC;IACjF,CAAC;IACD,4FAA4F;IAC5E,mBAAmB,CAAC,OAAsB;QACxD,MAAM,UAAU,GAAG,IAAI,CAAC,iBAAkB,CAAC;QAC3C,OAAO,CAAC,UAAU,CAAC,IAAI,CAAC,OAAO,CAAC,CAAC;QACjC,OAAO,CAAC,QAAQ,CAAC,IAAI,CAAC,OAAO,CAAC,CAAC;QAC/B,UAAU,CAAC,2BAA2B,CAAC,SAAS,EAAE,IAAI,CAAC,OAAO,EAAE,IAAI,CAAC,OAAO,CAAC,CAAC;IAChF,CAAC;IACD,uEAAuE;IACvD,UAAU,CAAC,IAAU;QACnC,MAAM,UAAU,GAAG,IAAI,CAAC,iBAAiB,GAAG,UAAU,CAAC,MAAM,EAAE,CAAC;QAChE,UAAU,CAAC,UAAU,GAAG,KAAK,CAAC;QAC9B,KAAK,MAAM,KAAK,IAAI,IAAI,CAAC,QAAQ;YAC/B,KAAK,CAAC,yBAAyB,CAAC,IAAI,CAAC,CAAC;QACxC,IAAI,CAAC,iBAAiB,GAAG,SAAS,CAAC;QACnC,OAAO,UAAU,CAAC;IACpB,CAAC;IACD,+EAA+E;IAC/D,kBAAkB,CAAC,MAAoB;QACrD,MAAM,eAAe,GAAiB,EAAE,CAAC;QACzC,IAAI,UAAU,GAAG,GAAG,CAAC;QACrB,IAAI,mBAA2C,CAAC;QAChD,KAAK,MAAM,KAAK,IAAI,MAAM,CAAC,QAAQ,EAAE,CAAC;YACpC,IAAI,KAAK,YAAY,IAAI,EAAE,CAAC;gBAC1B,MAAM,YAAY,GAAG,IAAI,CAAC,UAAU,CAAC,KAAK,CAAC,CAAC;gBAC5C,IAAI,YAAY,EAAE,CAAC;oBACjB,eAAe,CAAC,IAAI,CAAC,YAAY,CAAC,CAAC;oBACnC,MAAM,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,YAAY,CAAC,WAAW,CAAC,CAAC;oBAC7C,IAAI,CAAC,GAAG,UAAU,EAAE,CAAC;wBACnB,UAAU,GAAG,CAAC,CAAC;wBACf,mBAAmB,GAAG,YAAY,CAAC;oBACrC,CAAC;gBACH,CAAC;YACH,CAAC;QACH,CAAC;QACD,IAAI,mBAAmB,EAAE,CAAC;YACxB,MAAM,aAAa,GAAG,UAAU,CAAC,MAAM,EAAE,CAAC;YAC1C,MAAM,KAAK,GAAG,mBAAmB,CAAC,UAAU,CAAC,GAAG,CAAC,CAAC;YAClD,aAAa,CAAC,kBAAkB,CAAC,mBAAmB,EAAE,KAAK,CAAC,CAAC;YAC7D,KAAK,MAAM,YAAY,IAAI,eAAe,EAAE,CAAC;gBAC3C,IAAI,YAAY,KAAK,mBAAmB,EAAE,CAAC;oBACzC,MAAM,KAAK,GAAG,YAAY,CAAC,UAAU,CAAC,CAAC,GAAG,CAAC,CAAC;oBAC5C,aAAa,CAAC,kBAAkB,CAAC,YAAY,EAAE,KAAK,CAAC,CAAC;gBACxD,CAAC;YACH,CAAC;YACD,OAAO,aAAa,CAAC;QACvB,CAAC;QACD,OAAO,SAAS,CAAC;IACnB,CAAC;IACD,8EAA8E;IAC9D,iBAAiB,CAAC,MAAmB;QACnD,MAAM,aAAa,GAAG,UAAU,CAAC,MAAM,EAAE,CAAC;QAC1C,KAAK,MAAM,KAAK,IAAI,MAAM,CAAC,QAAQ,EAAE,CAAC;YACpC,MAAM,YAAY,GAAG,KAAK,CAAC,yBAAyB,CAAC,IAAI,CAAC,CAAC;YAC3D,IAAI,YAAY,EAAE,CAAC;gBACjB,MAAM,KAAK,GAAG,YAAY,CAAC,UAAU,CAAC,GAAG,CAAC,CAAC;gBAC3C,aAAa,CAAC,kBAAkB,CAAC,YAAY,EAAE,KAAK,CAAC,CAAC;YACxD,CAAC;QACH,CAAC;QACD,OAAO,aAAa,CAAC;IACvB,CAAC;IACO,cAAc,CAAiB;IAC/B,gBAAgB;QACtB,IAAI,IAAI,CAAC,cAAc;YACrB,OAAO,IAAI,CAAC,cAAc,CAAC;QAC7B,MAAM,OAAO,GAAG,aAAa,CAAC,eAAe,EAAE,CAAC;QAChD,gEAAgE;QAChE,OAAO,CAAC,QAAQ,GAAG,KAAK,CAAC,aAAa,CAAC,GAAG,CAAC,CAAC;QAC5C,IAAI,CAAC,cAAc,GAAG,OAAO,CAAC;QAC9B,OAAO,OAAO,CAAC;IACjB,CAAC;IACD;;;OAGG;IACI,oBAAoB,CAAC,EAAkB;QAC5C,MAAM,OAAO,GAAG,YAAY,CAAC,MAAM,EAAE,CAAC;QACtC,MAAM,OAAO,GAAG,IAAI,CAAC,gBAAgB,EAAE,CAAC;QACxC,EAAE,CAAC,WAAW,CAAC,OAAO,EAAE,OAAO,CAAC,CAAC;QACjC,IAAI,CAAC,kBAAkB,CAAC,OAAO,CAAC,CAAC;IACnC,CAAC;IACD,yEAAyE;IACzD,oBAAoB,CAAC,CAAiB;QACpD,OAAO,IAAI,CAAC,oBAAoB,CAAC,CAAC,CAAC,CAAC;IACtC,CAAC;IACD,0EAA0E;IAC1D,qBAAqB,CAAC,CAAkB;QACtD,OAAO,IAAI,CAAC,oBAAoB,CAAC,CAAC,CAAC,CAAC;IACtC,CAAC;IACD,0EAA0E;IAC1D,sBAAsB,CAAC,CAAqB;QAC1D,OAAO,IAAI,CAAC,oBAAoB,CAAC,CAAC,CAAC,CAAC;IACtC,CAAC;CACF","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 Curve\n */\n\n// import { Geometry, Angle, AngleSweep } from \"../Geometry\";\n\nimport { BSplineCurve3d } from \"../bspline/BSplineCurve\";\nimport { BSplineCurve3dH } from \"../bspline/BSplineCurve3dH\";\nimport { Angle } from \"../geometry3d/Angle\";\nimport { NullGeometryHandler } from \"../geometry3d/GeometryHandler\";\nimport { Point3d } from \"../geometry3d/Point3dVector3d\";\nimport { MomentData } from \"../geometry4d/MomentData\";\nimport { Arc3d } from \"./Arc3d\";\nimport { CurvePrimitive } from \"./CurvePrimitive\";\nimport { LineSegment3d } from \"./LineSegment3d\";\nimport { LineString3d } from \"./LineString3d\";\nimport { Loop } from \"./Loop\";\nimport { ParityRegion } from \"./ParityRegion\";\nimport { StrokeOptions } from \"./StrokeOptions\";\nimport { TransitionSpiral3d } from \"./spiral/TransitionSpiral3d\";\nimport { UnionRegion } from \"./UnionRegion\";\n\n/**\n * Implementation class for computing XY area moments.\n * @internal\n */\nexport class RegionMomentsXY extends NullGeometryHandler {\n  private _activeMomentData?: MomentData;\n  private _point0 = Point3d.create();\n  private _point1 = Point3d.create();\n  /**\n   * Accumulate (independent) integrations over:\n   * * The area between the arc and the chord connecting its endpoints.\n   * * The triangle with vertices: origin, arc start, arc end.\n   */\n  public override handleArc3d(arc: Arc3d): void {\n    const momentData = this._activeMomentData!;\n    const sweepRadians = arc.sweep.sweepRadians;\n    const alphaRadians = sweepRadians * 0.5;\n    let s = Math.sin(alphaRadians);\n    let c = Math.cos(alphaRadians);\n    let s1 = Math.sin(sweepRadians);\n    if (Angle.isFullCircleRadians(sweepRadians)) {\n      s = 0.0;\n      c = -1.0;\n      s1 = 0.0;\n    }\n    const q = 2 * s * s * s * c / (alphaRadians - s * c);\n    const s3 = s * s * s;\n    const s6 = s3 * s3;\n    const area = 0.5 * (sweepRadians - s1);\n    const inertiaXX = 0.25 * area * (1.0 - q / 3.0);\n    const inertiaYY1 = 0.25 * area * (1.0 + q);\n    const inertiaYY = inertiaYY1 - 4.0 * s6 / (9.0 * area);\n    const productXX = inertiaYY;\n    const productYY = inertiaXX;\n    const centerToCentroid = 4.0 * s * s * s / (3.0 * (sweepRadians - s1));\n    const midRadians = arc.sweep.fractionToRadians(0.5);\n    const centralPlane = arc.radiansToRotatedBasis(midRadians);\n    const centroid = centralPlane.origin.plusScaled(centralPlane.vectorU, centerToCentroid);\n    momentData.accumulateXYProductsInCentroidalFrame(\n      productXX, 0.0, productYY, area, centroid, centralPlane.vectorU, centralPlane.vectorV,\n    );\n    const pointB = arc.fractionToPoint(0.0);\n    const pointC = arc.fractionToPoint(1.0);\n    momentData.accumulateTriangleMomentsXY(undefined, pointB, pointC);\n  }\n  /** Accumulate integrals over the (triangular) areas from the origin to each line segment. */\n  public override handleLineString3d(ls: LineString3d): void {\n    const momentData = this._activeMomentData!;\n    momentData.accumulateTriangleToLineStringMomentsXY(undefined, ls.packedPoints);\n  }\n  /** Accumulate integrals over the (triangular) area from the origin to this line segment. */\n  public override handleLineSegment3d(segment: LineSegment3d): void {\n    const momentData = this._activeMomentData!;\n    segment.startPoint(this._point0);\n    segment.endPoint(this._point1);\n    momentData.accumulateTriangleMomentsXY(undefined, this._point0, this._point1);\n  }\n  /** Accumulate integrals from origin to all primitives in the chain. */\n  public override handleLoop(loop: Loop): MomentData | undefined {\n    const momentData = this._activeMomentData = MomentData.create();\n    momentData.needOrigin = false;\n    for (const child of loop.children)\n      child.dispatchToGeometryHandler(this);\n    this._activeMomentData = undefined;\n    return momentData;\n  }\n  /** Accumulate integrals from origin to the components of the parity region. */\n  public override handleParityRegion(region: ParityRegion): MomentData | undefined {\n    const allChildMoments: MomentData[] = [];\n    let maxAbsArea = 0.0;\n    let largestChildMoments: MomentData | undefined;\n    for (const child of region.children) {\n      if (child instanceof Loop) {\n        const childMoments = this.handleLoop(child);\n        if (childMoments) {\n          allChildMoments.push(childMoments);\n          const q = Math.abs(childMoments.quantitySum);\n          if (q > maxAbsArea) {\n            maxAbsArea = q;\n            largestChildMoments = childMoments;\n          }\n        }\n      }\n    }\n    if (largestChildMoments) {\n      const summedMoments = MomentData.create();\n      const sign0 = largestChildMoments.signFactor(1.0);\n      summedMoments.accumulateProducts(largestChildMoments, sign0);\n      for (const childMoments of allChildMoments) {\n        if (childMoments !== largestChildMoments) {\n          const sign1 = childMoments.signFactor(-1.0);\n          summedMoments.accumulateProducts(childMoments, sign1);\n        }\n      }\n      return summedMoments;\n    }\n    return undefined;\n  }\n  /** Accumulate integrals from origin to the components of the union region. */\n  public override handleUnionRegion(region: UnionRegion): MomentData | undefined {\n    const summedMoments = MomentData.create();\n    for (const child of region.children) {\n      const childMoments = child.dispatchToGeometryHandler(this);\n      if (childMoments) {\n        const sign0 = childMoments.signFactor(1.0);\n        summedMoments.accumulateProducts(childMoments, sign0);\n      }\n    }\n    return summedMoments;\n  }\n  private _strokeOptions?: StrokeOptions;\n  private getStrokeOptions(): StrokeOptions {\n    if (this._strokeOptions)\n      return this._strokeOptions;\n    const options = StrokeOptions.createForCurves();\n    // this is unusually fine for stroking, but appropriate for sum.\n    options.angleTol = Angle.createDegrees(5.0);\n    this._strokeOptions = options;\n    return options;\n  }\n  /**\n   * Handle a single curve primitive (not loop).\n   * * Stroke the curve and accumulate stroke array.\n   */\n  public handleCurvePrimitive(cp: CurvePrimitive) {\n    const strokes = LineString3d.create();\n    const options = this.getStrokeOptions();\n    cp.emitStrokes(strokes, options);\n    this.handleLineString3d(strokes);\n  }\n  /** Handle strongly typed  BSplineCurve3d  as generic curve primitive. */\n  public override handleBSplineCurve3d(g: BSplineCurve3d) {\n    return this.handleCurvePrimitive(g);\n  }\n  /** Handle strongly typed  BSplineCurve3dH  as generic curve primitive. */\n  public override handleBSplineCurve3dH(g: BSplineCurve3dH) {\n    return this.handleCurvePrimitive(g);\n  }\n  /** Handle strongly typed  TransitionSpiral as generic curve primitive. */\n  public override handleTransitionSpiral(g: TransitionSpiral3d) {\n    return this.handleCurvePrimitive(g);\n  }\n}\n"]}