@itwin/core-geometry 3.0.0-dev.76 → 3.0.0-dev.80

package/lib/esm/geometry4d/MomentData.js

@@ -0,0 +1,384 @@
+/*---------------------------------------------------------------------------------------------
+* Copyright (c) Bentley Systems, Incorporated. All rights reserved.
+* See LICENSE.md in the project root for license terms and full copyright notice.
+*--------------------------------------------------------------------------------------------*/
+/** @packageDocumentation
+ * @module Numerics
+ */
+// import { Point2d } from "./Geometry2d";
+/* eslint-disable @typescript-eslint/naming-convention */
+import { Geometry } from "../Geometry";
+import { Matrix3d } from "../geometry3d/Matrix3d";
+import { Point3d, Vector3d } from "../geometry3d/Point3dVector3d";
+import { Transform } from "../geometry3d/Transform";
+import { Matrix4d } from "./Matrix4d";
+import { Point4d } from "./Point4d";
+/**
+ * A MomentData structure exists in several levels:
+ * * First level: as a carrier of sums of inertial products that determine moments.
+ *   * 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
+ *         * 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,
+ *        * 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.
+ * @public
+ */
+export class MomentData {
+    constructor() {
+        this._point0 = Point3d.create();
+        this._point1 = Point3d.create();
+        this.origin = Point3d.createZero();
+        this.sums = Matrix4d.createZero();
+        this.localToWorldMap = Transform.createIdentity();
+        this.radiusOfGyration = Vector3d.create();
+        this.needOrigin = false;
+        this.absoluteQuantity = 0.1; // so optimizer sees its type
+        this.absoluteQuantity = undefined;
+    }
+    /** Return the lower-right (3,3) entry in the sums.
+     * * This is the quantity (i.e. length, area, or volume) summed
+     */
+    get quantitySum() { return this.sums.atIJ(3, 3); }
+    /** Return 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.
+     *
+     */
+    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.
+     *
+     */
+    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.
+     *
+     */
+    setOriginXYZIfNeeded(x, y, z) {
+        if (this.needOrigin) {
+            this.origin.set(x, y, z);
+            this.needOrigin = false;
+        }
+    }
+    /** 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
+    */
+    static create(origin, needOrigin = false) {
+        const data = new MomentData();
+        data.needOrigin = needOrigin;
+        if (origin) {
+            data.origin.setFromPoint3d(origin);
+            data.needOrigin = false;
+        }
+        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]`
+     */
+    static momentTensorFromInertiaProducts(products) {
+        const rr = products.sumDiagonal();
+        const result = Matrix3d.createScale(rr, rr, rr);
+        result.addScaledInPlace(products, -1.0);
+        return result;
+    }
+    /** Sort the columns of the matrix for increasing moments. */
+    static sortColumnsForIncreasingMoments(axes, moments) {
+        const points = [
+            axes.indexedColumnWithWeight(0, moments.x),
+            axes.indexedColumnWithWeight(1, moments.y),
+            axes.indexedColumnWithWeight(2, moments.z),
+        ].sort((dataA, dataB) => {
+            if (dataA.w < dataB.w)
+                return -1;
+            if (dataA.w > dataB.w)
+                return 1;
+            return 0;
+        });
+        axes.setColumnsPoint4dXYZ(points[0], points[1], points[2]);
+        if (axes.determinant() < 0)
+            axes.scaleColumnsInPlace(-1.0, -1.0, -1.0);
+        // prefer x and z positive -- y falls wherever . ..
+        if (axes.at(0, 0) < 0.0)
+            axes.scaleColumnsInPlace(-1.0, -1.0, 1.0);
+        if (axes.at(2, 2) < 0.0)
+            axes.scaleColumnsInPlace(1.0, -1.0, -1.0);
+        moments.set(points[0].w, points[1].w, points[2].w);
+    }
+    /**
+     * Return the principal moment data for an array of points.
+     * @param points array of points
+     */
+    static pointsToPrincipalAxes(points) {
+        const moments = new MomentData();
+        if (points.length === 0)
+            return moments;
+        moments.clearSums(points[0]);
+        moments.accumulatePointMomentsFromOrigin(points);
+        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
+     * @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]
+     */
+    static inertiaProductsToPrincipalAxes(origin, inertiaProducts) {
+        const moments = new MomentData();
+        moments.sums.setFrom(inertiaProducts);
+        moments.origin.setFrom(origin);
+        if (!moments.shiftOriginAndSumsToCentroidOfSums())
+            return undefined;
+        const products = moments.sums.matrixPart();
+        const w = moments.sums.weight();
+        if (w < 0.0)
+            products.scaleColumnsInPlace(-1, -1, -1);
+        const tensor = MomentData.momentTensorFromInertiaProducts(products);
+        const moment2 = Vector3d.create();
+        const axisVectors = Matrix3d.createZero();
+        tensor.fastSymmetricEigenvalues(axisVectors, moment2);
+        if (moment2.x < 0.0)
+            return undefined;
+        MomentData.sortColumnsForIncreasingMoments(axisVectors, moment2);
+        if (w < 0.0)
+            axisVectors.scaleColumnsInPlace(1, -1, -1);
+        moments.localToWorldMap = Transform.createOriginAndMatrix(moments.origin, axisVectors);
+        moments.radiusOfGyration.set(Math.sqrt(Math.abs(moment2.x)), Math.sqrt(Math.abs(moment2.y)), Math.sqrt(Math.abs(moment2.z)));
+        moments.radiusOfGyration.scaleInPlace(1.0 / Math.sqrt(Math.abs(w)));
+        moments.absoluteQuantity = Math.abs(w);
+        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
+     *   * 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.
+     * * If either or both are undefined, returns false.
+     * @param dataA first set of moments
+     * @param dataB second set of moments
+     */
+    static areEquivalentPrincipalAxes(dataA, dataB) {
+        if (dataA && dataB
+            && Geometry.isSameCoordinate(dataA.quantitySum, dataB.quantitySum)) { // um.. need different tolerance for area, volume?)
+            if (dataA.localToWorldMap.getOrigin().isAlmostEqual(dataB.localToWorldMap.getOrigin())
+                && dataA.radiusOfGyration.isAlmostEqual(dataB.radiusOfGyration)) {
+                if (Geometry.isSameCoordinate(dataA.radiusOfGyration.x, dataA.radiusOfGyration.y)) {
+                    // We have at least xy symmetry ....
+                    if (Geometry.isSameCoordinate(dataA.radiusOfGyration.x, dataA.radiusOfGyration.z))
+                        return true;
+                    // just xy.
+                    // allow opposite z directions.
+                    // If the z's are aligned, x an dy can spin freely.
+                    const zA = dataA.localToWorldMap.matrix.columnZ();
+                    const zB = dataB.localToWorldMap.matrix.columnZ();
+                    if (zA.isParallelTo(zB, true))
+                        return true;
+                    return false;
+                }
+                // no symmetry.  Test all three axes.
+                const vectorA = Vector3d.create();
+                const vectorB = Vector3d.create();
+                for (let i = 0; i < 3; i++) {
+                    dataA.localToWorldMap.matrix.getColumn(i, vectorA);
+                    dataB.localToWorldMap.matrix.getColumn(i, vectorB);
+                    if (!vectorA.isParallelTo(vectorB, true))
+                        return false;
+                }
+                return true;
+            }
+        }
+        return false;
+    }
+    /** Clear the MomentData sums to zero, and establish a new origin. */
+    clearSums(origin) {
+        this.sums.setZero();
+        if (origin)
+            this.origin.setFrom(origin);
+        else
+            this.origin.setZero();
+    }
+    /** Accumulate products-of-components for given points. */
+    accumulatePointMomentsFromOrigin(points) {
+        for (const p of points) {
+            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" */
+    shiftOriginAndSumsToCentroidOfSums() {
+        const xyz = this.sums.columnW().realPoint();
+        if (xyz) {
+            this.shiftOriginAndSumsByXYZ(xyz.x, xyz.y, xyz.z);
+            return true;
+        }
+        return false;
+    }
+    /** revise the accumulated sums
+     * * add ax,ay,ax to the origin coordinates.
+     * * apply the negative translation to the sums.
+    */
+    shiftOriginAndSumsByXYZ(ax, ay, az) {
+        this.origin.addXYZInPlace(ax, ay, az);
+        this.sums.multiplyTranslationSandwichInPlace(-ax, -ay, -az);
+    }
+    /** 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.
+     * Accumulate them to this.sums.
+     * * If `pointA` is undefined, use `this.origin` as pointA.
+     * * If `this.needOrigin` is set, pointB is used
+    */
+    accumulateTriangleMomentsXY(pointA, pointB, pointC) {
+        this.setOriginXYZIfNeeded(pointB.x, pointB.y, 0.0);
+        const x0 = this.origin.x;
+        const y0 = this.origin.y;
+        const vectorA = MomentData._vectorA =
+            pointA !== undefined ? Point4d.create(pointA.x - x0, pointA.y - y0, 0.0, 1.0, MomentData._vectorA)
+                : Point4d.create(this.origin.x, this.origin.y, 0.0, 1.0, MomentData._vectorA);
+        const vectorB = MomentData._vectorB = Point4d.create(pointB.x - x0, pointB.y - y0, 0.0, 1.0, MomentData._vectorB);
+        const vectorC = MomentData._vectorC = Point4d.create(pointC.x - x0, pointC.y - y0, 0.0, 1.0, MomentData._vectorC);
+        // accumulate Return product integrals I(0<=u<=1) I (0<=v<= u)  (w*W + u *U + v * V)(w*W + u *U + v * V)^  du dv
+        //  where w = 1-u-v
+        //  W = column vector (point00.x, point00.y, point00.z, 1.0) etc.
+        const detJ = Geometry.crossProductXYXY(vectorB.x - vectorA.x, vectorB.y - vectorA.y, vectorC.x - vectorA.x, vectorC.y - vectorA.y);
+        if (detJ !== 0.0) {
+            const r1_12 = detJ / 12.0;
+            const r1_24 = detJ / 24.0;
+            this.sums.addScaledOuterProductInPlace(vectorA, vectorA, r1_12);
+            this.sums.addScaledOuterProductInPlace(vectorA, vectorB, r1_24);
+            this.sums.addScaledOuterProductInPlace(vectorA, vectorC, r1_24);
+            this.sums.addScaledOuterProductInPlace(vectorB, vectorA, r1_24);
+            this.sums.addScaledOuterProductInPlace(vectorB, vectorB, r1_12);
+            this.sums.addScaledOuterProductInPlace(vectorB, vectorC, r1_24);
+            this.sums.addScaledOuterProductInPlace(vectorC, vectorA, r1_24);
+            this.sums.addScaledOuterProductInPlace(vectorC, vectorB, r1_24);
+            this.sums.addScaledOuterProductInPlace(vectorC, vectorC, r1_12);
+        }
+    }
+    /** 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);
+        this.sums.addScaledOuterProductInPlace(vectorA, vectorA, scaleFactor);
+    }
+    /** Accumulate wire moment integral from pointA to pointB */
+    accumulateLineMomentsXYZ(pointA, pointB) {
+        this.setOriginXYZIfNeeded(pointA.x, pointA.y, pointA.z);
+        const x0 = this.origin.x;
+        const y0 = this.origin.y;
+        const z0 = this.origin.z;
+        const vectorA = MomentData._vectorA = Point4d.create(pointA.x - x0, pointA.y - y0, pointA.z - z0, 1.0, MomentData._vectorA);
+        const vectorB = MomentData._vectorB = Point4d.create(pointB.x - x0, pointB.y - y0, pointB.z - z0, 1.0, MomentData._vectorB);
+        const detJ = pointA.distance(pointB);
+        const r1_3 = detJ / 3.0;
+        const r1_6 = detJ / 6.0;
+        this.sums.addScaledOuterProductInPlace(vectorA, vectorA, r1_3);
+        this.sums.addScaledOuterProductInPlace(vectorA, vectorB, r1_6);
+        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.
+     * 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.
+     *
+     */
+    accumulateTriangleToLineStringMomentsXY(sweepBase, points) {
+        const n = points.length;
+        if (n > 1) {
+            points.getPoint3dAtUncheckedPointIndex(0, this._point0);
+            for (let i = 1; i < n; i++) {
+                points.getPoint3dAtUncheckedPointIndex(i, this._point1);
+                this.accumulateTriangleMomentsXY(sweepBase, this._point0, this._point1);
+                this._point0.setFromPoint3d(this._point1);
+            }
+        }
+    }
+    // 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)
+     * * 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)
+     */
+    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);
+        const detJ = Geometry.crossProductXYXY(vectorU.x, vectorV.x, vectorU.y, vectorV.y);
+        const placement = Matrix4d.createRowValues(vectorU.x, vectorV.x, 0, origin.x - this.origin.x, vectorU.y, vectorV.y, 0, origin.y - this.origin.y, 0, 0, 0, 0, 0, 0, 0, 1);
+        const AB = placement.multiplyMatrixMatrix(centroidalProducts);
+        const ABAT = AB.multiplyMatrixMatrixTranspose(placement);
+        this.sums.addScaledInPlace(ABAT, detJ);
+    }
+    /**
+     * 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.
+     * *
+     */
+    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.
+   * *
+   */
+    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:
+     */
+    toJSON() {
+        return {
+            origin: this.origin,
+            sums: this.sums.toJSON(),
+            radiusOfGyration: this.radiusOfGyration.toJSON(),
+            localToWorld: this.localToWorldMap.toJSON(),
+        };
+    }
+}
+//# sourceMappingURL=MomentData.js.map

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

@@ -0,0 +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;IAkErB;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;IA5DD;;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;YACnB,IAAI,CAAC,MAAM,CAAC,cAAc,CAAC,MAAM,CAAC,CAAC;YACnC,IAAI,CAAC,UAAU,GAAG,KAAK,CAAC;SACzB;IACH,CAAC;IACD;;;OAGG;IACI,qCAAqC,CAAC,MAAwB;QACnE,IAAI,IAAI,CAAC,UAAU,IAAI,MAAM,CAAC,MAAM,GAAG,CAAC,EAAE;YACxC,MAAM,CAAC,6BAA6B,CAAC,CAAC,EAAE,IAAI,CAAC,MAAM,CAAC,CAAC;YACrD,IAAI,CAAC,UAAU,GAAG,KAAK,CAAC;SACzB;IACH,CAAC;IAED;;;OAGG;IACI,oBAAoB,CAAC,CAAS,EAAE,CAAS,EAAE,CAAS;QACzD,IAAI,IAAI,CAAC,UAAU,EAAE;YACnB,IAAI,CAAC,MAAM,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;YACzB,IAAI,CAAC,UAAU,GAAG,KAAK,CAAC;SACzB;IACH,CAAC;IAmBD;;;;;;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;YACV,IAAI,CAAC,MAAM,CAAC,cAAc,CAAC,MAAM,CAAC,CAAC;YACnC,IAAI,CAAC,UAAU,GAAG,KAAK,CAAC;SACzB;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,EAAG,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;gBACjE,IAAI,QAAQ,CAAC,gBAAgB,CAAC,KAAK,CAAC,gBAAgB,CAAC,CAAC,EAAE,KAAK,CAAC,gBAAgB,CAAC,CAAC,CAAC,EAAE;oBACjF,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;iBACd;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;oBAC1B,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;iBAChB;gBACD,OAAO,IAAI,CAAC;aACb;SACF;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;YACtB,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;SACR;IACH,CAAC;IACD,8DAA8D;IACvD,kCAAkC;QACvC,MAAM,GAAG,GAAG,IAAI,CAAC,IAAI,CAAC,OAAO,EAAE,CAAC,SAAS,EAAE,CAAC;QAC5C,IAAI,GAAG,EAAE;YACP,IAAI,CAAC,uBAAuB,CAAC,GAAG,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC;YAClD,OAAO,IAAI,CAAC;SACb;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;YAChB,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;SACjE;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;YACT,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;gBAC1B,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;aAC3C;SACF;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"]}

package/lib/esm/geometry4d/PlaneByOriginAndVectors4d.d.ts

@@ -0,0 +1,43 @@
+/** @packageDocumentation
+ * @module Numerics
+ */
+import { Point3d } from "../geometry3d/Point3dVector3d";
+import { Point4d } from "./Point4d";
+/**
+ * A Plane4dByOriginAndVectors is a 4d origin and pair of 4d "vectors" defining a 4d plane.
+ * * The parameterization of the plane is    `X = origin + vectorU*u + vectorV * v`
+ * * With particular weight values `origin.w === 1, vectorU.w === 0, vectorV.w === 0` this is like `Plane3dByOriginAndVectors`
+ * * With other weights, the deweighted xyz coordinates of points on the 4d plane still form a 3d plane.
+ * @public
+ */
+export declare class PlaneByOriginAndVectors4d {
+    /** homogeneous origin */
+    origin: Point4d;
+    /** homogeneous u-direction vector */
+    vectorU: Point4d;
+    /** homogeneous v-direction vector */
+    vectorV: Point4d;
+    private constructor();
+    /** Return a clone of this plane */
+    clone(result?: PlaneByOriginAndVectors4d): PlaneByOriginAndVectors4d;
+    /** copy all content from other plane */
+    setFrom(other: PlaneByOriginAndVectors4d): void;
+    /** Return true if origin, vectorU, and vectorV pass isAlmostEqual. */
+    isAlmostEqual(other: PlaneByOriginAndVectors4d): boolean;
+    /** Create a plane with (copies of) origin, vectorU, vectorV parameters, all given as full 4d points.
+     */
+    static createOriginAndVectors(origin: Point4d, vectorU: Point4d, vectorV: Point4d, result?: PlaneByOriginAndVectors4d): PlaneByOriginAndVectors4d;
+    /** Set all numeric data from complete list of (x,y,z,w) in origin, vectorU, and vectorV */
+    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;
+    /** Copy the contents of origin, vectorU, vectorV parameters to respective member variables */
+    setOriginAndVectors(origin: Point4d, vectorU: Point4d, vectorV: Point4d): PlaneByOriginAndVectors4d;
+    /** Create from complete list of (x,y,z,w) in origin, vectorU, and vectorV */
+    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;
+    /** create from origin point, (u=1,v=0) point, and (u=0,v=1) point. */
+    static createOriginAndTargets3d(origin: Point3d, targetU: Point3d, targetV: Point3d, result?: PlaneByOriginAndVectors4d): PlaneByOriginAndVectors4d;
+    /** evaluate plane point (full 3d) at given (u,v) coordinate. */
+    fractionToPoint(u: number, v: number, result?: Point4d): Point4d;
+    /** create a new plane which maps to the cartesian xy plane. */
+    static createXYPlane(result?: PlaneByOriginAndVectors4d): PlaneByOriginAndVectors4d;
+}
+//# sourceMappingURL=PlaneByOriginAndVectors4d.d.ts.map

package/lib/esm/geometry4d/PlaneByOriginAndVectors4d.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"PlaneByOriginAndVectors4d.d.ts","sourceRoot":"","sources":["../../../src/geometry4d/PlaneByOriginAndVectors4d.ts"],"names":[],"mappings":"AAKA;;GAEG;AAEH,OAAO,EAAE,OAAO,EAAE,MAAM,+BAA+B,CAAC;AACxD,OAAO,EAAE,OAAO,EAAE,MAAM,WAAW,CAAC;AAEpC;;;;;;GAMG;AACH,qBAAa,yBAAyB;IACpC,yBAAyB;IAClB,MAAM,EAAE,OAAO,CAAC;IACvB,qCAAqC;IAC9B,OAAO,EAAE,OAAO,CAAC;IACxB,qCAAqC;IAC9B,OAAO,EAAE,OAAO,CAAC;IACxB,OAAO;IAKP,mCAAmC;IAC5B,KAAK,CAAC,MAAM,CAAC,EAAE,yBAAyB,GAAG,yBAAyB;IAO3E,wCAAwC;IACjC,OAAO,CAAC,KAAK,EAAE,yBAAyB,GAAG,IAAI;IAKtD,sEAAsE;IAC/D,aAAa,CAAC,KAAK,EAAE,yBAAyB,GAAG,OAAO;IAK/D;OACG;WACW,sBAAsB,CAAC,MAAM,EAAE,OAAO,EAAE,OAAO,EAAE,OAAO,EAAE,OAAO,EAAE,OAAO,EAAE,MAAM,CAAC,EAAE,yBAAyB,GAAG,yBAAyB;IAOxJ,2FAA2F;IACpF,uBAAuB,CAAC,EAAE,EAAE,MAAM,EAAE,EAAE,EAAE,MAAM,EAAE,EAAE,EAAE,MAAM,EAAE,EAAE,EAAE,MAAM,EAAE,EAAE,EAAE,MAAM,EAAE,EAAE,EAAE,MAAM,EAAE,EAAE,EAAE,MAAM,EAAE,EAAE,EAAE,MAAM,EAAE,EAAE,EAAE,MAAM,EAAE,EAAE,EAAE,MAAM,EAAE,EAAE,EAAE,MAAM,EAAE,EAAE,EAAE,MAAM,GAAG,yBAAyB;IAMzM,8FAA8F;IACvF,mBAAmB,CAAC,MAAM,EAAE,OAAO,EAAE,OAAO,EAAE,OAAO,EAAE,OAAO,EAAE,OAAO,GAAG,yBAAyB;IAM1G,6EAA6E;WAC/D,0BAA0B,CAAC,EAAE,EAAE,MAAM,EAAE,EAAE,EAAE,MAAM,EAAE,EAAE,EAAE,MAAM,EAAE,EAAE,EAAE,MAAM,EAAE,EAAE,EAAE,MAAM,EAAE,EAAE,EAAE,MAAM,EAAE,EAAE,EAAE,MAAM,EAAE,EAAE,EAAE,MAAM,EAAE,EAAE,EAAE,MAAM,EAAE,EAAE,EAAE,MAAM,EAAE,EAAE,EAAE,MAAM,EAAE,EAAE,EAAE,MAAM,EAAE,MAAM,CAAC,EAAE,yBAAyB,GAAG,yBAAyB;IAKvP,sEAAsE;WACxD,wBAAwB,CAAC,MAAM,EAAE,OAAO,EAAE,OAAO,EAAE,OAAO,EAAE,OAAO,EAAE,OAAO,EAAE,MAAM,CAAC,EAAE,yBAAyB,GAAG,yBAAyB;IAG1J,gEAAgE;IACzD,eAAe,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,MAAM,CAAC,EAAE,OAAO,GAAG,OAAO;IAGvE,+DAA+D;WACjD,aAAa,CAAC,MAAM,CAAC,EAAE,yBAAyB,GAAG,yBAAyB;CAG3F"}

package/lib/esm/geometry4d/PlaneByOriginAndVectors4d.js

@@ -0,0 +1,81 @@
+/*---------------------------------------------------------------------------------------------
+* Copyright (c) Bentley Systems, Incorporated. All rights reserved.
+* See LICENSE.md in the project root for license terms and full copyright notice.
+*--------------------------------------------------------------------------------------------*/
+import { Point4d } from "./Point4d";
+/**
+ * A Plane4dByOriginAndVectors is a 4d origin and pair of 4d "vectors" defining a 4d plane.
+ * * The parameterization of the plane is    `X = origin + vectorU*u + vectorV * v`
+ * * With particular weight values `origin.w === 1, vectorU.w === 0, vectorV.w === 0` this is like `Plane3dByOriginAndVectors`
+ * * With other weights, the deweighted xyz coordinates of points on the 4d plane still form a 3d plane.
+ * @public
+ */
+export class PlaneByOriginAndVectors4d {
+    constructor(origin, vectorU, vectorV) {
+        this.origin = origin;
+        this.vectorU = vectorU;
+        this.vectorV = vectorV;
+    }
+    /** Return a clone of this plane */
+    clone(result) {
+        if (result) {
+            result.setFrom(this);
+            return result;
+        }
+        return new PlaneByOriginAndVectors4d(this.origin.clone(), this.vectorU.clone(), this.vectorV.clone());
+    }
+    /** copy all content from other plane */
+    setFrom(other) {
+        this.origin.setFrom(other.origin);
+        this.vectorU.setFrom(other.vectorU);
+        this.vectorV.setFrom(other.vectorV);
+    }
+    /** Return true if origin, vectorU, and vectorV pass isAlmostEqual. */
+    isAlmostEqual(other) {
+        return this.origin.isAlmostEqual(other.origin)
+            && this.vectorU.isAlmostEqual(other.vectorU)
+            && this.vectorV.isAlmostEqual(other.vectorV);
+    }
+    /** Create a plane with (copies of) origin, vectorU, vectorV parameters, all given as full 4d points.
+     */
+    static createOriginAndVectors(origin, vectorU, vectorV, result) {
+        if (result) {
+            result.setOriginAndVectors(origin, vectorU, vectorV);
+            return result;
+        }
+        return new PlaneByOriginAndVectors4d(origin.clone(), vectorU.clone(), vectorV.clone());
+    }
+    /** Set all numeric data from complete list of (x,y,z,w) in origin, vectorU, and vectorV */
+    setOriginAndVectorsXYZW(x0, y0, z0, w0, ux, uy, uz, uw, vx, vy, vz, vw) {
+        this.origin.set(x0, y0, z0, w0);
+        this.vectorU.set(ux, uy, uz, uw);
+        this.vectorV.set(vx, vy, vz, vw);
+        return this;
+    }
+    /** Copy the contents of origin, vectorU, vectorV parameters to respective member variables */
+    setOriginAndVectors(origin, vectorU, vectorV) {
+        this.origin.setFrom(origin);
+        this.vectorU.setFrom(vectorU);
+        this.vectorV.setFrom(vectorV);
+        return this;
+    }
+    /** Create from complete list of (x,y,z,w) in origin, vectorU, and vectorV */
+    static createOriginAndVectorsXYZW(x0, y0, z0, w0, ux, uy, uz, uw, vx, vy, vz, vw, result) {
+        if (result)
+            return result.setOriginAndVectorsXYZW(x0, y0, z0, w0, ux, uy, uz, uw, vx, vy, vz, vw);
+        return new PlaneByOriginAndVectors4d(Point4d.create(x0, y0, z0, w0), Point4d.create(ux, uy, uz, uw), Point4d.create(vx, vy, vz, uw));
+    }
+    /** create from origin point, (u=1,v=0) point, and (u=0,v=1) point. */
+    static createOriginAndTargets3d(origin, targetU, targetV, result) {
+        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);
+    }
+    /** evaluate plane point (full 3d) at given (u,v) coordinate. */
+    fractionToPoint(u, v, result) {
+        return this.origin.plus2Scaled(this.vectorU, u, this.vectorV, v, result);
+    }
+    /** create a new plane which maps to the cartesian xy plane. */
+    static createXYPlane(result) {
+        return PlaneByOriginAndVectors4d.createOriginAndVectorsXYZW(0, 0, 0, 1, 1, 0, 0, 0, 0, 1, 0, 0, result);
+    }
+}
+//# sourceMappingURL=PlaneByOriginAndVectors4d.js.map

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

@@ -0,0 +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;YACV,MAAM,CAAC,OAAO,CAAC,IAAI,CAAC,CAAC;YACrB,OAAO,MAAM,CAAC;SACf;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;YACV,MAAM,CAAC,mBAAmB,CAAC,MAAM,EAAE,OAAO,EAAE,OAAO,CAAC,CAAC;YACrD,OAAO,MAAM,CAAC;SACf;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 Plane4dByOriginAndVectors 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"]}