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

package/lib/esm/geometry3d/Ray3d.js

@@ -0,0 +1,325 @@
+/*---------------------------------------------------------------------------------------------
+* Copyright (c) Bentley Systems, Incorporated. All rights reserved.
+* See LICENSE.md in the project root for license terms and full copyright notice.
+*--------------------------------------------------------------------------------------------*/
+/** @packageDocumentation
+ * @module CartesianGeometry
+ */
+import { CurveCurveApproachType, CurveLocationDetail, CurveLocationDetailPair } from "../curve/CurveLocationDetail";
+import { AxisOrder, Geometry } from "../Geometry";
+import { SmallSystem } from "../numerics/Polynomials";
+import { Matrix3d } from "./Matrix3d";
+import { Vector2d } from "./Point2dVector2d";
+import { Point3d, Vector3d } from "./Point3dVector3d";
+import { Range1d } from "./Range";
+import { Transform } from "./Transform";
+/** A Ray3d contains
+ * * an origin point.
+ * * a direction vector.  The vector is NOT required to be normalized.
+ *  * an optional weight (number).
+ * @public
+ */
+export class Ray3d {
+    // constructor captures references !!!
+    constructor(origin, direction) {
+        this.origin = origin;
+        this.direction = direction;
+        this.a = undefined;
+    }
+    static _create(x, y, z, u, v, w) {
+        return new Ray3d(Point3d.create(x, y, z), Vector3d.create(u, v, w));
+    }
+    /** Create a ray on the x axis. */
+    static createXAxis() { return Ray3d._create(0, 0, 0, 1, 0, 0); }
+    /** Create a ray on the y axis. */
+    static createYAxis() { return Ray3d._create(0, 0, 0, 0, 1, 0); }
+    /** Create a ray on the z axis. */
+    static createZAxis() { return Ray3d._create(0, 0, 0, 0, 0, 1); }
+    /** Create a ray with all zeros. */
+    static createZero(result) {
+        if (result) {
+            result.origin.setZero();
+            result.direction.setZero();
+            return result;
+        }
+        return new Ray3d(Point3d.createZero(), Vector3d.createZero());
+    }
+    /** Test for nearly equal rays. */
+    isAlmostEqual(other) {
+        return this.origin.isAlmostEqual(other.origin) && this.direction.isAlmostEqual(other.direction);
+    }
+    /** Create a ray from origin and direction. */
+    static create(origin, direction, result) {
+        if (result) {
+            result.set(origin, direction);
+            return result;
+        }
+        return new Ray3d(origin.clone(), direction.clone());
+    }
+    /**
+     * Given a homogeneous point and its derivative components, construct a Ray3d with cartesian coordinates and derivatives.
+     * @param weightedPoint `[x,y,z,w]` parts of weighted point.
+     * @param weightedDerivative `[x,y,z,w]` derivatives
+     * @param result
+     */
+    static createWeightedDerivative(weightedPoint, weightedDerivative, result) {
+        const w = weightedPoint[3];
+        const dw = weightedDerivative[3];
+        const x = weightedPoint[0];
+        const y = weightedPoint[1];
+        const z = weightedPoint[2];
+        const dx = weightedDerivative[0] * w - weightedPoint[0] * dw;
+        const dy = weightedDerivative[1] * w - weightedPoint[1] * dw;
+        const dz = weightedDerivative[2] * w - weightedPoint[2] * dw;
+        if (Geometry.isSmallMetricDistance(w))
+            return undefined;
+        const divW = 1.0 / w;
+        const divWW = divW * divW;
+        return Ray3d.createXYZUVW(x * divW, y * divW, z * divW, dx * divWW, dy * divWW, dz * divWW, result);
+    }
+    /** Create from coordinates of the origin and direction. */
+    static createXYZUVW(originX, originY, originZ, directionX, directionY, directionZ, result) {
+        if (result) {
+            result.getOriginRef().set(originX, originY, originZ);
+            result.getDirectionRef().set(directionX, directionY, directionZ);
+            return result;
+        }
+        return new Ray3d(Point3d.create(originX, originY, originZ), Vector3d.create(directionX, directionY, directionZ));
+    }
+    /** Capture origin and direction in a new Ray3d. */
+    static createCapture(origin, direction) {
+        return new Ray3d(origin, direction);
+    }
+    /** Create from (clones of) origin, direction, and numeric weight. */
+    static createPointVectorNumber(origin, direction, a, result) {
+        if (result) {
+            result.origin.setFrom(origin);
+            result.direction.setFrom(direction);
+            result.a = a;
+            return result;
+        }
+        result = new Ray3d(origin.clone(), direction.clone());
+        result.a = a;
+        return result;
+    }
+    /** Create from origin and target.  The direction vector is the full length (non-unit) vector from origin to target. */
+    static createStartEnd(origin, target, result) {
+        if (result) {
+            result.origin.setFrom(origin);
+            result.direction.setStartEnd(origin, target);
+            return result;
+        }
+        return new Ray3d(origin, Vector3d.createStartEnd(origin, target));
+    }
+    /** Return a reference to the ray's origin. */
+    getOriginRef() { return this.origin; }
+    /** Return a reference to the ray's direction vector. */
+    getDirectionRef() { return this.direction; }
+    /** copy coordinates from origin and direction. */
+    set(origin, direction) {
+        this.origin.setFrom(origin);
+        this.direction.setFrom(direction);
+    }
+    /** Clone the ray. */
+    clone(result) {
+        if (result) {
+            result.set(this.origin.clone(), this.direction.clone());
+            return result;
+        }
+        return new Ray3d(this.origin.clone(), this.direction.clone());
+    }
+    /** Create a clone and return the transform of the clone. */
+    cloneTransformed(transform) {
+        return new Ray3d(transform.multiplyPoint3d(this.origin), transform.multiplyVector(this.direction));
+    }
+    /** Create a clone and return the inverse transform of the clone. */
+    cloneInverseTransformed(transform) {
+        const origin = transform.multiplyInversePoint3d(this.origin);
+        const direction = transform.matrix.multiplyInverseXYZAsVector3d(this.direction.x, this.direction.y, this.direction.z);
+        if (undefined !== origin && undefined !== direction)
+            return new Ray3d(origin, direction);
+        return undefined;
+    }
+    /** Apply a transform in place. */
+    transformInPlace(transform) {
+        transform.multiplyPoint3d(this.origin, this.origin);
+        transform.multiplyVector(this.direction, this.direction);
+    }
+    /** Copy data from another ray. */
+    setFrom(source) { this.set(source.origin, source.direction); }
+    /** * fraction 0 is the ray origin.
+     * * fraction 1 is at the end of the direction vector when placed at the origin.
+     * @returns Return a point at fractional position along the ray.
+     */
+    fractionToPoint(fraction, result) { return this.origin.plusScaled(this.direction, fraction, result); }
+    /** Return the dot product of the ray's direction vector with a vector from the ray origin to the space point. */
+    dotProductToPoint(spacePoint) { return this.direction.dotProductStartEnd(this.origin, spacePoint); }
+    /**
+     * Return the fractional coordinate (along the direction vector) of the spacePoint projected to the ray.
+     */
+    pointToFraction(spacePoint) {
+        return Geometry.safeDivideFraction(this.direction.dotProductStartEnd(this.origin, spacePoint), this.direction.magnitudeSquared(), 0);
+    }
+    /**
+     *
+     * Return the spacePoint projected onto the ray.
+     */
+    projectPointToRay(spacePoint) {
+        return this.origin.plusScaled(this.direction, this.pointToFraction(spacePoint));
+    }
+    /** Return a transform for rigid axes
+     * at ray origin with z in ray direction.  If the direction vector is zero, axes default to identity (from createHeadsUpTriad)
+     */
+    toRigidZFrame() {
+        const axes = Matrix3d.createRigidHeadsUp(this.direction, AxisOrder.ZXY);
+        return Transform.createOriginAndMatrix(this.origin, axes);
+    }
+    /**
+     * Convert {origin:[x,y,z], direction:[u,v,w]} to a Ray3d.
+     */
+    setFromJSON(json) {
+        if (!json) {
+            this.origin.set(0, 0, 0);
+            this.direction.set(0, 0, 1);
+            return;
+        }
+        this.origin.setFromJSON(json.origin);
+        this.direction.setFromJSON(json.direction);
+    }
+    /**
+     * try to scale the direction vector to a given magnitude.
+     * @returns Returns false if ray direction is a zero vector.
+     */
+    trySetDirectionMagnitudeInPlace(magnitude = 1.0) {
+        if (this.direction.tryNormalizeInPlace()) {
+            this.direction.scaleInPlace(magnitude);
+            return true;
+        }
+        this.direction.setZero();
+        this.a = 0.0;
+        return false;
+    }
+    /**
+     * * If parameter `a` is clearly nonzero and the direction vector can be normalized,
+     *    * save the parameter `a` as the optional `a` member of the ray.
+     *    * normalize the ray's direction vector
+     * * If parameter `a` is nearly zero,
+     *   * Set the `a` member to zero
+     *   * Set the ray's direction vector to zero.
+     * @param a area to be saved.
+     */
+    // input a ray and "a" understood as an area.
+    // if a is clearly nonzero metric squared and the vector can be normalized, install those and return true.
+    // otherwise set ray.z to zero and zero the vector of the ray and return false.
+    tryNormalizeInPlaceWithAreaWeight(a) {
+        const tolerance = Geometry.smallMetricDistanceSquared;
+        this.a = a;
+        if (Math.abs(a) > tolerance && this.direction.tryNormalizeInPlace(tolerance))
+            return true;
+        this.direction.setZero();
+        this.a = 0.0;
+        return false;
+    }
+    /**
+     * Convert an Angle to a JSON object.
+     * @return {*} [origin,normal]
+     */
+    toJSON() { return { origin: this.origin.toJSON(), direction: this.direction.toJSON() }; }
+    /** Create a new ray from json object.  See `setFromJSON` for json structure; */
+    static fromJSON(json) {
+        const result = Ray3d.createXAxis();
+        result.setFromJSON(json);
+        return result;
+    }
+    /** return distance from the ray to point in space */
+    distance(spacePoint) {
+        const uu = this.direction.magnitudeSquared();
+        const uv = this.dotProductToPoint(spacePoint);
+        const aa = Geometry.inverseMetricDistanceSquared(uu);
+        if (aa)
+            return Math.sqrt(this.origin.distanceSquared(spacePoint) - uv * uv * aa);
+        else
+            return Math.sqrt(this.origin.distanceSquared(spacePoint));
+    }
+    /**
+     * Return the intersection of the unbounded ray with a plane.
+     * Stores the point of intersection in the result point given as a parameter,
+     * and returns the parameter along the ray where the intersection occurs.
+     * Returns undefined if the ray and plane are parallel or coplanar.
+     */
+    intersectionWithPlane(plane, result) {
+        const vectorA = Vector3d.createStartEnd(plane.getOriginRef(), this.origin);
+        const uDotN = this.direction.dotProduct(plane.getNormalRef());
+        const nDotN = this.direction.magnitudeSquared();
+        const aDotN = vectorA.dotProduct(plane.getNormalRef());
+        const division = Geometry.conditionalDivideFraction(-aDotN, uDotN);
+        if (undefined === division)
+            return undefined;
+        const division1 = Geometry.conditionalDivideFraction(nDotN, uDotN);
+        if (undefined === division1)
+            return undefined;
+        if (result) {
+            this.origin.plusScaled(this.direction, division, result);
+        }
+        return division;
+    }
+    /**
+     * * Find intersection of the ray with a Range3d.
+     * * return the range of fractions (on the ray) which are "inside" the range.
+     * * Note that a range is always returned;  if there is no intersection it is indicated by the test `result.sNull`
+     */
+    intersectionWithRange3d(range, result) {
+        if (range.isNull)
+            return Range1d.createNull(result);
+        const interval = Range1d.createXX(-Geometry.largeCoordinateResult, Geometry.largeCoordinateResult, result);
+        if (interval.clipLinearMapToInterval(this.origin.x, this.direction.x, range.low.x, range.high.x)
+            && interval.clipLinearMapToInterval(this.origin.y, this.direction.y, range.low.y, range.high.y)
+            && interval.clipLinearMapToInterval(this.origin.z, this.direction.z, range.low.z, range.high.z))
+            return interval;
+        return interval;
+    }
+    /** Construct a vector from `ray.origin` to target point.
+     * * return the part of the vector that is perpendicular to `ray.direction`.
+     *  * i.e. return the shortest vector from the ray to the point.
+     */
+    perpendicularPartOfVectorToTarget(targetPoint, result) {
+        const vectorV = Vector3d.createStartEnd(this.origin, targetPoint);
+        const uu = this.direction.magnitudeSquared();
+        const uv = this.direction.dotProductStartEnd(this.origin, targetPoint);
+        const fraction = Geometry.safeDivideFraction(uv, uu, 0.0);
+        return vectorV.plusScaled(this.direction, -fraction, result);
+    }
+    /** Determine if two rays intersect, are fully overlapped, parallel but no coincident, or skew
+     * * Return a CurveLocationDetailPair which
+     * * contains fraction and point on each ray.
+     * * has (in the CurveLocationDetailPair structure, as member approachType) annotation indicating one of these relationships
+     *   * CurveCurveApproachType.Intersection -- the rays have a simple intersection, at fractions indicated in detailA and detailB
+     *   * CurveCurveApproachType.PerpendicularChord -- there is pair of where the rays have closest approach.  The rays are skew in space.
+     *   * CurveCurveApproachType.CoincidentGeometry -- the rays are the same unbounded line in space. The fractions and points are a representative single common point.
+     *   * CurveCurveApproachType.Parallel -- the rays are parallel (and not coincident).   The two points are at the minimum distance
+     */
+    static closestApproachRay3dRay3d(rayA, rayB) {
+        const intersectionFractions = Vector2d.create();
+        let fractionA, fractionB;
+        let pointA, pointB;
+        let pairType;
+        if (SmallSystem.ray3dXYZUVWClosestApproachUnbounded(rayA.origin.x, rayA.origin.y, rayA.origin.z, rayA.direction.x, rayA.direction.y, rayA.direction.z, rayB.origin.x, rayB.origin.y, rayB.origin.z, rayB.direction.x, rayB.direction.y, rayB.direction.z, intersectionFractions)) {
+            fractionA = intersectionFractions.x;
+            fractionB = intersectionFractions.y;
+            pointA = rayA.fractionToPoint(fractionA);
+            pointB = rayB.fractionToPoint(fractionB);
+            pairType = pointA.isAlmostEqualMetric(pointB) ? CurveCurveApproachType.Intersection : CurveCurveApproachType.PerpendicularChord;
+        }
+        else {
+            fractionB = 0.0;
+            fractionA = rayA.pointToFraction(rayB.origin);
+            pointA = rayA.fractionToPoint(fractionA);
+            pointB = rayB.fractionToPoint(fractionB);
+            pairType = pointA.isAlmostEqualMetric(pointB) ? CurveCurveApproachType.CoincidentGeometry : CurveCurveApproachType.ParallelGeometry;
+        }
+        const pair = CurveLocationDetailPair.createCapture(CurveLocationDetail.createRayFractionPoint(rayA, fractionA, rayA.fractionToPoint(fractionA)), CurveLocationDetail.createRayFractionPoint(rayB, fractionB, rayB.fractionToPoint(fractionB)));
+        pair.approachType = pairType;
+        return pair;
+    }
+}
+//# sourceMappingURL=Ray3d.js.map

package/lib/esm/geometry3d/Ray3d.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"Ray3d.js","sourceRoot":"","sources":["../../../src/geometry3d/Ray3d.ts"],"names":[],"mappings":"AAAA;;;+FAG+F;AAE/F;;GAEG;AACH,OAAO,EAAE,sBAAsB,EAAE,mBAAmB,EAAE,uBAAuB,EAAE,MAAM,8BAA8B,CAAC;AACpH,OAAO,EAAE,SAAS,EAAmB,QAAQ,EAAE,MAAM,aAAa,CAAC;AACnE,OAAO,EAAE,WAAW,EAAE,MAAM,yBAAyB,CAAC;AACtD,OAAO,EAAE,QAAQ,EAAE,MAAM,YAAY,CAAC;AAEtC,OAAO,EAAE,QAAQ,EAAE,MAAM,mBAAmB,CAAC;AAC7C,OAAO,EAAE,OAAO,EAAE,QAAQ,EAAE,MAAM,mBAAmB,CAAC;AACtD,OAAO,EAAE,OAAO,EAAW,MAAM,SAAS,CAAC;AAC3C,OAAO,EAAE,SAAS,EAAE,MAAM,aAAa,CAAC;AAGxC;;;;;GAKG;AACH,MAAM,OAAO,KAAK;IAOhB,sCAAsC;IACtC,YAAoB,MAAe,EAAE,SAAmB;QACtD,IAAI,CAAC,MAAM,GAAG,MAAM,CAAC;QACrB,IAAI,CAAC,SAAS,GAAG,SAAS,CAAC;QAC3B,IAAI,CAAC,CAAC,GAAG,SAAS,CAAC;IACrB,CAAC;IACO,MAAM,CAAC,OAAO,CAAC,CAAS,EAAE,CAAS,EAAE,CAAS,EAAE,CAAS,EAAE,CAAS,EAAE,CAAS;QACrF,OAAO,IAAI,KAAK,CAAC,OAAO,CAAC,MAAM,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,EAAE,QAAQ,CAAC,MAAM,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;IACtE,CAAC;IACD,kCAAkC;IAC3B,MAAM,CAAC,WAAW,KAAY,OAAO,KAAK,CAAC,OAAO,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;IAC9E,kCAAkC;IAC3B,MAAM,CAAC,WAAW,KAAY,OAAO,KAAK,CAAC,OAAO,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;IAC9E,kCAAkC;IAC3B,MAAM,CAAC,WAAW,KAAY,OAAO,KAAK,CAAC,OAAO,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;IAC9E,mCAAmC;IAC5B,MAAM,CAAC,UAAU,CAAC,MAAc;QACrC,IAAI,MAAM,EAAE;YACV,MAAM,CAAC,MAAM,CAAC,OAAO,EAAE,CAAC;YACxB,MAAM,CAAC,SAAS,CAAC,OAAO,EAAE,CAAC;YAC3B,OAAO,MAAM,CAAC;SACf;QACD,OAAO,IAAI,KAAK,CAAC,OAAO,CAAC,UAAU,EAAE,EAAE,QAAQ,CAAC,UAAU,EAAE,CAAC,CAAC;IAChE,CAAC;IACD,kCAAkC;IAC3B,aAAa,CAAC,KAAY;QAC/B,OAAO,IAAI,CAAC,MAAM,CAAC,aAAa,CAAC,KAAK,CAAC,MAAM,CAAC,IAAI,IAAI,CAAC,SAAS,CAAC,aAAa,CAAC,KAAK,CAAC,SAAS,CAAC,CAAC;IAClG,CAAC;IACD,8CAA8C;IACvC,MAAM,CAAC,MAAM,CAAC,MAAe,EAAE,SAAmB,EAAE,MAAc;QACvE,IAAI,MAAM,EAAE;YACV,MAAM,CAAC,GAAG,CAAC,MAAM,EAAE,SAAS,CAAC,CAAC;YAC9B,OAAO,MAAM,CAAC;SACf;QACD,OAAO,IAAI,KAAK,CAAC,MAAM,CAAC,KAAK,EAAE,EAAE,SAAS,CAAC,KAAK,EAAE,CAAC,CAAC;IACtD,CAAC;IACD;;;;;OAKG;IACI,MAAM,CAAC,wBAAwB,CAAC,aAA2B,EAAE,kBAAgC,EAAE,MAAc;QAClH,MAAM,CAAC,GAAG,aAAa,CAAC,CAAC,CAAC,CAAC;QAC3B,MAAM,EAAE,GAAG,kBAAkB,CAAC,CAAC,CAAC,CAAC;QACjC,MAAM,CAAC,GAAG,aAAa,CAAC,CAAC,CAAC,CAAC;QAC3B,MAAM,CAAC,GAAG,aAAa,CAAC,CAAC,CAAC,CAAC;QAC3B,MAAM,CAAC,GAAG,aAAa,CAAC,CAAC,CAAC,CAAC;QAC3B,MAAM,EAAE,GAAG,kBAAkB,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,aAAa,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC;QAC7D,MAAM,EAAE,GAAG,kBAAkB,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,aAAa,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC;QAC7D,MAAM,EAAE,GAAG,kBAAkB,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,aAAa,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC;QAC7D,IAAI,QAAQ,CAAC,qBAAqB,CAAC,CAAC,CAAC;YACnC,OAAO,SAAS,CAAC;QACnB,MAAM,IAAI,GAAG,GAAG,GAAG,CAAC,CAAC;QACrB,MAAM,KAAK,GAAG,IAAI,GAAG,IAAI,CAAC;QAC1B,OAAO,KAAK,CAAC,YAAY,CAAC,CAAC,GAAG,IAAI,EAAE,CAAC,GAAG,IAAI,EAAE,CAAC,GAAG,IAAI,EAAE,EAAE,GAAG,KAAK,EAAE,EAAE,GAAG,KAAK,EAAE,EAAE,GAAG,KAAK,EAAE,MAAM,CAAC,CAAC;IACtG,CAAC;IACD,2DAA2D;IACpD,MAAM,CAAC,YAAY,CAAC,OAAe,EAAE,OAAe,EAAE,OAAe,EAAE,UAAkB,EAAE,UAAkB,EAAE,UAAkB,EAAE,MAAc;QACtJ,IAAI,MAAM,EAAE;YACV,MAAM,CAAC,YAAY,EAAE,CAAC,GAAG,CAAC,OAAO,EAAE,OAAO,EAAE,OAAO,CAAC,CAAC;YACrD,MAAM,CAAC,eAAe,EAAE,CAAC,GAAG,CAAC,UAAU,EAAE,UAAU,EAAE,UAAU,CAAC,CAAC;YACjE,OAAO,MAAM,CAAC;SACf;QACD,OAAO,IAAI,KAAK,CAAC,OAAO,CAAC,MAAM,CAAC,OAAO,EAAE,OAAO,EAAE,OAAO,CAAC,EAAE,QAAQ,CAAC,MAAM,CAAC,UAAU,EAAE,UAAU,EAAE,UAAU,CAAC,CAAC,CAAC;IACnH,CAAC;IACD,mDAAmD;IAC5C,MAAM,CAAC,aAAa,CAAC,MAAe,EAAE,SAAmB;QAC9D,OAAO,IAAI,KAAK,CAAC,MAAM,EAAE,SAAS,CAAC,CAAC;IACtC,CAAC;IACD,qEAAqE;IAC9D,MAAM,CAAC,uBAAuB,CAAC,MAAe,EAAE,SAAmB,EAAE,CAAS,EAAE,MAAc;QACnG,IAAI,MAAM,EAAE;YACV,MAAM,CAAC,MAAM,CAAC,OAAO,CAAC,MAAM,CAAC,CAAC;YAC9B,MAAM,CAAC,SAAS,CAAC,OAAO,CAAC,SAAS,CAAC,CAAC;YACpC,MAAM,CAAC,CAAC,GAAG,CAAC,CAAC;YACb,OAAO,MAAM,CAAC;SACf;QACD,MAAM,GAAG,IAAI,KAAK,CAAC,MAAM,CAAC,KAAK,EAAE,EAAE,SAAS,CAAC,KAAK,EAAE,CAAC,CAAC;QACtD,MAAM,CAAC,CAAC,GAAG,CAAC,CAAC;QACb,OAAO,MAAM,CAAC;IAChB,CAAC;IACD,uHAAuH;IAChH,MAAM,CAAC,cAAc,CAAC,MAAe,EAAE,MAAe,EAAE,MAAc;QAC3E,IAAI,MAAM,EAAE;YACV,MAAM,CAAC,MAAM,CAAC,OAAO,CAAC,MAAM,CAAC,CAAC;YAC9B,MAAM,CAAC,SAAS,CAAC,WAAW,CAAC,MAAM,EAAE,MAAM,CAAC,CAAC;YAC7C,OAAO,MAAM,CAAC;SACf;QACD,OAAO,IAAI,KAAK,CAAC,MAAM,EAAE,QAAQ,CAAC,cAAc,CAAC,MAAM,EAAE,MAAM,CAAC,CAAC,CAAC;IACpE,CAAC;IACD,8CAA8C;IACvC,YAAY,KAAc,OAAO,IAAI,CAAC,MAAM,CAAC,CAAC,CAAC;IACtD,wDAAwD;IACjD,eAAe,KAAe,OAAO,IAAI,CAAC,SAAS,CAAC,CAAC,CAAC;IAC7D,kDAAkD;IAC3C,GAAG,CAAC,MAAe,EAAE,SAAmB;QAC7C,IAAI,CAAC,MAAM,CAAC,OAAO,CAAC,MAAM,CAAC,CAAC;QAC5B,IAAI,CAAC,SAAS,CAAC,OAAO,CAAC,SAAS,CAAC,CAAC;IACpC,CAAC;IACD,qBAAqB;IACd,KAAK,CAAC,MAAc;QACzB,IAAI,MAAM,EAAE;YACV,MAAM,CAAC,GAAG,CAAC,IAAI,CAAC,MAAM,CAAC,KAAK,EAAE,EAAE,IAAI,CAAC,SAAS,CAAC,KAAK,EAAE,CAAC,CAAC;YACxD,OAAO,MAAM,CAAC;SACf;QACD,OAAO,IAAI,KAAK,CAAC,IAAI,CAAC,MAAM,CAAC,KAAK,EAAE,EAAE,IAAI,CAAC,SAAS,CAAC,KAAK,EAAE,CAAC,CAAC;IAChE,CAAC;IACD,4DAA4D;IACrD,gBAAgB,CAAC,SAAoB;QAC1C,OAAO,IAAI,KAAK,CAAC,SAAS,CAAC,eAAe,CAAC,IAAI,CAAC,MAAM,CAAC,EAAE,SAAS,CAAC,cAAc,CAAC,IAAI,CAAC,SAAS,CAAC,CAAC,CAAC;IACrG,CAAC;IAED,oEAAoE;IAC7D,uBAAuB,CAAC,SAAoB;QACjD,MAAM,MAAM,GAAG,SAAS,CAAC,sBAAsB,CAAC,IAAI,CAAC,MAAM,CAAC,CAAC;QAC7D,MAAM,SAAS,GAAG,SAAS,CAAC,MAAM,CAAC,4BAA4B,CAAC,IAAI,CAAC,SAAS,CAAC,CAAC,EAAE,IAAI,CAAC,SAAS,CAAC,CAAC,EAAE,IAAI,CAAC,SAAS,CAAC,CAAC,CAAC,CAAC;QACtH,IAAI,SAAS,KAAK,MAAM,IAAI,SAAS,KAAK,SAAS;YACjD,OAAO,IAAI,KAAK,CAAC,MAAM,EAAE,SAAS,CAAC,CAAC;QACtC,OAAO,SAAS,CAAC;IACnB,CAAC;IAED,kCAAkC;IAC3B,gBAAgB,CAAC,SAAoB;QAC1C,SAAS,CAAC,eAAe,CAAC,IAAI,CAAC,MAAM,EAAE,IAAI,CAAC,MAAM,CAAC,CAAC;QACpD,SAAS,CAAC,cAAc,CAAC,IAAI,CAAC,SAAS,EAAE,IAAI,CAAC,SAAS,CAAC,CAAC;IAC3D,CAAC;IACD,kCAAkC;IAC3B,OAAO,CAAC,MAAa,IAAU,IAAI,CAAC,GAAG,CAAC,MAAM,CAAC,MAAM,EAAE,MAAM,CAAC,SAAS,CAAC,CAAC,CAAC,CAAC;IAClF;;;OAGG;IACI,eAAe,CAAC,QAAgB,EAAE,MAAgB,IAAa,OAAO,IAAI,CAAC,MAAM,CAAC,UAAU,CAAC,IAAI,CAAC,SAAS,EAAE,QAAQ,EAAE,MAAM,CAAC,CAAC,CAAC,CAAC;IACxI,iHAAiH;IAC1G,iBAAiB,CAAC,UAAmB,IAAY,OAAO,IAAI,CAAC,SAAS,CAAC,kBAAkB,CAAC,IAAI,CAAC,MAAM,EAAE,UAAU,CAAC,CAAC,CAAC,CAAC;IAC5H;;OAEG;IACI,eAAe,CAAC,UAAmB;QACxC,OAAO,QAAQ,CAAC,kBAAkB,CAAC,IAAI,CAAC,SAAS,CAAC,kBAAkB,CAAC,IAAI,CAAC,MAAM,EAAE,UAAU,CAAC,EAAE,IAAI,CAAC,SAAS,CAAC,gBAAgB,EAAE,EAAE,CAAC,CAAC,CAAC;IACvI,CAAC;IACD;;;OAGG;IACI,iBAAiB,CAAC,UAAmB;QAC1C,OAAO,IAAI,CAAC,MAAM,CAAC,UAAU,CAAC,IAAI,CAAC,SAAS,EAAE,IAAI,CAAC,eAAe,CAAC,UAAU,CAAC,CAAC,CAAC;IAClF,CAAC;IACD;;OAEG;IACI,aAAa;QAClB,MAAM,IAAI,GAAG,QAAQ,CAAC,kBAAkB,CAAC,IAAI,CAAC,SAAS,EAAE,SAAS,CAAC,GAAG,CAAC,CAAC;QACxE,OAAO,SAAS,CAAC,qBAAqB,CAAC,IAAI,CAAC,MAAM,EAAE,IAAI,CAAC,CAAC;IAC5D,CAAC;IACD;;OAEG;IACI,WAAW,CAAC,IAAU;QAC3B,IAAI,CAAC,IAAI,EAAE;YACT,IAAI,CAAC,MAAM,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;YACzB,IAAI,CAAC,SAAS,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;YAC5B,OAAO;SACR;QACD,IAAI,CAAC,MAAM,CAAC,WAAW,CAAC,IAAI,CAAC,MAAM,CAAC,CAAC;QACrC,IAAI,CAAC,SAAS,CAAC,WAAW,CAAC,IAAI,CAAC,SAAS,CAAC,CAAC;IAC7C,CAAC;IACD;;;OAGG;IACI,+BAA+B,CAAC,YAAoB,GAAG;QAC5D,IAAI,IAAI,CAAC,SAAS,CAAC,mBAAmB,EAAE,EAAE;YACxC,IAAI,CAAC,SAAS,CAAC,YAAY,CAAC,SAAS,CAAC,CAAC;YACvC,OAAO,IAAI,CAAC;SACb;QACD,IAAI,CAAC,SAAS,CAAC,OAAO,EAAE,CAAC;QACzB,IAAI,CAAC,CAAC,GAAG,GAAG,CAAC;QACb,OAAO,KAAK,CAAC;IACf,CAAC;IACD;;;;;;;;OAQG;IACH,6CAA6C;IAC7C,0GAA0G;IAC1G,+EAA+E;IACxE,iCAAiC,CAAC,CAAS;QAChD,MAAM,SAAS,GAAG,QAAQ,CAAC,0BAA0B,CAAC;QACtD,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC;QACX,IAAI,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,SAAS,IAAI,IAAI,CAAC,SAAS,CAAC,mBAAmB,CAAC,SAAS,CAAC;YAC1E,OAAO,IAAI,CAAC;QACd,IAAI,CAAC,SAAS,CAAC,OAAO,EAAE,CAAC;QACzB,IAAI,CAAC,CAAC,GAAG,GAAG,CAAC;QACb,OAAO,KAAK,CAAC;IACf,CAAC;IACD;;;OAGG;IACI,MAAM,KAAU,OAAO,EAAE,MAAM,EAAE,IAAI,CAAC,MAAM,CAAC,MAAM,EAAE,EAAE,SAAS,EAAE,IAAI,CAAC,SAAS,CAAC,MAAM,EAAE,EAAE,CAAC,CAAC,CAAC;IACrG,gFAAgF;IACzE,MAAM,CAAC,QAAQ,CAAC,IAAU;QAC/B,MAAM,MAAM,GAAG,KAAK,CAAC,WAAW,EAAE,CAAC;QACnC,MAAM,CAAC,WAAW,CAAC,IAAI,CAAC,CAAC;QACzB,OAAO,MAAM,CAAC;IAChB,CAAC;IACD,qDAAqD;IAC9C,QAAQ,CAAC,UAAmB;QACjC,MAAM,EAAE,GAAG,IAAI,CAAC,SAAS,CAAC,gBAAgB,EAAE,CAAC;QAC7C,MAAM,EAAE,GAAG,IAAI,CAAC,iBAAiB,CAAC,UAAU,CAAC,CAAC;QAC9C,MAAM,EAAE,GAAG,QAAQ,CAAC,4BAA4B,CAAC,EAAE,CAAC,CAAC;QACrD,IAAI,EAAE;YACJ,OAAO,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,MAAM,CAAC,eAAe,CAAC,UAAU,CAAC,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC,CAAC;;YAEzE,OAAO,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,MAAM,CAAC,eAAe,CAAC,UAAU,CAAC,CAAC,CAAC;IAC9D,CAAC;IACD;;;;;OAKG;IACI,qBAAqB,CAAC,KAAmC,EAAE,MAAgB;QAChF,MAAM,OAAO,GAAG,QAAQ,CAAC,cAAc,CAAC,KAAK,CAAC,YAAY,EAAE,EAAE,IAAI,CAAC,MAAM,CAAC,CAAC;QAC3E,MAAM,KAAK,GAAG,IAAI,CAAC,SAAS,CAAC,UAAU,CAAC,KAAK,CAAC,YAAY,EAAE,CAAC,CAAC;QAC9D,MAAM,KAAK,GAAG,IAAI,CAAC,SAAS,CAAC,gBAAgB,EAAE,CAAC;QAChD,MAAM,KAAK,GAAG,OAAO,CAAC,UAAU,CAAC,KAAK,CAAC,YAAY,EAAE,CAAC,CAAC;QACvD,MAAM,QAAQ,GAAG,QAAQ,CAAC,yBAAyB,CAAC,CAAC,KAAK,EAAE,KAAK,CAAC,CAAC;QACnE,IAAI,SAAS,KAAK,QAAQ;YACxB,OAAO,SAAS,CAAC;QACnB,MAAM,SAAS,GAAG,QAAQ,CAAC,yBAAyB,CAAC,KAAK,EAAE,KAAK,CAAC,CAAC;QACnE,IAAI,SAAS,KAAK,SAAS;YACzB,OAAO,SAAS,CAAC;QACnB,IAAI,MAAM,EAAE;YACV,IAAI,CAAC,MAAM,CAAC,UAAU,CAAC,IAAI,CAAC,SAAS,EAAE,QAAQ,EAAE,MAAM,CAAC,CAAC;SAC1D;QACD,OAAO,QAAQ,CAAC;IAClB,CAAC;IAED;;;;OAIG;IACI,uBAAuB,CAAC,KAAc,EAAE,MAAgB;QAC7D,IAAI,KAAK,CAAC,MAAM;YACd,OAAO,OAAO,CAAC,UAAU,CAAC,MAAM,CAAC,CAAC;QACpC,MAAM,QAAQ,GAAG,OAAO,CAAC,QAAQ,CAAC,CAAC,QAAQ,CAAC,qBAAqB,EAAE,QAAQ,CAAC,qBAAqB,EAAE,MAAM,CAAC,CAAC;QAC3G,IAAI,QAAQ,CAAC,uBAAuB,CAAC,IAAI,CAAC,MAAM,CAAC,CAAC,EAAE,IAAI,CAAC,SAAS,CAAC,CAAC,EAAE,KAAK,CAAC,GAAG,CAAC,CAAC,EAAE,KAAK,CAAC,IAAI,CAAC,CAAC,CAAC;eAC3F,QAAQ,CAAC,uBAAuB,CAAC,IAAI,CAAC,MAAM,CAAC,CAAC,EAAE,IAAI,CAAC,SAAS,CAAC,CAAC,EAAE,KAAK,CAAC,GAAG,CAAC,CAAC,EAAE,KAAK,CAAC,IAAI,CAAC,CAAC,CAAC;eAC5F,QAAQ,CAAC,uBAAuB,CAAC,IAAI,CAAC,MAAM,CAAC,CAAC,EAAE,IAAI,CAAC,SAAS,CAAC,CAAC,EAAE,KAAK,CAAC,GAAG,CAAC,CAAC,EAAE,KAAK,CAAC,IAAI,CAAC,CAAC,CAAC;YAC/F,OAAO,QAAQ,CAAC;QAClB,OAAO,QAAQ,CAAC;IAClB,CAAC;IAED;;;OAGG;IACI,iCAAiC,CAAC,WAAmB,EAAE,MAAiB;QAC7E,MAAM,OAAO,GAAG,QAAQ,CAAC,cAAc,CAAC,IAAI,CAAC,MAAM,EAAE,WAAW,CAAC,CAAC;QAClE,MAAM,EAAE,GAAG,IAAI,CAAC,SAAS,CAAC,gBAAgB,EAAE,CAAC;QAC7C,MAAM,EAAE,GAAG,IAAI,CAAC,SAAS,CAAC,kBAAkB,CAAC,IAAI,CAAC,MAAM,EAAE,WAAW,CAAC,CAAC;QACvE,MAAM,QAAQ,GAAG,QAAQ,CAAC,kBAAkB,CAAC,EAAE,EAAE,EAAE,EAAE,GAAG,CAAC,CAAC;QAC1D,OAAO,OAAO,CAAC,UAAU,CAAC,IAAI,CAAC,SAAS,EAAE,CAAC,QAAQ,EAAE,MAAM,CAAC,CAAC;IAC/D,CAAC;IACD;;;;;;;;OAQG;IACI,MAAM,CAAC,yBAAyB,CAAC,IAAW,EAAE,IAAW;QAC9D,MAAM,qBAAqB,GAAG,QAAQ,CAAC,MAAM,EAAE,CAAC;QAChD,IAAI,SAAS,EAAE,SAAS,CAAC;QACzB,IAAI,MAAM,EAAE,MAAM,CAAC;QACnB,IAAI,QAAQ,CAAC;QACb,IAAI,WAAW,CAAC,mCAAmC,CACjD,IAAI,CAAC,MAAM,CAAC,CAAC,EAAE,IAAI,CAAC,MAAM,CAAC,CAAC,EAAE,IAAI,CAAC,MAAM,CAAC,CAAC,EAAE,IAAI,CAAC,SAAS,CAAC,CAAC,EAAE,IAAI,CAAC,SAAS,CAAC,CAAC,EAAE,IAAI,CAAC,SAAS,CAAC,CAAC,EACjG,IAAI,CAAC,MAAM,CAAC,CAAC,EAAE,IAAI,CAAC,MAAM,CAAC,CAAC,EAAE,IAAI,CAAC,MAAM,CAAC,CAAC,EAAE,IAAI,CAAC,SAAS,CAAC,CAAC,EAAE,IAAI,CAAC,SAAS,CAAC,CAAC,EAAE,IAAI,CAAC,SAAS,CAAC,CAAC,EAAE,qBAAqB,CAAC,EAAE;YAC3H,SAAS,GAAG,qBAAqB,CAAC,CAAC,CAAC;YACpC,SAAS,GAAG,qBAAqB,CAAC,CAAC,CAAC;YACpC,MAAM,GAAG,IAAI,CAAC,eAAe,CAAC,SAAS,CAAC,CAAC;YACzC,MAAM,GAAG,IAAI,CAAC,eAAe,CAAC,SAAS,CAAC,CAAC;YACzC,QAAQ,GAAG,MAAM,CAAC,mBAAmB,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,sBAAsB,CAAC,YAAY,CAAC,CAAC,CAAC,sBAAsB,CAAC,kBAAkB,CAAC;SACjI;aAAM;YACL,SAAS,GAAG,GAAG,CAAC;YAChB,SAAS,GAAG,IAAI,CAAC,eAAe,CAAC,IAAI,CAAC,MAAM,CAAC,CAAC;YAC9C,MAAM,GAAG,IAAI,CAAC,eAAe,CAAC,SAAS,CAAC,CAAC;YACzC,MAAM,GAAG,IAAI,CAAC,eAAe,CAAC,SAAS,CAAC,CAAC;YACzC,QAAQ,GAAG,MAAM,CAAC,mBAAmB,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,sBAAsB,CAAC,kBAAkB,CAAC,CAAC,CAAC,sBAAsB,CAAC,gBAAgB,CAAC;SACrI;QACD,MAAM,IAAI,GAAG,uBAAuB,CAAC,aAAa,CAChD,mBAAmB,CAAC,sBAAsB,CAAC,IAAI,EAAE,SAAS,EAAE,IAAI,CAAC,eAAe,CAAC,SAAS,CAAC,CAAC,EAC5F,mBAAmB,CAAC,sBAAsB,CAAC,IAAI,EAAE,SAAS,EAAE,IAAI,CAAC,eAAe,CAAC,SAAS,CAAC,CAAC,CAAC,CAAC;QAChG,IAAI,CAAC,YAAY,GAAG,QAAQ,CAAC;QAC7B,OAAO,IAAI,CAAC;IACd,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 CartesianGeometry\r\n */\r\nimport { CurveCurveApproachType, CurveLocationDetail, CurveLocationDetailPair } from \"../curve/CurveLocationDetail\";\r\nimport { AxisOrder, BeJSONFunctions, Geometry } from \"../Geometry\";\r\nimport { SmallSystem } from \"../numerics/Polynomials\";\r\nimport { Matrix3d } from \"./Matrix3d\";\r\nimport { Plane3dByOriginAndUnitNormal } from \"./Plane3dByOriginAndUnitNormal\";\r\nimport { Vector2d } from \"./Point2dVector2d\";\r\nimport { Point3d, Vector3d } from \"./Point3dVector3d\";\r\nimport { Range1d, Range3d } from \"./Range\";\r\nimport { Transform } from \"./Transform\";\r\nimport { XYAndZ } from \"./XYZProps\";\r\n\r\n/** A Ray3d contains\r\n * * an origin point.\r\n * * a direction vector.  The vector is NOT required to be normalized.\r\n *  * an optional weight (number).\r\n * @public\r\n */\r\nexport class Ray3d implements BeJSONFunctions {\r\n  /** The ray origin */\r\n  public origin: Point3d;\r\n  /** The ray direction.  This is commonly (but not always) a unit vector. */\r\n  public direction: Vector3d;\r\n  /** Numeric annotation. */\r\n  public a?: number; // optional, e.g. weight.\r\n  // constructor captures references !!!\r\n  private constructor(origin: Point3d, direction: Vector3d) {\r\n    this.origin = origin;\r\n    this.direction = direction;\r\n    this.a = undefined;\r\n  }\r\n  private static _create(x: number, y: number, z: number, u: number, v: number, w: number) {\r\n    return new Ray3d(Point3d.create(x, y, z), Vector3d.create(u, v, w));\r\n  }\r\n  /** Create a ray on the x axis. */\r\n  public static createXAxis(): Ray3d { return Ray3d._create(0, 0, 0, 1, 0, 0); }\r\n  /** Create a ray on the y axis. */\r\n  public static createYAxis(): Ray3d { return Ray3d._create(0, 0, 0, 0, 1, 0); }\r\n  /** Create a ray on the z axis. */\r\n  public static createZAxis(): Ray3d { return Ray3d._create(0, 0, 0, 0, 0, 1); }\r\n  /** Create a ray with all zeros. */\r\n  public static createZero(result?: Ray3d): Ray3d {\r\n    if (result) {\r\n      result.origin.setZero();\r\n      result.direction.setZero();\r\n      return result;\r\n    }\r\n    return new Ray3d(Point3d.createZero(), Vector3d.createZero());\r\n  }\r\n  /** Test for nearly equal rays. */\r\n  public isAlmostEqual(other: Ray3d): boolean {\r\n    return this.origin.isAlmostEqual(other.origin) && this.direction.isAlmostEqual(other.direction);\r\n  }\r\n  /** Create a ray from origin and direction. */\r\n  public static create(origin: Point3d, direction: Vector3d, result?: Ray3d): Ray3d {\r\n    if (result) {\r\n      result.set(origin, direction);\r\n      return result;\r\n    }\r\n    return new Ray3d(origin.clone(), direction.clone());\r\n  }\r\n  /**\r\n   * Given a homogeneous point and its derivative components, construct a Ray3d with cartesian coordinates and derivatives.\r\n   * @param weightedPoint `[x,y,z,w]` parts of weighted point.\r\n   * @param weightedDerivative `[x,y,z,w]` derivatives\r\n   * @param result\r\n   */\r\n  public static createWeightedDerivative(weightedPoint: Float64Array, weightedDerivative: Float64Array, result?: Ray3d): Ray3d | undefined {\r\n    const w = weightedPoint[3];\r\n    const dw = weightedDerivative[3];\r\n    const x = weightedPoint[0];\r\n    const y = weightedPoint[1];\r\n    const z = weightedPoint[2];\r\n    const dx = weightedDerivative[0] * w - weightedPoint[0] * dw;\r\n    const dy = weightedDerivative[1] * w - weightedPoint[1] * dw;\r\n    const dz = weightedDerivative[2] * w - weightedPoint[2] * dw;\r\n    if (Geometry.isSmallMetricDistance(w))\r\n      return undefined;\r\n    const divW = 1.0 / w;\r\n    const divWW = divW * divW;\r\n    return Ray3d.createXYZUVW(x * divW, y * divW, z * divW, dx * divWW, dy * divWW, dz * divWW, result);\r\n  }\r\n  /** Create from coordinates of the origin and direction. */\r\n  public static createXYZUVW(originX: number, originY: number, originZ: number, directionX: number, directionY: number, directionZ: number, result?: Ray3d): Ray3d {\r\n    if (result) {\r\n      result.getOriginRef().set(originX, originY, originZ);\r\n      result.getDirectionRef().set(directionX, directionY, directionZ);\r\n      return result;\r\n    }\r\n    return new Ray3d(Point3d.create(originX, originY, originZ), Vector3d.create(directionX, directionY, directionZ));\r\n  }\r\n  /** Capture origin and direction in a new Ray3d. */\r\n  public static createCapture(origin: Point3d, direction: Vector3d): Ray3d {\r\n    return new Ray3d(origin, direction);\r\n  }\r\n  /** Create from (clones of) origin, direction, and numeric weight. */\r\n  public static createPointVectorNumber(origin: Point3d, direction: Vector3d, a: number, result?: Ray3d): Ray3d {\r\n    if (result) {\r\n      result.origin.setFrom(origin);\r\n      result.direction.setFrom(direction);\r\n      result.a = a;\r\n      return result;\r\n    }\r\n    result = new Ray3d(origin.clone(), direction.clone());\r\n    result.a = a;\r\n    return result;\r\n  }\r\n  /** Create from origin and target.  The direction vector is the full length (non-unit) vector from origin to target. */\r\n  public static createStartEnd(origin: Point3d, target: Point3d, result?: Ray3d): Ray3d {\r\n    if (result) {\r\n      result.origin.setFrom(origin);\r\n      result.direction.setStartEnd(origin, target);\r\n      return result;\r\n    }\r\n    return new Ray3d(origin, Vector3d.createStartEnd(origin, target));\r\n  }\r\n  /** Return a reference to the ray's origin. */\r\n  public getOriginRef(): Point3d { return this.origin; }\r\n  /** Return a reference to the ray's direction vector. */\r\n  public getDirectionRef(): Vector3d { return this.direction; }\r\n  /** copy coordinates from origin and direction. */\r\n  public set(origin: Point3d, direction: Vector3d): void {\r\n    this.origin.setFrom(origin);\r\n    this.direction.setFrom(direction);\r\n  }\r\n  /** Clone the ray. */\r\n  public clone(result?: Ray3d): Ray3d {\r\n    if (result) {\r\n      result.set(this.origin.clone(), this.direction.clone());\r\n      return result;\r\n    }\r\n    return new Ray3d(this.origin.clone(), this.direction.clone());\r\n  }\r\n  /** Create a clone and return the transform of the clone. */\r\n  public cloneTransformed(transform: Transform): Ray3d {\r\n    return new Ray3d(transform.multiplyPoint3d(this.origin), transform.multiplyVector(this.direction));\r\n  }\r\n\r\n  /** Create a clone and return the inverse transform of the clone. */\r\n  public cloneInverseTransformed(transform: Transform): Ray3d | undefined {\r\n    const origin = transform.multiplyInversePoint3d(this.origin);\r\n    const direction = transform.matrix.multiplyInverseXYZAsVector3d(this.direction.x, this.direction.y, this.direction.z);\r\n    if (undefined !== origin && undefined !== direction)\r\n      return new Ray3d(origin, direction);\r\n    return undefined;\r\n  }\r\n\r\n  /** Apply a transform in place. */\r\n  public transformInPlace(transform: Transform) {\r\n    transform.multiplyPoint3d(this.origin, this.origin);\r\n    transform.multiplyVector(this.direction, this.direction);\r\n  }\r\n  /** Copy data from another ray. */\r\n  public setFrom(source: Ray3d): void { this.set(source.origin, source.direction); }\r\n  /** * fraction 0 is the ray origin.\r\n   * * fraction 1 is at the end of the direction vector when placed at the origin.\r\n   * @returns Return a point at fractional position along the ray.\r\n   */\r\n  public fractionToPoint(fraction: number, result?: Point3d): Point3d { return this.origin.plusScaled(this.direction, fraction, result); }\r\n  /** Return the dot product of the ray's direction vector with a vector from the ray origin to the space point. */\r\n  public dotProductToPoint(spacePoint: Point3d): number { return this.direction.dotProductStartEnd(this.origin, spacePoint); }\r\n  /**\r\n   * Return the fractional coordinate (along the direction vector) of the spacePoint projected to the ray.\r\n   */\r\n  public pointToFraction(spacePoint: Point3d): number {\r\n    return Geometry.safeDivideFraction(this.direction.dotProductStartEnd(this.origin, spacePoint), this.direction.magnitudeSquared(), 0);\r\n  }\r\n  /**\r\n   *\r\n   * Return the spacePoint projected onto the ray.\r\n   */\r\n  public projectPointToRay(spacePoint: Point3d): Point3d {\r\n    return this.origin.plusScaled(this.direction, this.pointToFraction(spacePoint));\r\n  }\r\n  /** Return a transform for rigid axes\r\n   * at ray origin with z in ray direction.  If the direction vector is zero, axes default to identity (from createHeadsUpTriad)\r\n   */\r\n  public toRigidZFrame(): Transform | undefined {\r\n    const axes = Matrix3d.createRigidHeadsUp(this.direction, AxisOrder.ZXY);\r\n    return Transform.createOriginAndMatrix(this.origin, axes);\r\n  }\r\n  /**\r\n   * Convert {origin:[x,y,z], direction:[u,v,w]} to a Ray3d.\r\n   */\r\n  public setFromJSON(json?: any) {\r\n    if (!json) {\r\n      this.origin.set(0, 0, 0);\r\n      this.direction.set(0, 0, 1);\r\n      return;\r\n    }\r\n    this.origin.setFromJSON(json.origin);\r\n    this.direction.setFromJSON(json.direction);\r\n  }\r\n  /**\r\n   * try to scale the direction vector to a given magnitude.\r\n   * @returns Returns false if ray direction is a zero vector.\r\n   */\r\n  public trySetDirectionMagnitudeInPlace(magnitude: number = 1.0): boolean {\r\n    if (this.direction.tryNormalizeInPlace()) {\r\n      this.direction.scaleInPlace(magnitude);\r\n      return true;\r\n    }\r\n    this.direction.setZero();\r\n    this.a = 0.0;\r\n    return false;\r\n  }\r\n  /**\r\n   * * If parameter `a` is clearly nonzero and the direction vector can be normalized,\r\n   *    * save the parameter `a` as the optional `a` member of the ray.\r\n   *    * normalize the ray's direction vector\r\n   * * If parameter `a` is nearly zero,\r\n   *   * Set the `a` member to zero\r\n   *   * Set the ray's direction vector to zero.\r\n   * @param a area to be saved.\r\n   */\r\n  // input a ray and \"a\" understood as an area.\r\n  // if a is clearly nonzero metric squared and the vector can be normalized, install those and return true.\r\n  // otherwise set ray.z to zero and zero the vector of the ray and return false.\r\n  public tryNormalizeInPlaceWithAreaWeight(a: number): boolean {\r\n    const tolerance = Geometry.smallMetricDistanceSquared;\r\n    this.a = a;\r\n    if (Math.abs(a) > tolerance && this.direction.tryNormalizeInPlace(tolerance))\r\n      return true;\r\n    this.direction.setZero();\r\n    this.a = 0.0;\r\n    return false;\r\n  }\r\n  /**\r\n   * Convert an Angle to a JSON object.\r\n   * @return {*} [origin,normal]\r\n   */\r\n  public toJSON(): any { return { origin: this.origin.toJSON(), direction: this.direction.toJSON() }; }\r\n  /** Create a new ray from json object.  See `setFromJSON` for json structure; */\r\n  public static fromJSON(json?: any) {\r\n    const result = Ray3d.createXAxis();\r\n    result.setFromJSON(json);\r\n    return result;\r\n  }\r\n  /** return distance from the ray to point in space */\r\n  public distance(spacePoint: Point3d): number {\r\n    const uu = this.direction.magnitudeSquared();\r\n    const uv = this.dotProductToPoint(spacePoint);\r\n    const aa = Geometry.inverseMetricDistanceSquared(uu);\r\n    if (aa)\r\n      return Math.sqrt(this.origin.distanceSquared(spacePoint) - uv * uv * aa);\r\n    else\r\n      return Math.sqrt(this.origin.distanceSquared(spacePoint));\r\n  }\r\n  /**\r\n   * Return the intersection of the unbounded ray with a plane.\r\n   * Stores the point of intersection in the result point given as a parameter,\r\n   * and returns the parameter along the ray where the intersection occurs.\r\n   * Returns undefined if the ray and plane are parallel or coplanar.\r\n   */\r\n  public intersectionWithPlane(plane: Plane3dByOriginAndUnitNormal, result?: Point3d): number | undefined {\r\n    const vectorA = Vector3d.createStartEnd(plane.getOriginRef(), this.origin);\r\n    const uDotN = this.direction.dotProduct(plane.getNormalRef());\r\n    const nDotN = this.direction.magnitudeSquared();\r\n    const aDotN = vectorA.dotProduct(plane.getNormalRef());\r\n    const division = Geometry.conditionalDivideFraction(-aDotN, uDotN);\r\n    if (undefined === division)\r\n      return undefined;\r\n    const division1 = Geometry.conditionalDivideFraction(nDotN, uDotN);\r\n    if (undefined === division1)\r\n      return undefined;\r\n    if (result) {\r\n      this.origin.plusScaled(this.direction, division, result);\r\n    }\r\n    return division;\r\n  }\r\n\r\n  /**\r\n   * * Find intersection of the ray with a Range3d.\r\n   * * return the range of fractions (on the ray) which are \"inside\" the range.\r\n   * * Note that a range is always returned;  if there is no intersection it is indicated by the test `result.sNull`\r\n   */\r\n  public intersectionWithRange3d(range: Range3d, result?: Range1d): Range1d {\r\n    if (range.isNull)\r\n      return Range1d.createNull(result);\r\n    const interval = Range1d.createXX(-Geometry.largeCoordinateResult, Geometry.largeCoordinateResult, result);\r\n    if (interval.clipLinearMapToInterval(this.origin.x, this.direction.x, range.low.x, range.high.x)\r\n      && interval.clipLinearMapToInterval(this.origin.y, this.direction.y, range.low.y, range.high.y)\r\n      && interval.clipLinearMapToInterval(this.origin.z, this.direction.z, range.low.z, range.high.z))\r\n      return interval;\r\n    return interval;\r\n  }\r\n\r\n  /** Construct a vector from `ray.origin` to target point.\r\n   * * return the part of the vector that is perpendicular to `ray.direction`.\r\n   *  * i.e. return the shortest vector from the ray to the point.\r\n   */\r\n  public perpendicularPartOfVectorToTarget(targetPoint: XYAndZ, result?: Vector3d): Vector3d {\r\n    const vectorV = Vector3d.createStartEnd(this.origin, targetPoint);\r\n    const uu = this.direction.magnitudeSquared();\r\n    const uv = this.direction.dotProductStartEnd(this.origin, targetPoint);\r\n    const fraction = Geometry.safeDivideFraction(uv, uu, 0.0);\r\n    return vectorV.plusScaled(this.direction, -fraction, result);\r\n  }\r\n  /** Determine if two rays intersect, are fully overlapped, parallel but no coincident, or skew\r\n   * * Return a CurveLocationDetailPair which\r\n   * * contains fraction and point on each ray.\r\n   * * has (in the CurveLocationDetailPair structure, as member approachType) annotation indicating one of these relationships\r\n   *   * CurveCurveApproachType.Intersection -- the rays have a simple intersection, at fractions indicated in detailA and detailB\r\n   *   * CurveCurveApproachType.PerpendicularChord -- there is pair of where the rays have closest approach.  The rays are skew in space.\r\n   *   * CurveCurveApproachType.CoincidentGeometry -- the rays are the same unbounded line in space. The fractions and points are a representative single common point.\r\n   *   * CurveCurveApproachType.Parallel -- the rays are parallel (and not coincident).   The two points are at the minimum distance\r\n   */\r\n  public static closestApproachRay3dRay3d(rayA: Ray3d, rayB: Ray3d): CurveLocationDetailPair {\r\n    const intersectionFractions = Vector2d.create();\r\n    let fractionA, fractionB;\r\n    let pointA, pointB;\r\n    let pairType;\r\n    if (SmallSystem.ray3dXYZUVWClosestApproachUnbounded(\r\n      rayA.origin.x, rayA.origin.y, rayA.origin.z, rayA.direction.x, rayA.direction.y, rayA.direction.z,\r\n      rayB.origin.x, rayB.origin.y, rayB.origin.z, rayB.direction.x, rayB.direction.y, rayB.direction.z, intersectionFractions)) {\r\n      fractionA = intersectionFractions.x;\r\n      fractionB = intersectionFractions.y;\r\n      pointA = rayA.fractionToPoint(fractionA);\r\n      pointB = rayB.fractionToPoint(fractionB);\r\n      pairType = pointA.isAlmostEqualMetric(pointB) ? CurveCurveApproachType.Intersection : CurveCurveApproachType.PerpendicularChord;\r\n    } else {\r\n      fractionB = 0.0;\r\n      fractionA = rayA.pointToFraction(rayB.origin);\r\n      pointA = rayA.fractionToPoint(fractionA);\r\n      pointB = rayB.fractionToPoint(fractionB);\r\n      pairType = pointA.isAlmostEqualMetric(pointB) ? CurveCurveApproachType.CoincidentGeometry : CurveCurveApproachType.ParallelGeometry;\r\n    }\r\n    const pair = CurveLocationDetailPair.createCapture(\r\n      CurveLocationDetail.createRayFractionPoint(rayA, fractionA, rayA.fractionToPoint(fractionA)),\r\n      CurveLocationDetail.createRayFractionPoint(rayB, fractionB, rayB.fractionToPoint(fractionB)));\r\n    pair.approachType = pairType;\r\n    return pair;\r\n  }\r\n}\r\n"]}

package/lib/esm/geometry3d/ReusableObjectCache.d.ts

@@ -0,0 +1,55 @@
+import { GrowableXYZArray } from "./GrowableXYZArray";
+/** @packageDocumentation
+ * @module ArraysAndInterfaces
+ */
+/**
+ * abstract class managing an array of objects of type T, available for reuse by trusted callers.
+ * * Derived class must implement these methods:
+ *   * `createForCache()` -- create a new, ready-to-use object
+ *   * `clearForCache (data: T)` -- tidy up `data` so it can be reused.
+ * @internal
+ */
+export declare abstract class ReusableObjectCache<T> {
+    protected abstract clearForCache(data: T): void;
+    protected abstract createForCache(): T;
+    private _cachedObjects;
+    numDrop: number;
+    numCreate: number;
+    numReuse: number;
+    /**
+     * create a new cache for objects of type T
+     */
+    protected constructor();
+    /** Present `data` for storage in the cache, and hence reuse by any subsequent `grabFromCache`
+     *   * `data` will be sent to `clearForCache`.
+     *   * caller should never refer to this instance again.
+     */
+    dropToCache(data: T | undefined): void;
+    /**
+     * grab an object from the cache.
+     *  * The returned object becomes property of the caller.
+     *  * That is, the cache does not remember it for any further management
+     * @param data
+     */
+    grabFromCache(): T;
+    /** Drop all entries of data[] to the cache.
+     * @param data on input, the data to drop. on output, data is an empty array.
+     */
+    dropAllToCache(data: T[]): void;
+}
+/**
+ * Cache of GrowableXYZArray.
+ * Intended for use by (for instance) clipping methods that can be structured to have disciplined reuse of a small number of arrays for a large number of steps.
+ * @internal
+ */
+export declare class GrowableXYZArrayCache extends ReusableObjectCache<GrowableXYZArray> {
+    protected clearForCache(data: GrowableXYZArray): void;
+    protected createForCache(): GrowableXYZArray;
+    constructor();
+    /**
+     * Grab an array from the cache and immediately fill from a source
+     * @param source
+     */
+    grabAndFill(source: GrowableXYZArray): GrowableXYZArray;
+}
+//# sourceMappingURL=ReusableObjectCache.d.ts.map

package/lib/esm/geometry3d/ReusableObjectCache.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"ReusableObjectCache.d.ts","sourceRoot":"","sources":["../../../src/geometry3d/ReusableObjectCache.ts"],"names":[],"mappings":"AAIA,OAAO,EAAE,gBAAgB,EAAE,MAAM,oBAAoB,CAAC;AAEtD;;GAEG;AACH;;;;;;GAMG;AACH,8BAAsB,mBAAmB,CAAC,CAAC;IACzC,SAAS,CAAC,QAAQ,CAAC,aAAa,CAAC,IAAI,EAAE,CAAC,GAAG,IAAI;IAC/C,SAAS,CAAC,QAAQ,CAAC,cAAc,IAAI,CAAC;IACtC,OAAO,CAAC,cAAc,CAAM;IACrB,OAAO,EAAE,MAAM,CAAC;IAChB,SAAS,EAAE,MAAM,CAAC;IAClB,QAAQ,EAAE,MAAM,CAAC;IACxB;;OAEG;IACH,SAAS;IAMT;;;OAGG;IACI,WAAW,CAAC,IAAI,EAAE,CAAC,GAAG,SAAS;IAOtC;;;;;OAKG;IACI,aAAa,IAAI,CAAC;IAUzB;;OAEG;IACI,cAAc,CAAC,IAAI,EAAE,CAAC,EAAE;CAKhC;AACD;;;;GAIG;AACH,qBAAa,qBAAsB,SAAQ,mBAAmB,CAAC,gBAAgB,CAAC;IAC9E,SAAS,CAAC,aAAa,CAAC,IAAI,EAAE,gBAAgB,GAAG,IAAI;IACrD,SAAS,CAAC,cAAc,IAAI,gBAAgB;;IAE5C;;;OAGG;IACI,WAAW,CAAC,MAAM,EAAE,gBAAgB,GAAG,gBAAgB;CAM/D"}

package/lib/esm/geometry3d/ReusableObjectCache.js

@@ -0,0 +1,82 @@
+/*---------------------------------------------------------------------------------------------
+* Copyright (c) Bentley Systems, Incorporated. All rights reserved.
+* See LICENSE.md in the project root for license terms and full copyright notice.
+*--------------------------------------------------------------------------------------------*/
+import { GrowableXYZArray } from "./GrowableXYZArray";
+/** @packageDocumentation
+ * @module ArraysAndInterfaces
+ */
+/**
+ * abstract class managing an array of objects of type T, available for reuse by trusted callers.
+ * * Derived class must implement these methods:
+ *   * `createForCache()` -- create a new, ready-to-use object
+ *   * `clearForCache (data: T)` -- tidy up `data` so it can be reused.
+ * @internal
+ */
+export class ReusableObjectCache {
+    /**
+     * create a new cache for objects of type T
+     */
+    constructor() {
+        this._cachedObjects = [];
+        this.numDrop = 0;
+        this.numCreate = 0;
+        this.numReuse = 0;
+    }
+    /** Present `data` for storage in the cache, and hence reuse by any subsequent `grabFromCache`
+     *   * `data` will be sent to `clearForCache`.
+     *   * caller should never refer to this instance again.
+     */
+    dropToCache(data) {
+        if (data) {
+            this.numDrop++;
+            this.clearForCache(data);
+            this._cachedObjects.push(data);
+        }
+    }
+    /**
+     * grab an object from the cache.
+     *  * The returned object becomes property of the caller.
+     *  * That is, the cache does not remember it for any further management
+     * @param data
+     */
+    grabFromCache() {
+        let data = this._cachedObjects.pop();
+        if (data === undefined) {
+            data = this.createForCache();
+            this.numCreate++;
+        }
+        else {
+            this.numReuse++;
+        }
+        return data;
+    }
+    /** Drop all entries of data[] to the cache.
+     * @param data on input, the data to drop. on output, data is an empty array.
+     */
+    dropAllToCache(data) {
+        while (data.length > 0) {
+            this.dropToCache(data.pop());
+        }
+    }
+}
+/**
+ * Cache of GrowableXYZArray.
+ * Intended for use by (for instance) clipping methods that can be structured to have disciplined reuse of a small number of arrays for a large number of steps.
+ * @internal
+ */
+export class GrowableXYZArrayCache extends ReusableObjectCache {
+    clearForCache(data) { data.length = 0; }
+    createForCache() { return new GrowableXYZArray(10); }
+    constructor() { super(); }
+    /**
+     * Grab an array from the cache and immediately fill from a source
+     * @param source
+     */
+    grabAndFill(source) {
+        const dest = this.grabFromCache();
+        dest.pushFromGrowableXYZArray(source);
+        return dest;
+    }
+}
+//# sourceMappingURL=ReusableObjectCache.js.map

package/lib/esm/geometry3d/ReusableObjectCache.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"ReusableObjectCache.js","sourceRoot":"","sources":["../../../src/geometry3d/ReusableObjectCache.ts"],"names":[],"mappings":"AAAA;;;+FAG+F;AAC/F,OAAO,EAAE,gBAAgB,EAAE,MAAM,oBAAoB,CAAC;AAEtD;;GAEG;AACH;;;;;;GAMG;AACH,MAAM,OAAgB,mBAAmB;IAOvC;;OAEG;IACH;QACE,IAAI,CAAC,cAAc,GAAG,EAAE,CAAC;QACzB,IAAI,CAAC,OAAO,GAAG,CAAC,CAAC;QACjB,IAAI,CAAC,SAAS,GAAG,CAAC,CAAC;QACnB,IAAI,CAAC,QAAQ,GAAG,CAAC,CAAC;IACpB,CAAC;IACD;;;OAGG;IACI,WAAW,CAAC,IAAmB;QACpC,IAAI,IAAI,EAAE;YACR,IAAI,CAAC,OAAO,EAAE,CAAC;YACf,IAAI,CAAC,aAAa,CAAC,IAAI,CAAC,CAAC;YACzB,IAAI,CAAC,cAAc,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC;SAChC;IACH,CAAC;IACD;;;;;OAKG;IACI,aAAa;QAClB,IAAI,IAAI,GAAG,IAAI,CAAC,cAAc,CAAC,GAAG,EAAE,CAAC;QACrC,IAAI,IAAI,KAAK,SAAS,EAAE;YACtB,IAAI,GAAG,IAAI,CAAC,cAAc,EAAE,CAAC;YAC7B,IAAI,CAAC,SAAS,EAAE,CAAC;SAClB;aAAM;YACL,IAAI,CAAC,QAAQ,EAAE,CAAC;SACjB;QACD,OAAO,IAAI,CAAC;IACd,CAAC;IACD;;OAEG;IACI,cAAc,CAAC,IAAS;QAC7B,OAAO,IAAI,CAAC,MAAM,GAAG,CAAC,EAAE;YACtB,IAAI,CAAC,WAAW,CAAC,IAAI,CAAC,GAAG,EAAE,CAAC,CAAC;SAC9B;IACH,CAAC;CACF;AACD;;;;GAIG;AACH,MAAM,OAAO,qBAAsB,SAAQ,mBAAqC;IACpE,aAAa,CAAC,IAAsB,IAAU,IAAI,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC,CAAC;IAChE,cAAc,KAAuB,OAAO,IAAI,gBAAgB,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;IACjF,gBAAuB,KAAK,EAAE,CAAC,CAAC,CAAC;IACjC;;;OAGG;IACI,WAAW,CAAC,MAAwB;QACzC,MAAM,IAAI,GAAG,IAAI,CAAC,aAAa,EAAE,CAAC;QAClC,IAAI,CAAC,wBAAwB,CAAC,MAAM,CAAC,CAAC;QACtC,OAAO,IAAI,CAAC;IAEd,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\nimport { GrowableXYZArray } from \"./GrowableXYZArray\";\r\n\r\n/** @packageDocumentation\r\n * @module ArraysAndInterfaces\r\n */\r\n/**\r\n * abstract class managing an array of objects of type T, available for reuse by trusted callers.\r\n * * Derived class must implement these methods:\r\n *   * `createForCache()` -- create a new, ready-to-use object\r\n *   * `clearForCache (data: T)` -- tidy up `data` so it can be reused.\r\n * @internal\r\n */\r\nexport abstract class ReusableObjectCache<T> {\r\n  protected abstract clearForCache(data: T): void;\r\n  protected abstract createForCache(): T;\r\n  private _cachedObjects: T[];\r\n  public numDrop: number;\r\n  public numCreate: number;\r\n  public numReuse: number;\r\n  /**\r\n   * create a new cache for objects of type T\r\n   */\r\n  protected constructor() {\r\n    this._cachedObjects = [];\r\n    this.numDrop = 0;\r\n    this.numCreate = 0;\r\n    this.numReuse = 0;\r\n  }\r\n  /** Present `data` for storage in the cache, and hence reuse by any subsequent `grabFromCache`\r\n   *   * `data` will be sent to `clearForCache`.\r\n   *   * caller should never refer to this instance again.\r\n   */\r\n  public dropToCache(data: T | undefined) {\r\n    if (data) {\r\n      this.numDrop++;\r\n      this.clearForCache(data);\r\n      this._cachedObjects.push(data);\r\n    }\r\n  }\r\n  /**\r\n   * grab an object from the cache.\r\n   *  * The returned object becomes property of the caller.\r\n   *  * That is, the cache does not remember it for any further management\r\n   * @param data\r\n   */\r\n  public grabFromCache(): T {\r\n    let data = this._cachedObjects.pop();\r\n    if (data === undefined) {\r\n      data = this.createForCache();\r\n      this.numCreate++;\r\n    } else {\r\n      this.numReuse++;\r\n    }\r\n    return data;\r\n  }\r\n  /** Drop all entries of data[] to the cache.\r\n   * @param data on input, the data to drop. on output, data is an empty array.\r\n   */\r\n  public dropAllToCache(data: T[]) {\r\n    while (data.length > 0) {\r\n      this.dropToCache(data.pop());\r\n    }\r\n  }\r\n}\r\n/**\r\n * Cache of GrowableXYZArray.\r\n * Intended for use by (for instance) clipping methods that can be structured to have disciplined reuse of a small number of arrays for a large number of steps.\r\n * @internal\r\n */\r\nexport class GrowableXYZArrayCache extends ReusableObjectCache<GrowableXYZArray> {\r\n  protected clearForCache(data: GrowableXYZArray): void { data.length = 0; }\r\n  protected createForCache(): GrowableXYZArray { return new GrowableXYZArray(10); }\r\n  public constructor() { super(); }\r\n  /**\r\n   * Grab an array from the cache and immediately fill from a source\r\n   * @param source\r\n   */\r\n  public grabAndFill(source: GrowableXYZArray): GrowableXYZArray {\r\n    const dest = this.grabFromCache();\r\n    dest.pushFromGrowableXYZArray(source);\r\n    return dest;\r\n\r\n  }\r\n}\r\n"]}

package/lib/esm/geometry3d/Segment1d.d.ts

@@ -0,0 +1,106 @@
+/**
+ * * A Segment1d is an interval of an axis named x.
+ * * The interval is defined by two values x0 and x1.
+ * * The x0 and x1 values can be in either order.
+ *   * if `x0 < x1` fractional coordinates within the segment move from left to right.
+ *   * if `x0 > x1` fractional coordinates within the segment move from right to left.
+ * * This differs from a Range1d in that:
+ * * For a Range1d the reversed ordering of its limit values means "empty interval".
+ * * For a Segment1d the reversed ordering is a real interval but fractional positions move backwards.
+ * * The segment is parameterized with a fraction
+ * * * Fraction 0 is the start (`x0`)
+ * * * Fraction 1 is the end (`x1`)
+ * * * The fraction equation is `x = x0 + fraction * (x1-x0)` or (equivalently) `x = (1-fraction) * x0 + fraction * x1`
+ * @public
+ */
+export declare class Segment1d {
+    /** start coordinate */
+    x0: number;
+    /** end coordinate */
+    x1: number;
+    private constructor();
+    /**
+     * replace both end values.
+     * @param x0 new x0 value
+     * @param x1 new y0 value
+     */
+    set(x0: number, x1: number): void;
+    /**
+     * shift (translate) the segment along its axis by adding `dx` to both `x0` and `x1`.
+     * @param dx value to add to both x0 and x1
+     */
+    shift(dx: number): void;
+    /**
+     * create segment1d with given end values
+     * @param x0 start value
+     * @param x1 end value
+     * @param result optional pre-existing result to be reinitialized.
+     */
+    static create(x0?: number, x1?: number, result?: Segment1d): Segment1d;
+    /**
+     * Copy both end values from other Segment1d
+     * @param other source Segment1d
+     */
+    setFrom(other: Segment1d): void;
+    /**
+     * clone this Segment1d, return as a separate object.
+     */
+    clone(): Segment1d;
+    /**
+     * Returns true if both coordinates (`x0` and `x1`) are in the 0..1 range.
+     */
+    get isIn01(): boolean;
+    /**
+     * Evaluate the segment at fractional position
+     * @returns position within the segment
+     * @param fraction fractional position within this segment
+     */
+    fractionToPoint(fraction: number): number;
+    /**
+     * Return the signed start-to-end shift (aka signed distance)
+     */
+    signedDelta(): number;
+    /**
+     * Return the absolute start-to-end shift (aka distance)
+     */
+    absoluteDelta(): number;
+    /**
+     * * swap the x0 and x1 member values.
+     * * This makes the fractionToPoint evaluates reverse direction.
+     */
+    reverseInPlace(): void;
+    /**
+     * * if `x1<x0` multiplied by the scale factor is (strictly) negative, swap the x0 and x1 member values.
+     * * This makes the fractionToPoint evaluates reverse direction.
+     */
+    reverseIfNeededForDeltaSign(sign?: number): void;
+    /**
+     * Near equality test, using Geometry.isSameCoordinate for tolerances.
+     */
+    isAlmostEqual(other: Segment1d): boolean;
+    /**
+     * Return true if the segment limits are (exactly) 0 and 1
+     */
+    get isExact01(): boolean;
+    /**
+     * Return true if the segment limits are (exactly) 1 and 0
+     */
+    get isExact01Reversed(): boolean;
+    /** On input, `this` is an interval of a line.  On output, the interval has been clipped to positive parts of a linear function
+     * * f0 and f1 are values at parameter values 0 and 1 (which are in general NOT x0 and x1)
+     * * From that determine where the segment crosses function value 0.
+     * * The segment contains some interval in the same parameter space.
+     * * Clip the segment to the positive part of the space.
+     * * Return true (and modify the segment) if any of the segment remains.
+     * * Return false (but without modifying the segment) if the active part is entirely out.
+     */
+    clipBy01FunctionValuesPositive(f0: number, f1: number): boolean;
+    /**
+     * * On input, (f0,f1) is a (directed) segment.
+     * * On output, it is restricted to (0,1) while maintaining direction
+     * * If the clip leaves nothing, leave this segment alone and return false.
+     * * If the clip leaves something, update this segment and return true.
+     */
+    clampDirectedTo01(): boolean;
+}
+//# sourceMappingURL=Segment1d.d.ts.map

package/lib/esm/geometry3d/Segment1d.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"Segment1d.d.ts","sourceRoot":"","sources":["../../../src/geometry3d/Segment1d.ts"],"names":[],"mappings":"AAWA;;;;;;;;;;;;;;GAcG;AACH,qBAAa,SAAS;IACpB,uBAAuB;IAChB,EAAE,EAAE,MAAM,CAAC;IAClB,qBAAqB;IACd,EAAE,EAAE,MAAM,CAAC;IAClB,OAAO;IAIP;;;;OAIG;IACI,GAAG,CAAC,EAAE,EAAE,MAAM,EAAE,EAAE,EAAE,MAAM;IACjC;;;OAGG;IACI,KAAK,CAAC,EAAE,EAAE,MAAM;IACvB;;;;;OAKG;WACW,MAAM,CAAC,EAAE,GAAE,MAAU,EAAE,EAAE,GAAE,MAAU,EAAE,MAAM,CAAC,EAAE,SAAS,GAAG,SAAS;IAMnF;;;OAGG;IACI,OAAO,CAAC,KAAK,EAAE,SAAS;IAC/B;;OAEG;IACI,KAAK,IAAI,SAAS;IACzB;;OAEG;IACH,IAAW,MAAM,YAGhB;IACD;;;;OAIG;IACI,eAAe,CAAC,QAAQ,EAAE,MAAM,GAAG,MAAM;IAChD;;OAEG;IACI,WAAW,IAAI,MAAM;IAC5B;;OAEG;IACI,aAAa,IAAI,MAAM;IAC9B;;;OAGG;IACI,cAAc,IAAI,IAAI;IAC7B;;;OAGG;IACI,2BAA2B,CAAC,IAAI,GAAE,MAAU,GAAG,IAAI;IAI1D;;OAEG;IACI,aAAa,CAAC,KAAK,EAAE,SAAS,GAAG,OAAO;IAG/C;;OAEG;IACH,IAAW,SAAS,IAAI,OAAO,CAA+C;IAC9E;;OAEG;IACH,IAAW,iBAAiB,IAAI,OAAO,CAA+C;IAEtF;;;;;;;OAOG;IACI,8BAA8B,CAAC,EAAE,EAAE,MAAM,EAAE,EAAE,EAAE,MAAM,GAAG,OAAO;IAqBtE;;;;;OAKG;IACI,iBAAiB,IAAI,OAAO;CAkBpC"}