@itwin/core-geometry 5.0.0-dev.6 → 5.0.0-dev.61

package/lib/esm/geometry4d/MomentData.js

@@ -5,8 +5,8 @@
 /** @packageDocumentation
  * @module Numerics
  */
-// import { Point2d } from "./Geometry2d";
 /* eslint-disable @typescript-eslint/naming-convention */
+// cspell:word ABAT
 import { Geometry } from "../Geometry";
 import { Matrix3d } from "../geometry3d/Matrix3d";
 import { Point3d, Vector3d } from "../geometry3d/Point3dVector3d";
@@ -19,61 +19,84 @@ import { Point4d } from "./Point4d";
  *   * origin = local origin used as moments are summed.
  *   * sums = array of summed moments.
  *     * The [i,j] entry of the sums is a summed or integrated moment for product of axis i and j.
- *       * axes 0,1,2 are x,y,z
+ *       * axes 0,1,2 are x,y,z.
  *         * e.g. entry [0,1] is summed product xy
  *       * axis 3 is "w", which is 1 in sums.
- *         * e.g. entry 03 is summed x
- *    * In this level,
+ *         * 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 `radiiOfGyration` are filled in.
+ *  * Second level: after a call to inertiaProductsToPrincipalAxes, the `localToWorldMap`, `absoluteQuantity` and
+ *    `radiiOfGyration` are filled in.
  * @public
  */
 export class MomentData {
-    /** Return the lower-right (3,3) entry in the sums.
-     * * This is the quantity (i.e. length, area, or volume) summed
+    /** Origin used for sums. */
+    origin;
+    /** Flag to request deferred origin setup. */
+    needOrigin;
+    /**
+     * Moment sums.
+     * * Set to zero at initialization and if requested later.
+     * * Accumulated during data entry phase.
+     */
+    sums;
+    /**
+     * The mapping between principal and world system.
+     * * This set up with its inverse already constructed.
      */
-    get quantitySum() { return this.sums.atIJ(3, 3); }
-    /** Return a scale factor to make these sums match the target orientation sign.
+    localToWorldMap;
+    // private variables
+    static _vectorA;
+    static _vectorB;
+    static _vectorC;
+    _point0 = Point3d.create();
+    _point1 = Point3d.create();
+    /**
+     * 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 a scale factor to make these sums match the target orientation sign.
      * * 1.0 if `this.quantitySum` has the same sign as `targetSign`.
      * * -1.0 if `this.quantitySum` has the opposite sign from `targetSign`
      */
     signFactor(targetSign) {
         return targetSign * this.quantitySum > 0 ? 1.0 : -1.0;
     }
-    /**
-     *  If `this.needOrigin` flag is set, copy `origin` to `this.origin` and clear the flag.
-     *
-     */
+    /** If `this.needOrigin` flag is set, copy `origin` to `this.origin` and clear the flag. */
     setOriginIfNeeded(origin) {
         if (this.needOrigin) {
             this.origin.setFromPoint3d(origin);
             this.needOrigin = false;
         }
     }
-    /**
-     *  If `this.needOrigin` flag is set, copy `origin` to `this.origin` and clear the flag.
-     *
-     */
+    /** If `this.needOrigin` flag is set, copy `origin` to `this.origin` and clear the flag. */
     setOriginFromGrowableXYZArrayIfNeeded(points) {
         if (this.needOrigin && points.length > 0) {
             points.getPoint3dAtCheckedPointIndex(0, this.origin);
             this.needOrigin = false;
         }
     }
-    /**
-     *  If `this.needOrigin` flag is set, copy `origin` to `this.origin` and clear the flag.
-     *
-     */
+    /** If `this.needOrigin` flag is set, copy `origin` to `this.origin` and clear the flag. */
     setOriginXYZIfNeeded(x, y, z) {
         if (this.needOrigin) {
             this.origin.set(x, y, z);
             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._point0 = Point3d.create();
-        this._point1 = Point3d.create();
         this.origin = Point3d.createZero();
         this.sums = Matrix4d.createZero();
         this.localToWorldMap = Transform.createIdentity();
@@ -82,12 +105,13 @@ export class MomentData {
         this.absoluteQuantity = 0.1; // so optimizer sees its type
         this.absoluteQuantity = undefined;
     }
-    /** Create moments with optional origin.
-     * * origin and needOrigin are quirky.
+    /**
+     * Create moments with optional origin.
+     * * Origin and needOrigin are quirky.
      *   * (undefined, true) sets up to use first incoming point as origin.
-     *   * (origin) sets up to use that durable origin, set needsOrigin flag false
-     *   * (origin, true) the "true" is meaningless
-     *   * (undefined, false) makes 000 the durable origin
+     *   * (origin) sets up to use that durable origin, set needsOrigin flag false.
+     *   * (origin, true) the "true" is meaningless.
+     *   * (undefined, false) makes 000 the durable origin.
     */
     static create(origin, needOrigin = false) {
         const data = new MomentData();
@@ -99,8 +123,8 @@ export class MomentData {
         return data;
     }
     /**
-     * Return the formal tensor of integrated values `[yy+zz,xy,xz][yx,xx+zz,yz][zx,xy,xx+yy]`
-     * @param products matrix of (integrated) `[xx,xy,xz][yx,yy,yz][zx,xy,zz]`
+     * Return the formal tensor of integrated values `[yy+zz,xy,xz][yx,xx+zz,yz][zx,xy,xx+yy]`.
+     * @param products matrix of (integrated) `[xx,xy,xz][yx,yy,yz][zx,xy,zz]`.
      */
     static momentTensorFromInertiaProducts(products) {
         const rr = products.sumDiagonal();
@@ -133,7 +157,7 @@ export class MomentData {
     }
     /**
      * Return the principal moment data for an array of points.
-     * @param points array of points
+     * @param points array of points.
      */
     static pointsToPrincipalAxes(points) {
         const moments = new MomentData();
@@ -144,12 +168,12 @@ export class MomentData {
         return this.inertiaProductsToPrincipalAxes(moments.origin, moments.sums);
     }
     /**
-     * Compute principal axes from inertial products
-     * * The radii of gyration are sorted smallest to largest
-     * * Hence x axis is long direction
-     * * Hence planar data generates large moment as Z
+     * Compute principal axes from inertial products.
+     * * The radii of gyration are sorted smallest to largest.
+     * * 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();
@@ -177,18 +201,19 @@ export class MomentData {
         return moments;
     }
     /**
-     * Test for match among selected members as they exist after `inertiaProductsToPrincipalAxes`
-     * * The members considered are
-     *   * origin of local to world map (i.e. centroid)
-     *   * radius of gyration
+     * Test for match among selected members as they exist after `inertiaProductsToPrincipalAxes`.
+     * * The members considered are:
+     *   * origin of local to world map (i.e. centroid),
+     *   * radius of gyration,
      *   * axes of localToWorldMap.
      * * Axis direction tests allow these quirks:
      *   * opposite orientation is considered matched.
      * * Full xyz symmetry: If x,y,z radii are matched, axes are not tested.
-     * * Symmetry in xy plane: If x and y radii are matched, the x and y axes area allowed to spin freely.  Only Z direction is tested.
+     * * Symmetry in xy plane: If x and y radii are matched, the x and y axes area allowed to spin freely. Only Z direction
+     *   is tested.
      * * If either or both are undefined, returns false.
-     * @param dataA first set of moments
-     * @param dataB second set of moments
+     * @param dataA first set of moments.
+     * @param dataB second set of moments.
      */
     static areEquivalentPrincipalAxes(dataA, dataB) {
         if (dataA && dataB
@@ -236,7 +261,7 @@ export class MomentData {
             this.sums.addMomentsInPlace(p.x - this.origin.x, p.y - this.origin.y, p.z - this.origin.z, 1.0);
         }
     }
-    /** revise the accumulated sums to be "around the centroid" */
+    /** Revise the accumulated sums to be "around the centroid". */
     shiftOriginAndSumsToCentroidOfSums() {
         const xyz = this.sums.columnW().realPoint();
         if (xyz) {
@@ -245,7 +270,8 @@ export class MomentData {
         }
         return false;
     }
-    /** revise the accumulated sums
+    /**
+     * Revise the accumulated sums.
      * * add ax,ay,ax to the origin coordinates.
      * * apply the negative translation to the sums.
     */
@@ -253,11 +279,12 @@ export class MomentData {
         this.origin.addXYZInPlace(ax, ay, az);
         this.sums.multiplyTranslationSandwichInPlace(-ax, -ay, -az);
     }
-    /** revise the accumulated sums so they are based at a specified origin. */
+    /** Revise the accumulated sums so they are based at a specified origin. */
     shiftOriginAndSumsToNewOrigin(newOrigin) {
         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.
+    /**
+     * 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
@@ -289,7 +316,7 @@ export 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.create(point.x - this.origin.x, point.y - this.origin.y, point.z - this.origin.z, 1.0, MomentData._vectorA);
@@ -311,7 +338,8 @@ export class MomentData {
         this.sums.addScaledOuterProductInPlace(vectorB, vectorA, r1_6);
         this.sums.addScaledOuterProductInPlace(vectorB, vectorB, r1_3);
     }
-    /** compute moments of triangles from a base point to the given linestring.
+    /**
+     * 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.
@@ -328,18 +356,18 @@ export class MomentData {
             }
         }
     }
-    // cspell:word ABAT
     /**
      * * 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)
+     * * 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.createRowValues(productXX, productXY, 0, 0, productXY, productYY, 0, 0, 0, 0, 0, 0, 0, 0, 0, area);
@@ -351,27 +379,23 @@ export class MomentData {
     }
     /**
      * Accumulate sums from other moments.
-     * * scale by given scaleFactor (e.g. sign to correct orientation)
-     * * pull the origin from `other` if `this` needs an origin.
-     * *
+     * * Scale by given scaleFactor (e.g. sign to correct orientation).
+     * * Pull the origin from `other` if `this` needs an origin.
      */
     accumulateProducts(other, scale) {
         this.setOriginIfNeeded(other.origin);
         this.sums.addTranslationSandwichInPlace(other.sums, this.origin.x - other.origin.x, this.origin.y - other.origin.y, this.origin.z - other.origin.z, scale);
     }
     /**
-   * Accumulate sums from Matrix4d and origin.
-   * * scale by given scaleFactor (e.g. sign to correct orientation)
-   * * trap the origin if `this` needs an origin.
-   * *
-   */
+     * Accumulate sums from Matrix4d and origin.
+     * * Scale by given scaleFactor (e.g. sign to correct orientation).
+     * * Trap the origin if `this` needs an origin.
+     */
     accumulateProductsFromOrigin(origin, products, scale) {
         this.setOriginIfNeeded(origin);
         this.sums.addTranslationSandwichInPlace(products, this.origin.x - origin.x, this.origin.y - origin.y, this.origin.z - origin.z, scale);
     }
-    /**
-     * Convert to a json data object with:
-     */
+    /** Convert to a json data object. */
     toJSON() {
         return {
             origin: this.origin,

package/lib/esm/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,0CAA0C;AAC1C,yDAAyD;AAEzD,OAAO,EAAE,QAAQ,EAAE,MAAM,aAAa,CAAC;AAEvC,OAAO,EAAE,QAAQ,EAAE,MAAM,wBAAwB,CAAC;AAClD,OAAO,EAAE,OAAO,EAAE,QAAQ,EAAO,MAAM,+BAA+B,CAAC;AACvE,OAAO,EAAE,SAAS,EAAE,MAAM,yBAAyB,CAAC;AAEpD,OAAO,EAAE,QAAQ,EAAE,MAAM,YAAY,CAAC;AACtC,OAAO,EAAE,OAAO,EAAE,MAAM,WAAW,CAAC;AAEpC;;;;;;;;;;;;;;;GAeG;AACH,MAAM,OAAO,UAAU;IAcrB;;OAEG;IACH,IAAW,WAAW,KAAa,OAAO,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;IACjE;;;OAGG;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;IAED;;;OAGG;IACI,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;;;OAGG;IACI,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;IAED;;;OAGG;IACI,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;IAUD;QA6PQ,YAAO,GAAG,OAAO,CAAC,MAAM,EAAE,CAAC;QAC3B,YAAO,GAAG,OAAO,CAAC,MAAM,EAAE,CAAC;QA7PjC,IAAI,CAAC,MAAM,GAAG,OAAO,CAAC,UAAU,EAAE,CAAC;QACnC,IAAI,CAAC,IAAI,GAAG,QAAQ,CAAC,UAAU,EAAE,CAAC;QAClC,IAAI,CAAC,eAAe,GAAG,SAAS,CAAC,cAAc,EAAE,CAAC;QAClD,IAAI,CAAC,gBAAgB,GAAG,QAAQ,CAAC,MAAM,EAAE,CAAC;QAC1C,IAAI,CAAC,UAAU,GAAG,KAAK,CAAC;QACxB,IAAI,CAAC,gBAAgB,GAAG,GAAG,CAAC,CAAG,6BAA6B;QAC5D,IAAI,CAAC,gBAAgB,GAAG,SAAS,CAAC;IACpC,CAAC;IACD;;;;;;MAME;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,QAAQ,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;IAED;;;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,QAAQ,CAAC,MAAM,EAAE,CAAC;QAClC,MAAM,WAAW,GAAG,QAAQ,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,SAAS,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;;;;;;;;;;;;;OAaG;IACI,MAAM,CAAC,0BAA0B,CAAC,KAA6B,EAAE,KAA6B;QACnG,IAAI,KAAK,IAAI,KAAK;eACb,QAAQ,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,QAAQ,CAAC,gBAAgB,CAAC,KAAK,CAAC,gBAAgB,CAAC,CAAC,EAAE,KAAK,CAAC,gBAAgB,CAAC,CAAC,CAAC,EAAE,CAAC;oBAClF,oCAAoC;oBACpC,IAAI,QAAQ,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,QAAQ,CAAC,MAAM,EAAE,CAAC;gBAClC,MAAM,OAAO,GAAG,QAAQ,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,CAAC,CAAC;QACT,CAAC;IACH,CAAC;IACD,8DAA8D;IACvD,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;;;MAGE;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;IAKD;;;;MAIE;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,OAAO,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,OAAO,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,OAAO,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,OAAO,CAAC,MAAM,CAAC,MAAM,CAAC,CAAC,GAAG,EAAE,EAAE,MAAM,CAAC,CAAC,GAAG,EAAE,EAAE,GAAG,EAAE,GAAG,EAAE,UAAU,CAAC,QAAQ,CAAC,CAAC;QAElH,gHAAgH;QAChH,mBAAmB;QACnB,iEAAiE;QACjE,MAAM,IAAI,GAAG,QAAQ,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,uEAAuE;IAChE,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,OAAO,CAAC,MAAM,CAAC,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,CAAC,CAAC;QAC1J,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,OAAO,CAAC,MAAM,CAAC,MAAM,CAAC,CAAC,GAAG,EAAE,EAAE,MAAM,CAAC,CAAC,GAAG,EAAE,EAAE,MAAM,CAAC,CAAC,GAAG,EAAE,EAAE,GAAG,EAAE,UAAU,CAAC,QAAQ,CAAC,CAAC;QAC5H,MAAM,OAAO,GAAG,UAAU,CAAC,QAAQ,GAAG,OAAO,CAAC,MAAM,CAAC,MAAM,CAAC,CAAC,GAAG,EAAE,EAAE,MAAM,CAAC,CAAC,GAAG,EAAE,EAAE,MAAM,CAAC,CAAC,GAAG,EAAE,EAAE,GAAG,EAAE,UAAU,CAAC,QAAQ,CAAC,CAAC;QAC5H,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;IAEjE,CAAC;IAID;;;;;OAKG;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,mBAAmB;IACnB;;;;;;;;;;;OAWG;IACI,qCAAqC,CAAC,SAAiB,EAAE,SAAiB,EAAE,SAAiB,EAAE,IAAY,EAChH,MAAa,EAAE,OAAc,EAAE,OAAc;QAC7C,MAAM,kBAAkB,GAAG,QAAQ,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,CAAC,CAAC;QACjB,MAAM,IAAI,GAAG,QAAQ,CAAC,gBAAgB,CAAC,OAAO,CAAC,CAAC,EAAE,OAAO,CAAC,CAAC,EAAE,OAAO,CAAC,CAAC,EAAE,OAAO,CAAC,CAAC,CAAC,CAAC;QACnF,MAAM,SAAS,GAAG,QAAQ,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,CAAC,CAAC;QACd,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;;;;;OAKG;IACI,kBAAkB,CAAC,KAAiB,EAAE,KAAa;QACxD,IAAI,CAAC,iBAAiB,CAAC,KAAK,CAAC,MAAM,CAAC,CAAC;QACrC,IAAI,CAAC,IAAI,CAAC,6BAA6B,CAAC,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,CAAC,CAAC;IAC7J,CAAC;IAED;;;;;KAKC;IACM,4BAA4B,CAAC,MAAe,EAAE,QAAkB,EAAE,KAAa;QACpF,IAAI,CAAC,iBAAiB,CAAC,MAAM,CAAC,CAAC;QAC/B,IAAI,CAAC,IAAI,CAAC,6BAA6B,CAAC,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,CAAC,CAAC;IACzI,CAAC;IACD;;OAEG;IACI,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","sourcesContent":["/*---------------------------------------------------------------------------------------------\r\n* Copyright (c) Bentley Systems, Incorporated. All rights reserved.\r\n* See LICENSE.md in the project root for license terms and full copyright notice.\r\n*--------------------------------------------------------------------------------------------*/\r\n\r\n/** @packageDocumentation\r\n * @module Numerics\r\n */\r\n\r\n// import { Point2d } from \"./Geometry2d\";\r\n/* eslint-disable @typescript-eslint/naming-convention */\r\n\r\nimport { Geometry } from \"../Geometry\";\r\nimport { GrowableXYZArray } from \"../geometry3d/GrowableXYZArray\";\r\nimport { Matrix3d } from \"../geometry3d/Matrix3d\";\r\nimport { Point3d, Vector3d, XYZ } from \"../geometry3d/Point3dVector3d\";\r\nimport { Transform } from \"../geometry3d/Transform\";\r\nimport { XAndY, XYAndZ } from \"../geometry3d/XYZProps\";\r\nimport { Matrix4d } from \"./Matrix4d\";\r\nimport { Point4d } from \"./Point4d\";\r\n\r\n/**\r\n * A MomentData structure exists in several levels:\r\n * * First level: as a carrier of sums of inertial products that determine moments.\r\n *   * origin = local origin used as moments are summed.\r\n *   * sums = array of summed moments.\r\n *     * The [i,j] entry of the sums is a summed or integrated moment for product of axis i and j.\r\n *       * axes 0,1,2 are x,y,z\r\n *         * e.g. entry [0,1] is summed product xy\r\n *       * axis 3 is \"w\", which is 1 in sums.\r\n *         * e.g. entry 03 is summed x\r\n *    * In this level,\r\n *        * the `absoluteQuantity` member is undefined.\r\n *        * the `localToWorldMap` and `radiiOfGyration` are created by have undefined contents.\r\n *  * Second level: after a call to inertiaProductsToPrincipalAxes, the `localToWorldMap`, `absoluteQuantity` and `radiiOfGyration` are filled in.\r\n * @public\r\n */\r\nexport class MomentData {\r\n  /** Origin used for sums. */\r\n  public origin: Point3d;\r\n  /** flag to request deferred origin setup. */\r\n  public needOrigin: boolean;\r\n  /** Moment sums.\r\n   * * Set to zero at initialization and if requested later.\r\n   * * Accumulated during data entry phase.\r\n   */\r\n  public sums: Matrix4d;\r\n  /** the mapping between principal and world system.\r\n   * * This set up with its inverse already constructed.\r\n   */\r\n  public localToWorldMap: Transform;\r\n  /** Return the lower-right (3,3) entry in the sums.\r\n   * * This is the quantity (i.e. length, area, or volume) summed\r\n   */\r\n  public get quantitySum(): number { return this.sums.atIJ(3, 3); }\r\n  /** Return a scale factor to make these sums match the target orientation sign.\r\n   * * 1.0 if `this.quantitySum` has the same sign as `targetSign`.\r\n   * * -1.0 if `this.quantitySum` has the opposite sign from `targetSign`\r\n   */\r\n  public signFactor(targetSign: number): number {\r\n    return targetSign * this.quantitySum > 0 ? 1.0 : -1.0;\r\n  }\r\n\r\n  /**\r\n   *  If `this.needOrigin` flag is set, copy `origin` to `this.origin` and clear the flag.\r\n   *\r\n   */\r\n  public setOriginIfNeeded(origin: Point3d) {\r\n    if (this.needOrigin) {\r\n      this.origin.setFromPoint3d(origin);\r\n      this.needOrigin = false;\r\n    }\r\n  }\r\n  /**\r\n   *  If `this.needOrigin` flag is set, copy `origin` to `this.origin` and clear the flag.\r\n   *\r\n   */\r\n  public setOriginFromGrowableXYZArrayIfNeeded(points: GrowableXYZArray) {\r\n    if (this.needOrigin && points.length > 0) {\r\n      points.getPoint3dAtCheckedPointIndex(0, this.origin);\r\n      this.needOrigin = false;\r\n    }\r\n  }\r\n\r\n  /**\r\n   *  If `this.needOrigin` flag is set, copy `origin` to `this.origin` and clear the flag.\r\n   *\r\n   */\r\n  public setOriginXYZIfNeeded(x: number, y: number, z: number) {\r\n    if (this.needOrigin) {\r\n      this.origin.set(x, y, z);\r\n      this.needOrigin = false;\r\n    }\r\n  }\r\n\r\n  /** radii of gyration (square roots of principal second moments)\r\n   */\r\n  public radiusOfGyration: Vector3d;\r\n  /** principal quantity (e.g. length, area, or volume).  This is undefined in raw moments, and becomes defined by\r\n   *\r\n   */\r\n  public absoluteQuantity?: number;\r\n\r\n  private constructor() {\r\n    this.origin = Point3d.createZero();\r\n    this.sums = Matrix4d.createZero();\r\n    this.localToWorldMap = Transform.createIdentity();\r\n    this.radiusOfGyration = Vector3d.create();\r\n    this.needOrigin = false;\r\n    this.absoluteQuantity = 0.1;   // so optimizer sees its type\r\n    this.absoluteQuantity = undefined;\r\n  }\r\n  /** Create moments with optional origin.\r\n   * * origin and needOrigin are quirky.\r\n   *   * (undefined, true) sets up to use first incoming point as origin.\r\n   *   * (origin) sets up to use that durable origin, set needsOrigin flag false\r\n   *   * (origin, true) the \"true\" is meaningless\r\n   *   * (undefined, false) makes 000 the durable origin\r\n  */\r\n  public static create(origin?: Point3d | undefined, needOrigin: boolean = false): MomentData {\r\n    const data = new MomentData();\r\n    data.needOrigin = needOrigin;\r\n    if (origin) {\r\n      data.origin.setFromPoint3d(origin);\r\n      data.needOrigin = false;\r\n    }\r\n    return data;\r\n  }\r\n  /**\r\n   * Return the formal tensor of integrated values `[yy+zz,xy,xz][yx,xx+zz,yz][zx,xy,xx+yy]`\r\n   * @param products matrix of (integrated) `[xx,xy,xz][yx,yy,yz][zx,xy,zz]`\r\n   */\r\n  public static momentTensorFromInertiaProducts(products: Matrix3d): Matrix3d {\r\n    const rr = products.sumDiagonal();\r\n    const result = Matrix3d.createScale(rr, rr, rr);\r\n    result.addScaledInPlace(products, -1.0);\r\n    return result;\r\n  }\r\n  /** Sort the columns of the matrix for increasing moments. */\r\n  public static sortColumnsForIncreasingMoments(axes: Matrix3d, moments: Vector3d) {\r\n    const points = [\r\n      axes.indexedColumnWithWeight(0, moments.x),\r\n      axes.indexedColumnWithWeight(1, moments.y),\r\n      axes.indexedColumnWithWeight(2, moments.z),\r\n    ].sort((dataA: Point4d, dataB: Point4d): number => {\r\n      if (dataA.w < dataB.w) return -1;\r\n      if (dataA.w > dataB.w) return 1;\r\n      return 0;\r\n    });\r\n    axes.setColumnsPoint4dXYZ(points[0], points[1], points[2]);\r\n    if (axes.determinant() < 0)\r\n      axes.scaleColumnsInPlace(-1.0, -1.0, -1.0);\r\n    // prefer x and z positive -- y falls wherever . ..\r\n    if (axes.at(0, 0) < 0.0)\r\n      axes.scaleColumnsInPlace(-1.0, -1.0, 1.0);\r\n    if (axes.at(2, 2) < 0.0)\r\n      axes.scaleColumnsInPlace(1.0, -1.0, -1.0);\r\n    moments.set(points[0].w, points[1].w, points[2].w);\r\n  }\r\n\r\n  /**\r\n   * Return the principal moment data for an array of points.\r\n   * @param points array of points\r\n   */\r\n  public static pointsToPrincipalAxes(points: Point3d[]): MomentData | undefined {\r\n    const moments = new MomentData();\r\n    if (points.length === 0)\r\n      return moments;\r\n    moments.clearSums(points[0]);\r\n    moments.accumulatePointMomentsFromOrigin(points);\r\n    return this.inertiaProductsToPrincipalAxes(moments.origin, moments.sums);\r\n  }\r\n  /**\r\n   * Compute principal axes from inertial products\r\n   * * The radii of gyration are sorted smallest to largest\r\n   * * Hence x axis is long direction\r\n   * * Hence planar data generates large moment as Z\r\n   * @param origin The origin used for the inertia products.\r\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]\r\n   */\r\n  public static inertiaProductsToPrincipalAxes(origin: XYZ, inertiaProducts: Matrix4d): MomentData | undefined {\r\n    const moments = new MomentData();\r\n    moments.sums.setFrom(inertiaProducts);\r\n    moments.origin.setFrom(origin);\r\n    if (!moments.shiftOriginAndSumsToCentroidOfSums())\r\n      return undefined;\r\n    const products = moments.sums.matrixPart();\r\n    const w = moments.sums.weight();\r\n    if (w < 0.0)\r\n      products.scaleColumnsInPlace(-1, -1, -1);\r\n    const tensor = MomentData.momentTensorFromInertiaProducts(products);\r\n    const moment2 = Vector3d.create();\r\n    const axisVectors = Matrix3d.createZero();\r\n    tensor.fastSymmetricEigenvalues(axisVectors, moment2);\r\n    if (moment2.x < 0.0)\r\n      return undefined;\r\n    MomentData.sortColumnsForIncreasingMoments(axisVectors, moment2);\r\n    if (w < 0.0)\r\n      axisVectors.scaleColumnsInPlace(1, -1, -1);\r\n    moments.localToWorldMap = Transform.createOriginAndMatrix(moments.origin, axisVectors);\r\n    moments.radiusOfGyration.set(\r\n      Math.sqrt(Math.abs(moment2.x)), Math.sqrt(Math.abs(moment2.y)), Math.sqrt(Math.abs(moment2.z)));\r\n    moments.radiusOfGyration.scaleInPlace(1.0 / Math.sqrt(Math.abs(w)));\r\n    moments.absoluteQuantity = Math.abs(w);\r\n    return moments;\r\n  }\r\n  /**\r\n   * Test for match among selected members as they exist after `inertiaProductsToPrincipalAxes`\r\n   * * The members considered are\r\n   *   * origin of local to world map (i.e. centroid)\r\n   *   * radius of gyration\r\n   *   * axes of localToWorldMap.\r\n   * * Axis direction tests allow these quirks:\r\n   *   * opposite orientation is considered matched.\r\n   * * Full xyz symmetry: If x,y,z radii are matched, axes are not tested.\r\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 is tested.\r\n   * * If either or both are undefined, returns false.\r\n   * @param dataA first set of moments\r\n   * @param dataB second set of moments\r\n   */\r\n  public static areEquivalentPrincipalAxes(dataA: MomentData | undefined, dataB: MomentData | undefined): boolean {\r\n    if (dataA && dataB\r\n      && Geometry.isSameCoordinate(dataA.quantitySum, dataB.quantitySum)) {  // um.. need different tolerance for area, volume?)\r\n      if (dataA.localToWorldMap.getOrigin().isAlmostEqual(dataB.localToWorldMap.getOrigin())\r\n        && dataA.radiusOfGyration.isAlmostEqual(dataB.radiusOfGyration)) {\r\n        if (Geometry.isSameCoordinate(dataA.radiusOfGyration.x, dataA.radiusOfGyration.y)) {\r\n          // We have at least xy symmetry ....\r\n          if (Geometry.isSameCoordinate(dataA.radiusOfGyration.x, dataA.radiusOfGyration.z))\r\n            return true;\r\n          // just xy.\r\n          // allow opposite z directions.\r\n          // If the z's are aligned, x an dy can spin freely.\r\n          const zA = dataA.localToWorldMap.matrix.columnZ();\r\n          const zB = dataB.localToWorldMap.matrix.columnZ();\r\n          if (zA.isParallelTo(zB, true))\r\n            return true;\r\n          return false;\r\n        }\r\n        // no symmetry.  Test all three axes.\r\n        const vectorA = Vector3d.create();\r\n        const vectorB = Vector3d.create();\r\n        for (let i = 0; i < 3; i++) {\r\n          dataA.localToWorldMap.matrix.getColumn(i, vectorA);\r\n          dataB.localToWorldMap.matrix.getColumn(i, vectorB);\r\n          if (!vectorA.isParallelTo(vectorB, true))\r\n            return false;\r\n        }\r\n        return true;\r\n      }\r\n    }\r\n    return false;\r\n  }\r\n  /** Clear the MomentData sums to zero, and establish a new origin. */\r\n  public clearSums(origin?: Point3d) {\r\n    this.sums.setZero();\r\n    if (origin)\r\n      this.origin.setFrom(origin);\r\n    else\r\n      this.origin.setZero();\r\n  }\r\n  /** Accumulate products-of-components for given points. */\r\n  public accumulatePointMomentsFromOrigin(points: Point3d[]) {\r\n    for (const p of points) {\r\n      this.sums.addMomentsInPlace(\r\n        p.x - this.origin.x,\r\n        p.y - this.origin.y,\r\n        p.z - this.origin.z,\r\n        1.0);\r\n    }\r\n  }\r\n  /** revise the accumulated sums to be \"around the centroid\" */\r\n  public shiftOriginAndSumsToCentroidOfSums(): boolean {\r\n    const xyz = this.sums.columnW().realPoint();\r\n    if (xyz) {\r\n      this.shiftOriginAndSumsByXYZ(xyz.x, xyz.y, xyz.z);\r\n      return true;\r\n    }\r\n    return false;\r\n  }\r\n\r\n  /** revise the accumulated sums\r\n   * * add ax,ay,ax to the origin coordinates.\r\n   * * apply the negative translation to the sums.\r\n  */\r\n  public shiftOriginAndSumsByXYZ(ax: number, ay: number, az: number) {\r\n    this.origin.addXYZInPlace(ax, ay, az);\r\n    this.sums.multiplyTranslationSandwichInPlace(-ax, -ay, -az);\r\n  }\r\n  /** revise the accumulated sums so they are based at a specified origin. */\r\n  public shiftOriginAndSumsToNewOrigin(newOrigin: XYAndZ) {\r\n    this.shiftOriginAndSumsByXYZ(newOrigin.x - this.origin.x, newOrigin.y - this.origin.y, newOrigin.z - this.origin.z);\r\n  }\r\n  private static _vectorA?: Point4d;\r\n  private static _vectorB?: Point4d;\r\n  private static _vectorC?: Point4d;\r\n\r\n  /** compute moments of a triangle from the origin to the given line.\r\n   * Accumulate them to this.sums.\r\n   * * If `pointA` is undefined, use `this.origin` as pointA.\r\n   * * If `this.needOrigin` is set, pointB is used\r\n  */\r\n  public accumulateTriangleMomentsXY(pointA: XAndY | undefined, pointB: XAndY, pointC: XAndY) {\r\n    this.setOriginXYZIfNeeded(pointB.x, pointB.y, 0.0);\r\n    const x0 = this.origin.x;\r\n    const y0 = this.origin.y;\r\n    const vectorA = MomentData._vectorA =\r\n      pointA !== undefined ? Point4d.create(pointA.x - x0, pointA.y - y0, 0.0, 1.0, MomentData._vectorA)\r\n        : Point4d.create(this.origin.x, this.origin.y, 0.0, 1.0, MomentData._vectorA);\r\n    const vectorB = MomentData._vectorB = Point4d.create(pointB.x - x0, pointB.y - y0, 0.0, 1.0, MomentData._vectorB);\r\n    const vectorC = MomentData._vectorC = Point4d.create(pointC.x - x0, pointC.y - y0, 0.0, 1.0, MomentData._vectorC);\r\n\r\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\r\n    //  where w = 1-u-v\r\n    //  W = column vector (point00.x, point00.y, point00.z, 1.0) etc.\r\n    const detJ = Geometry.crossProductXYXY(vectorB.x - vectorA.x, vectorB.y - vectorA.y, vectorC.x - vectorA.x, vectorC.y - vectorA.y);\r\n    if (detJ !== 0.0) {\r\n      const r1_12 = detJ / 12.0;\r\n      const r1_24 = detJ / 24.0;\r\n\r\n      this.sums.addScaledOuterProductInPlace(vectorA, vectorA, r1_12);\r\n      this.sums.addScaledOuterProductInPlace(vectorA, vectorB, r1_24);\r\n      this.sums.addScaledOuterProductInPlace(vectorA, vectorC, r1_24);\r\n\r\n      this.sums.addScaledOuterProductInPlace(vectorB, vectorA, r1_24);\r\n      this.sums.addScaledOuterProductInPlace(vectorB, vectorB, r1_12);\r\n      this.sums.addScaledOuterProductInPlace(vectorB, vectorC, r1_24);\r\n\r\n      this.sums.addScaledOuterProductInPlace(vectorC, vectorA, r1_24);\r\n      this.sums.addScaledOuterProductInPlace(vectorC, vectorB, r1_24);\r\n      this.sums.addScaledOuterProductInPlace(vectorC, vectorC, r1_12);\r\n    }\r\n  }\r\n  /** add scaled outer product of (4d, unit weight) point to this.sums */\r\n  public accumulateScaledOuterProduct(point: XYAndZ, scaleFactor: number) {\r\n    this.setOriginXYZIfNeeded(point.x, point.y, 0.0);\r\n    const vectorA = MomentData._vectorA = Point4d.create(point.x - this.origin.x, point.y - this.origin.y, point.z - this.origin.z, 1.0, MomentData._vectorA);\r\n    this.sums.addScaledOuterProductInPlace(vectorA, vectorA, scaleFactor);\r\n  }\r\n  /** Accumulate wire moment integral from pointA to pointB */\r\n  public accumulateLineMomentsXYZ(pointA: Point3d, pointB: Point3d) {\r\n    this.setOriginXYZIfNeeded(pointA.x, pointA.y, pointA.z);\r\n    const x0 = this.origin.x;\r\n    const y0 = this.origin.y;\r\n    const z0 = this.origin.z;\r\n    const vectorA = MomentData._vectorA = Point4d.create(pointA.x - x0, pointA.y - y0, pointA.z - z0, 1.0, MomentData._vectorA);\r\n    const vectorB = MomentData._vectorB = Point4d.create(pointB.x - x0, pointB.y - y0, pointB.z - z0, 1.0, MomentData._vectorB);\r\n    const detJ = pointA.distance(pointB);\r\n    const r1_3 = detJ / 3.0;\r\n    const r1_6 = detJ / 6.0;\r\n    this.sums.addScaledOuterProductInPlace(vectorA, vectorA, r1_3);\r\n    this.sums.addScaledOuterProductInPlace(vectorA, vectorB, r1_6);\r\n    this.sums.addScaledOuterProductInPlace(vectorB, vectorA, r1_6);\r\n    this.sums.addScaledOuterProductInPlace(vectorB, vectorB, r1_3);\r\n\r\n  }\r\n\r\n  private _point0 = Point3d.create();\r\n  private _point1 = Point3d.create();\r\n  /** compute moments of triangles from a base point to the given linestring.\r\n   * Accumulate them to this.sums.\r\n   * * If `pointA` is undefined, use `this.origin` as pointA.\r\n   * * If `this.needOrigin` is set, the first point of the array is captured as local origin for subsequent sums.\r\n   *\r\n   */\r\n  public accumulateTriangleToLineStringMomentsXY(sweepBase: XAndY | undefined, points: GrowableXYZArray) {\r\n    const n = points.length;\r\n    if (n > 1) {\r\n      points.getPoint3dAtUncheckedPointIndex(0, this._point0);\r\n      for (let i = 1; i < n; i++) {\r\n        points.getPoint3dAtUncheckedPointIndex(i, this._point1);\r\n        this.accumulateTriangleMomentsXY(sweepBase, this._point0, this._point1);\r\n        this._point0.setFromPoint3d(this._point1);\r\n      }\r\n    }\r\n  }\r\n  // cspell:word ABAT\r\n  /**\r\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].\r\n   * * Sandwich this between transforms with columns [vectorU, vectorV, 0000, origin].  (Column weights 0001) (only xy parts of vectors)\r\n   * * scale by detJ for the xy-only determinant of the vectors.\r\n   * @param productXX\r\n   * @param productXY\r\n   * @param productYY\r\n   * @param area Area in caller's system\r\n   * @param origin Caller's origin\r\n   * @param vectorU Caller's U axis (not necessarily unit)\r\n   * @param vectorV Caller's V axis (not necessarily unit)\r\n   */\r\n  public accumulateXYProductsInCentroidalFrame(productXX: number, productXY: number, productYY: number, area: number,\r\n    origin: XAndY, vectorU: XAndY, vectorV: XAndY) {\r\n    const centroidalProducts = Matrix4d.createRowValues(\r\n      productXX, productXY, 0, 0,\r\n      productXY, productYY, 0, 0,\r\n      0, 0, 0, 0,\r\n      0, 0, 0, area);\r\n    const detJ = Geometry.crossProductXYXY(vectorU.x, vectorV.x, vectorU.y, vectorV.y);\r\n    const placement = Matrix4d.createRowValues(\r\n      vectorU.x, vectorV.x, 0, origin.x - this.origin.x,\r\n      vectorU.y, vectorV.y, 0, origin.y - this.origin.y,\r\n      0, 0, 0, 0,\r\n      0, 0, 0, 1);\r\n    const AB = placement.multiplyMatrixMatrix(centroidalProducts);\r\n    const ABAT = AB.multiplyMatrixMatrixTranspose(placement);\r\n    this.sums.addScaledInPlace(ABAT, detJ);\r\n  }\r\n  /**\r\n   * Accumulate sums from other moments.\r\n   * * scale by given scaleFactor (e.g. sign to correct orientation)\r\n   * * pull the origin from `other` if `this` needs an origin.\r\n   * *\r\n   */\r\n  public accumulateProducts(other: MomentData, scale: number) {\r\n    this.setOriginIfNeeded(other.origin);\r\n    this.sums.addTranslationSandwichInPlace(other.sums, this.origin.x - other.origin.x, this.origin.y - other.origin.y, this.origin.z - other.origin.z, scale);\r\n  }\r\n\r\n  /**\r\n * Accumulate sums from Matrix4d and origin.\r\n * * scale by given scaleFactor (e.g. sign to correct orientation)\r\n * * trap the origin if `this` needs an origin.\r\n * *\r\n */\r\n  public accumulateProductsFromOrigin(origin: Point3d, products: Matrix4d, scale: number) {\r\n    this.setOriginIfNeeded(origin);\r\n    this.sums.addTranslationSandwichInPlace(products, this.origin.x - origin.x, this.origin.y - origin.y, this.origin.z - origin.z, scale);\r\n  }\r\n  /**\r\n   * Convert to a json data object with:\r\n   */\r\n  public toJSON(): any {\r\n    return {\r\n      origin: this.origin,\r\n      sums: this.sums.toJSON(),\r\n      radiusOfGyration: this.radiusOfGyration.toJSON(),\r\n      localToWorld: this.localToWorldMap.toJSON(),\r\n    };\r\n  }\r\n}\r\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,OAAO,EAAE,QAAQ,EAAE,MAAM,aAAa,CAAC;AAEvC,OAAO,EAAE,QAAQ,EAAE,MAAM,wBAAwB,CAAC;AAClD,OAAO,EAAE,OAAO,EAAE,QAAQ,EAAO,MAAM,+BAA+B,CAAC;AACvE,OAAO,EAAE,SAAS,EAAE,MAAM,yBAAyB,CAAC;AAEpD,OAAO,EAAE,QAAQ,EAAE,MAAM,YAAY,CAAC;AACtC,OAAO,EAAE,OAAO,EAAE,MAAM,WAAW,CAAC;AAEpC;;;;;;;;;;;;;;;;GAgBG;AACH,MAAM,OAAO,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,OAAO,CAAC,MAAM,EAAE,CAAC;IAC3B,OAAO,GAAG,OAAO,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,OAAO,CAAC,UAAU,EAAE,CAAC;QACnC,IAAI,CAAC,IAAI,GAAG,QAAQ,CAAC,UAAU,EAAE,CAAC;QAClC,IAAI,CAAC,eAAe,GAAG,SAAS,CAAC,cAAc,EAAE,CAAC;QAClD,IAAI,CAAC,gBAAgB,GAAG,QAAQ,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,QAAQ,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,QAAQ,CAAC,MAAM,EAAE,CAAC;QAClC,MAAM,WAAW,GAAG,QAAQ,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,SAAS,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,QAAQ,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,QAAQ,CAAC,gBAAgB,CAAC,KAAK,CAAC,gBAAgB,CAAC,CAAC,EAAE,KAAK,CAAC,gBAAgB,CAAC,CAAC,CAAC,EAAE,CAAC;oBAClF,oCAAoC;oBACpC,IAAI,QAAQ,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,QAAQ,CAAC,MAAM,EAAE,CAAC;gBAClC,MAAM,OAAO,GAAG,QAAQ,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,OAAO,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,OAAO,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,OAAO,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,OAAO,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,QAAQ,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,OAAO,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,OAAO,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,OAAO,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,QAAQ,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,QAAQ,CAAC,gBAAgB,CAAC,OAAO,CAAC,CAAC,EAAE,OAAO,CAAC,CAAC,EAAE,OAAO,CAAC,CAAC,EAAE,OAAO,CAAC,CAAC,CAAC,CAAC;QACnF,MAAM,SAAS,GAAG,QAAQ,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","sourcesContent":["/*---------------------------------------------------------------------------------------------\r\n* Copyright (c) Bentley Systems, Incorporated. All rights reserved.\r\n* See LICENSE.md in the project root for license terms and full copyright notice.\r\n*--------------------------------------------------------------------------------------------*/\r\n\r\n/** @packageDocumentation\r\n * @module Numerics\r\n */\r\n\r\n/* eslint-disable @typescript-eslint/naming-convention */\r\n// cspell:word ABAT\r\n\r\nimport { Geometry } from \"../Geometry\";\r\nimport { GrowableXYZArray } from \"../geometry3d/GrowableXYZArray\";\r\nimport { Matrix3d } from \"../geometry3d/Matrix3d\";\r\nimport { Point3d, Vector3d, XYZ } from \"../geometry3d/Point3dVector3d\";\r\nimport { Transform } from \"../geometry3d/Transform\";\r\nimport { XAndY, XYAndZ } from \"../geometry3d/XYZProps\";\r\nimport { Matrix4d } from \"./Matrix4d\";\r\nimport { Point4d } from \"./Point4d\";\r\n\r\n/**\r\n * A MomentData structure exists in several levels:\r\n * * First level: as a carrier of sums of inertial products that determine moments.\r\n *   * origin = local origin used as moments are summed.\r\n *   * sums = array of summed moments.\r\n *     * The [i,j] entry of the sums is a summed or integrated moment for product of axis i and j.\r\n *       * axes 0,1,2 are x,y,z.\r\n *         * e.g. entry [0,1] is summed product xy\r\n *       * axis 3 is \"w\", which is 1 in sums.\r\n *         * e.g. entry 03 is summed x.\r\n *    * In this level:\r\n *        * the `absoluteQuantity` member is undefined.\r\n *        * the `localToWorldMap` and `radiiOfGyration` are created by have undefined contents.\r\n *  * Second level: after a call to inertiaProductsToPrincipalAxes, the `localToWorldMap`, `absoluteQuantity` and\r\n *    `radiiOfGyration` are filled in.\r\n * @public\r\n */\r\nexport class MomentData {\r\n  /** Origin used for sums. */\r\n  public origin: Point3d;\r\n  /** Flag to request deferred origin setup. */\r\n  public needOrigin: boolean;\r\n  /**\r\n   * Moment sums.\r\n   * * Set to zero at initialization and if requested later.\r\n   * * Accumulated during data entry phase.\r\n   */\r\n  public sums: Matrix4d;\r\n  /**\r\n   * The mapping between principal and world system.\r\n   * * This set up with its inverse already constructed.\r\n   */\r\n  public localToWorldMap: Transform;\r\n  // private variables\r\n  private static _vectorA?: Point4d;\r\n  private static _vectorB?: Point4d;\r\n  private static _vectorC?: Point4d;\r\n  private _point0 = Point3d.create();\r\n  private _point1 = Point3d.create();\r\n  /**\r\n   * Return the lower-right (3,3) entry in the sums.\r\n   * * This is the quantity (i.e. length, area, or volume) summed.\r\n   */\r\n  public get quantitySum(): number {\r\n    return this.sums.atIJ(3, 3);\r\n  }\r\n  /**\r\n   * Return a scale factor to make these sums match the target orientation sign.\r\n   * * 1.0 if `this.quantitySum` has the same sign as `targetSign`.\r\n   * * -1.0 if `this.quantitySum` has the opposite sign from `targetSign`\r\n   */\r\n  public signFactor(targetSign: number): number {\r\n    return targetSign * this.quantitySum > 0 ? 1.0 : -1.0;\r\n  }\r\n  /** If `this.needOrigin` flag is set, copy `origin` to `this.origin` and clear the flag. */\r\n  public setOriginIfNeeded(origin: Point3d) {\r\n    if (this.needOrigin) {\r\n      this.origin.setFromPoint3d(origin);\r\n      this.needOrigin = false;\r\n    }\r\n  }\r\n  /** If `this.needOrigin` flag is set, copy `origin` to `this.origin` and clear the flag. */\r\n  public setOriginFromGrowableXYZArrayIfNeeded(points: GrowableXYZArray) {\r\n    if (this.needOrigin && points.length > 0) {\r\n      points.getPoint3dAtCheckedPointIndex(0, this.origin);\r\n      this.needOrigin = false;\r\n    }\r\n  }\r\n  /** If `this.needOrigin` flag is set, copy `origin` to `this.origin` and clear the flag. */\r\n  public setOriginXYZIfNeeded(x: number, y: number, z: number) {\r\n    if (this.needOrigin) {\r\n      this.origin.set(x, y, z);\r\n      this.needOrigin = false;\r\n    }\r\n  }\r\n  /** Radii of gyration (square roots of principal second moments). */\r\n  public radiusOfGyration: Vector3d;\r\n  /**\r\n   * Principal quantity (e.g. length, area, or volume). This is undefined in raw moments, and becomes defined by\r\n   * inertiaProductsToPrincipalAxes.\r\n   */\r\n  public absoluteQuantity?: number;\r\n  /** Constructor. */\r\n  private constructor() {\r\n    this.origin = Point3d.createZero();\r\n    this.sums = Matrix4d.createZero();\r\n    this.localToWorldMap = Transform.createIdentity();\r\n    this.radiusOfGyration = Vector3d.create();\r\n    this.needOrigin = false;\r\n    this.absoluteQuantity = 0.1; // so optimizer sees its type\r\n    this.absoluteQuantity = undefined;\r\n  }\r\n  /**\r\n   * Create moments with optional origin.\r\n   * * Origin and needOrigin are quirky.\r\n   *   * (undefined, true) sets up to use first incoming point as origin.\r\n   *   * (origin) sets up to use that durable origin, set needsOrigin flag false.\r\n   *   * (origin, true) the \"true\" is meaningless.\r\n   *   * (undefined, false) makes 000 the durable origin.\r\n  */\r\n  public static create(origin?: Point3d | undefined, needOrigin: boolean = false): MomentData {\r\n    const data = new MomentData();\r\n    data.needOrigin = needOrigin;\r\n    if (origin) {\r\n      data.origin.setFromPoint3d(origin);\r\n      data.needOrigin = false;\r\n    }\r\n    return data;\r\n  }\r\n  /**\r\n   * Return the formal tensor of integrated values `[yy+zz,xy,xz][yx,xx+zz,yz][zx,xy,xx+yy]`.\r\n   * @param products matrix of (integrated) `[xx,xy,xz][yx,yy,yz][zx,xy,zz]`.\r\n   */\r\n  public static momentTensorFromInertiaProducts(products: Matrix3d): Matrix3d {\r\n    const rr = products.sumDiagonal();\r\n    const result = Matrix3d.createScale(rr, rr, rr);\r\n    result.addScaledInPlace(products, -1.0);\r\n    return result;\r\n  }\r\n  /** Sort the columns of the matrix for increasing moments. */\r\n  public static sortColumnsForIncreasingMoments(axes: Matrix3d, moments: Vector3d) {\r\n    const points = [\r\n      axes.indexedColumnWithWeight(0, moments.x),\r\n      axes.indexedColumnWithWeight(1, moments.y),\r\n      axes.indexedColumnWithWeight(2, moments.z),\r\n    ].sort((dataA: Point4d, dataB: Point4d): number => {\r\n      if (dataA.w < dataB.w) return -1;\r\n      if (dataA.w > dataB.w) return 1;\r\n      return 0;\r\n    });\r\n    axes.setColumnsPoint4dXYZ(points[0], points[1], points[2]);\r\n    if (axes.determinant() < 0)\r\n      axes.scaleColumnsInPlace(-1.0, -1.0, -1.0);\r\n    // prefer x and z positive -- y falls wherever . ..\r\n    if (axes.at(0, 0) < 0.0)\r\n      axes.scaleColumnsInPlace(-1.0, -1.0, 1.0);\r\n    if (axes.at(2, 2) < 0.0)\r\n      axes.scaleColumnsInPlace(1.0, -1.0, -1.0);\r\n    moments.set(points[0].w, points[1].w, points[2].w);\r\n  }\r\n  /**\r\n   * Return the principal moment data for an array of points.\r\n   * @param points array of points.\r\n   */\r\n  public static pointsToPrincipalAxes(points: Point3d[]): MomentData | undefined {\r\n    const moments = new MomentData();\r\n    if (points.length === 0)\r\n      return moments;\r\n    moments.clearSums(points[0]);\r\n    moments.accumulatePointMomentsFromOrigin(points);\r\n    return this.inertiaProductsToPrincipalAxes(moments.origin, moments.sums);\r\n  }\r\n  /**\r\n   * Compute principal axes from inertial products.\r\n   * * The radii of gyration are sorted smallest to largest.\r\n   * * Hence x axis is long direction.\r\n   * * Hence planar data generates large moment as Z.\r\n   * @param origin The origin used for the inertia products.\r\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].\r\n   */\r\n  public static inertiaProductsToPrincipalAxes(origin: XYZ, inertiaProducts: Matrix4d): MomentData | undefined {\r\n    const moments = new MomentData();\r\n    moments.sums.setFrom(inertiaProducts);\r\n    moments.origin.setFrom(origin);\r\n    if (!moments.shiftOriginAndSumsToCentroidOfSums())\r\n      return undefined;\r\n    const products = moments.sums.matrixPart();\r\n    const w = moments.sums.weight();\r\n    if (w < 0.0)\r\n      products.scaleColumnsInPlace(-1, -1, -1);\r\n    const tensor = MomentData.momentTensorFromInertiaProducts(products);\r\n    const moment2 = Vector3d.create();\r\n    const axisVectors = Matrix3d.createZero();\r\n    tensor.fastSymmetricEigenvalues(axisVectors, moment2);\r\n    if (moment2.x < 0.0)\r\n      return undefined;\r\n    MomentData.sortColumnsForIncreasingMoments(axisVectors, moment2);\r\n    if (w < 0.0)\r\n      axisVectors.scaleColumnsInPlace(1, -1, -1);\r\n    moments.localToWorldMap = Transform.createOriginAndMatrix(moments.origin, axisVectors);\r\n    moments.radiusOfGyration.set(\r\n      Math.sqrt(Math.abs(moment2.x)), Math.sqrt(Math.abs(moment2.y)), Math.sqrt(Math.abs(moment2.z)));\r\n    moments.radiusOfGyration.scaleInPlace(1.0 / Math.sqrt(Math.abs(w)));\r\n    moments.absoluteQuantity = Math.abs(w);\r\n    return moments;\r\n  }\r\n  /**\r\n   * Test for match among selected members as they exist after `inertiaProductsToPrincipalAxes`.\r\n   * * The members considered are:\r\n   *   * origin of local to world map (i.e. centroid),\r\n   *   * radius of gyration,\r\n   *   * axes of localToWorldMap.\r\n   * * Axis direction tests allow these quirks:\r\n   *   * opposite orientation is considered matched.\r\n   * * Full xyz symmetry: If x,y,z radii are matched, axes are not tested.\r\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\r\n   *   is tested.\r\n   * * If either or both are undefined, returns false.\r\n   * @param dataA first set of moments.\r\n   * @param dataB second set of moments.\r\n   */\r\n  public static areEquivalentPrincipalAxes(dataA: MomentData | undefined, dataB: MomentData | undefined): boolean {\r\n    if (dataA && dataB\r\n      && Geometry.isSameCoordinate(dataA.quantitySum, dataB.quantitySum)) {  // um.. need different tolerance for area, volume?)\r\n      if (dataA.localToWorldMap.getOrigin().isAlmostEqual(dataB.localToWorldMap.getOrigin())\r\n        && dataA.radiusOfGyration.isAlmostEqual(dataB.radiusOfGyration)) {\r\n        if (Geometry.isSameCoordinate(dataA.radiusOfGyration.x, dataA.radiusOfGyration.y)) {\r\n          // We have at least xy symmetry ....\r\n          if (Geometry.isSameCoordinate(dataA.radiusOfGyration.x, dataA.radiusOfGyration.z))\r\n            return true;\r\n          // just xy.\r\n          // allow opposite z directions.\r\n          // If the z's are aligned, x an dy can spin freely.\r\n          const zA = dataA.localToWorldMap.matrix.columnZ();\r\n          const zB = dataB.localToWorldMap.matrix.columnZ();\r\n          if (zA.isParallelTo(zB, true))\r\n            return true;\r\n          return false;\r\n        }\r\n        // no symmetry.  Test all three axes.\r\n        const vectorA = Vector3d.create();\r\n        const vectorB = Vector3d.create();\r\n        for (let i = 0; i < 3; i++) {\r\n          dataA.localToWorldMap.matrix.getColumn(i, vectorA);\r\n          dataB.localToWorldMap.matrix.getColumn(i, vectorB);\r\n          if (!vectorA.isParallelTo(vectorB, true))\r\n            return false;\r\n        }\r\n        return true;\r\n      }\r\n    }\r\n    return false;\r\n  }\r\n  /** Clear the MomentData sums to zero, and establish a new origin. */\r\n  public clearSums(origin?: Point3d) {\r\n    this.sums.setZero();\r\n    if (origin)\r\n      this.origin.setFrom(origin);\r\n    else\r\n      this.origin.setZero();\r\n  }\r\n  /** Accumulate products-of-components for given points. */\r\n  public accumulatePointMomentsFromOrigin(points: Point3d[]) {\r\n    for (const p of points) {\r\n      this.sums.addMomentsInPlace(\r\n        p.x - this.origin.x,\r\n        p.y - this.origin.y,\r\n        p.z - this.origin.z,\r\n        1.0,\r\n      );\r\n    }\r\n  }\r\n  /** Revise the accumulated sums to be \"around the centroid\". */\r\n  public shiftOriginAndSumsToCentroidOfSums(): boolean {\r\n    const xyz = this.sums.columnW().realPoint();\r\n    if (xyz) {\r\n      this.shiftOriginAndSumsByXYZ(xyz.x, xyz.y, xyz.z);\r\n      return true;\r\n    }\r\n    return false;\r\n  }\r\n\r\n  /**\r\n   * Revise the accumulated sums.\r\n   * * add ax,ay,ax to the origin coordinates.\r\n   * * apply the negative translation to the sums.\r\n  */\r\n  public shiftOriginAndSumsByXYZ(ax: number, ay: number, az: number) {\r\n    this.origin.addXYZInPlace(ax, ay, az);\r\n    this.sums.multiplyTranslationSandwichInPlace(-ax, -ay, -az);\r\n  }\r\n  /** Revise the accumulated sums so they are based at a specified origin. */\r\n  public shiftOriginAndSumsToNewOrigin(newOrigin: XYAndZ) {\r\n    this.shiftOriginAndSumsByXYZ(newOrigin.x - this.origin.x, newOrigin.y - this.origin.y, newOrigin.z - this.origin.z);\r\n  }\r\n  /**\r\n   * Compute moments of a triangle from the origin to the given line.\r\n   * Accumulate them to this.sums.\r\n   * * If `pointA` is undefined, use `this.origin` as pointA.\r\n   * * If `this.needOrigin` is set, pointB is used\r\n  */\r\n  public accumulateTriangleMomentsXY(pointA: XAndY | undefined, pointB: XAndY, pointC: XAndY) {\r\n    this.setOriginXYZIfNeeded(pointB.x, pointB.y, 0.0);\r\n    const x0 = this.origin.x;\r\n    const y0 = this.origin.y;\r\n    const vectorA = MomentData._vectorA =\r\n      pointA !== undefined ? Point4d.create(pointA.x - x0, pointA.y - y0, 0.0, 1.0, MomentData._vectorA)\r\n        : Point4d.create(this.origin.x, this.origin.y, 0.0, 1.0, MomentData._vectorA);\r\n    const vectorB = MomentData._vectorB = Point4d.create(pointB.x - x0, pointB.y - y0, 0.0, 1.0, MomentData._vectorB);\r\n    const vectorC = MomentData._vectorC = Point4d.create(pointC.x - x0, pointC.y - y0, 0.0, 1.0, MomentData._vectorC);\r\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\r\n    //  where w = 1-u-v\r\n    //  W = column vector (point00.x, point00.y, point00.z, 1.0) etc.\r\n    const detJ = Geometry.crossProductXYXY(vectorB.x - vectorA.x, vectorB.y - vectorA.y, vectorC.x - vectorA.x, vectorC.y - vectorA.y);\r\n    if (detJ !== 0.0) {\r\n      const r1_12 = detJ / 12.0;\r\n      const r1_24 = detJ / 24.0;\r\n\r\n      this.sums.addScaledOuterProductInPlace(vectorA, vectorA, r1_12);\r\n      this.sums.addScaledOuterProductInPlace(vectorA, vectorB, r1_24);\r\n      this.sums.addScaledOuterProductInPlace(vectorA, vectorC, r1_24);\r\n\r\n      this.sums.addScaledOuterProductInPlace(vectorB, vectorA, r1_24);\r\n      this.sums.addScaledOuterProductInPlace(vectorB, vectorB, r1_12);\r\n      this.sums.addScaledOuterProductInPlace(vectorB, vectorC, r1_24);\r\n\r\n      this.sums.addScaledOuterProductInPlace(vectorC, vectorA, r1_24);\r\n      this.sums.addScaledOuterProductInPlace(vectorC, vectorB, r1_24);\r\n      this.sums.addScaledOuterProductInPlace(vectorC, vectorC, r1_12);\r\n    }\r\n  }\r\n  /** Add scaled outer product of (4d, unit weight) point to this.sums. */\r\n  public accumulateScaledOuterProduct(point: XYAndZ, scaleFactor: number) {\r\n    this.setOriginXYZIfNeeded(point.x, point.y, 0.0);\r\n    const vectorA = MomentData._vectorA = Point4d.create(\r\n      point.x - this.origin.x, point.y - this.origin.y, point.z - this.origin.z, 1.0, MomentData._vectorA,\r\n    );\r\n    this.sums.addScaledOuterProductInPlace(vectorA, vectorA, scaleFactor);\r\n  }\r\n  /** Accumulate wire moment integral from pointA to pointB */\r\n  public accumulateLineMomentsXYZ(pointA: Point3d, pointB: Point3d) {\r\n    this.setOriginXYZIfNeeded(pointA.x, pointA.y, pointA.z);\r\n    const x0 = this.origin.x;\r\n    const y0 = this.origin.y;\r\n    const z0 = this.origin.z;\r\n    const vectorA = MomentData._vectorA = Point4d.create(\r\n      pointA.x - x0, pointA.y - y0, pointA.z - z0, 1.0, MomentData._vectorA,\r\n    );\r\n    const vectorB = MomentData._vectorB = Point4d.create(\r\n      pointB.x - x0, pointB.y - y0, pointB.z - z0, 1.0, MomentData._vectorB,\r\n    );\r\n    const detJ = pointA.distance(pointB);\r\n    const r1_3 = detJ / 3.0;\r\n    const r1_6 = detJ / 6.0;\r\n    this.sums.addScaledOuterProductInPlace(vectorA, vectorA, r1_3);\r\n    this.sums.addScaledOuterProductInPlace(vectorA, vectorB, r1_6);\r\n    this.sums.addScaledOuterProductInPlace(vectorB, vectorA, r1_6);\r\n    this.sums.addScaledOuterProductInPlace(vectorB, vectorB, r1_3);\r\n  }\r\n  /**\r\n   * Compute moments of triangles from a base point to the given linestring.\r\n   * Accumulate them to this.sums.\r\n   * * If `pointA` is undefined, use `this.origin` as pointA.\r\n   * * If `this.needOrigin` is set, the first point of the array is captured as local origin for subsequent sums.\r\n   *\r\n   */\r\n  public accumulateTriangleToLineStringMomentsXY(sweepBase: XAndY | undefined, points: GrowableXYZArray) {\r\n    const n = points.length;\r\n    if (n > 1) {\r\n      points.getPoint3dAtUncheckedPointIndex(0, this._point0);\r\n      for (let i = 1; i < n; i++) {\r\n        points.getPoint3dAtUncheckedPointIndex(i, this._point1);\r\n        this.accumulateTriangleMomentsXY(sweepBase, this._point0, this._point1);\r\n        this._point0.setFromPoint3d(this._point1);\r\n      }\r\n    }\r\n  }\r\n  /**\r\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].\r\n   * * Sandwich this between transforms with columns [vectorU, vectorV, 0000, origin].  (Column weights 0001) (only xy\r\n   *   parts of vectors).\r\n   * * scale by detJ for the xy-only determinant of the vectors.\r\n   * @param productXX\r\n   * @param productXY\r\n   * @param productYY\r\n   * @param area Area in caller's system.\r\n   * @param origin Caller's origin.\r\n   * @param vectorU Caller's U axis (not necessarily unit).\r\n   * @param vectorV Caller's V axis (not necessarily unit).\r\n   */\r\n  public accumulateXYProductsInCentroidalFrame(\r\n    productXX: number, productXY: number, productYY: number,\r\n    area: number, origin: XAndY, vectorU: XAndY, vectorV: XAndY,\r\n  ) {\r\n    const centroidalProducts = Matrix4d.createRowValues(\r\n      productXX, productXY, 0, 0,\r\n      productXY, productYY, 0, 0,\r\n      0, 0, 0, 0,\r\n      0, 0, 0, area,\r\n    );\r\n    const detJ = Geometry.crossProductXYXY(vectorU.x, vectorV.x, vectorU.y, vectorV.y);\r\n    const placement = Matrix4d.createRowValues(\r\n      vectorU.x, vectorV.x, 0, origin.x - this.origin.x,\r\n      vectorU.y, vectorV.y, 0, origin.y - this.origin.y,\r\n      0, 0, 0, 0,\r\n      0, 0, 0, 1,\r\n    );\r\n    const AB = placement.multiplyMatrixMatrix(centroidalProducts);\r\n    const ABAT = AB.multiplyMatrixMatrixTranspose(placement);\r\n    this.sums.addScaledInPlace(ABAT, detJ);\r\n  }\r\n  /**\r\n   * Accumulate sums from other moments.\r\n   * * Scale by given scaleFactor (e.g. sign to correct orientation).\r\n   * * Pull the origin from `other` if `this` needs an origin.\r\n   */\r\n  public accumulateProducts(other: MomentData, scale: number) {\r\n    this.setOriginIfNeeded(other.origin);\r\n    this.sums.addTranslationSandwichInPlace(\r\n      other.sums, this.origin.x - other.origin.x, this.origin.y - other.origin.y, this.origin.z - other.origin.z, scale,\r\n    );\r\n  }\r\n\r\n  /**\r\n   * Accumulate sums from Matrix4d and origin.\r\n   * * Scale by given scaleFactor (e.g. sign to correct orientation).\r\n   * * Trap the origin if `this` needs an origin.\r\n   */\r\n  public accumulateProductsFromOrigin(origin: Point3d, products: Matrix4d, scale: number) {\r\n    this.setOriginIfNeeded(origin);\r\n    this.sums.addTranslationSandwichInPlace(\r\n      products, this.origin.x - origin.x, this.origin.y - origin.y, this.origin.z - origin.z, scale,\r\n    );\r\n  }\r\n  /** Convert to a json data object. */\r\n  public toJSON(): any {\r\n    return {\r\n      origin: this.origin,\r\n      sums: this.sums.toJSON(),\r\n      radiusOfGyration: this.radiusOfGyration.toJSON(),\r\n      localToWorld: this.localToWorldMap.toJSON(),\r\n    };\r\n  }\r\n}\r\n"]}

package/lib/esm/geometry4d/PlaneByOriginAndVectors4d.js

@@ -11,6 +11,12 @@ import { Point4d } from "./Point4d";
  * @public
  */
 export class PlaneByOriginAndVectors4d {
+    /** homogeneous origin */
+    origin;
+    /** homogeneous u-direction vector */
+    vectorU;
+    /** homogeneous v-direction vector */
+    vectorV;
     constructor(origin, vectorU, vectorV) {
         this.origin = origin;
         this.vectorU = vectorU;

package/lib/esm/geometry4d/PlaneByOriginAndVectors4d.js.map

@@ -1 +1 @@
-{"version":3,"file":"PlaneByOriginAndVectors4d.js","sourceRoot":"","sources":["../../../src/geometry4d/PlaneByOriginAndVectors4d.ts"],"names":[],"mappings":"AAAA;;;+FAG+F;AAO/F,OAAO,EAAE,OAAO,EAAE,MAAM,WAAW,CAAC;AAEpC;;;;;;GAMG;AACH,MAAM,OAAO,yBAAyB;IAOpC,YAAoB,MAAe,EAAE,OAAgB,EAAE,OAAgB;QACrE,IAAI,CAAC,MAAM,GAAG,MAAM,CAAC;QACrB,IAAI,CAAC,OAAO,GAAG,OAAO,CAAC;QACvB,IAAI,CAAC,OAAO,GAAG,OAAO,CAAC;IACzB,CAAC;IACD,mCAAmC;IAC5B,KAAK,CAAC,MAAkC;QAC7C,IAAI,MAAM,EAAE,CAAC;YACX,MAAM,CAAC,OAAO,CAAC,IAAI,CAAC,CAAC;YACrB,OAAO,MAAM,CAAC;QAChB,CAAC;QACD,OAAO,IAAI,yBAAyB,CAAC,IAAI,CAAC,MAAM,CAAC,KAAK,EAAE,EAAE,IAAI,CAAC,OAAO,CAAC,KAAK,EAAE,EAAE,IAAI,CAAC,OAAO,CAAC,KAAK,EAAE,CAAC,CAAC;IACxG,CAAC;IACD,wCAAwC;IACjC,OAAO,CAAC,KAAgC;QAC7C,IAAI,CAAC,MAAM,CAAC,OAAO,CAAC,KAAK,CAAC,MAAM,CAAC,CAAC;QAClC,IAAI,CAAC,OAAO,CAAC,OAAO,CAAC,KAAK,CAAC,OAAO,CAAC,CAAC;QACpC,IAAI,CAAC,OAAO,CAAC,OAAO,CAAC,KAAK,CAAC,OAAO,CAAC,CAAC;IACtC,CAAC;IACD,sEAAsE;IAC/D,aAAa,CAAC,KAAgC;QACnD,OAAO,IAAI,CAAC,MAAM,CAAC,aAAa,CAAC,KAAK,CAAC,MAAM,CAAC;eACzC,IAAI,CAAC,OAAO,CAAC,aAAa,CAAC,KAAK,CAAC,OAAO,CAAC;eACzC,IAAI,CAAC,OAAO,CAAC,aAAa,CAAC,KAAK,CAAC,OAAO,CAAC,CAAC;IACjD,CAAC;IACD;OACG;IACI,MAAM,CAAC,sBAAsB,CAAC,MAAe,EAAE,OAAgB,EAAE,OAAgB,EAAE,MAAkC;QAC1H,IAAI,MAAM,EAAE,CAAC;YACX,MAAM,CAAC,mBAAmB,CAAC,MAAM,EAAE,OAAO,EAAE,OAAO,CAAC,CAAC;YACrD,OAAO,MAAM,CAAC;QAChB,CAAC;QACD,OAAO,IAAI,yBAAyB,CAAC,MAAM,CAAC,KAAK,EAAE,EAAE,OAAO,CAAC,KAAK,EAAE,EAAE,OAAO,CAAC,KAAK,EAAE,CAAC,CAAC;IACzF,CAAC;IACD,2FAA2F;IACpF,uBAAuB,CAAC,EAAU,EAAE,EAAU,EAAE,EAAU,EAAE,EAAU,EAAE,EAAU,EAAE,EAAU,EAAE,EAAU,EAAE,EAAU,EAAE,EAAU,EAAE,EAAU,EAAE,EAAU,EAAE,EAAU;QAC3K,IAAI,CAAC,MAAM,CAAC,GAAG,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC;QAChC,IAAI,CAAC,OAAO,CAAC,GAAG,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC;QACjC,IAAI,CAAC,OAAO,CAAC,GAAG,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC;QACjC,OAAO,IAAI,CAAC;IACd,CAAC;IACD,8FAA8F;IACvF,mBAAmB,CAAC,MAAe,EAAE,OAAgB,EAAE,OAAgB;QAC5E,IAAI,CAAC,MAAM,CAAC,OAAO,CAAC,MAAM,CAAC,CAAC;QAC5B,IAAI,CAAC,OAAO,CAAC,OAAO,CAAC,OAAO,CAAC,CAAC;QAC9B,IAAI,CAAC,OAAO,CAAC,OAAO,CAAC,OAAO,CAAC,CAAC;QAC9B,OAAO,IAAI,CAAC;IACd,CAAC;IACD,6EAA6E;IACtE,MAAM,CAAC,0BAA0B,CAAC,EAAU,EAAE,EAAU,EAAE,EAAU,EAAE,EAAU,EAAE,EAAU,EAAE,EAAU,EAAE,EAAU,EAAE,EAAU,EAAE,EAAU,EAAE,EAAU,EAAE,EAAU,EAAE,EAAU,EAAE,MAAkC;QACzN,IAAI,MAAM;YACR,OAAO,MAAM,CAAC,uBAAuB,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC;QACxF,OAAO,IAAI,yBAAyB,CAAC,OAAO,CAAC,MAAM,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,EAAE,OAAO,CAAC,MAAM,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,EAAE,OAAO,CAAC,MAAM,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC,CAAC;IACvI,CAAC;IACD,sEAAsE;IAC/D,MAAM,CAAC,wBAAwB,CAAC,MAAe,EAAE,OAAgB,EAAE,OAAgB,EAAE,MAAkC;QAC5H,OAAO,yBAAyB,CAAC,0BAA0B,CAAC,MAAM,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC,EAAE,GAAG,EAAE,OAAO,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,EAAE,OAAO,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,EAAE,OAAO,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,EAAE,GAAG,EAAE,OAAO,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,EAAE,OAAO,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,EAAE,OAAO,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,EAAE,GAAG,EAAE,MAAM,CAAC,CAAC;IACvP,CAAC;IACD,gEAAgE;IACzD,eAAe,CAAC,CAAS,EAAE,CAAS,EAAE,MAAgB;QAC3D,OAAO,IAAI,CAAC,MAAM,CAAC,WAAW,CAAC,IAAI,CAAC,OAAO,EAAE,CAAC,EAAE,IAAI,CAAC,OAAO,EAAE,CAAC,EAAE,MAAM,CAAC,CAAC;IAC3E,CAAC;IACD,+DAA+D;IACxD,MAAM,CAAC,aAAa,CAAC,MAAkC;QAC5D,OAAO,yBAAyB,CAAC,0BAA0B,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,MAAM,CAAC,CAAC;IAC1G,CAAC;CACF","sourcesContent":["/*---------------------------------------------------------------------------------------------\r\n* Copyright (c) Bentley Systems, Incorporated. All rights reserved.\r\n* See LICENSE.md in the project root for license terms and full copyright notice.\r\n*--------------------------------------------------------------------------------------------*/\r\n\r\n/** @packageDocumentation\r\n * @module Numerics\r\n */\r\n\r\nimport { Point3d } from \"../geometry3d/Point3dVector3d\";\r\nimport { Point4d } from \"./Point4d\";\r\n\r\n/**\r\n * A PlaneByOriginAndVectors4d is a 4d origin and pair of 4d \"vectors\" defining a 4d plane.\r\n * * The parameterization of the plane is    `X = origin + vectorU * u + vectorV * v`\r\n * * With particular weight values `origin.w === 1, vectorU.w === 0, vectorV.w === 0` this is like `Plane3dByOriginAndVectors`\r\n * * With other weights, the deweighted xyz coordinates of points on the 4d plane still form a 3d plane.\r\n * @public\r\n */\r\nexport class PlaneByOriginAndVectors4d {\r\n  /** homogeneous origin */\r\n  public origin: Point4d;\r\n  /** homogeneous u-direction vector */\r\n  public vectorU: Point4d;\r\n  /** homogeneous v-direction vector */\r\n  public vectorV: Point4d;\r\n  private constructor(origin: Point4d, vectorU: Point4d, vectorV: Point4d) {\r\n    this.origin = origin;\r\n    this.vectorU = vectorU;\r\n    this.vectorV = vectorV;\r\n  }\r\n  /** Return a clone of this plane */\r\n  public clone(result?: PlaneByOriginAndVectors4d): PlaneByOriginAndVectors4d {\r\n    if (result) {\r\n      result.setFrom(this);\r\n      return result;\r\n    }\r\n    return new PlaneByOriginAndVectors4d(this.origin.clone(), this.vectorU.clone(), this.vectorV.clone());\r\n  }\r\n  /** copy all content from other plane */\r\n  public setFrom(other: PlaneByOriginAndVectors4d): void {\r\n    this.origin.setFrom(other.origin);\r\n    this.vectorU.setFrom(other.vectorU);\r\n    this.vectorV.setFrom(other.vectorV);\r\n  }\r\n  /** Return true if origin, vectorU, and vectorV pass isAlmostEqual. */\r\n  public isAlmostEqual(other: PlaneByOriginAndVectors4d): boolean {\r\n    return this.origin.isAlmostEqual(other.origin)\r\n      && this.vectorU.isAlmostEqual(other.vectorU)\r\n      && this.vectorV.isAlmostEqual(other.vectorV);\r\n  }\r\n  /** Create a plane with (copies of) origin, vectorU, vectorV parameters, all given as full 4d points.\r\n   */\r\n  public static createOriginAndVectors(origin: Point4d, vectorU: Point4d, vectorV: Point4d, result?: PlaneByOriginAndVectors4d): PlaneByOriginAndVectors4d {\r\n    if (result) {\r\n      result.setOriginAndVectors(origin, vectorU, vectorV);\r\n      return result;\r\n    }\r\n    return new PlaneByOriginAndVectors4d(origin.clone(), vectorU.clone(), vectorV.clone());\r\n  }\r\n  /** Set all numeric data from complete list of (x,y,z,w) in origin, vectorU, and vectorV */\r\n  public setOriginAndVectorsXYZW(x0: number, y0: number, z0: number, w0: number, ux: number, uy: number, uz: number, uw: number, vx: number, vy: number, vz: number, vw: number): PlaneByOriginAndVectors4d {\r\n    this.origin.set(x0, y0, z0, w0);\r\n    this.vectorU.set(ux, uy, uz, uw);\r\n    this.vectorV.set(vx, vy, vz, vw);\r\n    return this;\r\n  }\r\n  /** Copy the contents of origin, vectorU, vectorV parameters to respective member variables */\r\n  public setOriginAndVectors(origin: Point4d, vectorU: Point4d, vectorV: Point4d): PlaneByOriginAndVectors4d {\r\n    this.origin.setFrom(origin);\r\n    this.vectorU.setFrom(vectorU);\r\n    this.vectorV.setFrom(vectorV);\r\n    return this;\r\n  }\r\n  /** Create from complete list of (x,y,z,w) in origin, vectorU, and vectorV */\r\n  public static createOriginAndVectorsXYZW(x0: number, y0: number, z0: number, w0: number, ux: number, uy: number, uz: number, uw: number, vx: number, vy: number, vz: number, vw: number, result?: PlaneByOriginAndVectors4d): PlaneByOriginAndVectors4d {\r\n    if (result)\r\n      return result.setOriginAndVectorsXYZW(x0, y0, z0, w0, ux, uy, uz, uw, vx, vy, vz, vw);\r\n    return new PlaneByOriginAndVectors4d(Point4d.create(x0, y0, z0, w0), Point4d.create(ux, uy, uz, uw), Point4d.create(vx, vy, vz, uw));\r\n  }\r\n  /** create from origin point, (u=1,v=0) point, and (u=0,v=1) point. */\r\n  public static createOriginAndTargets3d(origin: Point3d, targetU: Point3d, targetV: Point3d, result?: PlaneByOriginAndVectors4d): PlaneByOriginAndVectors4d {\r\n    return PlaneByOriginAndVectors4d.createOriginAndVectorsXYZW(origin.x, origin.y, origin.z, 1.0, targetU.x - origin.x, targetU.y - origin.y, targetU.z - origin.z, 0.0, targetV.x - origin.x, targetV.y - origin.y, targetV.z - origin.z, 0.0, result);\r\n  }\r\n  /** evaluate plane point (full 3d) at given (u,v) coordinate. */\r\n  public fractionToPoint(u: number, v: number, result?: Point4d): Point4d {\r\n    return this.origin.plus2Scaled(this.vectorU, u, this.vectorV, v, result);\r\n  }\r\n  /** create a new plane which maps to the cartesian xy plane. */\r\n  public static createXYPlane(result?: PlaneByOriginAndVectors4d): PlaneByOriginAndVectors4d {\r\n    return PlaneByOriginAndVectors4d.createOriginAndVectorsXYZW(0, 0, 0, 1, 1, 0, 0, 0, 0, 1, 0, 0, result);\r\n  }\r\n}\r\n"]}
+{"version":3,"file":"PlaneByOriginAndVectors4d.js","sourceRoot":"","sources":["../../../src/geometry4d/PlaneByOriginAndVectors4d.ts"],"names":[],"mappings":"AAAA;;;+FAG+F;AAO/F,OAAO,EAAE,OAAO,EAAE,MAAM,WAAW,CAAC;AAEpC;;;;;;GAMG;AACH,MAAM,OAAO,yBAAyB;IACpC,yBAAyB;IAClB,MAAM,CAAU;IACvB,qCAAqC;IAC9B,OAAO,CAAU;IACxB,qCAAqC;IAC9B,OAAO,CAAU;IACxB,YAAoB,MAAe,EAAE,OAAgB,EAAE,OAAgB;QACrE,IAAI,CAAC,MAAM,GAAG,MAAM,CAAC;QACrB,IAAI,CAAC,OAAO,GAAG,OAAO,CAAC;QACvB,IAAI,CAAC,OAAO,GAAG,OAAO,CAAC;IACzB,CAAC;IACD,mCAAmC;IAC5B,KAAK,CAAC,MAAkC;QAC7C,IAAI,MAAM,EAAE,CAAC;YACX,MAAM,CAAC,OAAO,CAAC,IAAI,CAAC,CAAC;YACrB,OAAO,MAAM,CAAC;QAChB,CAAC;QACD,OAAO,IAAI,yBAAyB,CAAC,IAAI,CAAC,MAAM,CAAC,KAAK,EAAE,EAAE,IAAI,CAAC,OAAO,CAAC,KAAK,EAAE,EAAE,IAAI,CAAC,OAAO,CAAC,KAAK,EAAE,CAAC,CAAC;IACxG,CAAC;IACD,wCAAwC;IACjC,OAAO,CAAC,KAAgC;QAC7C,IAAI,CAAC,MAAM,CAAC,OAAO,CAAC,KAAK,CAAC,MAAM,CAAC,CAAC;QAClC,IAAI,CAAC,OAAO,CAAC,OAAO,CAAC,KAAK,CAAC,OAAO,CAAC,CAAC;QACpC,IAAI,CAAC,OAAO,CAAC,OAAO,CAAC,KAAK,CAAC,OAAO,CAAC,CAAC;IACtC,CAAC;IACD,sEAAsE;IAC/D,aAAa,CAAC,KAAgC;QACnD,OAAO,IAAI,CAAC,MAAM,CAAC,aAAa,CAAC,KAAK,CAAC,MAAM,CAAC;eACzC,IAAI,CAAC,OAAO,CAAC,aAAa,CAAC,KAAK,CAAC,OAAO,CAAC;eACzC,IAAI,CAAC,OAAO,CAAC,aAAa,CAAC,KAAK,CAAC,OAAO,CAAC,CAAC;IACjD,CAAC;IACD;OACG;IACI,MAAM,CAAC,sBAAsB,CAAC,MAAe,EAAE,OAAgB,EAAE,OAAgB,EAAE,MAAkC;QAC1H,IAAI,MAAM,EAAE,CAAC;YACX,MAAM,CAAC,mBAAmB,CAAC,MAAM,EAAE,OAAO,EAAE,OAAO,CAAC,CAAC;YACrD,OAAO,MAAM,CAAC;QAChB,CAAC;QACD,OAAO,IAAI,yBAAyB,CAAC,MAAM,CAAC,KAAK,EAAE,EAAE,OAAO,CAAC,KAAK,EAAE,EAAE,OAAO,CAAC,KAAK,EAAE,CAAC,CAAC;IACzF,CAAC;IACD,2FAA2F;IACpF,uBAAuB,CAAC,EAAU,EAAE,EAAU,EAAE,EAAU,EAAE,EAAU,EAAE,EAAU,EAAE,EAAU,EAAE,EAAU,EAAE,EAAU,EAAE,EAAU,EAAE,EAAU,EAAE,EAAU,EAAE,EAAU;QAC3K,IAAI,CAAC,MAAM,CAAC,GAAG,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC;QAChC,IAAI,CAAC,OAAO,CAAC,GAAG,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC;QACjC,IAAI,CAAC,OAAO,CAAC,GAAG,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC;QACjC,OAAO,IAAI,CAAC;IACd,CAAC;IACD,8FAA8F;IACvF,mBAAmB,CAAC,MAAe,EAAE,OAAgB,EAAE,OAAgB;QAC5E,IAAI,CAAC,MAAM,CAAC,OAAO,CAAC,MAAM,CAAC,CAAC;QAC5B,IAAI,CAAC,OAAO,CAAC,OAAO,CAAC,OAAO,CAAC,CAAC;QAC9B,IAAI,CAAC,OAAO,CAAC,OAAO,CAAC,OAAO,CAAC,CAAC;QAC9B,OAAO,IAAI,CAAC;IACd,CAAC;IACD,6EAA6E;IACtE,MAAM,CAAC,0BAA0B,CAAC,EAAU,EAAE,EAAU,EAAE,EAAU,EAAE,EAAU,EAAE,EAAU,EAAE,EAAU,EAAE,EAAU,EAAE,EAAU,EAAE,EAAU,EAAE,EAAU,EAAE,EAAU,EAAE,EAAU,EAAE,MAAkC;QACzN,IAAI,MAAM;YACR,OAAO,MAAM,CAAC,uBAAuB,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC;QACxF,OAAO,IAAI,yBAAyB,CAAC,OAAO,CAAC,MAAM,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,EAAE,OAAO,CAAC,MAAM,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,EAAE,OAAO,CAAC,MAAM,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC,CAAC;IACvI,CAAC;IACD,sEAAsE;IAC/D,MAAM,CAAC,wBAAwB,CAAC,MAAe,EAAE,OAAgB,EAAE,OAAgB,EAAE,MAAkC;QAC5H,OAAO,yBAAyB,CAAC,0BAA0B,CAAC,MAAM,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC,EAAE,GAAG,EAAE,OAAO,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,EAAE,OAAO,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,EAAE,OAAO,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,EAAE,GAAG,EAAE,OAAO,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,EAAE,OAAO,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,EAAE,OAAO,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,EAAE,GAAG,EAAE,MAAM,CAAC,CAAC;IACvP,CAAC;IACD,gEAAgE;IACzD,eAAe,CAAC,CAAS,EAAE,CAAS,EAAE,MAAgB;QAC3D,OAAO,IAAI,CAAC,MAAM,CAAC,WAAW,CAAC,IAAI,CAAC,OAAO,EAAE,CAAC,EAAE,IAAI,CAAC,OAAO,EAAE,CAAC,EAAE,MAAM,CAAC,CAAC;IAC3E,CAAC;IACD,+DAA+D;IACxD,MAAM,CAAC,aAAa,CAAC,MAAkC;QAC5D,OAAO,yBAAyB,CAAC,0BAA0B,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,MAAM,CAAC,CAAC;IAC1G,CAAC;CACF","sourcesContent":["/*---------------------------------------------------------------------------------------------\r\n* Copyright (c) Bentley Systems, Incorporated. All rights reserved.\r\n* See LICENSE.md in the project root for license terms and full copyright notice.\r\n*--------------------------------------------------------------------------------------------*/\r\n\r\n/** @packageDocumentation\r\n * @module Numerics\r\n */\r\n\r\nimport { Point3d } from \"../geometry3d/Point3dVector3d\";\r\nimport { Point4d } from \"./Point4d\";\r\n\r\n/**\r\n * A PlaneByOriginAndVectors4d is a 4d origin and pair of 4d \"vectors\" defining a 4d plane.\r\n * * The parameterization of the plane is    `X = origin + vectorU * u + vectorV * v`\r\n * * With particular weight values `origin.w === 1, vectorU.w === 0, vectorV.w === 0` this is like `Plane3dByOriginAndVectors`\r\n * * With other weights, the deweighted xyz coordinates of points on the 4d plane still form a 3d plane.\r\n * @public\r\n */\r\nexport class PlaneByOriginAndVectors4d {\r\n  /** homogeneous origin */\r\n  public origin: Point4d;\r\n  /** homogeneous u-direction vector */\r\n  public vectorU: Point4d;\r\n  /** homogeneous v-direction vector */\r\n  public vectorV: Point4d;\r\n  private constructor(origin: Point4d, vectorU: Point4d, vectorV: Point4d) {\r\n    this.origin = origin;\r\n    this.vectorU = vectorU;\r\n    this.vectorV = vectorV;\r\n  }\r\n  /** Return a clone of this plane */\r\n  public clone(result?: PlaneByOriginAndVectors4d): PlaneByOriginAndVectors4d {\r\n    if (result) {\r\n      result.setFrom(this);\r\n      return result;\r\n    }\r\n    return new PlaneByOriginAndVectors4d(this.origin.clone(), this.vectorU.clone(), this.vectorV.clone());\r\n  }\r\n  /** copy all content from other plane */\r\n  public setFrom(other: PlaneByOriginAndVectors4d): void {\r\n    this.origin.setFrom(other.origin);\r\n    this.vectorU.setFrom(other.vectorU);\r\n    this.vectorV.setFrom(other.vectorV);\r\n  }\r\n  /** Return true if origin, vectorU, and vectorV pass isAlmostEqual. */\r\n  public isAlmostEqual(other: PlaneByOriginAndVectors4d): boolean {\r\n    return this.origin.isAlmostEqual(other.origin)\r\n      && this.vectorU.isAlmostEqual(other.vectorU)\r\n      && this.vectorV.isAlmostEqual(other.vectorV);\r\n  }\r\n  /** Create a plane with (copies of) origin, vectorU, vectorV parameters, all given as full 4d points.\r\n   */\r\n  public static createOriginAndVectors(origin: Point4d, vectorU: Point4d, vectorV: Point4d, result?: PlaneByOriginAndVectors4d): PlaneByOriginAndVectors4d {\r\n    if (result) {\r\n      result.setOriginAndVectors(origin, vectorU, vectorV);\r\n      return result;\r\n    }\r\n    return new PlaneByOriginAndVectors4d(origin.clone(), vectorU.clone(), vectorV.clone());\r\n  }\r\n  /** Set all numeric data from complete list of (x,y,z,w) in origin, vectorU, and vectorV */\r\n  public setOriginAndVectorsXYZW(x0: number, y0: number, z0: number, w0: number, ux: number, uy: number, uz: number, uw: number, vx: number, vy: number, vz: number, vw: number): PlaneByOriginAndVectors4d {\r\n    this.origin.set(x0, y0, z0, w0);\r\n    this.vectorU.set(ux, uy, uz, uw);\r\n    this.vectorV.set(vx, vy, vz, vw);\r\n    return this;\r\n  }\r\n  /** Copy the contents of origin, vectorU, vectorV parameters to respective member variables */\r\n  public setOriginAndVectors(origin: Point4d, vectorU: Point4d, vectorV: Point4d): PlaneByOriginAndVectors4d {\r\n    this.origin.setFrom(origin);\r\n    this.vectorU.setFrom(vectorU);\r\n    this.vectorV.setFrom(vectorV);\r\n    return this;\r\n  }\r\n  /** Create from complete list of (x,y,z,w) in origin, vectorU, and vectorV */\r\n  public static createOriginAndVectorsXYZW(x0: number, y0: number, z0: number, w0: number, ux: number, uy: number, uz: number, uw: number, vx: number, vy: number, vz: number, vw: number, result?: PlaneByOriginAndVectors4d): PlaneByOriginAndVectors4d {\r\n    if (result)\r\n      return result.setOriginAndVectorsXYZW(x0, y0, z0, w0, ux, uy, uz, uw, vx, vy, vz, vw);\r\n    return new PlaneByOriginAndVectors4d(Point4d.create(x0, y0, z0, w0), Point4d.create(ux, uy, uz, uw), Point4d.create(vx, vy, vz, uw));\r\n  }\r\n  /** create from origin point, (u=1,v=0) point, and (u=0,v=1) point. */\r\n  public static createOriginAndTargets3d(origin: Point3d, targetU: Point3d, targetV: Point3d, result?: PlaneByOriginAndVectors4d): PlaneByOriginAndVectors4d {\r\n    return PlaneByOriginAndVectors4d.createOriginAndVectorsXYZW(origin.x, origin.y, origin.z, 1.0, targetU.x - origin.x, targetU.y - origin.y, targetU.z - origin.z, 0.0, targetV.x - origin.x, targetV.y - origin.y, targetV.z - origin.z, 0.0, result);\r\n  }\r\n  /** evaluate plane point (full 3d) at given (u,v) coordinate. */\r\n  public fractionToPoint(u: number, v: number, result?: Point4d): Point4d {\r\n    return this.origin.plus2Scaled(this.vectorU, u, this.vectorV, v, result);\r\n  }\r\n  /** create a new plane which maps to the cartesian xy plane. */\r\n  public static createXYPlane(result?: PlaneByOriginAndVectors4d): PlaneByOriginAndVectors4d {\r\n    return PlaneByOriginAndVectors4d.createOriginAndVectorsXYZW(0, 0, 0, 1, 1, 0, 0, 0, 0, 1, 0, 0, result);\r\n  }\r\n}\r\n"]}

package/lib/esm/geometry4d/Point4d.js

@@ -34,6 +34,8 @@ function quotientDerivative2(ddg, dh, ddh, f, df, divH) {
  * @public
  */
 export class Point4d extends Plane3d {
+    /** x,y,z,w are packed into a Float64Array */
+    xyzw;
     /** Set x,y,z,w of this point.  */
     set(x = 0, y = 0, z = 0, w = 0) {
         this.xyzw[0] = x;