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

package/lib/esm/solid/TorusPipe.js

@@ -0,0 +1,296 @@
+/*---------------------------------------------------------------------------------------------
+* Copyright (c) Bentley Systems, Incorporated. All rights reserved.
+* See LICENSE.md in the project root for license terms and full copyright notice.
+*--------------------------------------------------------------------------------------------*/
+/** @packageDocumentation
+ * @module Solid
+ */
+import { Arc3d } from "../curve/Arc3d";
+import { Loop } from "../curve/Loop";
+import { Path } from "../curve/Path";
+import { Geometry } from "../Geometry";
+import { Angle } from "../geometry3d/Angle";
+import { AngleSweep } from "../geometry3d/AngleSweep";
+import { Plane3dByOriginAndVectors } from "../geometry3d/Plane3dByOriginAndVectors";
+import { Vector2d } from "../geometry3d/Point2dVector2d";
+import { Transform } from "../geometry3d/Transform";
+import { SolidPrimitive } from "./SolidPrimitive";
+/**
+ * A torus pipe is a partial torus (donut).  In a local coordinate system
+ * * The z axis passes through the hole.
+ * * The "major hoop" arc has
+ *   * vectorTheta0 = (radiusA,0,0)
+ *   * vectorTheta90 = (0, radiusA,0)
+ *   * The major arc point at angle theta is `C(theta) = vectorTheta0 * cos(theta) + vectorTheta90 * sin(theta)
+ * * The minor hoop at theta various with phi "around the minor hoop"
+ *    * (x,y,z) = C(theta) + (radiusB *cos(theta), radiusB * sin(theta)) * cos(phi) + (0,radiusB,0) * sin(phi)
+ * * The stored form of the torus pipe is oriented for positive volume:
+ *   * Both radii are positive, with r0 >= r1 > 0
+ *   * The sweep is positive
+ *   * The coordinate system has positive determinant.
+ * * For uv parameterization,
+ *   * u is around the minor hoop, with (0..1) mapping to phi of (0 degrees ..360 degrees)
+ *   * v is along the major hoop with (0..1) mapping to theta of (0 .. sweep)
+ *   * a constant v section is a full circle
+ *   * a constant u section is an arc with sweep angle matching the torusPipe sweep angle.
+ * @public
+ */
+export class TorusPipe extends SolidPrimitive {
+    constructor(map, radiusA, radiusB, sweep, capped) {
+        super(capped);
+        /** String name for schema properties */
+        this.solidPrimitiveType = "torusPipe";
+        this._localToWorld = map;
+        this._radiusA = radiusA;
+        this._radiusB = radiusB;
+        this._sweep = sweep;
+        this._isReversed = false;
+    }
+    /** return a copy of the TorusPipe */
+    clone() {
+        const result = new TorusPipe(this._localToWorld.clone(), this._radiusA, this._radiusB, this._sweep.clone(), this.capped);
+        result._isReversed = this._isReversed;
+        return result;
+    }
+    /** Apply `transform` to the local coordinate system. */
+    tryTransformInPlace(transform) {
+        if (transform.matrix.isSingular())
+            return false;
+        transform.multiplyTransformTransform(this._localToWorld, this._localToWorld);
+        return true;
+    }
+    /** Clone this TorusPipe and transform the clone */
+    cloneTransformed(transform) {
+        const result = this.clone();
+        transform.multiplyTransformTransform(result._localToWorld, result._localToWorld);
+        return result;
+    }
+    /** Create a new `TorusPipe`
+     * @param frame local to world transformation
+     * @param majorRadius major hoop radius
+     * @param minorRadius minor hoop radius
+     * @param sweep sweep angle for major circle, with positive sweep from x axis towards y axis.
+     * @param capped true for circular caps
+     */
+    static createInFrame(frame, majorRadius, minorRadius, sweep, capped) {
+        // force near-zero radii to true zero
+        majorRadius = Math.abs(Geometry.correctSmallMetricDistance(majorRadius));
+        minorRadius = Math.abs(Geometry.correctSmallMetricDistance(minorRadius));
+        if (majorRadius < minorRadius)
+            return undefined;
+        if (majorRadius === 0.0)
+            return undefined;
+        if (minorRadius === 0.0)
+            return undefined;
+        if (sweep.isAlmostZero)
+            return undefined;
+        const xScale = 1.0;
+        let yScale = 1.0;
+        let zScale = 1.0;
+        if (frame.matrix.determinant() < 0.0)
+            zScale *= -1.0;
+        let isReversed = false;
+        const sweep1 = sweep.clone();
+        if (sweep.radians < 0.0) {
+            sweep1.setRadians(-sweep.radians);
+            zScale *= -1.0;
+            yScale *= -1.0;
+            isReversed = true;
+        }
+        const frame1 = frame.clone();
+        frame1.matrix.scaleColumns(xScale, yScale, zScale, frame1.matrix);
+        const result = new TorusPipe(frame1, majorRadius, minorRadius, sweep1, capped);
+        result._isReversed = isReversed;
+        return result;
+    }
+    /** Create a TorusPipe from the typical parameters of the Dgn file */
+    static createDgnTorusPipe(center, vectorX, vectorY, majorRadius, minorRadius, sweep, capped) {
+        const vectorZ = vectorX.crossProduct(vectorY);
+        vectorZ.scaleToLength(vectorX.magnitude(), vectorZ);
+        const frame = Transform.createOriginAndMatrixColumns(center, vectorX, vectorY, vectorZ);
+        return TorusPipe.createInFrame(frame, majorRadius, minorRadius, sweep, capped);
+    }
+    /** Create a TorusPipe from its primary arc and minor radius */
+    static createAlongArc(arc, minorRadius, capped) {
+        if (!Angle.isAlmostEqualRadiansAllowPeriodShift(0.0, arc.sweep.startRadians))
+            arc = arc.cloneInRotatedBasis(arc.sweep.startAngle);
+        const sweepRadians = arc.sweep.sweepRadians;
+        const data = arc.toScaledMatrix3d();
+        const frame = Transform.createOriginAndMatrix(data.center, data.axes);
+        return TorusPipe.createInFrame(frame, data.r0, minorRadius, Angle.createRadians(sweepRadians), capped);
+    }
+    /** Return a coordinate frame (right handed, unit axes)
+     * * origin at center of major circle
+     * * major circle in xy plane
+     * * z axis perpendicular
+     */
+    getConstructiveFrame() {
+        return this._localToWorld.cloneRigid();
+    }
+    /** Return the center of the torus pipe (inside the donut hole) */
+    cloneCenter() { return this._localToWorld.getOrigin(); }
+    /** return the vector along the x axis (in the major hoops plane) */
+    cloneVectorX() { return this._localToWorld.matrix.columnX(); }
+    /** return the vector along the y axis (in the major hoop plane) */
+    cloneVectorY() { return this._localToWorld.matrix.columnY(); }
+    /** get the minor hoop radius (`radiusA`) */
+    getMinorRadius() { return this._radiusB; }
+    /** get the major hoop radius (`radiusB`) */
+    getMajorRadius() { return this._radiusA; }
+    /** get the sweep angle along the major circle. */
+    getSweepAngle() { return this._sweep.clone(); }
+    /** Ask if this TorusPipe is labeled as reversed */
+    getIsReversed() { return this._isReversed; }
+    /** Return the sweep angle as a fraction of full 360 degrees (2PI radians) */
+    getThetaFraction() { return this._sweep.radians / (Math.PI * 2.0); }
+    /** ask if `other` is an instance of `TorusPipe` */
+    isSameGeometryClass(other) { return other instanceof TorusPipe; }
+    /** test if `this` and `other` have nearly equal geometry */
+    isAlmostEqual(other) {
+        if (other instanceof TorusPipe) {
+            if ((!this._sweep.isFullCircle) && this.capped !== other.capped)
+                return false;
+            if (!this._localToWorld.isAlmostEqual(other._localToWorld))
+                return false;
+            return Geometry.isSameCoordinate(this._radiusA, other._radiusA)
+                && Geometry.isSameCoordinate(this._radiusB, other._radiusB)
+                && this._sweep.isAlmostEqualNoPeriodShift(other._sweep);
+        }
+        return false;
+    }
+    /** Return the angle (in radians) for given fractional position around the major hoop.
+     */
+    vFractionToRadians(v) { return this._sweep.radians * v; }
+    /** Second step of double dispatch:  call `handler.handleTorusPipe(this)` */
+    dispatchToGeometryHandler(handler) {
+        return handler.handleTorusPipe(this);
+    }
+    /**
+     * Return the Arc3d section at vFraction.  For the TorusPipe, this is a minor circle.
+     * @param vFraction fractional position along the sweep direction
+     */
+    constantVSection(v) {
+        const thetaRadians = this.vFractionToRadians(v);
+        const c0 = Math.cos(thetaRadians);
+        const s0 = Math.sin(thetaRadians);
+        const majorRadius = this.getMajorRadius();
+        const minorRadius = this.getMinorRadius();
+        const center = this._localToWorld.multiplyXYZ(majorRadius * c0, majorRadius * s0, 0);
+        const vector0 = this._localToWorld.multiplyVectorXYZ(minorRadius * c0, minorRadius * s0, 0);
+        const vector90 = this._localToWorld.multiplyVectorXYZ(0, 0, minorRadius);
+        return Loop.create(Arc3d.create(center, vector0, vector90));
+    }
+    /** Return an arc at constant u, and arc sweep  matching this TorusPipe sweep. */
+    constantUSection(uFraction) {
+        const theta1Radians = this._sweep.radians;
+        const phiRadians = uFraction * Math.PI;
+        const majorRadius = this.getMajorRadius();
+        const minorRadius = this.getMinorRadius();
+        const transform = this._localToWorld;
+        const axes = transform.matrix;
+        const center = this._localToWorld.multiplyXYZ(0, 0, minorRadius * Math.sin(phiRadians));
+        const rxy = majorRadius + minorRadius * Math.cos(phiRadians);
+        const vector0 = axes.multiplyXYZ(rxy, 0, 0);
+        const vector90 = axes.multiplyXYZ(0, rxy, 0);
+        return Path.create(Arc3d.create(center, vector0, vector90, AngleSweep.createStartEndRadians(0.0, theta1Radians)));
+    }
+    /** extend `rangeToExtend` to include this `TorusPipe` */
+    extendRange(rangeToExtend, transform) {
+        const theta1Radians = this._sweep.radians;
+        const majorRadius = this.getMajorRadius();
+        const minorRadius = this.getMinorRadius();
+        const transform0 = this._localToWorld;
+        const numThetaSample = Math.ceil(theta1Radians / (Math.PI / 16.0));
+        const numHalfPhiSample = 16;
+        let phi0 = 0;
+        let dPhi = 0;
+        let numPhiSample = 0;
+        let theta = 0;
+        let cosTheta = 0;
+        let sinTheta = 0;
+        let rxy = 0;
+        let phi = 0;
+        let j = 0;
+        const dTheta = theta1Radians / numThetaSample;
+        for (let i = 0; i <= numThetaSample; i++) {
+            theta = i * dTheta;
+            cosTheta = Math.cos(theta);
+            sinTheta = Math.sin(theta);
+            // At the ends, do the entire phi circle.
+            // Otherwise only do the outer half
+            if (i === 0 || i === numThetaSample) {
+                phi0 = -Math.PI;
+                dPhi = 2.0 * Math.PI / numHalfPhiSample;
+                numPhiSample = numHalfPhiSample;
+            }
+            else {
+                phi0 = -0.5 * Math.PI;
+                dPhi = Math.PI / numHalfPhiSample;
+                numPhiSample = 2 * numHalfPhiSample - 1;
+            }
+            if (transform) {
+                for (j = 0; j <= numPhiSample; j++) {
+                    phi = phi0 + j * dPhi;
+                    rxy = majorRadius + minorRadius * Math.cos(phi);
+                    rangeToExtend.extendTransformTransformedXYZ(transform, transform0, cosTheta * rxy, sinTheta * rxy, Math.sin(phi) * minorRadius);
+                }
+            }
+            else {
+                for (j = 0; j <= numPhiSample; j++) {
+                    phi = phi0 + j * dPhi;
+                    rxy = majorRadius + minorRadius * Math.sin(phi);
+                    rangeToExtend.extendTransformedXYZ(transform0, cosTheta * rxy, sinTheta * rxy, Math.sin(phi) * minorRadius);
+                }
+            }
+        }
+    }
+    /** Evaluate as a uv surface
+     * @param u fractional position in minor (phi)
+     * @param v fractional position on major (theta) arc
+     */
+    uvFractionToPoint(u, v, result) {
+        const thetaRadians = v * this._sweep.radians;
+        const phiRadians = u * Math.PI * 2.0;
+        const cosTheta = Math.cos(thetaRadians);
+        const sinTheta = Math.sin(thetaRadians);
+        const minorRadius = this.getMinorRadius();
+        const rxy = this.getMajorRadius() + Math.cos(phiRadians) * minorRadius;
+        return this._localToWorld.multiplyXYZ(rxy * cosTheta, rxy * sinTheta, minorRadius * Math.sin(phiRadians), result);
+    }
+    /** Evaluate as a uv surface, returning point and two vectors.
+     * @param u fractional position in minor (phi)
+     * @param v fractional position on major (theta) arc
+     */
+    uvFractionToPointAndTangents(u, v, result) {
+        const thetaRadians = v * this._sweep.radians;
+        const phiRadians = u * Math.PI * 2.0;
+        const fTheta = this._sweep.radians;
+        const fPhi = Math.PI * 2.0;
+        const cosTheta = Math.cos(thetaRadians);
+        const sinTheta = Math.sin(thetaRadians);
+        const sinPhi = Math.sin(phiRadians);
+        const cosPhi = Math.cos(phiRadians);
+        const minorRadius = this.getMinorRadius();
+        const rxy = this.getMajorRadius() + Math.cos(phiRadians) * minorRadius;
+        const rSinPhi = minorRadius * sinPhi;
+        const rCosPhi = minorRadius * cosPhi; // appears only as derivative of rSinPhi.
+        return Plane3dByOriginAndVectors.createOriginAndVectors(this._localToWorld.multiplyXYZ(cosTheta * rxy, sinTheta * rxy, rSinPhi), this._localToWorld.multiplyVectorXYZ(-cosTheta * rSinPhi * fPhi, -sinTheta * rSinPhi * fPhi, rCosPhi * fPhi), this._localToWorld.multiplyVectorXYZ(-rxy * sinTheta * fTheta, rxy * cosTheta * fTheta, 0), result);
+    }
+    /**
+     * Directional distance query
+     * * u direction is around the (full) minor hoop
+     * * v direction is around the outer radius, sum of (absolute values of) major and minor radii.
+     */
+    maxIsoParametricDistance() {
+        const a = Math.abs(this.getMajorRadius());
+        const b = Math.abs(this.getMinorRadius());
+        return Vector2d.create(b * Math.PI * 2.0, (a + b) * this._sweep.radians);
+    }
+    /**
+     * @return true if this is a closed volume.
+     */
+    get isClosedVolume() {
+        return this.capped || this._sweep.isFullCircle;
+    }
+}
+//# sourceMappingURL=TorusPipe.js.map

package/lib/esm/solid/TorusPipe.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"TorusPipe.js","sourceRoot":"","sources":["../../../src/solid/TorusPipe.ts"],"names":[],"mappings":"AAAA;;;+FAG+F;AAE/F;;GAEG;AAEH,OAAO,EAAE,KAAK,EAAE,MAAM,gBAAgB,CAAC;AAGvC,OAAO,EAAE,IAAI,EAAE,MAAM,eAAe,CAAC;AACrC,OAAO,EAAE,IAAI,EAAE,MAAM,eAAe,CAAC;AACrC,OAAO,EAAE,QAAQ,EAAE,MAAM,aAAa,CAAC;AACvC,OAAO,EAAE,KAAK,EAAE,MAAM,qBAAqB,CAAC;AAC5C,OAAO,EAAE,UAAU,EAAE,MAAM,0BAA0B,CAAC;AAEtD,OAAO,EAAE,yBAAyB,EAAE,MAAM,yCAAyC,CAAC;AACpF,OAAO,EAAE,QAAQ,EAAE,MAAM,+BAA+B,CAAC;AAGzD,OAAO,EAAE,SAAS,EAAE,MAAM,yBAAyB,CAAC;AACpD,OAAO,EAAE,cAAc,EAAE,MAAM,kBAAkB,CAAC;AAElD;;;;;;;;;;;;;;;;;;;GAmBG;AACH,MAAM,OAAO,SAAU,SAAQ,cAAc;IAU3C,YAAsB,GAAc,EAAE,OAAe,EAAE,OAAe,EAAE,KAAY,EAAE,MAAe;QACnG,KAAK,CAAC,MAAM,CAAC,CAAC;QAVhB,wCAAwC;QACxB,uBAAkB,GAAG,WAAW,CAAC;QAU/C,IAAI,CAAC,aAAa,GAAG,GAAG,CAAC;QACzB,IAAI,CAAC,QAAQ,GAAG,OAAO,CAAC;QACxB,IAAI,CAAC,QAAQ,GAAG,OAAO,CAAC;QACxB,IAAI,CAAC,MAAM,GAAG,KAAK,CAAC;QACpB,IAAI,CAAC,WAAW,GAAG,KAAK,CAAC;IAC3B,CAAC;IACD,qCAAqC;IAC9B,KAAK;QACV,MAAM,MAAM,GAAG,IAAI,SAAS,CAAC,IAAI,CAAC,aAAa,CAAC,KAAK,EAAE,EAAE,IAAI,CAAC,QAAQ,EAAE,IAAI,CAAC,QAAQ,EAAE,IAAI,CAAC,MAAM,CAAC,KAAK,EAAE,EAAE,IAAI,CAAC,MAAM,CAAC,CAAC;QACzH,MAAM,CAAC,WAAW,GAAG,IAAI,CAAC,WAAW,CAAC;QACtC,OAAO,MAAM,CAAC;IAChB,CAAC;IACD,wDAAwD;IACjD,mBAAmB,CAAC,SAAoB;QAC7C,IAAI,SAAS,CAAC,MAAM,CAAC,UAAU,EAAE;YAC/B,OAAO,KAAK,CAAC;QACf,SAAS,CAAC,0BAA0B,CAAC,IAAI,CAAC,aAAa,EAAE,IAAI,CAAC,aAAa,CAAC,CAAC;QAC7E,OAAO,IAAI,CAAC;IACd,CAAC;IACD,mDAAmD;IAC5C,gBAAgB,CAAC,SAAoB;QAC1C,MAAM,MAAM,GAAG,IAAI,CAAC,KAAK,EAAE,CAAC;QAC5B,SAAS,CAAC,0BAA0B,CAAC,MAAM,CAAC,aAAa,EAAE,MAAM,CAAC,aAAa,CAAC,CAAC;QACjF,OAAO,MAAM,CAAC;IAChB,CAAC;IACD;;;;;;OAMG;IACI,MAAM,CAAC,aAAa,CAAC,KAAgB,EAAE,WAAmB,EAAE,WAAmB,EAAE,KAAY,EAAE,MAAe;QACnH,qCAAqC;QACrC,WAAW,GAAG,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,0BAA0B,CAAC,WAAW,CAAC,CAAC,CAAC;QACzE,WAAW,GAAG,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,0BAA0B,CAAC,WAAW,CAAC,CAAC,CAAC;QACzE,IAAI,WAAW,GAAG,WAAW;YAAE,OAAO,SAAS,CAAC;QAChD,IAAI,WAAW,KAAK,GAAG;YAAE,OAAO,SAAS,CAAC;QAC1C,IAAI,WAAW,KAAK,GAAG;YAAE,OAAO,SAAS,CAAC;QAE1C,IAAI,KAAK,CAAC,YAAY;YAAE,OAAO,SAAS,CAAC;QACzC,MAAM,MAAM,GAAG,GAAG,CAAC;QACnB,IAAI,MAAM,GAAG,GAAG,CAAC;QACjB,IAAI,MAAM,GAAG,GAAG,CAAC;QACjB,IAAI,KAAK,CAAC,MAAM,CAAC,WAAW,EAAE,GAAG,GAAG;YAAE,MAAM,IAAI,CAAC,GAAG,CAAC;QACrD,IAAI,UAAU,GAAG,KAAK,CAAC;QACvB,MAAM,MAAM,GAAG,KAAK,CAAC,KAAK,EAAE,CAAC;QAC7B,IAAI,KAAK,CAAC,OAAO,GAAG,GAAG,EAAE;YACvB,MAAM,CAAC,UAAU,CAAC,CAAC,KAAK,CAAC,OAAO,CAAC,CAAC;YAClC,MAAM,IAAI,CAAC,GAAG,CAAC;YACf,MAAM,IAAI,CAAC,GAAG,CAAC;YACf,UAAU,GAAG,IAAI,CAAC;SACnB;QACD,MAAM,MAAM,GAAG,KAAK,CAAC,KAAK,EAAE,CAAC;QAC7B,MAAM,CAAC,MAAM,CAAC,YAAY,CAAC,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,CAAC,MAAM,CAAC,CAAC;QAClE,MAAM,MAAM,GAAG,IAAI,SAAS,CAAC,MAAM,EAAE,WAAW,EAAE,WAAW,EAAE,MAAM,EAAE,MAAM,CAAC,CAAC;QAC/E,MAAM,CAAC,WAAW,GAAG,UAAU,CAAC;QAChC,OAAO,MAAM,CAAC;IAChB,CAAC;IAED,qEAAqE;IAC9D,MAAM,CAAC,kBAAkB,CAAC,MAAe,EAAE,OAAiB,EAAE,OAAiB,EAAE,WAAmB,EAAE,WAAmB,EAC9H,KAAY,EAAE,MAAe;QAC7B,MAAM,OAAO,GAAG,OAAO,CAAC,YAAY,CAAC,OAAO,CAAC,CAAC;QAC9C,OAAO,CAAC,aAAa,CAAC,OAAO,CAAC,SAAS,EAAE,EAAE,OAAO,CAAC,CAAC;QACpD,MAAM,KAAK,GAAG,SAAS,CAAC,4BAA4B,CAAC,MAAM,EAAE,OAAO,EAAE,OAAO,EAAE,OAAO,CAAC,CAAC;QACxF,OAAO,SAAS,CAAC,aAAa,CAAC,KAAK,EAAE,WAAW,EAAE,WAAW,EAAE,KAAK,EAAE,MAAM,CAAC,CAAC;IACjF,CAAC;IACD,+DAA+D;IACxD,MAAM,CAAC,cAAc,CAAC,GAAU,EAAE,WAAmB,EAAE,MAAe;QAC3E,IAAI,CAAC,KAAK,CAAC,oCAAoC,CAAC,GAAG,EAAE,GAAG,CAAC,KAAK,CAAC,YAAY,CAAC;YAC1E,GAAG,GAAG,GAAG,CAAC,mBAAmB,CAAC,GAAG,CAAC,KAAK,CAAC,UAAU,CAAC,CAAC;QACtD,MAAM,YAAY,GAAG,GAAG,CAAC,KAAK,CAAC,YAAY,CAAC;QAC5C,MAAM,IAAI,GAAG,GAAG,CAAC,gBAAgB,EAAE,CAAC;QACpC,MAAM,KAAK,GAAG,SAAS,CAAC,qBAAqB,CAAC,IAAI,CAAC,MAAM,EAAE,IAAI,CAAC,IAAI,CAAC,CAAC;QACtE,OAAO,SAAS,CAAC,aAAa,CAAC,KAAK,EAAE,IAAI,CAAC,EAAE,EAAE,WAAW,EAAE,KAAK,CAAC,aAAa,CAAC,YAAY,CAAC,EAAE,MAAM,CAAC,CAAC;IACzG,CAAC;IAED;;;;OAIG;IACI,oBAAoB;QACzB,OAAO,IAAI,CAAC,aAAa,CAAC,UAAU,EAAE,CAAC;IACzC,CAAC;IACD,kEAAkE;IAC3D,WAAW,KAAc,OAAO,IAAI,CAAC,aAAa,CAAC,SAAS,EAAE,CAAC,CAAC,CAAC;IACxE,oEAAoE;IAC7D,YAAY,KAAe,OAAO,IAAI,CAAC,aAAa,CAAC,MAAM,CAAC,OAAO,EAAE,CAAC,CAAC,CAAC;IAC/E,mEAAmE;IAC5D,YAAY,KAAe,OAAO,IAAI,CAAC,aAAa,CAAC,MAAM,CAAC,OAAO,EAAE,CAAC,CAAC,CAAC;IAC/E,4CAA4C;IACrC,cAAc,KAAa,OAAO,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC;IACzD,4CAA4C;IACrC,cAAc,KAAa,OAAO,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC;IACzD,kDAAkD;IAC3C,aAAa,KAAY,OAAO,IAAI,CAAC,MAAM,CAAC,KAAK,EAAE,CAAC,CAAC,CAAC;IAC7D,mDAAmD;IAC5C,aAAa,KAAc,OAAO,IAAI,CAAC,WAAW,CAAC,CAAC,CAAC;IAC5D,6EAA6E;IACtE,gBAAgB,KAAa,OAAO,IAAI,CAAC,MAAM,CAAC,OAAO,GAAG,CAAC,IAAI,CAAC,EAAE,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC;IACnF,mDAAmD;IAC5C,mBAAmB,CAAC,KAAU,IAAa,OAAO,KAAK,YAAY,SAAS,CAAC,CAAC,CAAC;IACtF,4DAA4D;IAC5C,aAAa,CAAC,KAAoB;QAChD,IAAI,KAAK,YAAY,SAAS,EAAE;YAC9B,IAAI,CAAC,CAAC,IAAI,CAAC,MAAM,CAAC,YAAY,CAAC,IAAI,IAAI,CAAC,MAAM,KAAK,KAAK,CAAC,MAAM;gBAAE,OAAO,KAAK,CAAC;YAC9E,IAAI,CAAC,IAAI,CAAC,aAAa,CAAC,aAAa,CAAC,KAAK,CAAC,aAAa,CAAC;gBAAE,OAAO,KAAK,CAAC;YACzE,OAAO,QAAQ,CAAC,gBAAgB,CAAC,IAAI,CAAC,QAAQ,EAAE,KAAK,CAAC,QAAQ,CAAC;mBAC1D,QAAQ,CAAC,gBAAgB,CAAC,IAAI,CAAC,QAAQ,EAAE,KAAK,CAAC,QAAQ,CAAC;mBACxD,IAAI,CAAC,MAAM,CAAC,0BAA0B,CAAC,KAAK,CAAC,MAAM,CAAC,CAAC;SAC3D;QACD,OAAO,KAAK,CAAC;IACf,CAAC;IACD;OACG;IACI,kBAAkB,CAAC,CAAS,IAAY,OAAO,IAAI,CAAC,MAAM,CAAC,OAAO,GAAG,CAAC,CAAC,CAAC,CAAC;IAChF,4EAA4E;IACrE,yBAAyB,CAAC,OAAwB;QACvD,OAAO,OAAO,CAAC,eAAe,CAAC,IAAI,CAAC,CAAC;IACvC,CAAC;IAED;;;OAGG;IACI,gBAAgB,CAAC,CAAS;QAC/B,MAAM,YAAY,GAAG,IAAI,CAAC,kBAAkB,CAAC,CAAC,CAAC,CAAC;QAChD,MAAM,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,YAAY,CAAC,CAAC;QAClC,MAAM,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,YAAY,CAAC,CAAC;QAClC,MAAM,WAAW,GAAG,IAAI,CAAC,cAAc,EAAE,CAAC;QAC1C,MAAM,WAAW,GAAG,IAAI,CAAC,cAAc,EAAE,CAAC;QAE1C,MAAM,MAAM,GAAG,IAAI,CAAC,aAAa,CAAC,WAAW,CAAC,WAAW,GAAG,EAAE,EAAE,WAAW,GAAG,EAAE,EAAE,CAAC,CAAC,CAAC;QACrF,MAAM,OAAO,GAAG,IAAI,CAAC,aAAa,CAAC,iBAAiB,CAAC,WAAW,GAAG,EAAE,EAAE,WAAW,GAAG,EAAE,EAAE,CAAC,CAAC,CAAC;QAC5F,MAAM,QAAQ,GAAG,IAAI,CAAC,aAAa,CAAC,iBAAiB,CAAC,CAAC,EAAE,CAAC,EAAE,WAAW,CAAC,CAAC;QACzE,OAAO,IAAI,CAAC,MAAM,CAAC,KAAK,CAAC,MAAM,CAAC,MAAM,EAAE,OAAO,EAAE,QAAQ,CAAC,CAAC,CAAC;IAC9D,CAAC;IACD,iFAAiF;IAC1E,gBAAgB,CAAC,SAAiB;QACvC,MAAM,aAAa,GAAG,IAAI,CAAC,MAAM,CAAC,OAAO,CAAC;QAC1C,MAAM,UAAU,GAAG,SAAS,GAAG,IAAI,CAAC,EAAE,CAAC;QAEvC,MAAM,WAAW,GAAG,IAAI,CAAC,cAAc,EAAE,CAAC;QAC1C,MAAM,WAAW,GAAG,IAAI,CAAC,cAAc,EAAE,CAAC;QAC1C,MAAM,SAAS,GAAG,IAAI,CAAC,aAAa,CAAC;QACrC,MAAM,IAAI,GAAG,SAAS,CAAC,MAAM,CAAC;QAC9B,MAAM,MAAM,GAAG,IAAI,CAAC,aAAa,CAAC,WAAW,CAAC,CAAC,EAAE,CAAC,EAAE,WAAW,GAAG,IAAI,CAAC,GAAG,CAAC,UAAU,CAAC,CAAC,CAAC;QACxF,MAAM,GAAG,GAAG,WAAW,GAAG,WAAW,GAAG,IAAI,CAAC,GAAG,CAAC,UAAU,CAAC,CAAC;QAC7D,MAAM,OAAO,GAAG,IAAI,CAAC,WAAW,CAAC,GAAG,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;QAC5C,MAAM,QAAQ,GAAG,IAAI,CAAC,WAAW,CAAC,CAAC,EAAE,GAAG,EAAE,CAAC,CAAC,CAAC;QAC7C,OAAO,IAAI,CAAC,MAAM,CAAC,KAAK,CAAC,MAAM,CAAC,MAAM,EAAE,OAAO,EAAE,QAAQ,EAAE,UAAU,CAAC,qBAAqB,CAAC,GAAG,EAAE,aAAa,CAAC,CAAC,CAAC,CAAC;IACpH,CAAC;IACD,yDAAyD;IAClD,WAAW,CAAC,aAAsB,EAAE,SAAqB;QAC9D,MAAM,aAAa,GAAG,IAAI,CAAC,MAAM,CAAC,OAAO,CAAC;QAC1C,MAAM,WAAW,GAAG,IAAI,CAAC,cAAc,EAAE,CAAC;QAC1C,MAAM,WAAW,GAAG,IAAI,CAAC,cAAc,EAAE,CAAC;QAC1C,MAAM,UAAU,GAAG,IAAI,CAAC,aAAa,CAAC;QACtC,MAAM,cAAc,GAAG,IAAI,CAAC,IAAI,CAAC,aAAa,GAAG,CAAC,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,CAAC,CAAC;QACnE,MAAM,gBAAgB,GAAG,EAAE,CAAC;QAC5B,IAAI,IAAI,GAAG,CAAC,CAAC;QACb,IAAI,IAAI,GAAG,CAAC,CAAC;QACb,IAAI,YAAY,GAAG,CAAC,CAAC;QACrB,IAAI,KAAK,GAAG,CAAC,CAAC;QACd,IAAI,QAAQ,GAAG,CAAC,CAAC;QACjB,IAAI,QAAQ,GAAG,CAAC,CAAC;QACjB,IAAI,GAAG,GAAG,CAAC,CAAC;QACZ,IAAI,GAAG,GAAG,CAAC,CAAC;QACZ,IAAI,CAAC,GAAG,CAAC,CAAC;QACV,MAAM,MAAM,GAAG,aAAa,GAAG,cAAc,CAAC;QAC9C,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,cAAc,EAAE,CAAC,EAAE,EAAE;YACxC,KAAK,GAAG,CAAC,GAAG,MAAM,CAAC;YACnB,QAAQ,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC;YAC3B,QAAQ,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC;YAC3B,yCAAyC;YACzC,mCAAmC;YACnC,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,KAAK,cAAc,EAAE;gBACnC,IAAI,GAAG,CAAC,IAAI,CAAC,EAAE,CAAC;gBAChB,IAAI,GAAG,GAAG,GAAG,IAAI,CAAC,EAAE,GAAG,gBAAgB,CAAC;gBACxC,YAAY,GAAG,gBAAgB,CAAC;aACjC;iBAAM;gBACL,IAAI,GAAG,CAAC,GAAG,GAAG,IAAI,CAAC,EAAE,CAAC;gBACtB,IAAI,GAAG,IAAI,CAAC,EAAE,GAAG,gBAAgB,CAAC;gBAClC,YAAY,GAAG,CAAC,GAAG,gBAAgB,GAAG,CAAC,CAAC;aACzC;YACD,IAAI,SAAS,EAAE;gBACb,KAAK,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,YAAY,EAAE,CAAC,EAAE,EAAE;oBAClC,GAAG,GAAG,IAAI,GAAG,CAAC,GAAG,IAAI,CAAC;oBACtB,GAAG,GAAG,WAAW,GAAG,WAAW,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;oBAChD,aAAa,CAAC,6BAA6B,CAAC,SAAS,EAAE,UAAU,EAC/D,QAAQ,GAAG,GAAG,EAAE,QAAQ,GAAG,GAAG,EAC9B,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,WAAW,CAAC,CAAC;iBAChC;aACF;iBAAM;gBACL,KAAK,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,YAAY,EAAE,CAAC,EAAE,EAAE;oBAClC,GAAG,GAAG,IAAI,GAAG,CAAC,GAAG,IAAI,CAAC;oBACtB,GAAG,GAAG,WAAW,GAAG,WAAW,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;oBAChD,aAAa,CAAC,oBAAoB,CAAC,UAAU,EAC3C,QAAQ,GAAG,GAAG,EAAE,QAAQ,GAAG,GAAG,EAC9B,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,WAAW,CAAC,CAAC;iBAChC;aACF;SACF;IACH,CAAC;IACD;;;OAGG;IACI,iBAAiB,CAAC,CAAS,EAAE,CAAS,EAAE,MAAgB;QAC7D,MAAM,YAAY,GAAG,CAAC,GAAG,IAAI,CAAC,MAAM,CAAC,OAAO,CAAC;QAC7C,MAAM,UAAU,GAAG,CAAC,GAAG,IAAI,CAAC,EAAE,GAAG,GAAG,CAAC;QACrC,MAAM,QAAQ,GAAG,IAAI,CAAC,GAAG,CAAC,YAAY,CAAC,CAAC;QACxC,MAAM,QAAQ,GAAG,IAAI,CAAC,GAAG,CAAC,YAAY,CAAC,CAAC;QACxC,MAAM,WAAW,GAAG,IAAI,CAAC,cAAc,EAAE,CAAC;QAC1C,MAAM,GAAG,GAAG,IAAI,CAAC,cAAc,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,UAAU,CAAC,GAAG,WAAW,CAAC;QACvE,OAAO,IAAI,CAAC,aAAa,CAAC,WAAW,CAAC,GAAG,GAAG,QAAQ,EAAE,GAAG,GAAG,QAAQ,EAAE,WAAW,GAAG,IAAI,CAAC,GAAG,CAAC,UAAU,CAAC,EAAE,MAAM,CAAC,CAAC;IACpH,CAAC;IACD;;;OAGG;IACI,4BAA4B,CAAC,CAAS,EAAE,CAAS,EAAE,MAAkC;QAC1F,MAAM,YAAY,GAAG,CAAC,GAAG,IAAI,CAAC,MAAM,CAAC,OAAO,CAAC;QAC7C,MAAM,UAAU,GAAG,CAAC,GAAG,IAAI,CAAC,EAAE,GAAG,GAAG,CAAC;QACrC,MAAM,MAAM,GAAG,IAAI,CAAC,MAAM,CAAC,OAAO,CAAC;QACnC,MAAM,IAAI,GAAG,IAAI,CAAC,EAAE,GAAG,GAAG,CAAC;QAC3B,MAAM,QAAQ,GAAG,IAAI,CAAC,GAAG,CAAC,YAAY,CAAC,CAAC;QACxC,MAAM,QAAQ,GAAG,IAAI,CAAC,GAAG,CAAC,YAAY,CAAC,CAAC;QACxC,MAAM,MAAM,GAAG,IAAI,CAAC,GAAG,CAAC,UAAU,CAAC,CAAC;QACpC,MAAM,MAAM,GAAG,IAAI,CAAC,GAAG,CAAC,UAAU,CAAC,CAAC;QACpC,MAAM,WAAW,GAAG,IAAI,CAAC,cAAc,EAAE,CAAC;QAC1C,MAAM,GAAG,GAAG,IAAI,CAAC,cAAc,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,UAAU,CAAC,GAAG,WAAW,CAAC;QACvE,MAAM,OAAO,GAAG,WAAW,GAAG,MAAM,CAAC;QACrC,MAAM,OAAO,GAAG,WAAW,GAAG,MAAM,CAAC,CAAG,yCAAyC;QACjF,OAAO,yBAAyB,CAAC,sBAAsB,CACrD,IAAI,CAAC,aAAa,CAAC,WAAW,CAAC,QAAQ,GAAG,GAAG,EAAE,QAAQ,GAAG,GAAG,EAAE,OAAO,CAAC,EACvE,IAAI,CAAC,aAAa,CAAC,iBAAiB,CAAC,CAAC,QAAQ,GAAG,OAAO,GAAG,IAAI,EAAE,CAAC,QAAQ,GAAG,OAAO,GAAG,IAAI,EAAE,OAAO,GAAG,IAAI,CAAC,EAC5G,IAAI,CAAC,aAAa,CAAC,iBAAiB,CAAC,CAAC,GAAG,GAAG,QAAQ,GAAG,MAAM,EAAE,GAAG,GAAG,QAAQ,GAAG,MAAM,EAAE,CAAC,CAAC,EAC1F,MAAM,CAAC,CAAC;IACZ,CAAC;IACD;;;;OAIG;IACI,wBAAwB;QAC7B,MAAM,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,cAAc,EAAE,CAAC,CAAC;QAC1C,MAAM,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,cAAc,EAAE,CAAC,CAAC;QAC1C,OAAO,QAAQ,CAAC,MAAM,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,GAAG,GAAG,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,MAAM,CAAC,OAAO,CAAC,CAAC;IAC3E,CAAC;IACD;;OAEG;IACH,IAAW,cAAc;QACvB,OAAO,IAAI,CAAC,MAAM,IAAI,IAAI,CAAC,MAAM,CAAC,YAAY,CAAC;IACjD,CAAC;CAEF","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 Solid\r\n */\r\n\r\nimport { Arc3d } from \"../curve/Arc3d\";\r\nimport { CurveCollection } from \"../curve/CurveCollection\";\r\nimport { GeometryQuery } from \"../curve/GeometryQuery\";\r\nimport { Loop } from \"../curve/Loop\";\r\nimport { Path } from \"../curve/Path\";\r\nimport { Geometry } from \"../Geometry\";\r\nimport { Angle } from \"../geometry3d/Angle\";\r\nimport { AngleSweep } from \"../geometry3d/AngleSweep\";\r\nimport { GeometryHandler, UVSurface, UVSurfaceIsoParametricDistance } from \"../geometry3d/GeometryHandler\";\r\nimport { Plane3dByOriginAndVectors } from \"../geometry3d/Plane3dByOriginAndVectors\";\r\nimport { Vector2d } from \"../geometry3d/Point2dVector2d\";\r\nimport { Point3d, Vector3d } from \"../geometry3d/Point3dVector3d\";\r\nimport { Range3d } from \"../geometry3d/Range\";\r\nimport { Transform } from \"../geometry3d/Transform\";\r\nimport { SolidPrimitive } from \"./SolidPrimitive\";\r\n\r\n/**\r\n * A torus pipe is a partial torus (donut).  In a local coordinate system\r\n * * The z axis passes through the hole.\r\n * * The \"major hoop\" arc has\r\n *   * vectorTheta0 = (radiusA,0,0)\r\n *   * vectorTheta90 = (0, radiusA,0)\r\n *   * The major arc point at angle theta is `C(theta) = vectorTheta0 * cos(theta) + vectorTheta90 * sin(theta)\r\n * * The minor hoop at theta various with phi \"around the minor hoop\"\r\n *    * (x,y,z) = C(theta) + (radiusB *cos(theta), radiusB * sin(theta)) * cos(phi) + (0,radiusB,0) * sin(phi)\r\n * * The stored form of the torus pipe is oriented for positive volume:\r\n *   * Both radii are positive, with r0 >= r1 > 0\r\n *   * The sweep is positive\r\n *   * The coordinate system has positive determinant.\r\n * * For uv parameterization,\r\n *   * u is around the minor hoop, with (0..1) mapping to phi of (0 degrees ..360 degrees)\r\n *   * v is along the major hoop with (0..1) mapping to theta of (0 .. sweep)\r\n *   * a constant v section is a full circle\r\n *   * a constant u section is an arc with sweep angle matching the torusPipe sweep angle.\r\n * @public\r\n */\r\nexport class TorusPipe extends SolidPrimitive implements UVSurface, UVSurfaceIsoParametricDistance {\r\n  /** String name for schema properties */\r\n  public readonly solidPrimitiveType = \"torusPipe\";\r\n\r\n  private _localToWorld: Transform;\r\n  private _radiusA: number;  // radius of (large) circle in xy plane\r\n  private _radiusB: number;  // radius of (small) circle in xz plane.\r\n  private _sweep: Angle;\r\n  private _isReversed: boolean;\r\n\r\n  protected constructor(map: Transform, radiusA: number, radiusB: number, sweep: Angle, capped: boolean) {\r\n    super(capped);\r\n    this._localToWorld = map;\r\n    this._radiusA = radiusA;\r\n    this._radiusB = radiusB;\r\n    this._sweep = sweep;\r\n    this._isReversed = false;\r\n  }\r\n  /** return a copy of the TorusPipe */\r\n  public clone(): TorusPipe {\r\n    const result = new TorusPipe(this._localToWorld.clone(), this._radiusA, this._radiusB, this._sweep.clone(), this.capped);\r\n    result._isReversed = this._isReversed;\r\n    return result;\r\n  }\r\n  /** Apply `transform` to the local coordinate system. */\r\n  public tryTransformInPlace(transform: Transform): boolean {\r\n    if (transform.matrix.isSingular())\r\n      return false;\r\n    transform.multiplyTransformTransform(this._localToWorld, this._localToWorld);\r\n    return true;\r\n  }\r\n  /** Clone this TorusPipe and transform the clone */\r\n  public cloneTransformed(transform: Transform): TorusPipe | undefined {\r\n    const result = this.clone();\r\n    transform.multiplyTransformTransform(result._localToWorld, result._localToWorld);\r\n    return result;\r\n  }\r\n  /** Create a new `TorusPipe`\r\n   * @param frame local to world transformation\r\n   * @param majorRadius major hoop radius\r\n   * @param minorRadius minor hoop radius\r\n   * @param sweep sweep angle for major circle, with positive sweep from x axis towards y axis.\r\n   * @param capped true for circular caps\r\n   */\r\n  public static createInFrame(frame: Transform, majorRadius: number, minorRadius: number, sweep: Angle, capped: boolean): TorusPipe | undefined {\r\n    // force near-zero radii to true zero\r\n    majorRadius = Math.abs(Geometry.correctSmallMetricDistance(majorRadius));\r\n    minorRadius = Math.abs(Geometry.correctSmallMetricDistance(minorRadius));\r\n    if (majorRadius < minorRadius) return undefined;\r\n    if (majorRadius === 0.0) return undefined;\r\n    if (minorRadius === 0.0) return undefined;\r\n\r\n    if (sweep.isAlmostZero) return undefined;\r\n    const xScale = 1.0;\r\n    let yScale = 1.0;\r\n    let zScale = 1.0;\r\n    if (frame.matrix.determinant() < 0.0) zScale *= -1.0;\r\n    let isReversed = false;\r\n    const sweep1 = sweep.clone();\r\n    if (sweep.radians < 0.0) {\r\n      sweep1.setRadians(-sweep.radians);\r\n      zScale *= -1.0;\r\n      yScale *= -1.0;\r\n      isReversed = true;\r\n    }\r\n    const frame1 = frame.clone();\r\n    frame1.matrix.scaleColumns(xScale, yScale, zScale, frame1.matrix);\r\n    const result = new TorusPipe(frame1, majorRadius, minorRadius, sweep1, capped);\r\n    result._isReversed = isReversed;\r\n    return result;\r\n  }\r\n\r\n  /** Create a TorusPipe from the typical parameters of the Dgn file */\r\n  public static createDgnTorusPipe(center: Point3d, vectorX: Vector3d, vectorY: Vector3d, majorRadius: number, minorRadius: number,\r\n    sweep: Angle, capped: boolean) {\r\n    const vectorZ = vectorX.crossProduct(vectorY);\r\n    vectorZ.scaleToLength(vectorX.magnitude(), vectorZ);\r\n    const frame = Transform.createOriginAndMatrixColumns(center, vectorX, vectorY, vectorZ);\r\n    return TorusPipe.createInFrame(frame, majorRadius, minorRadius, sweep, capped);\r\n  }\r\n  /** Create a TorusPipe from its primary arc and minor radius */\r\n  public static createAlongArc(arc: Arc3d, minorRadius: number, capped: boolean) {\r\n    if (!Angle.isAlmostEqualRadiansAllowPeriodShift(0.0, arc.sweep.startRadians))\r\n      arc = arc.cloneInRotatedBasis(arc.sweep.startAngle);\r\n    const sweepRadians = arc.sweep.sweepRadians;\r\n    const data = arc.toScaledMatrix3d();\r\n    const frame = Transform.createOriginAndMatrix(data.center, data.axes);\r\n    return TorusPipe.createInFrame(frame, data.r0, minorRadius, Angle.createRadians(sweepRadians), capped);\r\n  }\r\n\r\n  /** Return a coordinate frame (right handed, unit axes)\r\n   * * origin at center of major circle\r\n   * * major circle in xy plane\r\n   * * z axis perpendicular\r\n   */\r\n  public getConstructiveFrame(): Transform | undefined {\r\n    return this._localToWorld.cloneRigid();\r\n  }\r\n  /** Return the center of the torus pipe (inside the donut hole) */\r\n  public cloneCenter(): Point3d { return this._localToWorld.getOrigin(); }\r\n  /** return the vector along the x axis (in the major hoops plane) */\r\n  public cloneVectorX(): Vector3d { return this._localToWorld.matrix.columnX(); }\r\n  /** return the vector along the y axis (in the major hoop plane) */\r\n  public cloneVectorY(): Vector3d { return this._localToWorld.matrix.columnY(); }\r\n  /** get the minor hoop radius (`radiusA`) */\r\n  public getMinorRadius(): number { return this._radiusB; }\r\n  /** get the major hoop radius (`radiusB`) */\r\n  public getMajorRadius(): number { return this._radiusA; }\r\n  /** get the sweep angle along the major circle. */\r\n  public getSweepAngle(): Angle { return this._sweep.clone(); }\r\n  /** Ask if this TorusPipe is labeled as reversed */\r\n  public getIsReversed(): boolean { return this._isReversed; }\r\n  /** Return the sweep angle as a fraction of full 360 degrees (2PI radians) */\r\n  public getThetaFraction(): number { return this._sweep.radians / (Math.PI * 2.0); }\r\n  /** ask if `other` is an instance of `TorusPipe` */\r\n  public isSameGeometryClass(other: any): boolean { return other instanceof TorusPipe; }\r\n  /** test if `this` and `other` have nearly equal geometry */\r\n  public override isAlmostEqual(other: GeometryQuery): boolean {\r\n    if (other instanceof TorusPipe) {\r\n      if ((!this._sweep.isFullCircle) && this.capped !== other.capped) return false;\r\n      if (!this._localToWorld.isAlmostEqual(other._localToWorld)) return false;\r\n      return Geometry.isSameCoordinate(this._radiusA, other._radiusA)\r\n        && Geometry.isSameCoordinate(this._radiusB, other._radiusB)\r\n        && this._sweep.isAlmostEqualNoPeriodShift(other._sweep);\r\n    }\r\n    return false;\r\n  }\r\n  /** Return the angle (in radians) for given fractional position around the major hoop.\r\n   */\r\n  public vFractionToRadians(v: number): number { return this._sweep.radians * v; }\r\n  /** Second step of double dispatch:  call `handler.handleTorusPipe(this)` */\r\n  public dispatchToGeometryHandler(handler: GeometryHandler): any {\r\n    return handler.handleTorusPipe(this);\r\n  }\r\n\r\n  /**\r\n   * Return the Arc3d section at vFraction.  For the TorusPipe, this is a minor circle.\r\n   * @param vFraction fractional position along the sweep direction\r\n   */\r\n  public constantVSection(v: number): CurveCollection | undefined {\r\n    const thetaRadians = this.vFractionToRadians(v);\r\n    const c0 = Math.cos(thetaRadians);\r\n    const s0 = Math.sin(thetaRadians);\r\n    const majorRadius = this.getMajorRadius();\r\n    const minorRadius = this.getMinorRadius();\r\n\r\n    const center = this._localToWorld.multiplyXYZ(majorRadius * c0, majorRadius * s0, 0);\r\n    const vector0 = this._localToWorld.multiplyVectorXYZ(minorRadius * c0, minorRadius * s0, 0);\r\n    const vector90 = this._localToWorld.multiplyVectorXYZ(0, 0, minorRadius);\r\n    return Loop.create(Arc3d.create(center, vector0, vector90));\r\n  }\r\n  /** Return an arc at constant u, and arc sweep  matching this TorusPipe sweep. */\r\n  public constantUSection(uFraction: number): CurveCollection | undefined {\r\n    const theta1Radians = this._sweep.radians;\r\n    const phiRadians = uFraction * Math.PI;\r\n\r\n    const majorRadius = this.getMajorRadius();\r\n    const minorRadius = this.getMinorRadius();\r\n    const transform = this._localToWorld;\r\n    const axes = transform.matrix;\r\n    const center = this._localToWorld.multiplyXYZ(0, 0, minorRadius * Math.sin(phiRadians));\r\n    const rxy = majorRadius + minorRadius * Math.cos(phiRadians);\r\n    const vector0 = axes.multiplyXYZ(rxy, 0, 0);\r\n    const vector90 = axes.multiplyXYZ(0, rxy, 0);\r\n    return Path.create(Arc3d.create(center, vector0, vector90, AngleSweep.createStartEndRadians(0.0, theta1Radians)));\r\n  }\r\n  /** extend `rangeToExtend` to include this `TorusPipe` */\r\n  public extendRange(rangeToExtend: Range3d, transform?: Transform) {\r\n    const theta1Radians = this._sweep.radians;\r\n    const majorRadius = this.getMajorRadius();\r\n    const minorRadius = this.getMinorRadius();\r\n    const transform0 = this._localToWorld;\r\n    const numThetaSample = Math.ceil(theta1Radians / (Math.PI / 16.0));\r\n    const numHalfPhiSample = 16;\r\n    let phi0 = 0;\r\n    let dPhi = 0;\r\n    let numPhiSample = 0;\r\n    let theta = 0;\r\n    let cosTheta = 0;\r\n    let sinTheta = 0;\r\n    let rxy = 0;\r\n    let phi = 0;\r\n    let j = 0;\r\n    const dTheta = theta1Radians / numThetaSample;\r\n    for (let i = 0; i <= numThetaSample; i++) {\r\n      theta = i * dTheta;\r\n      cosTheta = Math.cos(theta);\r\n      sinTheta = Math.sin(theta);\r\n      // At the ends, do the entire phi circle.\r\n      // Otherwise only do the outer half\r\n      if (i === 0 || i === numThetaSample) {\r\n        phi0 = -Math.PI;\r\n        dPhi = 2.0 * Math.PI / numHalfPhiSample;\r\n        numPhiSample = numHalfPhiSample;\r\n      } else {\r\n        phi0 = -0.5 * Math.PI;\r\n        dPhi = Math.PI / numHalfPhiSample;\r\n        numPhiSample = 2 * numHalfPhiSample - 1;\r\n      }\r\n      if (transform) {\r\n        for (j = 0; j <= numPhiSample; j++) {\r\n          phi = phi0 + j * dPhi;\r\n          rxy = majorRadius + minorRadius * Math.cos(phi);\r\n          rangeToExtend.extendTransformTransformedXYZ(transform, transform0,\r\n            cosTheta * rxy, sinTheta * rxy,\r\n            Math.sin(phi) * minorRadius);\r\n        }\r\n      } else {\r\n        for (j = 0; j <= numPhiSample; j++) {\r\n          phi = phi0 + j * dPhi;\r\n          rxy = majorRadius + minorRadius * Math.sin(phi);\r\n          rangeToExtend.extendTransformedXYZ(transform0,\r\n            cosTheta * rxy, sinTheta * rxy,\r\n            Math.sin(phi) * minorRadius);\r\n        }\r\n      }\r\n    }\r\n  }\r\n  /** Evaluate as a uv surface\r\n   * @param u fractional position in minor (phi)\r\n   * @param v fractional position on major (theta) arc\r\n   */\r\n  public uvFractionToPoint(u: number, v: number, result?: Point3d): Point3d {\r\n    const thetaRadians = v * this._sweep.radians;\r\n    const phiRadians = u * Math.PI * 2.0;\r\n    const cosTheta = Math.cos(thetaRadians);\r\n    const sinTheta = Math.sin(thetaRadians);\r\n    const minorRadius = this.getMinorRadius();\r\n    const rxy = this.getMajorRadius() + Math.cos(phiRadians) * minorRadius;\r\n    return this._localToWorld.multiplyXYZ(rxy * cosTheta, rxy * sinTheta, minorRadius * Math.sin(phiRadians), result);\r\n  }\r\n  /** Evaluate as a uv surface, returning point and two vectors.\r\n   * @param u fractional position in minor (phi)\r\n   * @param v fractional position on major (theta) arc\r\n   */\r\n  public uvFractionToPointAndTangents(u: number, v: number, result?: Plane3dByOriginAndVectors): Plane3dByOriginAndVectors {\r\n    const thetaRadians = v * this._sweep.radians;\r\n    const phiRadians = u * Math.PI * 2.0;\r\n    const fTheta = this._sweep.radians;\r\n    const fPhi = Math.PI * 2.0;\r\n    const cosTheta = Math.cos(thetaRadians);\r\n    const sinTheta = Math.sin(thetaRadians);\r\n    const sinPhi = Math.sin(phiRadians);\r\n    const cosPhi = Math.cos(phiRadians);\r\n    const minorRadius = this.getMinorRadius();\r\n    const rxy = this.getMajorRadius() + Math.cos(phiRadians) * minorRadius;\r\n    const rSinPhi = minorRadius * sinPhi;\r\n    const rCosPhi = minorRadius * cosPhi;   // appears only as derivative of rSinPhi.\r\n    return Plane3dByOriginAndVectors.createOriginAndVectors(\r\n      this._localToWorld.multiplyXYZ(cosTheta * rxy, sinTheta * rxy, rSinPhi),\r\n      this._localToWorld.multiplyVectorXYZ(-cosTheta * rSinPhi * fPhi, -sinTheta * rSinPhi * fPhi, rCosPhi * fPhi),\r\n      this._localToWorld.multiplyVectorXYZ(-rxy * sinTheta * fTheta, rxy * cosTheta * fTheta, 0),\r\n      result);\r\n  }\r\n  /**\r\n   * Directional distance query\r\n   * * u direction is around the (full) minor hoop\r\n   * * v direction is around the outer radius, sum of (absolute values of) major and minor radii.\r\n   */\r\n  public maxIsoParametricDistance(): Vector2d {\r\n    const a = Math.abs(this.getMajorRadius());\r\n    const b = Math.abs(this.getMinorRadius());\r\n    return Vector2d.create(b * Math.PI * 2.0, (a + b) * this._sweep.radians);\r\n  }\r\n  /**\r\n   * @return true if this is a closed volume.\r\n   */\r\n  public get isClosedVolume(): boolean {\r\n    return this.capped || this._sweep.isFullCircle;\r\n  }\r\n\r\n}\r\n"]}

package/lib/esm/topology/ChainMerge.d.ts

@@ -0,0 +1,119 @@
+/** @packageDocumentation
+ * @module Topology
+ */
+import { ConvexClipPlaneSet } from "../clipping/ConvexClipPlaneSet";
+import { LineSegment3d } from "../curve/LineSegment3d";
+import { LineString3d } from "../curve/LineString3d";
+import { PlaneAltitudeEvaluator } from "../Geometry";
+import { GrowableXYZArray } from "../geometry3d/GrowableXYZArray";
+import { Point3d, Vector3d } from "../geometry3d/Point3dVector3d";
+/**
+ * interface containing various options appropriate to merging lines segments into chains.
+ * @internal
+ */
+export interface ChainMergeContextOptions {
+    /**
+     * Tolerance for declaring points equal
+     * * Default is `Geometry.smallMetricDistance`
+     */
+    tolerance?: number;
+    /**
+     * Direction for primary sort.
+     * * Default is `Vector3d.create(0.294234298, 0.72391399, 0.45234328798)`
+     * * this vector should NOT be along a principal x,y,z direction.
+     * * The internal form will be normalized.
+     */
+    primarySortDirection?: Vector3d;
+}
+/**
+ * * Context for assembling line segments into chains.
+ * * Use the context in phases:
+ *   * Create the context:   `context = ChainMergeContext.create ()`
+ *   * Add line with any combination of:
+ *      * `context.addSegment(pointA, pointB)`
+ *      * `context.addLineSegment3dArray (segments)`
+ *   * Scan all coordinate data for common coordinates.  Twist nodes together to form chains:
+ *      * `context.clusterAndMergeVerticesXYZ ()`
+ *   * Collect the chains:
+ *      * myLinestringArray = context.collectMaximalChains();
+ * * The context carries an optional plane which is used by addSegmentsOnPlane
+ *
+ * @internal
+ */
+export declare class ChainMergeContext {
+    private _graph;
+    private _options;
+    private _plane?;
+    private _convexClipper?;
+    private constructor();
+    /** Save a reference plane for later use, e.g. in addSegmentsOnPlane */
+    set plane(value: PlaneAltitudeEvaluator | undefined);
+    /** Property access for the reference plane. */
+    get plane(): PlaneAltitudeEvaluator | undefined;
+    /** Save a reference plane for later use, e.g. in addSegmentsOnPlane */
+    set convexClipper(value: ConvexClipPlaneSet | undefined);
+    /** Property access for the reference plane. */
+    get convexClipper(): ConvexClipPlaneSet | undefined;
+    /**
+     * * Construct an empty chain merge graph.
+     * * The options parameter may contain any combination of the options values.
+     *   * tolerance = absolute tolerance for declaring points equal.
+     *     * Default is `Geometry.smallMetricDistance`
+     *   * primarySortDirection = direction for first sort.
+     *     * To minimize clash among points on primary sort, this should NOT be perpendicular to any principal plane.
+     *     * The default points into the first octant with non-obvious components.
+     */
+    static create(options?: ChainMergeContextOptions): ChainMergeContext;
+    /** Add a segment to the evolving graph. */
+    addSegment(pointA: Point3d, pointB: Point3d): void;
+    /** Add all segments from an array to the graph. */
+    addLineSegment3dArray(data: LineSegment3d[]): void;
+    /** Add edges for all segments that are "on" the plane.
+     * * No action if `this.plane` is undefined.
+     */
+    addSegmentsOnPlane(points: GrowableXYZArray, addClosure?: boolean): void;
+    /**
+     * Return a numeric value to be used for sorting, with little chance widely separated nodes will have identical key.
+     * * Any single x,y,z component is a poor choice because horizontal and vertical alignment is common.
+     * * Hence take dot product of x,y,z with non-trivial fraction values.
+     * @param node node with x,y,z coordinates
+     */
+    private primarySortKey;
+    /** Return difference of sortData members as sort comparison */
+    private static nodeCompareSortData;
+    /** test if nodeA is a dangling edge end (i.e. edges around vertex equal 1, but detect it without walking all the way around. */
+    private static isIsolatedEnd;
+    /** test if nodeA is at a vertex with exactly 2 edges (i.e. edges around vertex equal w, but detect it without walking all the way around. */
+    private static isChainInteriorVertex;
+    /**
+     * * isolate all edge ends.
+     * * perform cluster analysis to determine nearly coincident points.
+     * * pinch all edges together at each vertex.
+     */
+    clusterAndMergeVerticesXYZ(): void;
+    /**
+     * If node0 is not visited, creating a linestring with that starting edge and all successive edges along a chain.
+     * @param chains growing array of chains.
+     * @param node0 start node for search.
+     */
+    private collectMaximalLineString3dFromStartNode;
+    /**
+     * If node0 is not visited, creating a linestring with that starting edge and all successive edges along a chain.
+     * @param chains growing array of chains.
+     * @param node0 start node for search.
+     */
+    private collectMaximalGrowableXYXArrayFromStartNode;
+    /**
+     * * find edges with start, end in same vertex loop.
+     * * pinch them away from the loop
+     * * set mask on both sides.
+     * * Return the number of excisions.
+     */
+    private exciseAndMarkSlingEdges;
+    /** Collect chains which have maximum edge count, broken at an vertex with other than 2 edges.
+     * * This is assumed to be preceded by a call to a vertex-cluster step such as `clusterAndMergeVerticesYXZ`
+     */
+    collectMaximalChains(): LineString3d[];
+    collectMaximalGrowableXYZArrays(): GrowableXYZArray[];
+}
+//# sourceMappingURL=ChainMerge.d.ts.map

package/lib/esm/topology/ChainMerge.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"ChainMerge.d.ts","sourceRoot":"","sources":["../../../src/topology/ChainMerge.ts"],"names":[],"mappings":"AAKA;;GAEG;AAEH,OAAO,EAAE,kBAAkB,EAAE,MAAM,gCAAgC,CAAC;AACpE,OAAO,EAAE,aAAa,EAAE,MAAM,wBAAwB,CAAC;AACvD,OAAO,EAAE,YAAY,EAAE,MAAM,uBAAuB,CAAC;AACrD,OAAO,EAAY,sBAAsB,EAAE,MAAM,aAAa,CAAC;AAC/D,OAAO,EAAE,gBAAgB,EAAE,MAAM,gCAAgC,CAAC;AAClE,OAAO,EAAE,OAAO,EAAE,QAAQ,EAAE,MAAM,+BAA+B,CAAC;AAIlE;;;GAGG;AACH,MAAM,WAAW,wBAAwB;IACvC;;;OAGG;IACH,SAAS,CAAC,EAAE,MAAM,CAAC;IACnB;;;;;OAKG;IACH,oBAAoB,CAAC,EAAE,QAAQ,CAAC;CACjC;AAoDD;;;;;;;;;;;;;;GAcG;AACH,qBAAa,iBAAiB;IAE5B,OAAO,CAAC,MAAM,CAAgB;IAC9B,OAAO,CAAC,QAAQ,CAAoC;IACpD,OAAO,CAAC,MAAM,CAAC,CAAyB;IACxC,OAAO,CAAC,cAAc,CAAC,CAAqB;IAC5C,OAAO;IAIP,uEAAuE;IACvE,IAAW,KAAK,CAAC,KAAK,EAAE,sBAAsB,GAAG,SAAS,EAA0B;IACpF,+CAA+C;IAC/C,IAAW,KAAK,IAAI,sBAAsB,GAAG,SAAS,CAAwB;IAE9E,uEAAuE;IACvE,IAAW,aAAa,CAAC,KAAK,EAAE,kBAAkB,GAAG,SAAS,EAAkC;IAChG,+CAA+C;IAC/C,IAAW,aAAa,IAAI,kBAAkB,GAAG,SAAS,CAAgC;IAE1F;;;;;;;;OAQG;WACW,MAAM,CAAC,OAAO,CAAC,EAAE,wBAAwB,GAAG,iBAAiB;IAI3E,2CAA2C;IACpC,UAAU,CAAC,MAAM,EAAE,OAAO,EAAE,MAAM,EAAE,OAAO;IAGlD,mDAAmD;IAC5C,qBAAqB,CAAC,IAAI,EAAE,aAAa,EAAE;IAKlD;;OAEG;IACI,kBAAkB,CAAC,MAAM,EAAE,gBAAgB,EAAE,UAAU,GAAE,OAAe;IAgB/E;;;;;OAKG;IACH,OAAO,CAAC,cAAc;IAGtB,+DAA+D;IAC/D,OAAO,CAAC,MAAM,CAAC,mBAAmB;IAGlC,gIAAgI;IAChI,OAAO,CAAC,MAAM,CAAC,aAAa;IAI5B,6IAA6I;IAC7I,OAAO,CAAC,MAAM,CAAC,qBAAqB;IAKpC;;;;OAIG;IACI,0BAA0B;IA8BjC;;;;OAIG;IACH,OAAO,CAAC,uCAAuC;IAe/C;;;;OAIG;IACH,OAAO,CAAC,2CAA2C;IAiBnD;;;;;OAKG;IACH,OAAO,CAAC,uBAAuB;IAe/B;;OAEG;IACI,oBAAoB,IAAI,YAAY,EAAE;IAmBtC,+BAA+B,IAAI,gBAAgB,EAAE;CAmB7D"}