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

package/lib/esm/clipping/ClipPlane.js

@@ -0,0 +1,505 @@
+/*---------------------------------------------------------------------------------------------
+* Copyright (c) Bentley Systems, Incorporated. All rights reserved.
+* See LICENSE.md in the project root for license terms and full copyright notice.
+*--------------------------------------------------------------------------------------------*/
+import { AxisOrder, Geometry } from "../Geometry";
+import { GrowableFloat64Array } from "../geometry3d/GrowableFloat64Array";
+import { GrowableXYZArray } from "../geometry3d/GrowableXYZArray";
+import { Matrix3d } from "../geometry3d/Matrix3d";
+import { Plane3dByOriginAndUnitNormal } from "../geometry3d/Plane3dByOriginAndUnitNormal";
+import { Point3d, Vector3d } from "../geometry3d/Point3dVector3d";
+import { IndexedXYZCollectionPolygonOps } from "../geometry3d/PolygonOps";
+import { Range1d, Range3d } from "../geometry3d/Range";
+import { Transform } from "../geometry3d/Transform";
+import { Point4d } from "../geometry4d/Point4d";
+import { AnalyticRoots } from "../numerics/Polynomials";
+import { ClipUtilities } from "./ClipUtils";
+/** A ClipPlane is a single plane represented as
+ * * An inward unit normal (u,v,w)
+ * * A signedDistance
+ *
+ * Hence
+ * * The halfspace function evaluation for "point" [x,y,z,] is: ([x,y,z] DOT (u,v,w)l - signedDistance)
+ * * POSITIVE values of the halfspace function are "inside"
+ * * ZERO value of the halfspace function is "on"
+ * * NEGATIVE value of the halfspace function is "outside"
+ * * A representative point on the plane is (signedDistance*u, signedDistance * v, signedDistance *w)
+ * * Given a point and inward normal, the signedDistance is (point DOT normal)
+ * @public
+ */
+export class ClipPlane {
+    constructor(normal, distance, invisible, interior) {
+        this._invisible = invisible;
+        this._interior = interior;
+        this._inwardNormal = normal;
+        this._distanceFromOrigin = distance;
+    }
+    /*
+      private safeSetXYZDistance(nx: number, ny: number, nz: number, d: number) {
+        this._inwardNormal.set(nx, ny, nz);
+        this._distanceFromOrigin = d;
+      }
+    */
+    /**
+     * Return true if all members are almostEqual to corresponding members of other.
+     * @param other clip plane to compare
+     */
+    isAlmostEqual(other) {
+        return Geometry.isSameCoordinate(this._distanceFromOrigin, other._distanceFromOrigin)
+            && this._inwardNormal.isAlmostEqual(other._inwardNormal)
+            && this._interior === other._interior
+            && this._invisible === other._invisible;
+    }
+    /** return a cloned plane */
+    clone() {
+        const result = new ClipPlane(this._inwardNormal.clone(), this._distanceFromOrigin, this._invisible, this._interior);
+        return result;
+    }
+    /** return Return a cloned plane with coordinate data negated. */
+    cloneNegated() {
+        const plane = new ClipPlane(this._inwardNormal.clone(), this._distanceFromOrigin, this._invisible, this._interior);
+        plane.negateInPlace();
+        return plane;
+    }
+    /** Create a ClipPlane from Plane3dByOriginAndUnitNormal. */
+    static createPlane(plane, invisible = false, interior = false, result) {
+        const distance = plane.getNormalRef().dotProduct(plane.getOriginRef());
+        if (result) {
+            result._invisible = invisible;
+            result._interior = interior;
+            result._inwardNormal = plane.getNormalRef().clone();
+            result._distanceFromOrigin = distance;
+            return result;
+        }
+        return new ClipPlane(plane.getNormalRef().clone(), distance, invisible, interior);
+    }
+    /**
+     * * Create a ClipPlane with direct normal and signedDistance.
+     * * The vector is normalized for storage.
+     */
+    static createNormalAndDistance(normal, distance, invisible = false, interior = false, result) {
+        const normalized = normal.normalize();
+        if (normalized) {
+            if (result) {
+                result._invisible = invisible;
+                result._interior = interior;
+                result._inwardNormal = normalized;
+                result._distanceFromOrigin = distance;
+            }
+            return new ClipPlane(normalized, distance, invisible, interior);
+        }
+        return undefined;
+    }
+    /** Create a ClipPlane
+     * * "normal" is the inward normal of the plane. (It is internally normalized)
+     * * "point" is any point of the plane.
+     * * The stored distance for the plane is the dot product of the point with the normal (i.e. treat the point's xyz as a vector from the origin.)
+     */
+    static createNormalAndPoint(normal, point, invisible = false, interior = false, result) {
+        const normalized = normal.normalize();
+        if (normalized) {
+            const distance = normalized.dotProduct(point);
+            if (result) {
+                result._invisible = invisible;
+                result._interior = interior;
+                result._inwardNormal = normalized;
+                result._distanceFromOrigin = distance;
+            }
+            return new ClipPlane(normalized, distance, invisible, interior);
+        }
+        return undefined;
+    }
+    /** Create a ClipPlane
+     * * "normal" (normalX, normalY, nz) is the inward normal of the plane.
+     * * The given (normalX,normalY,normalZ)
+     * * "point" is any point of the plane.
+     * * The stored distance for the plane is the dot product of the point with the normal (i.e. treat the point's xyz as a vector from the origin.)
+     */
+    static createNormalAndPointXYZXYZ(normalX, normalY, normalZ, originX, originY, originZ, invisible = false, interior = false, result) {
+        const q = Geometry.hypotenuseXYZ(normalX, normalY, normalZ);
+        const r = Geometry.conditionalDivideFraction(1, q);
+        if (r !== undefined) {
+            if (result) {
+                result._inwardNormal.set(normalX * r, normalY * r, normalZ * r);
+                result._distanceFromOrigin = result._inwardNormal.dotProductXYZ(originX, originY, originZ);
+                result._invisible = invisible;
+                result._interior = interior;
+                return result;
+            }
+            const normal = Vector3d.create(normalX * r, normalY * r, normalZ * r);
+            return new ClipPlane(normal, normal.dotProductXYZ(originX, originY, originZ), invisible, interior);
+        }
+        return undefined;
+    }
+    /**
+     * return a json object of the form
+     * `{"normal":[u,v,w],"dist":signedDistanceValue,"interior":true,"invisible":true}`
+     */
+    toJSON() {
+        const props = {
+            normal: this.inwardNormalRef.toJSON(),
+            dist: this.distance,
+        };
+        if (this.interior)
+            props.interior = true;
+        if (this.invisible)
+            props.invisible = true;
+        return props;
+    }
+    /** parse json object to ClipPlane instance */
+    static fromJSON(json, result) {
+        if (json && json.normal && undefined !== json.dist && Number.isFinite(json.dist))
+            return ClipPlane.createNormalAndDistance(Vector3d.fromJSON(json.normal), json.dist, !!json.invisible, !!json.interior);
+        return ClipPlane.createNormalAndDistance(Vector3d.unitZ(), 0, false, false, result);
+    }
+    /** Set both the invisible and interior flags. */
+    setFlags(invisible, interior) {
+        this._invisible = invisible;
+        this._interior = interior;
+    }
+    /**
+     * Return the stored distanceFromOrigin property.
+     */
+    get distance() { return this._distanceFromOrigin; }
+    /**
+     * Return the stored inward normal property.
+     */
+    get inwardNormalRef() { return this._inwardNormal; }
+    /**
+     * Return the "interior" property bit
+     */
+    get interior() { return this._interior; }
+    /**
+     * Return the "invisible" property bit.
+     */
+    get invisible() { return this._invisible; }
+    /**
+     * Create a plane defined by two points, an up vector, and a tilt angle relative to the up vector.
+     * @param point0 start point of the edge
+     * @param point1 end point of the edge
+     * @param upVector vector perpendicular to the plane
+     * @param tiltAngle angle to tilt the plane around the edge in the direction of the up vector.
+     * @param result optional preallocated plane
+     */
+    static createEdgeAndUpVector(point0, point1, upVector, tiltAngle, result) {
+        const edgeVector = Vector3d.createFrom(point1.minus(point0));
+        let normal = (upVector.crossProduct(edgeVector)).normalize();
+        if (normal) {
+            if (tiltAngle !== undefined && !tiltAngle.isAlmostZero) {
+                const tiltNormal = Vector3d.createRotateVectorAroundVector(normal, edgeVector, tiltAngle);
+                if (tiltNormal) {
+                    normal = tiltNormal.clone();
+                }
+            }
+            normal.negate(normal);
+            return ClipPlane.createNormalAndPoint(normal, point0, false, false, result);
+        }
+        return undefined;
+    }
+    /**
+     * Create a plane perpendicular to the edge between the xy parts of point0 and point1
+     */
+    static createEdgeXY(point0, point1, result) {
+        const normal = Vector3d.create(point0.y - point1.y, point1.x - point0.x);
+        if (normal.normalizeInPlace())
+            return ClipPlane.createNormalAndPoint(normal, point0, false, false, result);
+        return undefined;
+    }
+    /**
+     * Return the Plane3d form of the plane.
+     * * The plane origin is the point `distance * inwardNormal`
+     * * The plane normal is the inward normal of the ClipPlane.
+     */
+    getPlane3d() {
+        const d = this._distanceFromOrigin;
+        // Normal should be normalized, will not return undefined
+        return Plane3dByOriginAndUnitNormal.create(Point3d.create(this._inwardNormal.x * d, this._inwardNormal.y * d, this._inwardNormal.z * d), this._inwardNormal);
+    }
+    /**
+     * Return the Point4d d form of the plane.
+     * * The homogeneous xyz are the inward normal xyz.
+     * * The homogeneous weight is the negated ClipPlane distance.
+     */
+    getPlane4d() {
+        return Point4d.create(this._inwardNormal.x, this._inwardNormal.y, this._inwardNormal.z, -this._distanceFromOrigin);
+    }
+    /**
+     * Set the plane from DPoint4d style plane.
+     * * The saved plane has its direction normalized.
+     * * This preserves the plane itself as a zero set but make plane evaluations act as true distances (even if the plane coefficients are scaled otherwise)
+     * @param plane
+     */
+    setPlane4d(plane) {
+        const a = Math.sqrt(plane.x * plane.x + plane.y * plane.y + plane.z * plane.z);
+        const r = a === 0.0 ? 1.0 : 1.0 / a;
+        this._inwardNormal.x = r * plane.x;
+        this._inwardNormal.y = r * plane.y;
+        this._inwardNormal.z = r * plane.z;
+        this._distanceFromOrigin = -r * plane.w;
+    }
+    /**
+     * Evaluate the altitude in weighted space, i.e. (dot product with inward normal) minus distance, with point.w scale applied to distance)
+     * @param point space point to test
+     */
+    weightedAltitude(point) {
+        return point.x * this._inwardNormal.x + point.y * this._inwardNormal.y + point.z * this._inwardNormal.z - point.w * this._distanceFromOrigin;
+    }
+    /**
+     * Evaluate the distance from the plane to a point in space, i.e. (dot product with inward normal) minus distance
+     * @param point space point to test
+     */
+    altitude(point) {
+        return point.x * this._inwardNormal.x + point.y * this._inwardNormal.y + point.z * this._inwardNormal.z - this._distanceFromOrigin;
+    }
+    /**
+     * Evaluate the distance from the plane to a point in space with point given as x,y,z, i.e. (dot product with inward normal) minus distance
+     * @param point space point to test
+     */
+    altitudeXYZ(x, y, z) {
+        return x * this._inwardNormal.x + y * this._inwardNormal.y + z * this._inwardNormal.z - this._distanceFromOrigin;
+    }
+    /** Return the dot product of the plane normal with the vector (NOT using the plane's distanceFromOrigin).
+     */
+    velocity(vector) {
+        return vector.x * this._inwardNormal.x + vector.y * this._inwardNormal.y + vector.z * this._inwardNormal.z;
+    }
+    /** Return the dot product of the plane normal with the x,yz, vector components (NOT using the plane's distanceFromOrigin).
+     */
+    velocityXYZ(x, y, z) {
+        return x * this._inwardNormal.x + y * this._inwardNormal.y + z * this._inwardNormal.z;
+    }
+    /** Return the dot product of the plane normal with the point (treating the point xyz as a vector, and NOT using the plane's distanceFromOrigin).
+     */
+    dotProductPlaneNormalPoint(point) {
+        return point.x * this._inwardNormal.x + point.y * this._inwardNormal.y + point.z * this._inwardNormal.z;
+    }
+    /**
+     * Return true if spacePoint is inside or on the plane, with tolerance applied to "on".
+     * @param spacePoint point to test.
+     * @param tolerance tolerance for considering "near plane" to be "on plane"
+     */
+    isPointOnOrInside(spacePoint, tolerance = Geometry.smallMetricDistance) {
+        let value = this.altitude(spacePoint);
+        if (tolerance) {
+            value += tolerance;
+        }
+        return value >= 0.0;
+    }
+    /**
+     * Return true if spacePoint is strictly inside the halfspace, with tolerance applied to "on".
+     * @param spacePoint point to test.
+     * @param tolerance tolerance for considering "near plane" to be "on plane"
+     */
+    isPointInside(point, tolerance = Geometry.smallMetricDistance) {
+        let value = this.altitude(point);
+        if (tolerance) {
+            value -= tolerance;
+        }
+        return value > 0.0;
+    }
+    /**
+     * Return true if spacePoint is strictly on the plane, within tolerance
+     * @param spacePoint point to test.
+     * @param tolerance tolerance for considering "near plane" to be "on plane"
+     */
+    isPointOn(point, tolerance = Geometry.smallMetricDistance) {
+        return Math.abs(this.altitude(point)) <= tolerance;
+    }
+    /**
+     * Compute intersections of an (UNBOUNDED) arc with the plane.  Append them (as radians) to a growing array.
+     * @param arc arc to test.  The angle limits of the arc are NOT considered.
+     * @param intersectionRadians array to receive results
+     */
+    appendIntersectionRadians(arc, intersectionRadians) {
+        const arcVectors = arc.toVectors();
+        const alpha = this.altitude(arc.center);
+        const beta = this.velocity(arcVectors.vector0);
+        const gamma = this.velocity(arcVectors.vector90);
+        AnalyticRoots.appendImplicitLineUnitCircleIntersections(alpha, beta, gamma, undefined, undefined, intersectionRadians);
+    }
+    /** Announce fractional intervals of arc clip.
+     * * Each call to `announce(fraction0, fraction1, arc)` announces one interval that is inside the clip plane.
+     */
+    announceClippedArcIntervals(arc, announce) {
+        const breaks = ClipPlane._clipArcFractionArray;
+        breaks.clear();
+        this.appendIntersectionRadians(arc, breaks);
+        arc.sweep.radiansArrayToPositivePeriodicFractions(breaks);
+        return ClipUtilities.selectIntervals01(arc, breaks, this, announce);
+    }
+    /**
+     * * Compute intersection of (unbounded) segment with the plane.
+     * * If the ends are on the same side of the plane, return undefined.
+     * * If the intersection is an endpoint or interior to the segment return the fraction.
+     * * If both ends are on, return undefined.
+     */
+    getBoundedSegmentSimpleIntersection(pointA, pointB) {
+        const h0 = this.altitude(pointA);
+        const h1 = this.altitude(pointB);
+        if (h0 * h1 > 0.0)
+            return undefined;
+        if (h0 === 0.0 && h1 === 0.0) {
+            return undefined;
+        }
+        return -h0 / (h1 - h0);
+    }
+    /** Apply transform to the origin.  Apply inverse transpose of the matrix part to th normal vector. */
+    transformInPlace(transform) {
+        const plane = this.getPlane3d();
+        const matrix = transform.matrix;
+        const newPoint = transform.multiplyPoint3d(plane.getOriginRef());
+        // Normal transforms as the inverse transpose of the matrix part
+        // BTW: If the matrix is orthogonal, this is a long way to multiply by the matrix part (mumble grumble)
+        const newNormal = matrix.multiplyInverseTranspose(plane.getNormalRef());
+        if (!newNormal)
+            return false;
+        plane.set(newPoint, newNormal);
+        const normalized = (plane.getNormalRef()).normalize();
+        if (!normalized)
+            return false;
+        this._inwardNormal = normalized;
+        this._distanceFromOrigin = this._inwardNormal.dotProduct(plane.getOriginRef());
+        return true;
+    }
+    /** Set the invisible flag.   Interpretation of this is up to the use code algorithms. */
+    setInvisible(invisible) {
+        this._invisible = invisible;
+    }
+    /**  reverse the sign of all coefficients, so outside and inside reverse */
+    negateInPlace() {
+        this._inwardNormal = this._inwardNormal.negate();
+        this._distanceFromOrigin = -this._distanceFromOrigin;
+    }
+    /**
+     * Move the plane INWARD by given distance
+     * @param offset distance of shift inwards
+     */
+    offsetDistance(offset) {
+        this._distanceFromOrigin += offset;
+    }
+    /**
+     * Clip a polygon to the inside or outside of the plane.
+     * * Results with 2 or fewer points are ignored.
+     * * Other than ensuring capacity in the arrays, there are no object allocations during execution of this function.
+     * @param xyz input points.
+     * @param work work buffer
+     * @param tolerance tolerance for "on plane" decision.
+     */
+    clipConvexPolygonInPlace(xyz, work, inside = true, tolerance = Geometry.smallMetricDistance) {
+        return IndexedXYZCollectionPolygonOps.clipConvexPolygonInPlace(this, xyz, work, inside, tolerance);
+    }
+    /**
+     * Multiply the ClipPlane's DPoint4d by matrix.
+     * @param matrix matrix to apply.
+     * @param invert if true, use in verse of the matrix.
+     * @param transpose if true, use the transpose of the matrix (or inverse, per invert parameter)
+     * * Note that if matrixA is applied to all of space, the matrix to send to this method to get a corresponding effect on the plane is the inverse transpose of matrixA
+     * * Callers that will apply the same matrix to many planes should pre-invert the matrix for efficiency.
+     * * Both params default to true to get the full effect of transforming space.
+     * @param matrix matrix to apply
+     * @return false if unable to invert
+     */
+    multiplyPlaneByMatrix4d(matrix, invert = true, transpose = true) {
+        const plane = this.getPlane4d();
+        if (invert) {
+            const inverse = matrix.createInverse();
+            if (inverse)
+                return this.multiplyPlaneByMatrix4d(inverse, false, transpose);
+            return false;
+        }
+        if (transpose)
+            matrix.multiplyTransposePoint4d(plane, plane);
+        else
+            matrix.multiplyPoint4d(plane, plane);
+        this.setPlane4d(plane);
+        return true;
+    }
+    /** announce the interval (if any) where a line is within the clip plane half space. */
+    announceClippedSegmentIntervals(f0, f1, pointA, pointB, announce) {
+        if (f1 < f0)
+            return false;
+        const h0 = -this.altitude(pointA);
+        const h1 = -this.altitude(pointB);
+        const delta = h1 - h0;
+        const f = Geometry.conditionalDivideFraction(-h0, delta);
+        if (f === undefined) { // The segment is parallel to the plane.
+            if (h0 <= 0.0) {
+                if (announce)
+                    announce(f0, f1);
+                return true;
+            }
+            return false;
+        }
+        if (delta > 0) { // segment aims OUT
+            if (f < f1)
+                f1 = f;
+        }
+        else {
+            // segment aims IN
+            if (f > f0)
+                f0 = f;
+        }
+        if (f1 < f0)
+            return false;
+        if (announce)
+            announce(f0, f1);
+        return true;
+    }
+    /**
+     * Return a coordinate frame with
+     * * origin at closest point to global origin
+     * * z axis points in
+     * x and y are "in plane"
+     */
+    getFrame() {
+        const d = this._distanceFromOrigin;
+        const origin = Point3d.create(this._inwardNormal.x * d, this._inwardNormal.y * d, this._inwardNormal.z * d);
+        const matrix = Matrix3d.createRigidHeadsUp(this._inwardNormal, AxisOrder.ZXY);
+        return Transform.createOriginAndMatrix(origin, matrix);
+    }
+    /**
+     * Return the intersection of the plane with a range cube.
+     * @param range
+     * @param xyzOut intersection polygon.  This is convex.
+     */
+    intersectRange(range, addClosurePoint = false) {
+        if (range.isNull)
+            return undefined;
+        const corners = range.corners();
+        const frameOnPlane = this.getFrame();
+        frameOnPlane.multiplyInversePoint3dArrayInPlace(corners);
+        const localRange = Range3d.createArray(corners);
+        if (localRange.low.z * localRange.high.z > 0.0)
+            return undefined;
+        // oversized polygon on local z= 0
+        const xyzOut = new GrowableXYZArray();
+        xyzOut.pushXYZ(localRange.low.x, localRange.low.y, 0);
+        xyzOut.pushXYZ(localRange.high.x, localRange.low.y, 0);
+        xyzOut.pushXYZ(localRange.high.x, localRange.high.y, 0);
+        xyzOut.pushXYZ(localRange.low.x, localRange.high.y, 0);
+        xyzOut.multiplyTransformInPlace(frameOnPlane);
+        IndexedXYZCollectionPolygonOps.intersectRangeConvexPolygonInPlace(range, xyzOut);
+        if (xyzOut.length === 0)
+            return undefined;
+        if (addClosurePoint)
+            xyzOut.pushWrap(1);
+        return xyzOut;
+    }
+    /** Implement appendPolygonClip, as defined in interface PolygonClipper.  /**
+     *
+     * @param xyz input polygon.  This is not changed.
+     * @param insideFragments Array to receive "inside" fragments.  Each fragment is a GrowableXYZArray grabbed from the cache.  This is NOT cleared.
+     * @param outsideFragments Array to receive "outside" fragments.  Each fragment is a GrowableXYZArray grabbed from the cache.  This is NOT cleared.
+     * @param arrayCache cache for reusable GrowableXYZArray.
+     */
+    appendPolygonClip(xyz, insideFragments, outsideFragments, arrayCache) {
+        const perpendicularRange = Range1d.createNull();
+        const newInside = arrayCache.grabFromCache();
+        const newOutside = arrayCache.grabFromCache();
+        IndexedXYZCollectionPolygonOps.splitConvexPolygonInsideOutsidePlane(this, xyz, newInside, newOutside, perpendicularRange);
+        ClipUtilities.captureOrDrop(newInside, 3, insideFragments, arrayCache);
+        ClipUtilities.captureOrDrop(newOutside, 3, outsideFragments, arrayCache);
+    }
+}
+ClipPlane._clipArcFractionArray = new GrowableFloat64Array();
+//# sourceMappingURL=ClipPlane.js.map

package/lib/esm/clipping/ClipPlane.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"ClipPlane.js","sourceRoot":"","sources":["../../../src/clipping/ClipPlane.ts"],"names":[],"mappings":"AAAA;;;+FAG+F;AAQ/F,OAAO,EAAE,SAAS,EAAE,QAAQ,EAA0B,MAAM,aAAa,CAAC;AAG1E,OAAO,EAAE,oBAAoB,EAAE,MAAM,oCAAoC,CAAC;AAC1E,OAAO,EAAE,gBAAgB,EAAE,MAAM,gCAAgC,CAAC;AAClE,OAAO,EAAE,QAAQ,EAAE,MAAM,wBAAwB,CAAC;AAClD,OAAO,EAAE,4BAA4B,EAAE,MAAM,4CAA4C,CAAC;AAC1F,OAAO,EAAE,OAAO,EAAE,QAAQ,EAAE,MAAM,+BAA+B,CAAC;AAClE,OAAO,EAAE,8BAA8B,EAAE,MAAM,0BAA0B,CAAC;AAC1E,OAAO,EAAE,OAAO,EAAE,OAAO,EAAE,MAAM,qBAAqB,CAAC;AACvD,OAAO,EAAE,SAAS,EAAE,MAAM,yBAAyB,CAAC;AAEpD,OAAO,EAAE,OAAO,EAAE,MAAM,uBAAuB,CAAC;AAChD,OAAO,EAAE,aAAa,EAAE,MAAM,yBAAyB,CAAC;AACxD,OAAO,EAAW,aAAa,EAAkB,MAAM,aAAa,CAAC;AAkBrE;;;;;;;;;;;;GAYG;AACH,MAAM,OAAO,SAAS;IAUpB,YAAoB,MAAgB,EAAE,QAAgB,EAAE,SAAkB,EAAE,QAAiB;QAC3F,IAAI,CAAC,UAAU,GAAG,SAAS,CAAC;QAC5B,IAAI,CAAC,SAAS,GAAG,QAAQ,CAAC;QAC1B,IAAI,CAAC,aAAa,GAAG,MAAM,CAAC;QAC5B,IAAI,CAAC,mBAAmB,GAAG,QAAQ,CAAC;IACtC,CAAC;IACD;;;;;MAKE;IACF;;;OAGG;IACI,aAAa,CAAC,KAAgB;QACnC,OAAO,QAAQ,CAAC,gBAAgB,CAAC,IAAI,CAAC,mBAAmB,EAAE,KAAK,CAAC,mBAAmB,CAAC;eAChF,IAAI,CAAC,aAAa,CAAC,aAAa,CAAC,KAAK,CAAC,aAAa,CAAC;eACrD,IAAI,CAAC,SAAS,KAAK,KAAK,CAAC,SAAS;eAClC,IAAI,CAAC,UAAU,KAAK,KAAK,CAAC,UAAU,CAAC;IAC5C,CAAC;IACD,4BAA4B;IACrB,KAAK;QACV,MAAM,MAAM,GAAG,IAAI,SAAS,CAAC,IAAI,CAAC,aAAa,CAAC,KAAK,EAAE,EAAE,IAAI,CAAC,mBAAmB,EAAE,IAAI,CAAC,UAAU,EAAE,IAAI,CAAC,SAAS,CAAC,CAAC;QACpH,OAAO,MAAM,CAAC;IAChB,CAAC;IACD,iEAAiE;IAC1D,YAAY;QACjB,MAAM,KAAK,GAAG,IAAI,SAAS,CAAC,IAAI,CAAC,aAAa,CAAC,KAAK,EAAE,EAAE,IAAI,CAAC,mBAAmB,EAAE,IAAI,CAAC,UAAU,EAAE,IAAI,CAAC,SAAS,CAAC,CAAC;QACnH,KAAK,CAAC,aAAa,EAAE,CAAC;QACtB,OAAO,KAAK,CAAC;IACf,CAAC;IACD,4DAA4D;IACrD,MAAM,CAAC,WAAW,CAAC,KAAmC,EAAE,YAAqB,KAAK,EAAE,WAAoB,KAAK,EAAE,MAAkB;QACtI,MAAM,QAAQ,GAAG,KAAK,CAAC,YAAY,EAAE,CAAC,UAAU,CAAC,KAAK,CAAC,YAAY,EAAE,CAAC,CAAC;QACvE,IAAI,MAAM,EAAE;YACV,MAAM,CAAC,UAAU,GAAG,SAAS,CAAC;YAC9B,MAAM,CAAC,SAAS,GAAG,QAAQ,CAAC;YAC5B,MAAM,CAAC,aAAa,GAAG,KAAK,CAAC,YAAY,EAAE,CAAC,KAAK,EAAE,CAAC;YACpD,MAAM,CAAC,mBAAmB,GAAG,QAAQ,CAAC;YACtC,OAAO,MAAM,CAAC;SACf;QACD,OAAO,IAAI,SAAS,CAAC,KAAK,CAAC,YAAY,EAAE,CAAC,KAAK,EAAE,EAAE,QAAQ,EAAE,SAAS,EAAE,QAAQ,CAAC,CAAC;IACpF,CAAC;IACD;;;OAGG;IACI,MAAM,CAAC,uBAAuB,CAAC,MAAgB,EAAE,QAAgB,EAAE,YAAqB,KAAK,EAAE,WAAoB,KAAK,EAAE,MAAkB;QACjJ,MAAM,UAAU,GAAG,MAAM,CAAC,SAAS,EAAE,CAAC;QACtC,IAAI,UAAU,EAAE;YACd,IAAI,MAAM,EAAE;gBACV,MAAM,CAAC,UAAU,GAAG,SAAS,CAAC;gBAC9B,MAAM,CAAC,SAAS,GAAG,QAAQ,CAAC;gBAC5B,MAAM,CAAC,aAAa,GAAG,UAAU,CAAC;gBAClC,MAAM,CAAC,mBAAmB,GAAG,QAAQ,CAAC;aACvC;YACD,OAAO,IAAI,SAAS,CAAC,UAAU,EAAE,QAAQ,EAAE,SAAS,EAAE,QAAQ,CAAC,CAAC;SACjE;QACD,OAAO,SAAS,CAAC;IACnB,CAAC;IACD;;;;OAIG;IACI,MAAM,CAAC,oBAAoB,CAAC,MAAgB,EAAE,KAAc,EAAE,YAAqB,KAAK,EAAE,WAAoB,KAAK,EAAE,MAAkB;QAC5I,MAAM,UAAU,GAAG,MAAM,CAAC,SAAS,EAAE,CAAC;QACtC,IAAI,UAAU,EAAE;YACd,MAAM,QAAQ,GAAG,UAAU,CAAC,UAAU,CAAC,KAAK,CAAC,CAAC;YAC9C,IAAI,MAAM,EAAE;gBACV,MAAM,CAAC,UAAU,GAAG,SAAS,CAAC;gBAC9B,MAAM,CAAC,SAAS,GAAG,QAAQ,CAAC;gBAC5B,MAAM,CAAC,aAAa,GAAG,UAAU,CAAC;gBAClC,MAAM,CAAC,mBAAmB,GAAG,QAAQ,CAAC;aACvC;YACD,OAAO,IAAI,SAAS,CAAC,UAAU,EAAE,QAAQ,EAAE,SAAS,EAAE,QAAQ,CAAC,CAAC;SACjE;QACD,OAAO,SAAS,CAAC;IACnB,CAAC;IAED;;;;;OAKG;IACI,MAAM,CAAC,0BAA0B,CAAC,OAAe,EAAE,OAAe,EAAE,OAAe,EACxF,OAAe,EAAE,OAAe,EAAE,OAAe,EACjD,YAAqB,KAAK,EAAE,WAAoB,KAAK,EAAE,MAAkB;QACzE,MAAM,CAAC,GAAG,QAAQ,CAAC,aAAa,CAAC,OAAO,EAAE,OAAO,EAAE,OAAO,CAAC,CAAC;QAC5D,MAAM,CAAC,GAAG,QAAQ,CAAC,yBAAyB,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;QACnD,IAAI,CAAC,KAAK,SAAS,EAAE;YACnB,IAAI,MAAM,EAAE;gBACV,MAAM,CAAC,aAAa,CAAC,GAAG,CAAC,OAAO,GAAG,CAAC,EAAE,OAAO,GAAG,CAAC,EAAE,OAAO,GAAG,CAAC,CAAC,CAAC;gBAChE,MAAM,CAAC,mBAAmB,GAAG,MAAM,CAAC,aAAa,CAAC,aAAa,CAAC,OAAO,EAAE,OAAO,EAAE,OAAO,CAAC,CAAC;gBAC3F,MAAM,CAAC,UAAU,GAAG,SAAS,CAAC;gBAC9B,MAAM,CAAC,SAAS,GAAG,QAAQ,CAAC;gBAC5B,OAAO,MAAM,CAAC;aACf;YACD,MAAM,MAAM,GAAG,QAAQ,CAAC,MAAM,CAAC,OAAO,GAAG,CAAC,EAAE,OAAO,GAAG,CAAC,EAAE,OAAO,GAAG,CAAC,CAAC,CAAC;YACtE,OAAO,IAAI,SAAS,CAAC,MAAM,EAAE,MAAM,CAAC,aAAa,CAAC,OAAO,EAAE,OAAO,EAAE,OAAO,CAAC,EAAE,SAAS,EAAE,QAAQ,CAAC,CAAC;SACpG;QACD,OAAO,SAAS,CAAC;IACnB,CAAC;IAED;;;OAGG;IACI,MAAM;QACX,MAAM,KAAK,GAAmB;YAC5B,MAAM,EAAE,IAAI,CAAC,eAAe,CAAC,MAAM,EAAE;YACrC,IAAI,EAAE,IAAI,CAAC,QAAQ;SACpB,CAAC;QAEF,IAAI,IAAI,CAAC,QAAQ;YACf,KAAK,CAAC,QAAQ,GAAG,IAAI,CAAC;QAExB,IAAI,IAAI,CAAC,SAAS;YAChB,KAAK,CAAC,SAAS,GAAG,IAAI,CAAC;QAEzB,OAAO,KAAK,CAAC;IACf,CAAC;IAED,8CAA8C;IACvC,MAAM,CAAC,QAAQ,CAAC,IAAoB,EAAE,MAAkB;QAC7D,IAAI,IAAI,IAAI,IAAI,CAAC,MAAM,IAAI,SAAS,KAAK,IAAI,CAAC,IAAI,IAAI,MAAM,CAAC,QAAQ,CAAC,IAAI,CAAC,IAAI,CAAC;YAC9E,OAAO,SAAS,CAAC,uBAAuB,CAAC,QAAQ,CAAC,QAAQ,CAAC,IAAI,CAAC,MAAM,CAAC,EAAE,IAAI,CAAC,IAAI,EAAE,CAAC,CAAC,IAAI,CAAC,SAAS,EAAE,CAAC,CAAC,IAAI,CAAC,QAAQ,CAAC,CAAC;QAEzH,OAAO,SAAS,CAAC,uBAAuB,CAAC,QAAQ,CAAC,KAAK,EAAE,EAAE,CAAC,EAAE,KAAK,EAAE,KAAK,EAAE,MAAM,CAAC,CAAC;IACtF,CAAC;IAED,iDAAiD;IAC1C,QAAQ,CAAC,SAAkB,EAAE,QAAiB;QACnD,IAAI,CAAC,UAAU,GAAG,SAAS,CAAC;QAC5B,IAAI,CAAC,SAAS,GAAG,QAAQ,CAAC;IAC5B,CAAC;IAED;;OAEG;IACH,IAAW,QAAQ,KAAK,OAAO,IAAI,CAAC,mBAAmB,CAAC,CAAC,CAAC;IAC1D;;OAEG;IACH,IAAW,eAAe,KAAe,OAAO,IAAI,CAAC,aAAa,CAAC,CAAC,CAAC;IACrE;;OAEG;IACH,IAAW,QAAQ,KAAK,OAAO,IAAI,CAAC,SAAS,CAAC,CAAC,CAAC;IAChD;;OAEG;IACH,IAAW,SAAS,KAAK,OAAO,IAAI,CAAC,UAAU,CAAC,CAAC,CAAC;IAElD;;;;;;;OAOG;IACI,MAAM,CAAC,qBAAqB,CAAC,MAAe,EAAE,MAAe,EAAE,QAAkB,EAAE,SAAiB,EAAE,MAAkB;QAC7H,MAAM,UAAU,GAAG,QAAQ,CAAC,UAAU,CAAC,MAAM,CAAC,KAAK,CAAC,MAAM,CAAC,CAAC,CAAC;QAC7D,IAAI,MAAM,GAAG,CAAC,QAAQ,CAAC,YAAY,CAAC,UAAU,CAAC,CAAC,CAAC,SAAS,EAAE,CAAC;QAE7D,IAAI,MAAM,EAAE;YACV,IAAI,SAAS,KAAK,SAAS,IAAI,CAAC,SAAS,CAAC,YAAY,EAAE;gBACtD,MAAM,UAAU,GAAG,QAAQ,CAAC,8BAA8B,CAAC,MAAM,EAAE,UAAU,EAAE,SAAS,CAAC,CAAC;gBAC1F,IAAI,UAAU,EAAE;oBACd,MAAM,GAAG,UAAU,CAAC,KAAK,EAAE,CAAC;iBAC7B;aACF;YACD,MAAM,CAAC,MAAM,CAAC,MAAM,CAAC,CAAC;YACtB,OAAO,SAAS,CAAC,oBAAoB,CAAC,MAAM,EAAE,MAAM,EAAE,KAAK,EAAE,KAAK,EAAE,MAAM,CAAC,CAAC;SAC7E;QACD,OAAO,SAAS,CAAC;IACnB,CAAC;IACD;;OAEG;IACI,MAAM,CAAC,YAAY,CAAC,MAAe,EAAE,MAAe,EAAE,MAAkB;QAC7E,MAAM,MAAM,GAAG,QAAQ,CAAC,MAAM,CAAC,MAAM,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC;QACzE,IAAI,MAAM,CAAC,gBAAgB,EAAE;YAC3B,OAAO,SAAS,CAAC,oBAAoB,CAAC,MAAM,EAAE,MAAM,EAAE,KAAK,EAAE,KAAK,EAAE,MAAM,CAAC,CAAC;QAC9E,OAAO,SAAS,CAAC;IACnB,CAAC;IACD;;;;OAIG;IACI,UAAU;QACf,MAAM,CAAC,GAAG,IAAI,CAAC,mBAAmB,CAAC;QACnC,yDAAyD;QACzD,OAAO,4BAA4B,CAAC,MAAM,CAAC,OAAO,CAAC,MAAM,CAAC,IAAI,CAAC,aAAa,CAAC,CAAC,GAAG,CAAC,EAAE,IAAI,CAAC,aAAa,CAAC,CAAC,GAAG,CAAC,EAAE,IAAI,CAAC,aAAa,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,IAAI,CAAC,aAAa,CAAE,CAAC;IAChK,CAAC;IAED;;;;OAIG;IACI,UAAU;QACf,OAAO,OAAO,CAAC,MAAM,CAAC,IAAI,CAAC,aAAa,CAAC,CAAC,EAAE,IAAI,CAAC,aAAa,CAAC,CAAC,EAAE,IAAI,CAAC,aAAa,CAAC,CAAC,EAAE,CAAE,IAAI,CAAC,mBAAmB,CAAC,CAAC;IACtH,CAAC;IACD;;;;;OAKG;IACI,UAAU,CAAC,KAAc;QAC9B,MAAM,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC;QAC/E,MAAM,CAAC,GAAG,CAAC,KAAK,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC;QACpC,IAAI,CAAC,aAAa,CAAC,CAAC,GAAG,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC;QACnC,IAAI,CAAC,aAAa,CAAC,CAAC,GAAG,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC;QACnC,IAAI,CAAC,aAAa,CAAC,CAAC,GAAG,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC;QACnC,IAAI,CAAC,mBAAmB,GAAG,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC;IAC1C,CAAC;IAED;;;OAGG;IACI,gBAAgB,CAAC,KAAc;QACpC,OAAO,KAAK,CAAC,CAAC,GAAG,IAAI,CAAC,aAAa,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,GAAG,IAAI,CAAC,aAAa,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,GAAG,IAAI,CAAC,aAAa,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,GAAG,IAAI,CAAC,mBAAmB,CAAC;IAC/I,CAAC;IAED;;;OAGG;IACI,QAAQ,CAAC,KAAc;QAC5B,OAAO,KAAK,CAAC,CAAC,GAAG,IAAI,CAAC,aAAa,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,GAAG,IAAI,CAAC,aAAa,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,GAAG,IAAI,CAAC,aAAa,CAAC,CAAC,GAAG,IAAI,CAAC,mBAAmB,CAAC;IACrI,CAAC;IAED;;;OAGG;IACI,WAAW,CAAC,CAAS,EAAE,CAAS,EAAE,CAAS;QAChD,OAAO,CAAC,GAAG,IAAI,CAAC,aAAa,CAAC,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,aAAa,CAAC,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,aAAa,CAAC,CAAC,GAAG,IAAI,CAAC,mBAAmB,CAAC;IACnH,CAAC;IAED;OACG;IACI,QAAQ,CAAC,MAAgB;QAC9B,OAAO,MAAM,CAAC,CAAC,GAAG,IAAI,CAAC,aAAa,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,GAAG,IAAI,CAAC,aAAa,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,GAAG,IAAI,CAAC,aAAa,CAAC,CAAC,CAAC;IAC7G,CAAC;IAED;OACG;IACI,WAAW,CAAC,CAAS,EAAE,CAAS,EAAE,CAAS;QAChD,OAAO,CAAC,GAAG,IAAI,CAAC,aAAa,CAAC,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,aAAa,CAAC,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,aAAa,CAAC,CAAC,CAAC;IACxF,CAAC;IAED;OACG;IACI,0BAA0B,CAAC,KAAc;QAC9C,OAAO,KAAK,CAAC,CAAC,GAAG,IAAI,CAAC,aAAa,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,GAAG,IAAI,CAAC,aAAa,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,GAAG,IAAI,CAAC,aAAa,CAAC,CAAC,CAAC;IAC1G,CAAC;IACD;;;;OAIG;IACI,iBAAiB,CAAC,UAAmB,EAAE,YAAoB,QAAQ,CAAC,mBAAmB;QAC5F,IAAI,KAAK,GAAG,IAAI,CAAC,QAAQ,CAAC,UAAU,CAAC,CAAC;QACtC,IAAI,SAAS,EAAE;YAAE,KAAK,IAAI,SAAS,CAAC;SAAE;QACtC,OAAO,KAAK,IAAI,GAAG,CAAC;IACtB,CAAC;IAED;;;;OAIG;IACI,aAAa,CAAC,KAAc,EAAE,YAAoB,QAAQ,CAAC,mBAAmB;QACnF,IAAI,KAAK,GAAG,IAAI,CAAC,QAAQ,CAAC,KAAK,CAAC,CAAC;QACjC,IAAI,SAAS,EAAE;YAAE,KAAK,IAAI,SAAS,CAAC;SAAE;QACtC,OAAO,KAAK,GAAG,GAAG,CAAC;IACrB,CAAC;IAED;;;;OAIG;IACI,SAAS,CAAC,KAAc,EAAE,YAAoB,QAAQ,CAAC,mBAAmB;QAC/E,OAAO,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,QAAQ,CAAC,KAAK,CAAC,CAAC,IAAI,SAAS,CAAC;IACrD,CAAC;IACD;;;;OAIG;IACI,yBAAyB,CAAC,GAAU,EAAE,mBAAyC;QACpF,MAAM,UAAU,GAAG,GAAG,CAAC,SAAS,EAAE,CAAC;QACnC,MAAM,KAAK,GAAG,IAAI,CAAC,QAAQ,CAAC,GAAG,CAAC,MAAM,CAAC,CAAC;QACxC,MAAM,IAAI,GAAG,IAAI,CAAC,QAAQ,CAAC,UAAU,CAAC,OAAO,CAAC,CAAC;QAC/C,MAAM,KAAK,GAAG,IAAI,CAAC,QAAQ,CAAC,UAAU,CAAC,QAAQ,CAAC,CAAC;QACjD,aAAa,CAAC,yCAAyC,CAAC,KAAK,EAAE,IAAI,EAAE,KAAK,EAAE,SAAS,EAAE,SAAS,EAAE,mBAAmB,CAAC,CAAC;IACzH,CAAC;IAGD;;OAEG;IACI,2BAA2B,CAAC,GAAU,EAAE,QAA6C;QAC1F,MAAM,MAAM,GAAG,SAAS,CAAC,qBAAqB,CAAC;QAC/C,MAAM,CAAC,KAAK,EAAE,CAAC;QACf,IAAI,CAAC,yBAAyB,CAAC,GAAG,EAAE,MAAM,CAAC,CAAC;QAC5C,GAAG,CAAC,KAAK,CAAC,uCAAuC,CAAC,MAAM,CAAC,CAAC;QAC1D,OAAO,aAAa,CAAC,iBAAiB,CAAC,GAAG,EAAE,MAAM,EAAE,IAAI,EAAE,QAAQ,CAAC,CAAC;IACtE,CAAC;IACD;;;;;OAKG;IACI,mCAAmC,CAAC,MAAe,EAAE,MAAe;QACzE,MAAM,EAAE,GAAG,IAAI,CAAC,QAAQ,CAAC,MAAM,CAAC,CAAC;QACjC,MAAM,EAAE,GAAG,IAAI,CAAC,QAAQ,CAAC,MAAM,CAAC,CAAC;QACjC,IAAI,EAAE,GAAG,EAAE,GAAG,GAAG;YACf,OAAO,SAAS,CAAC;QACnB,IAAI,EAAE,KAAK,GAAG,IAAI,EAAE,KAAK,GAAG,EAAE;YAC5B,OAAO,SAAS,CAAC;SAClB;QACD,OAAO,CAAE,EAAE,GAAG,CAAC,EAAE,GAAG,EAAE,CAAC,CAAC;IAC1B,CAAC;IAED,sGAAsG;IAC/F,gBAAgB,CAAC,SAAoB;QAC1C,MAAM,KAAK,GAAiC,IAAI,CAAC,UAAU,EAAE,CAAC;QAC9D,MAAM,MAAM,GAAa,SAAS,CAAC,MAAM,CAAC;QAC1C,MAAM,QAAQ,GAAG,SAAS,CAAC,eAAe,CAAC,KAAK,CAAC,YAAY,EAAE,CAAC,CAAC;QACjE,gEAAgE;QAChE,uGAAuG;QACvG,MAAM,SAAS,GAAG,MAAM,CAAC,wBAAwB,CAAC,KAAK,CAAC,YAAY,EAAE,CAAC,CAAC;QACxE,IAAI,CAAC,SAAS;YACZ,OAAO,KAAK,CAAC;QAEf,KAAK,CAAC,GAAG,CAAC,QAAQ,EAAE,SAAS,CAAC,CAAC;QAC/B,MAAM,UAAU,GAAG,CAAC,KAAK,CAAC,YAAY,EAAE,CAAC,CAAC,SAAS,EAAE,CAAC;QACtD,IAAI,CAAC,UAAU;YACb,OAAO,KAAK,CAAC;QACf,IAAI,CAAC,aAAa,GAAG,UAAU,CAAC;QAChC,IAAI,CAAC,mBAAmB,GAAG,IAAI,CAAC,aAAa,CAAC,UAAU,CAAC,KAAK,CAAC,YAAY,EAAE,CAAC,CAAC;QAC/E,OAAO,IAAI,CAAC;IACd,CAAC;IACD,yFAAyF;IAClF,YAAY,CAAC,SAAkB;QACpC,IAAI,CAAC,UAAU,GAAG,SAAS,CAAC;IAC9B,CAAC;IAED,2EAA2E;IACpE,aAAa;QAClB,IAAI,CAAC,aAAa,GAAG,IAAI,CAAC,aAAa,CAAC,MAAM,EAAE,CAAC;QACjD,IAAI,CAAC,mBAAmB,GAAG,CAAE,IAAI,CAAC,mBAAmB,CAAC;IACxD,CAAC;IACD;;;OAGG;IACI,cAAc,CAAC,MAAc;QAClC,IAAI,CAAC,mBAAmB,IAAI,MAAM,CAAC;IACrC,CAAC;IAED;;;;;;;OAOG;IACI,wBAAwB,CAAC,GAAqB,EAAE,IAAsB,EAAE,SAAkB,IAAI,EAAE,YAAoB,QAAQ,CAAC,mBAAmB;QACrJ,OAAO,8BAA8B,CAAC,wBAAwB,CAAC,IAAI,EAAE,GAAG,EAAE,IAAI,EAAE,MAAM,EAAE,SAAS,CAAC,CAAC;IACrG,CAAC;IAED;;;;;;;;;;OAUG;IACI,uBAAuB,CAAC,MAAgB,EAAE,SAAkB,IAAI,EAAE,YAAqB,IAAI;QAChG,MAAM,KAAK,GAAY,IAAI,CAAC,UAAU,EAAE,CAAC;QACzC,IAAI,MAAM,EAAE;YACV,MAAM,OAAO,GAAG,MAAM,CAAC,aAAa,EAAE,CAAC;YACvC,IAAI,OAAO;gBACT,OAAO,IAAI,CAAC,uBAAuB,CAAC,OAAO,EAAE,KAAK,EAAE,SAAS,CAAC,CAAC;YACjE,OAAO,KAAK,CAAC;SACd;QACD,IAAI,SAAS;YACX,MAAM,CAAC,wBAAwB,CAAC,KAAK,EAAE,KAAK,CAAC,CAAC;;YAE9C,MAAM,CAAC,eAAe,CAAC,KAAK,EAAE,KAAK,CAAC,CAAC;QACvC,IAAI,CAAC,UAAU,CAAC,KAAK,CAAC,CAAC;QACvB,OAAO,IAAI,CAAC;IACd,CAAC;IAED,uFAAuF;IAChF,+BAA+B,CAAC,EAAU,EAAE,EAAU,EAAE,MAAe,EAAE,MAAe,EAAE,QAAyD;QACxJ,IAAI,EAAE,GAAG,EAAE;YACT,OAAO,KAAK,CAAC;QACf,MAAM,EAAE,GAAG,CAAE,IAAI,CAAC,QAAQ,CAAC,MAAM,CAAC,CAAC;QACnC,MAAM,EAAE,GAAG,CAAE,IAAI,CAAC,QAAQ,CAAC,MAAM,CAAC,CAAC;QACnC,MAAM,KAAK,GAAG,EAAE,GAAG,EAAE,CAAC;QACtB,MAAM,CAAC,GAAG,QAAQ,CAAC,yBAAyB,CAAC,CAAC,EAAE,EAAE,KAAK,CAAC,CAAC;QACzD,IAAI,CAAC,KAAK,SAAS,EAAE,EAAE,wCAAwC;YAC7D,IAAI,EAAE,IAAI,GAAG,EAAE;gBAAE,IAAI,QAAQ;oBAAE,QAAQ,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC;gBAAC,OAAO,IAAI,CAAC;aAAE;YAC/D,OAAO,KAAK,CAAC;SACd;QACD,IAAI,KAAK,GAAG,CAAC,EAAE,EAAE,mBAAmB;YAClC,IAAI,CAAC,GAAG,EAAE;gBAAE,EAAE,GAAG,CAAC,CAAC;SACpB;aAAM;YACL,kBAAkB;YAClB,IAAI,CAAC,GAAG,EAAE;gBACR,EAAE,GAAG,CAAC,CAAC;SACV;QACD,IAAI,EAAE,GAAG,EAAE;YACT,OAAO,KAAK,CAAC;QACf,IAAI,QAAQ;YAAE,QAAQ,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC;QAC/B,OAAO,IAAI,CAAC;IACd,CAAC;IACD;;;;;OAKG;IACI,QAAQ;QACb,MAAM,CAAC,GAAG,IAAI,CAAC,mBAAmB,CAAC;QACnC,MAAM,MAAM,GAAG,OAAO,CAAC,MAAM,CAAC,IAAI,CAAC,aAAa,CAAC,CAAC,GAAG,CAAC,EAAE,IAAI,CAAC,aAAa,CAAC,CAAC,GAAG,CAAC,EAAE,IAAI,CAAC,aAAa,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;QAC5G,MAAM,MAAM,GAAG,QAAQ,CAAC,kBAAkB,CAAC,IAAI,CAAC,aAAa,EAAE,SAAS,CAAC,GAAG,CAAC,CAAC;QAC9E,OAAO,SAAS,CAAC,qBAAqB,CAAC,MAAM,EAAE,MAAM,CAAC,CAAC;IACzD,CAAC;IACD;;;;OAIG;IACI,cAAc,CAAC,KAAc,EAAE,kBAA2B,KAAK;QACpE,IAAI,KAAK,CAAC,MAAM;YACd,OAAO,SAAS,CAAC;QACnB,MAAM,OAAO,GAAG,KAAK,CAAC,OAAO,EAAE,CAAC;QAChC,MAAM,YAAY,GAAG,IAAI,CAAC,QAAQ,EAAE,CAAC;QACrC,YAAY,CAAC,kCAAkC,CAAC,OAAO,CAAC,CAAC;QACzD,MAAM,UAAU,GAAG,OAAO,CAAC,WAAW,CAAC,OAAO,CAAC,CAAC;QAChD,IAAI,UAAU,CAAC,GAAG,CAAC,CAAC,GAAG,UAAU,CAAC,IAAI,CAAC,CAAC,GAAG,GAAG;YAC5C,OAAO,SAAS,CAAC;QACnB,kCAAkC;QAClC,MAAM,MAAM,GAAG,IAAI,gBAAgB,EAAE,CAAC;QACtC,MAAM,CAAC,OAAO,CAAC,UAAU,CAAC,GAAG,CAAC,CAAC,EAAE,UAAU,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;QACtD,MAAM,CAAC,OAAO,CAAC,UAAU,CAAC,IAAI,CAAC,CAAC,EAAE,UAAU,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;QACvD,MAAM,CAAC,OAAO,CAAC,UAAU,CAAC,IAAI,CAAC,CAAC,EAAE,UAAU,CAAC,IAAI,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;QACxD,MAAM,CAAC,OAAO,CAAC,UAAU,CAAC,GAAG,CAAC,CAAC,EAAE,UAAU,CAAC,IAAI,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;QACvD,MAAM,CAAC,wBAAwB,CAAC,YAAY,CAAC,CAAC;QAC9C,8BAA8B,CAAC,kCAAkC,CAAC,KAAK,EAAE,MAAM,CAAC,CAAC;QACjF,IAAI,MAAM,CAAC,MAAM,KAAK,CAAC;YACrB,OAAO,SAAS,CAAC;QACnB,IAAI,eAAe;YACjB,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC;QACrB,OAAO,MAAM,CAAC;IAChB,CAAC;IACD;;;;;;OAMG;IACI,iBAAiB,CACtB,GAAqB,EACrB,eAAmC,EACnC,gBAAoC,EACpC,UAAiC;QACjC,MAAM,kBAAkB,GAAG,OAAO,CAAC,UAAU,EAAE,CAAC;QAChD,MAAM,SAAS,GAAG,UAAU,CAAC,aAAa,EAAE,CAAC;QAC7C,MAAM,UAAU,GAAG,UAAU,CAAC,aAAa,EAAE,CAAC;QAC9C,8BAA8B,CAAC,oCAAoC,CAAC,IAAI,EAAE,GAAG,EAAE,SAAS,EAAE,UAAU,EAAE,kBAAkB,CAAC,CAAC;QAC1H,aAAa,CAAC,aAAa,CAAC,SAAS,EAAE,CAAC,EAAE,eAAe,EAAE,UAAU,CAAC,CAAC;QACvE,aAAa,CAAC,aAAa,CAAC,UAAU,EAAE,CAAC,EAAE,gBAAgB,EAAE,UAAU,CAAC,CAAC;IAC3E,CAAC;;AA1Lc,+BAAqB,GAAG,IAAI,oBAAoB,EAAE,CAAC","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 CartesianGeometry\r\n */\r\n\r\nimport { Arc3d } from \"../curve/Arc3d\";\r\nimport { AnnounceNumberNumberCurvePrimitive } from \"../curve/CurvePrimitive\";\r\nimport { AxisOrder, Geometry, PlaneAltitudeEvaluator } from \"../Geometry\";\r\nimport { Angle } from \"../geometry3d/Angle\";\r\nimport { XYZProps } from \"../geometry3d/XYZProps\";\r\nimport { GrowableFloat64Array } from \"../geometry3d/GrowableFloat64Array\";\r\nimport { GrowableXYZArray } from \"../geometry3d/GrowableXYZArray\";\r\nimport { Matrix3d } from \"../geometry3d/Matrix3d\";\r\nimport { Plane3dByOriginAndUnitNormal } from \"../geometry3d/Plane3dByOriginAndUnitNormal\";\r\nimport { Point3d, Vector3d } from \"../geometry3d/Point3dVector3d\";\r\nimport { IndexedXYZCollectionPolygonOps } from \"../geometry3d/PolygonOps\";\r\nimport { Range1d, Range3d } from \"../geometry3d/Range\";\r\nimport { Transform } from \"../geometry3d/Transform\";\r\nimport { Matrix4d } from \"../geometry4d/Matrix4d\";\r\nimport { Point4d } from \"../geometry4d/Point4d\";\r\nimport { AnalyticRoots } from \"../numerics/Polynomials\";\r\nimport { Clipper, ClipUtilities, PolygonClipper } from \"./ClipUtils\";\r\nimport { GrowableXYZArrayCache } from \"../geometry3d/ReusableObjectCache\";\r\n\r\n/** Wire format describing a [[ClipPlane]].\r\n * If either [[normal]] or [[dist]] are omitted, defaults to a normal of [[Vector3d.unitZ]] and a distance of zero.\r\n * @public\r\n */\r\nexport interface ClipPlaneProps {\r\n  /** The plane's inward normal. */\r\n  normal?: XYZProps;\r\n  /** The plane's distance from the origin. */\r\n  dist?: number;\r\n  /** Defaults to `false`. */\r\n  invisible?: boolean;\r\n  /** Defaults to `false`. */\r\n  interior?: boolean;\r\n}\r\n\r\n/** A ClipPlane is a single plane represented as\r\n * * An inward unit normal (u,v,w)\r\n * * A signedDistance\r\n *\r\n * Hence\r\n * * The halfspace function evaluation for \"point\" [x,y,z,] is: ([x,y,z] DOT (u,v,w)l - signedDistance)\r\n * * POSITIVE values of the halfspace function are \"inside\"\r\n * * ZERO value of the halfspace function is \"on\"\r\n * * NEGATIVE value of the halfspace function is \"outside\"\r\n * * A representative point on the plane is (signedDistance*u, signedDistance * v, signedDistance *w)\r\n * * Given a point and inward normal, the signedDistance is (point DOT normal)\r\n * @public\r\n */\r\nexport class ClipPlane implements Clipper, PlaneAltitudeEvaluator, PolygonClipper {\r\n  private _inwardNormal: Vector3d;\r\n  /** Construct a parallel plane through the origin.\r\n   * * Move it to the actual position.\r\n   * * _distanceFromOrigin is the distance it moved, with the (inward) normal direction as positive\r\n   */\r\n  private _distanceFromOrigin: number;\r\n  private _invisible: boolean;\r\n  private _interior: boolean;\r\n\r\n  private constructor(normal: Vector3d, distance: number, invisible: boolean, interior: boolean) {\r\n    this._invisible = invisible;\r\n    this._interior = interior;\r\n    this._inwardNormal = normal;\r\n    this._distanceFromOrigin = distance;\r\n  }\r\n  /*\r\n    private safeSetXYZDistance(nx: number, ny: number, nz: number, d: number) {\r\n      this._inwardNormal.set(nx, ny, nz);\r\n      this._distanceFromOrigin = d;\r\n    }\r\n  */\r\n  /**\r\n   * Return true if all members are almostEqual to corresponding members of other.\r\n   * @param other clip plane to compare\r\n   */\r\n  public isAlmostEqual(other: ClipPlane): boolean {\r\n    return Geometry.isSameCoordinate(this._distanceFromOrigin, other._distanceFromOrigin)\r\n      && this._inwardNormal.isAlmostEqual(other._inwardNormal)\r\n      && this._interior === other._interior\r\n      && this._invisible === other._invisible;\r\n  }\r\n  /** return a cloned plane */\r\n  public clone(): ClipPlane {\r\n    const result = new ClipPlane(this._inwardNormal.clone(), this._distanceFromOrigin, this._invisible, this._interior);\r\n    return result;\r\n  }\r\n  /** return Return a cloned plane with coordinate data negated. */\r\n  public cloneNegated(): ClipPlane {\r\n    const plane = new ClipPlane(this._inwardNormal.clone(), this._distanceFromOrigin, this._invisible, this._interior);\r\n    plane.negateInPlace();\r\n    return plane;\r\n  }\r\n  /** Create a ClipPlane from Plane3dByOriginAndUnitNormal. */\r\n  public static createPlane(plane: Plane3dByOriginAndUnitNormal, invisible: boolean = false, interior: boolean = false, result?: ClipPlane): ClipPlane {\r\n    const distance = plane.getNormalRef().dotProduct(plane.getOriginRef());\r\n    if (result) {\r\n      result._invisible = invisible;\r\n      result._interior = interior;\r\n      result._inwardNormal = plane.getNormalRef().clone();\r\n      result._distanceFromOrigin = distance;\r\n      return result;\r\n    }\r\n    return new ClipPlane(plane.getNormalRef().clone(), distance, invisible, interior);\r\n  }\r\n  /**\r\n   * * Create a ClipPlane with direct normal and signedDistance.\r\n   * * The vector is normalized for storage.\r\n   */\r\n  public static createNormalAndDistance(normal: Vector3d, distance: number, invisible: boolean = false, interior: boolean = false, result?: ClipPlane): ClipPlane | undefined {\r\n    const normalized = normal.normalize();\r\n    if (normalized) {\r\n      if (result) {\r\n        result._invisible = invisible;\r\n        result._interior = interior;\r\n        result._inwardNormal = normalized;\r\n        result._distanceFromOrigin = distance;\r\n      }\r\n      return new ClipPlane(normalized, distance, invisible, interior);\r\n    }\r\n    return undefined;\r\n  }\r\n  /** Create a ClipPlane\r\n   * * \"normal\" is the inward normal of the plane. (It is internally normalized)\r\n   * * \"point\" is any point of the plane.\r\n   * * The stored distance for the plane is the dot product of the point with the normal (i.e. treat the point's xyz as a vector from the origin.)\r\n   */\r\n  public static createNormalAndPoint(normal: Vector3d, point: Point3d, invisible: boolean = false, interior: boolean = false, result?: ClipPlane): ClipPlane | undefined {\r\n    const normalized = normal.normalize();\r\n    if (normalized) {\r\n      const distance = normalized.dotProduct(point);\r\n      if (result) {\r\n        result._invisible = invisible;\r\n        result._interior = interior;\r\n        result._inwardNormal = normalized;\r\n        result._distanceFromOrigin = distance;\r\n      }\r\n      return new ClipPlane(normalized, distance, invisible, interior);\r\n    }\r\n    return undefined;\r\n  }\r\n\r\n  /** Create a ClipPlane\r\n   * * \"normal\" (normalX, normalY, nz) is the inward normal of the plane.\r\n   * * The given (normalX,normalY,normalZ)\r\n   * * \"point\" is any point of the plane.\r\n   * * The stored distance for the plane is the dot product of the point with the normal (i.e. treat the point's xyz as a vector from the origin.)\r\n   */\r\n  public static createNormalAndPointXYZXYZ(normalX: number, normalY: number, normalZ: number,\r\n    originX: number, originY: number, originZ: number,\r\n    invisible: boolean = false, interior: boolean = false, result?: ClipPlane): ClipPlane | undefined {\r\n    const q = Geometry.hypotenuseXYZ(normalX, normalY, normalZ);\r\n    const r = Geometry.conditionalDivideFraction(1, q);\r\n    if (r !== undefined) {\r\n      if (result) {\r\n        result._inwardNormal.set(normalX * r, normalY * r, normalZ * r);\r\n        result._distanceFromOrigin = result._inwardNormal.dotProductXYZ(originX, originY, originZ);\r\n        result._invisible = invisible;\r\n        result._interior = interior;\r\n        return result;\r\n      }\r\n      const normal = Vector3d.create(normalX * r, normalY * r, normalZ * r);\r\n      return new ClipPlane(normal, normal.dotProductXYZ(originX, originY, originZ), invisible, interior);\r\n    }\r\n    return undefined;\r\n  }\r\n\r\n  /**\r\n   * return a json object of the form\r\n   * `{\"normal\":[u,v,w],\"dist\":signedDistanceValue,\"interior\":true,\"invisible\":true}`\r\n   */\r\n  public toJSON(): ClipPlaneProps {\r\n    const props: ClipPlaneProps = {\r\n      normal: this.inwardNormalRef.toJSON(),\r\n      dist: this.distance,\r\n    };\r\n\r\n    if (this.interior)\r\n      props.interior = true;\r\n\r\n    if (this.invisible)\r\n      props.invisible = true;\r\n\r\n    return props;\r\n  }\r\n\r\n  /** parse json object to ClipPlane instance */\r\n  public static fromJSON(json: ClipPlaneProps, result?: ClipPlane): ClipPlane | undefined {\r\n    if (json && json.normal && undefined !== json.dist && Number.isFinite(json.dist))\r\n      return ClipPlane.createNormalAndDistance(Vector3d.fromJSON(json.normal), json.dist, !!json.invisible, !!json.interior);\r\n\r\n    return ClipPlane.createNormalAndDistance(Vector3d.unitZ(), 0, false, false, result);\r\n  }\r\n\r\n  /** Set both the invisible and interior flags. */\r\n  public setFlags(invisible: boolean, interior: boolean) {\r\n    this._invisible = invisible;\r\n    this._interior = interior;\r\n  }\r\n\r\n  /**\r\n   * Return the stored distanceFromOrigin property.\r\n   */\r\n  public get distance() { return this._distanceFromOrigin; }\r\n  /**\r\n   * Return the stored inward normal property.\r\n   */\r\n  public get inwardNormalRef(): Vector3d { return this._inwardNormal; }\r\n  /**\r\n   * Return the \"interior\" property bit\r\n   */\r\n  public get interior() { return this._interior; }\r\n  /**\r\n   * Return the \"invisible\" property bit.\r\n   */\r\n  public get invisible() { return this._invisible; }\r\n\r\n  /**\r\n   * Create a plane defined by two points, an up vector, and a tilt angle relative to the up vector.\r\n   * @param point0 start point of the edge\r\n   * @param point1 end point of the edge\r\n   * @param upVector vector perpendicular to the plane\r\n   * @param tiltAngle angle to tilt the plane around the edge in the direction of the up vector.\r\n   * @param result optional preallocated plane\r\n   */\r\n  public static createEdgeAndUpVector(point0: Point3d, point1: Point3d, upVector: Vector3d, tiltAngle?: Angle, result?: ClipPlane): ClipPlane | undefined {\r\n    const edgeVector = Vector3d.createFrom(point1.minus(point0));\r\n    let normal = (upVector.crossProduct(edgeVector)).normalize();\r\n\r\n    if (normal) {\r\n      if (tiltAngle !== undefined && !tiltAngle.isAlmostZero) {\r\n        const tiltNormal = Vector3d.createRotateVectorAroundVector(normal, edgeVector, tiltAngle);\r\n        if (tiltNormal) {\r\n          normal = tiltNormal.clone();\r\n        }\r\n      }\r\n      normal.negate(normal);\r\n      return ClipPlane.createNormalAndPoint(normal, point0, false, false, result);\r\n    }\r\n    return undefined;\r\n  }\r\n  /**\r\n   * Create a plane perpendicular to the edge between the xy parts of point0 and point1\r\n   */\r\n  public static createEdgeXY(point0: Point3d, point1: Point3d, result?: ClipPlane): ClipPlane | undefined {\r\n    const normal = Vector3d.create(point0.y - point1.y, point1.x - point0.x);\r\n    if (normal.normalizeInPlace())\r\n      return ClipPlane.createNormalAndPoint(normal, point0, false, false, result);\r\n    return undefined;\r\n  }\r\n  /**\r\n   * Return the Plane3d form of the plane.\r\n   * * The plane origin is the point `distance * inwardNormal`\r\n   * * The plane normal is the inward normal of the ClipPlane.\r\n   */\r\n  public getPlane3d(): Plane3dByOriginAndUnitNormal {\r\n    const d = this._distanceFromOrigin;\r\n    // Normal should be normalized, will not return undefined\r\n    return Plane3dByOriginAndUnitNormal.create(Point3d.create(this._inwardNormal.x * d, this._inwardNormal.y * d, this._inwardNormal.z * d), this._inwardNormal)!;\r\n  }\r\n\r\n  /**\r\n   * Return the Point4d d form of the plane.\r\n   * * The homogeneous xyz are the inward normal xyz.\r\n   * * The homogeneous weight is the negated ClipPlane distance.\r\n   */\r\n  public getPlane4d(): Point4d {\r\n    return Point4d.create(this._inwardNormal.x, this._inwardNormal.y, this._inwardNormal.z, - this._distanceFromOrigin);\r\n  }\r\n  /**\r\n   * Set the plane from DPoint4d style plane.\r\n   * * The saved plane has its direction normalized.\r\n   * * This preserves the plane itself as a zero set but make plane evaluations act as true distances (even if the plane coefficients are scaled otherwise)\r\n   * @param plane\r\n   */\r\n  public setPlane4d(plane: Point4d) {\r\n    const a = Math.sqrt(plane.x * plane.x + plane.y * plane.y + plane.z * plane.z);\r\n    const r = a === 0.0 ? 1.0 : 1.0 / a;\r\n    this._inwardNormal.x = r * plane.x;\r\n    this._inwardNormal.y = r * plane.y;\r\n    this._inwardNormal.z = r * plane.z;\r\n    this._distanceFromOrigin = -r * plane.w;\r\n  }\r\n\r\n  /**\r\n   * Evaluate the altitude in weighted space, i.e. (dot product with inward normal) minus distance, with point.w scale applied to distance)\r\n   * @param point space point to test\r\n   */\r\n  public weightedAltitude(point: Point4d): number {\r\n    return point.x * this._inwardNormal.x + point.y * this._inwardNormal.y + point.z * this._inwardNormal.z - point.w * this._distanceFromOrigin;\r\n  }\r\n\r\n  /**\r\n   * Evaluate the distance from the plane to a point in space, i.e. (dot product with inward normal) minus distance\r\n   * @param point space point to test\r\n   */\r\n  public altitude(point: Point3d): number {\r\n    return point.x * this._inwardNormal.x + point.y * this._inwardNormal.y + point.z * this._inwardNormal.z - this._distanceFromOrigin;\r\n  }\r\n\r\n  /**\r\n   * Evaluate the distance from the plane to a point in space with point given as x,y,z, i.e. (dot product with inward normal) minus distance\r\n   * @param point space point to test\r\n   */\r\n  public altitudeXYZ(x: number, y: number, z: number): number {\r\n    return x * this._inwardNormal.x + y * this._inwardNormal.y + z * this._inwardNormal.z - this._distanceFromOrigin;\r\n  }\r\n\r\n  /** Return the dot product of the plane normal with the vector (NOT using the plane's distanceFromOrigin).\r\n   */\r\n  public velocity(vector: Vector3d): number {\r\n    return vector.x * this._inwardNormal.x + vector.y * this._inwardNormal.y + vector.z * this._inwardNormal.z;\r\n  }\r\n\r\n  /** Return the dot product of the plane normal with the x,yz, vector components (NOT using the plane's distanceFromOrigin).\r\n   */\r\n  public velocityXYZ(x: number, y: number, z: number): number {\r\n    return x * this._inwardNormal.x + y * this._inwardNormal.y + z * this._inwardNormal.z;\r\n  }\r\n\r\n  /** Return the dot product of the plane normal with the point (treating the point xyz as a vector, and NOT using the plane's distanceFromOrigin).\r\n   */\r\n  public dotProductPlaneNormalPoint(point: Point3d): number {\r\n    return point.x * this._inwardNormal.x + point.y * this._inwardNormal.y + point.z * this._inwardNormal.z;\r\n  }\r\n  /**\r\n   * Return true if spacePoint is inside or on the plane, with tolerance applied to \"on\".\r\n   * @param spacePoint point to test.\r\n   * @param tolerance tolerance for considering \"near plane\" to be \"on plane\"\r\n   */\r\n  public isPointOnOrInside(spacePoint: Point3d, tolerance: number = Geometry.smallMetricDistance): boolean {\r\n    let value = this.altitude(spacePoint);\r\n    if (tolerance) { value += tolerance; }\r\n    return value >= 0.0;\r\n  }\r\n\r\n  /**\r\n   * Return true if spacePoint is strictly inside the halfspace, with tolerance applied to \"on\".\r\n   * @param spacePoint point to test.\r\n   * @param tolerance tolerance for considering \"near plane\" to be \"on plane\"\r\n   */\r\n  public isPointInside(point: Point3d, tolerance: number = Geometry.smallMetricDistance): boolean {\r\n    let value = this.altitude(point);\r\n    if (tolerance) { value -= tolerance; }\r\n    return value > 0.0;\r\n  }\r\n\r\n  /**\r\n   * Return true if spacePoint is strictly on the plane, within tolerance\r\n   * @param spacePoint point to test.\r\n   * @param tolerance tolerance for considering \"near plane\" to be \"on plane\"\r\n   */\r\n  public isPointOn(point: Point3d, tolerance: number = Geometry.smallMetricDistance): boolean {\r\n    return Math.abs(this.altitude(point)) <= tolerance;\r\n  }\r\n  /**\r\n   * Compute intersections of an (UNBOUNDED) arc with the plane.  Append them (as radians) to a growing array.\r\n   * @param arc arc to test.  The angle limits of the arc are NOT considered.\r\n   * @param intersectionRadians array to receive results\r\n   */\r\n  public appendIntersectionRadians(arc: Arc3d, intersectionRadians: GrowableFloat64Array) {\r\n    const arcVectors = arc.toVectors();\r\n    const alpha = this.altitude(arc.center);\r\n    const beta = this.velocity(arcVectors.vector0);\r\n    const gamma = this.velocity(arcVectors.vector90);\r\n    AnalyticRoots.appendImplicitLineUnitCircleIntersections(alpha, beta, gamma, undefined, undefined, intersectionRadians);\r\n  }\r\n\r\n  private static _clipArcFractionArray = new GrowableFloat64Array();\r\n  /** Announce fractional intervals of arc clip.\r\n   * * Each call to `announce(fraction0, fraction1, arc)` announces one interval that is inside the clip plane.\r\n   */\r\n  public announceClippedArcIntervals(arc: Arc3d, announce?: AnnounceNumberNumberCurvePrimitive): boolean {\r\n    const breaks = ClipPlane._clipArcFractionArray;\r\n    breaks.clear();\r\n    this.appendIntersectionRadians(arc, breaks);\r\n    arc.sweep.radiansArrayToPositivePeriodicFractions(breaks);\r\n    return ClipUtilities.selectIntervals01(arc, breaks, this, announce);\r\n  }\r\n  /**\r\n   * * Compute intersection of (unbounded) segment with the plane.\r\n   * * If the ends are on the same side of the plane, return undefined.\r\n   * * If the intersection is an endpoint or interior to the segment return the fraction.\r\n   * * If both ends are on, return undefined.\r\n   */\r\n  public getBoundedSegmentSimpleIntersection(pointA: Point3d, pointB: Point3d): number | undefined {\r\n    const h0 = this.altitude(pointA);\r\n    const h1 = this.altitude(pointB);\r\n    if (h0 * h1 > 0.0)\r\n      return undefined;\r\n    if (h0 === 0.0 && h1 === 0.0) {\r\n      return undefined;\r\n    }\r\n    return - h0 / (h1 - h0);\r\n  }\r\n\r\n  /** Apply transform to the origin.  Apply inverse transpose of the matrix part to th normal vector. */\r\n  public transformInPlace(transform: Transform): boolean {\r\n    const plane: Plane3dByOriginAndUnitNormal = this.getPlane3d();\r\n    const matrix: Matrix3d = transform.matrix;\r\n    const newPoint = transform.multiplyPoint3d(plane.getOriginRef());\r\n    // Normal transforms as the inverse transpose of the matrix part\r\n    // BTW: If the matrix is orthogonal, this is a long way to multiply by the matrix part (mumble grumble)\r\n    const newNormal = matrix.multiplyInverseTranspose(plane.getNormalRef());\r\n    if (!newNormal)\r\n      return false;\r\n\r\n    plane.set(newPoint, newNormal);\r\n    const normalized = (plane.getNormalRef()).normalize();\r\n    if (!normalized)\r\n      return false;\r\n    this._inwardNormal = normalized;\r\n    this._distanceFromOrigin = this._inwardNormal.dotProduct(plane.getOriginRef());\r\n    return true;\r\n  }\r\n  /** Set the invisible flag.   Interpretation of this is up to the use code algorithms. */\r\n  public setInvisible(invisible: boolean) {\r\n    this._invisible = invisible;\r\n  }\r\n\r\n  /**  reverse the sign of all coefficients, so outside and inside reverse */\r\n  public negateInPlace() {\r\n    this._inwardNormal = this._inwardNormal.negate();\r\n    this._distanceFromOrigin = - this._distanceFromOrigin;\r\n  }\r\n  /**\r\n   * Move the plane INWARD by given distance\r\n   * @param offset distance of shift inwards\r\n   */\r\n  public offsetDistance(offset: number) {\r\n    this._distanceFromOrigin += offset;\r\n  }\r\n\r\n  /**\r\n   * Clip a polygon to the inside or outside of the plane.\r\n   * * Results with 2 or fewer points are ignored.\r\n   * * Other than ensuring capacity in the arrays, there are no object allocations during execution of this function.\r\n   * @param xyz input points.\r\n   * @param work work buffer\r\n   * @param tolerance tolerance for \"on plane\" decision.\r\n   */\r\n  public clipConvexPolygonInPlace(xyz: GrowableXYZArray, work: GrowableXYZArray, inside: boolean = true, tolerance: number = Geometry.smallMetricDistance) {\r\n    return IndexedXYZCollectionPolygonOps.clipConvexPolygonInPlace(this, xyz, work, inside, tolerance);\r\n  }\r\n\r\n  /**\r\n   * Multiply the ClipPlane's DPoint4d by matrix.\r\n   * @param matrix matrix to apply.\r\n   * @param invert if true, use in verse of the matrix.\r\n   * @param transpose if true, use the transpose of the matrix (or inverse, per invert parameter)\r\n   * * Note that if matrixA is applied to all of space, the matrix to send to this method to get a corresponding effect on the plane is the inverse transpose of matrixA\r\n   * * Callers that will apply the same matrix to many planes should pre-invert the matrix for efficiency.\r\n   * * Both params default to true to get the full effect of transforming space.\r\n   * @param matrix matrix to apply\r\n   * @return false if unable to invert\r\n   */\r\n  public multiplyPlaneByMatrix4d(matrix: Matrix4d, invert: boolean = true, transpose: boolean = true): boolean {\r\n    const plane: Point4d = this.getPlane4d();\r\n    if (invert) {\r\n      const inverse = matrix.createInverse();\r\n      if (inverse)\r\n        return this.multiplyPlaneByMatrix4d(inverse, false, transpose);\r\n      return false;\r\n    }\r\n    if (transpose)\r\n      matrix.multiplyTransposePoint4d(plane, plane);\r\n    else\r\n      matrix.multiplyPoint4d(plane, plane);\r\n    this.setPlane4d(plane);\r\n    return true;\r\n  }\r\n\r\n  /** announce the interval (if any) where a line is within the clip plane half space. */\r\n  public announceClippedSegmentIntervals(f0: number, f1: number, pointA: Point3d, pointB: Point3d, announce?: (fraction0: number, fraction1: number) => void): boolean {\r\n    if (f1 < f0)\r\n      return false;\r\n    const h0 = - this.altitude(pointA);\r\n    const h1 = - this.altitude(pointB);\r\n    const delta = h1 - h0;\r\n    const f = Geometry.conditionalDivideFraction(-h0, delta);\r\n    if (f === undefined) { // The segment is parallel to the plane.\r\n      if (h0 <= 0.0) { if (announce) announce(f0, f1); return true; }\r\n      return false;\r\n    }\r\n    if (delta > 0) { // segment aims OUT\r\n      if (f < f1) f1 = f;\r\n    } else {\r\n      // segment aims IN\r\n      if (f > f0)\r\n        f0 = f;\r\n    }\r\n    if (f1 < f0)\r\n      return false;\r\n    if (announce) announce(f0, f1);\r\n    return true;\r\n  }\r\n  /**\r\n   * Return a coordinate frame with\r\n   * * origin at closest point to global origin\r\n   * * z axis points in\r\n   * x and y are \"in plane\"\r\n   */\r\n  public getFrame(): Transform {\r\n    const d = this._distanceFromOrigin;\r\n    const origin = Point3d.create(this._inwardNormal.x * d, this._inwardNormal.y * d, this._inwardNormal.z * d);\r\n    const matrix = Matrix3d.createRigidHeadsUp(this._inwardNormal, AxisOrder.ZXY);\r\n    return Transform.createOriginAndMatrix(origin, matrix);\r\n  }\r\n  /**\r\n   * Return the intersection of the plane with a range cube.\r\n   * @param range\r\n   * @param xyzOut intersection polygon.  This is convex.\r\n   */\r\n  public intersectRange(range: Range3d, addClosurePoint: boolean = false): GrowableXYZArray | undefined {\r\n    if (range.isNull)\r\n      return undefined;\r\n    const corners = range.corners();\r\n    const frameOnPlane = this.getFrame();\r\n    frameOnPlane.multiplyInversePoint3dArrayInPlace(corners);\r\n    const localRange = Range3d.createArray(corners);\r\n    if (localRange.low.z * localRange.high.z > 0.0)\r\n      return undefined;\r\n    // oversized polygon on local z= 0\r\n    const xyzOut = new GrowableXYZArray();\r\n    xyzOut.pushXYZ(localRange.low.x, localRange.low.y, 0);\r\n    xyzOut.pushXYZ(localRange.high.x, localRange.low.y, 0);\r\n    xyzOut.pushXYZ(localRange.high.x, localRange.high.y, 0);\r\n    xyzOut.pushXYZ(localRange.low.x, localRange.high.y, 0);\r\n    xyzOut.multiplyTransformInPlace(frameOnPlane);\r\n    IndexedXYZCollectionPolygonOps.intersectRangeConvexPolygonInPlace(range, xyzOut);\r\n    if (xyzOut.length === 0)\r\n      return undefined;\r\n    if (addClosurePoint)\r\n      xyzOut.pushWrap(1);\r\n    return xyzOut;\r\n  }\r\n  /** Implement appendPolygonClip, as defined in interface PolygonClipper.  /**\r\n   *\r\n   * @param xyz input polygon.  This is not changed.\r\n   * @param insideFragments Array to receive \"inside\" fragments.  Each fragment is a GrowableXYZArray grabbed from the cache.  This is NOT cleared.\r\n   * @param outsideFragments Array to receive \"outside\" fragments.  Each fragment is a GrowableXYZArray grabbed from the cache.  This is NOT cleared.\r\n   * @param arrayCache cache for reusable GrowableXYZArray.\r\n   */\r\n  public appendPolygonClip(\r\n    xyz: GrowableXYZArray,\r\n    insideFragments: GrowableXYZArray[],\r\n    outsideFragments: GrowableXYZArray[],\r\n    arrayCache: GrowableXYZArrayCache): void {\r\n    const perpendicularRange = Range1d.createNull();\r\n    const newInside = arrayCache.grabFromCache();\r\n    const newOutside = arrayCache.grabFromCache();\r\n    IndexedXYZCollectionPolygonOps.splitConvexPolygonInsideOutsidePlane(this, xyz, newInside, newOutside, perpendicularRange);\r\n    ClipUtilities.captureOrDrop(newInside, 3, insideFragments, arrayCache);\r\n    ClipUtilities.captureOrDrop(newOutside, 3, outsideFragments, arrayCache);\r\n  }\r\n}\r\n"]}