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

package/lib/esm/solid/Sphere.js

@@ -0,0 +1,294 @@
+/*---------------------------------------------------------------------------------------------
+* 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 { LineString3d } from "../curve/LineString3d";
+import { Loop } from "../curve/Loop";
+import { StrokeOptions } from "../curve/StrokeOptions";
+import { Geometry } from "../Geometry";
+import { AngleSweep } from "../geometry3d/AngleSweep";
+import { Matrix3d } from "../geometry3d/Matrix3d";
+import { Plane3dByOriginAndVectors } from "../geometry3d/Plane3dByOriginAndVectors";
+import { Vector2d } from "../geometry3d/Point2dVector2d";
+import { Vector3d } from "../geometry3d/Point3dVector3d";
+import { Transform } from "../geometry3d/Transform";
+import { SolidPrimitive } from "./SolidPrimitive";
+/**
+ * A Sphere is
+ *
+ * * A unit sphere (but read on ....)
+ * * mapped by an arbitrary (possibly skewed, non-uniform scaled) transform
+ * * hence possibly the final geometry is ellipsoidal
+ * @public
+ */
+export class Sphere extends SolidPrimitive {
+    constructor(localToWorld, latitudeSweep, capped) {
+        super(capped);
+        /** String name for schema properties */
+        this.solidPrimitiveType = "sphere";
+        this._localToWorld = localToWorld;
+        this._latitudeSweep = latitudeSweep ? latitudeSweep : AngleSweep.createFullLatitude();
+    }
+    /** Return the latitude (in radians) all fractional v. */
+    vFractionToRadians(v) {
+        return this._latitudeSweep.fractionToRadians(v);
+    }
+    /** Return the longitude (in radians) all fractional u. */
+    uFractionToRadians(u) {
+        return u * Math.PI * 2.0;
+    }
+    /** return a deep clone */
+    clone() {
+        return new Sphere(this._localToWorld.clone(), this._latitudeSweep.clone(), this.capped);
+    }
+    /** Transform the sphere in place.
+     * * Fails if the transform is singular.
+     */
+    tryTransformInPlace(transform) {
+        if (transform.matrix.isSingular())
+            return false;
+        transform.multiplyTransformTransform(this._localToWorld, this._localToWorld);
+        return true;
+    }
+    /** Return a transformed clone. */
+    cloneTransformed(transform) {
+        const sphere1 = this.clone();
+        transform.multiplyTransformTransform(sphere1._localToWorld, sphere1._localToWorld);
+        if (transform.matrix.determinant() < 0.0) {
+            if (sphere1._latitudeSweep !== undefined) {
+                sphere1._latitudeSweep.reverseInPlace();
+            }
+        }
+        return sphere1;
+    }
+    /** Return a coordinate frame (right handed, unit axes)
+     * * origin at sphere center
+     * * equator in xy plane
+     * * z axis perpendicular
+     */
+    getConstructiveFrame() {
+        return this._localToWorld.cloneRigid();
+    }
+    /** Return the latitude sweep as fraction of south pole to north pole. */
+    get latitudeSweepFraction() { return this._latitudeSweep.sweepRadians / Math.PI; }
+    /** Create from center and radius, with optional restricted latitudes. */
+    static createCenterRadius(center, radius, latitudeSweep) {
+        const localToWorld = Transform.createOriginAndMatrix(center, Matrix3d.createUniformScale(radius));
+        return new Sphere(localToWorld, latitudeSweep ? latitudeSweep : AngleSweep.createFullLatitude(), false);
+    }
+    /** Create an ellipsoid which is a unit sphere mapped to position by an (arbitrary, possibly skewed and scaled) transform. */
+    static createEllipsoid(localToWorld, latitudeSweep, capped) {
+        return new Sphere(localToWorld, latitudeSweep, capped);
+    }
+    /** Create a sphere from the typical parameters of the Dgn file */
+    static createDgnSphere(center, vectorX, vectorZ, radiusXY, radiusZ, latitudeSweep, capped) {
+        const vectorY = vectorX.rotate90Around(vectorZ);
+        if (vectorY && !vectorX.isParallelTo(vectorZ)) {
+            const matrix = Matrix3d.createColumns(vectorX, vectorY, vectorZ);
+            matrix.scaleColumns(radiusXY, radiusXY, radiusZ, matrix);
+            const frame = Transform.createOriginAndMatrix(center, matrix);
+            return new Sphere(frame, latitudeSweep.clone(), capped);
+        }
+        return undefined;
+    }
+    /** Create a sphere from the typical parameters of the Dgn file */
+    static createFromAxesAndScales(center, axes, radiusX, radiusY, radiusZ, latitudeSweep, capped) {
+        const localToWorld = Transform.createOriginAndMatrix(center, axes);
+        localToWorld.matrix.scaleColumnsInPlace(radiusX, radiusY, radiusZ);
+        return new Sphere(localToWorld, latitudeSweep ? latitudeSweep.clone() : AngleSweep.createFullLatitude(), capped);
+    }
+    /** return (copy of) sphere center */
+    cloneCenter() { return this._localToWorld.getOrigin(); }
+    /** return the (full length, i.e. scaled by radius) X vector from the sphere transform */
+    cloneVectorX() { return this._localToWorld.matrix.columnX(); }
+    /** return the (full length, i.e. scaled by radius) Y vector from the sphere transform */
+    cloneVectorY() { return this._localToWorld.matrix.columnY(); }
+    /** return the (full length, i.e. scaled by radius) Z vector from the sphere transform */
+    cloneVectorZ() { return this._localToWorld.matrix.columnZ(); }
+    /** return (a copy of) the sphere's angle sweep. */
+    cloneLatitudeSweep() { return this._latitudeSweep.clone(); }
+    /** Test if the geometry is a true sphere taking the transform (which might have nonuniform scaling) is applied. */
+    trueSphereRadius() {
+        const factors = this._localToWorld.matrix.factorRigidWithSignedScale();
+        if (!factors)
+            return undefined;
+        if (factors && factors.scale > 0)
+            return factors.scale;
+        return undefined;
+    }
+    /**
+     * Return the larger of the primary xyz axis radii
+     */
+    maxAxisRadius() {
+        const matrix = this._localToWorld.matrix;
+        return Geometry.maxXYZ(matrix.columnXMagnitude(), matrix.columnYMagnitude(), matrix.columnYMagnitude());
+    }
+    /**
+     * Return a (clone of) the sphere's local to world transformation.
+     */
+    cloneLocalToWorld() { return this._localToWorld.clone(); }
+    /** Test if `other` is a `Sphere` */
+    isSameGeometryClass(other) { return other instanceof Sphere; }
+    /** Test for same geometry in `other` */
+    isAlmostEqual(other) {
+        if (other instanceof Sphere) {
+            if (this.capped !== other.capped)
+                return false;
+            if (!this._localToWorld.isAlmostEqual(other._localToWorld))
+                return false;
+            return true;
+        }
+        return false;
+    }
+    /**
+     *  return strokes for a cross-section (elliptic arc) at specified fraction v along the axis.
+     * * if strokeOptions is supplied, it is applied to the equator radii.
+     * @param v fractional position along the cone axis
+     * @param strokes stroke count or options.
+     */
+    strokeConstantVSection(v, fixedStrokeCount, options) {
+        let strokeCount = 16;
+        if (fixedStrokeCount !== undefined && Number.isFinite(fixedStrokeCount)) {
+            strokeCount = fixedStrokeCount;
+        }
+        else if (options instanceof StrokeOptions) {
+            strokeCount = options.applyTolerancesToArc(Geometry.maxXY(this._localToWorld.matrix.columnXMagnitude(), this._localToWorld.matrix.columnYMagnitude()));
+        }
+        strokeCount = Geometry.clampToStartEnd(strokeCount, 4, 64);
+        const transform = this._localToWorld;
+        const phi = this.vFractionToRadians(v);
+        const c1 = Math.cos(phi);
+        const s1 = Math.sin(phi);
+        let c0, s0;
+        const result = LineString3d.createForStrokes(fixedStrokeCount, options);
+        const deltaRadians = Math.PI * 2.0 / strokeCount;
+        const fractions = result.fractions; // possibly undefined !!!
+        const derivatives = result.packedDerivatives; // possibly undefined !!!
+        const uvParams = result.packedUVParams; // possibly undefined !!
+        const surfaceNormals = result.packedSurfaceNormals;
+        const dXdu = Vector3d.create();
+        const dXdv = Vector3d.create();
+        const normal = Vector3d.create();
+        let radians = 0;
+        for (let i = 0; i <= strokeCount; i++) {
+            if (i * 2 <= strokeCount)
+                radians = i * deltaRadians;
+            else
+                radians = (i - strokeCount) * deltaRadians;
+            c0 = Math.cos(radians);
+            s0 = Math.sin(radians);
+            const xyz = transform.multiplyXYZ(c1 * c0, c1 * s0, s1);
+            result.addPoint(xyz);
+            if (fractions)
+                fractions.push(i / strokeCount);
+            if (derivatives) {
+                transform.matrix.multiplyXYZ(-c1 * s0, c1 * c0, 0.0, dXdu);
+                derivatives.push(dXdu);
+            }
+            if (uvParams) {
+                uvParams.pushXY(i / strokeCount, v);
+            }
+            if (surfaceNormals) {
+                transform.matrix.multiplyXYZ(-s0, c0, 0, dXdu);
+                transform.matrix.multiplyXYZ(-s1 * c0, -s1 * s0, c1, dXdv);
+                dXdu.unitCrossProduct(dXdv, normal);
+                surfaceNormals.push(normal);
+            }
+        }
+        return result;
+    }
+    /** Second step of double dispatch:  call `handler.handleSphere(this)` */
+    dispatchToGeometryHandler(handler) {
+        return handler.handleSphere(this);
+    }
+    /**
+     * Return the Arc3d section at vFraction.  For the sphere, this is a latitude circle.
+     * @param vFraction fractional position along the sweep direction
+     */
+    constantVSection(vFraction) {
+        const phi = this._latitudeSweep.fractionToRadians(vFraction);
+        const s1 = Math.sin(phi);
+        const c1 = Math.cos(phi);
+        const transform = this._localToWorld;
+        const center = transform.multiplyXYZ(0, 0, s1);
+        const vector0 = transform.matrix.multiplyXYZ(c1, 0, 0);
+        const vector90 = transform.matrix.multiplyXYZ(0, c1, 0);
+        return Loop.create(Arc3d.create(center, vector0, vector90));
+    }
+    /** Extend a range to contain this sphere. */
+    extendRange(range, transform) {
+        let placement = this._localToWorld;
+        if (transform) {
+            placement = transform.multiplyTransformTransform(placement);
+        }
+        range.extendTransformedXYZ(placement, -1, -1, -1);
+        range.extendTransformedXYZ(placement, 1, -1, -1);
+        range.extendTransformedXYZ(placement, -1, 1, -1);
+        range.extendTransformedXYZ(placement, 1, 1, -1);
+        range.extendTransformedXYZ(placement, -1, -1, 1);
+        range.extendTransformedXYZ(placement, 1, -1, 1);
+        range.extendTransformedXYZ(placement, -1, 1, 1);
+        range.extendTransformedXYZ(placement, 1, 1, 1);
+    }
+    /** Evaluate as a uv surface
+     * @param uFraction fractional position in minor (phi)
+     * @param vFraction fractional position on major (theta) arc
+     */
+    uvFractionToPoint(uFraction, vFraction, result) {
+        // sphere with radius 1 . . .
+        const thetaRadians = this.uFractionToRadians(uFraction);
+        const phiRadians = this.vFractionToRadians(vFraction);
+        const cosTheta = Math.cos(thetaRadians);
+        const sinTheta = Math.sin(thetaRadians);
+        const sinPhi = Math.sin(phiRadians);
+        const cosPhi = Math.cos(phiRadians);
+        return this._localToWorld.multiplyXYZ(cosTheta * cosPhi, sinTheta * cosPhi, sinPhi, 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(uFraction, vFraction, result) {
+        const thetaRadians = this.uFractionToRadians(uFraction);
+        const phiRadians = this.vFractionToRadians(vFraction);
+        const fTheta = Math.PI * 2.0;
+        const fPhi = this._latitudeSweep.sweepRadians;
+        const cosTheta = Math.cos(thetaRadians);
+        const sinTheta = Math.sin(thetaRadians);
+        const sinPhi = Math.sin(phiRadians);
+        const cosPhi = Math.cos(phiRadians);
+        return Plane3dByOriginAndVectors.createOriginAndVectors(this._localToWorld.multiplyXYZ(cosTheta * cosPhi, sinTheta * cosPhi, sinPhi), this._localToWorld.matrix.multiplyXYZ(-fTheta * sinTheta, fTheta * cosTheta, 0), // !!! note cosTheta term is omitted -- scale is wrong, but remains non-zero at poles.
+        this._localToWorld.matrix.multiplyXYZ(-fPhi * cosTheta * sinPhi, -fPhi * sinTheta * sinPhi, fPhi * cosPhi), result);
+    }
+    /**
+     * * A sphere is can be closed two ways:
+     *   * full sphere (no caps needed for closure)
+     *   * incomplete but with caps
+     * @return true if this is a closed volume.
+     */
+    get isClosedVolume() {
+        return this.capped || this._latitudeSweep.isFullLatitudeSweep;
+    }
+    /**
+     * Directional distance query
+     * * u direction is around longitude circle at maximum distance from axis.
+     * * v direction is on a line of longitude between the latitude limits.
+     */
+    maxIsoParametricDistance() {
+        // approximate radius at equator .. if elliptic, this is not exact . . .
+        const rX = this._localToWorld.matrix.columnXMagnitude();
+        const rY = this._localToWorld.matrix.columnYMagnitude();
+        const rZ = this._localToWorld.matrix.columnZMagnitude();
+        const rMaxU = Math.max(rX, rY);
+        let dMaxU = Math.PI * 2.0 * rMaxU;
+        if (!this._latitudeSweep.isRadiansInSweep(0.0))
+            dMaxU *= Math.max(Math.cos(Math.abs(this._latitudeSweep.startRadians)), Math.cos(Math.abs(this._latitudeSweep.endRadians)));
+        const dMaxV = Math.max(rMaxU, rZ) * Math.abs(this._latitudeSweep.sweepRadians);
+        return Vector2d.create(dMaxU, dMaxV);
+    }
+}
+//# sourceMappingURL=Sphere.js.map

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

@@ -0,0 +1 @@
+{"version":3,"file":"Sphere.js","sourceRoot":"","sources":["../../../src/solid/Sphere.ts"],"names":[],"mappings":"AAAA;;;+FAG+F;AAE/F;;GAEG;AAEH,OAAO,EAAE,KAAK,EAAE,MAAM,gBAAgB,CAAC;AAGvC,OAAO,EAAE,YAAY,EAAE,MAAM,uBAAuB,CAAC;AACrD,OAAO,EAAE,IAAI,EAAE,MAAM,eAAe,CAAC;AACrC,OAAO,EAAE,aAAa,EAAE,MAAM,wBAAwB,CAAC;AACvD,OAAO,EAAE,QAAQ,EAAE,MAAM,aAAa,CAAC;AACvC,OAAO,EAAE,UAAU,EAAE,MAAM,0BAA0B,CAAC;AAEtD,OAAO,EAAE,QAAQ,EAAE,MAAM,wBAAwB,CAAC;AAClD,OAAO,EAAE,yBAAyB,EAAE,MAAM,yCAAyC,CAAC;AACpF,OAAO,EAAE,QAAQ,EAAE,MAAM,+BAA+B,CAAC;AACzD,OAAO,EAAW,QAAQ,EAAE,MAAM,+BAA+B,CAAC;AAElE,OAAO,EAAE,SAAS,EAAE,MAAM,yBAAyB,CAAC;AACpD,OAAO,EAAE,cAAc,EAAE,MAAM,kBAAkB,CAAC;AAElD;;;;;;;GAOG;AACH,MAAM,OAAO,MAAO,SAAQ,cAAc;IAexC,YAAoB,YAAuB,EAAE,aAAyB,EAAE,MAAe;QACrF,KAAK,CAAC,MAAM,CAAC,CAAC;QAfhB,wCAAwC;QACxB,uBAAkB,GAAG,QAAQ,CAAC;QAe5C,IAAI,CAAC,aAAa,GAAG,YAAY,CAAC;QAClC,IAAI,CAAC,cAAc,GAAG,aAAa,CAAC,CAAC,CAAC,aAAa,CAAC,CAAC,CAAC,UAAU,CAAC,kBAAkB,EAAE,CAAC;IACxF,CAAC;IAbD,yDAAyD;IAClD,kBAAkB,CAAC,CAAS;QACjC,OAAO,IAAI,CAAC,cAAc,CAAC,iBAAiB,CAAC,CAAC,CAAC,CAAC;IAClD,CAAC;IACD,0DAA0D;IACnD,kBAAkB,CAAC,CAAS;QACjC,OAAO,CAAC,GAAG,IAAI,CAAC,EAAE,GAAG,GAAG,CAAC;IAC3B,CAAC;IAOD,0BAA0B;IACnB,KAAK;QACV,OAAO,IAAI,MAAM,CAAC,IAAI,CAAC,aAAa,CAAC,KAAK,EAAE,EAAE,IAAI,CAAC,cAAc,CAAC,KAAK,EAAE,EAAE,IAAI,CAAC,MAAM,CAAC,CAAC;IAC1F,CAAC;IACD;;OAEG;IACI,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,kCAAkC;IAC3B,gBAAgB,CAAC,SAAoB;QAC1C,MAAM,OAAO,GAAG,IAAI,CAAC,KAAK,EAAE,CAAC;QAC7B,SAAS,CAAC,0BAA0B,CAAC,OAAO,CAAC,aAAa,EAAE,OAAO,CAAC,aAAa,CAAC,CAAC;QACnF,IAAI,SAAS,CAAC,MAAM,CAAC,WAAW,EAAE,GAAG,GAAG,EAAE;YACxC,IAAI,OAAO,CAAC,cAAc,KAAK,SAAS,EAAE;gBACxC,OAAO,CAAC,cAAc,CAAC,cAAc,EAAE,CAAC;aACzC;SACF;QACD,OAAO,OAAO,CAAC;IACjB,CAAC;IACD;;;;OAIG;IACI,oBAAoB;QACzB,OAAO,IAAI,CAAC,aAAa,CAAC,UAAU,EAAE,CAAC;IACzC,CAAC;IACD,yEAAyE;IACzE,IAAW,qBAAqB,KAAa,OAAO,IAAI,CAAC,cAAc,CAAC,YAAY,GAAG,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC;IACjG,yEAAyE;IAClE,MAAM,CAAC,kBAAkB,CAAC,MAAe,EAAE,MAAc,EAAE,aAA0B;QAC1F,MAAM,YAAY,GAAG,SAAS,CAAC,qBAAqB,CAAC,MAAM,EAAE,QAAQ,CAAC,kBAAkB,CAAC,MAAM,CAAC,CAAC,CAAC;QAClG,OAAO,IAAI,MAAM,CAAC,YAAY,EAC5B,aAAa,CAAC,CAAC,CAAC,aAAa,CAAC,CAAC,CAAC,UAAU,CAAC,kBAAkB,EAAE,EAAE,KAAK,CAAC,CAAC;IAC5E,CAAC;IACD,6HAA6H;IACtH,MAAM,CAAC,eAAe,CAAC,YAAuB,EAAE,aAAyB,EAAE,MAAe;QAC/F,OAAO,IAAI,MAAM,CAAC,YAAY,EAAE,aAAa,EAAE,MAAM,CAAC,CAAC;IACzD,CAAC;IAED,kEAAkE;IAC3D,MAAM,CAAC,eAAe,CAAC,MAAe,EAAE,OAAiB,EAAE,OAAiB,EAAE,QAAgB,EAAE,OAAe,EACpH,aAAyB,EACzB,MAAe;QACf,MAAM,OAAO,GAAG,OAAO,CAAC,cAAc,CAAC,OAAO,CAAC,CAAC;QAChD,IAAI,OAAO,IAAI,CAAC,OAAO,CAAC,YAAY,CAAC,OAAO,CAAC,EAAE;YAC7C,MAAM,MAAM,GAAG,QAAQ,CAAC,aAAa,CAAC,OAAO,EAAE,OAAO,EAAE,OAAO,CAAC,CAAC;YACjE,MAAM,CAAC,YAAY,CAAC,QAAQ,EAAE,QAAQ,EAAE,OAAO,EAAE,MAAM,CAAC,CAAC;YACzD,MAAM,KAAK,GAAG,SAAS,CAAC,qBAAqB,CAAC,MAAM,EAAE,MAAM,CAAC,CAAC;YAC9D,OAAO,IAAI,MAAM,CAAC,KAAK,EAAE,aAAa,CAAC,KAAK,EAAE,EAAE,MAAM,CAAC,CAAC;SACzD;QACD,OAAO,SAAS,CAAC;IACnB,CAAC;IAED,kEAAkE;IAC3D,MAAM,CAAC,uBAAuB,CAAC,MAAe,EAAE,IAA0B,EAAE,OAAe,EAAE,OAAe,EAAE,OAAe,EAClI,aAAqC,EACrC,MAAe;QACf,MAAM,YAAY,GAAG,SAAS,CAAC,qBAAqB,CAAC,MAAM,EAAE,IAAI,CAAC,CAAC;QACnE,YAAY,CAAC,MAAM,CAAC,mBAAmB,CAAC,OAAO,EAAE,OAAO,EAAE,OAAO,CAAC,CAAC;QACnE,OAAO,IAAI,MAAM,CAAC,YAAY,EAAE,aAAa,CAAC,CAAC,CAAC,aAAa,CAAC,KAAK,EAAE,CAAC,CAAC,CAAC,UAAU,CAAC,kBAAkB,EAAE,EAAE,MAAM,CAAC,CAAC;IACnH,CAAC;IAED,qCAAqC;IAC9B,WAAW,KAAc,OAAO,IAAI,CAAC,aAAa,CAAC,SAAS,EAAE,CAAC,CAAC,CAAC;IACxE,yFAAyF;IAClF,YAAY,KAAe,OAAO,IAAI,CAAC,aAAa,CAAC,MAAM,CAAC,OAAO,EAAE,CAAC,CAAC,CAAC;IAC/E,yFAAyF;IAClF,YAAY,KAAe,OAAO,IAAI,CAAC,aAAa,CAAC,MAAM,CAAC,OAAO,EAAE,CAAC,CAAC,CAAC;IAC/E,yFAAyF;IAClF,YAAY,KAAe,OAAO,IAAI,CAAC,aAAa,CAAC,MAAM,CAAC,OAAO,EAAE,CAAC,CAAC,CAAC;IAC/E,mDAAmD;IAC5C,kBAAkB,KAAiB,OAAO,IAAI,CAAC,cAAc,CAAC,KAAK,EAAE,CAAC,CAAC,CAAC;IAC/E,mHAAmH;IAC5G,gBAAgB;QACrB,MAAM,OAAO,GAAG,IAAI,CAAC,aAAa,CAAC,MAAM,CAAC,0BAA0B,EAAE,CAAC;QACvE,IAAI,CAAC,OAAO;YAAE,OAAO,SAAS,CAAC;QAC/B,IAAI,OAAO,IAAI,OAAO,CAAC,KAAK,GAAG,CAAC;YAC9B,OAAO,OAAO,CAAC,KAAK,CAAC;QACvB,OAAO,SAAS,CAAC;IACnB,CAAC;IACD;;OAEG;IACI,aAAa;QAClB,MAAM,MAAM,GAAG,IAAI,CAAC,aAAa,CAAC,MAAM,CAAC;QACzC,OAAO,QAAQ,CAAC,MAAM,CAAC,MAAM,CAAC,gBAAgB,EAAE,EAAE,MAAM,CAAC,gBAAgB,EAAE,EAAE,MAAM,CAAC,gBAAgB,EAAE,CAAC,CAAC;IAC1G,CAAC;IACD;;OAEG;IACI,iBAAiB,KAAgB,OAAO,IAAI,CAAC,aAAa,CAAC,KAAK,EAAE,CAAC,CAAC,CAAC;IAC5E,oCAAoC;IAC7B,mBAAmB,CAAC,KAAU,IAAa,OAAO,KAAK,YAAY,MAAM,CAAC,CAAC,CAAC;IACnF,wCAAwC;IACxB,aAAa,CAAC,KAAoB;QAChD,IAAI,KAAK,YAAY,MAAM,EAAE;YAC3B,IAAI,IAAI,CAAC,MAAM,KAAK,KAAK,CAAC,MAAM;gBAAE,OAAO,KAAK,CAAC;YAC/C,IAAI,CAAC,IAAI,CAAC,aAAa,CAAC,aAAa,CAAC,KAAK,CAAC,aAAa,CAAC;gBAAE,OAAO,KAAK,CAAC;YACzE,OAAO,IAAI,CAAC;SACb;QACD,OAAO,KAAK,CAAC;IACf,CAAC;IACD;;;;;OAKG;IACI,sBAAsB,CAAC,CAAS,EAAE,gBAAoC,EAC3E,OAAuB;QACvB,IAAI,WAAW,GAAG,EAAE,CAAC;QACrB,IAAI,gBAAgB,KAAK,SAAS,IAAI,MAAM,CAAC,QAAQ,CAAC,gBAAgB,CAAC,EAAE;YACvE,WAAW,GAAG,gBAAgB,CAAC;SAChC;aAAM,IAAI,OAAO,YAAY,aAAa,EAAE;YAC3C,WAAW,GAAG,OAAO,CAAC,oBAAoB,CAAC,QAAQ,CAAC,KAAK,CAAC,IAAI,CAAC,aAAa,CAAC,MAAM,CAAC,gBAAgB,EAAE,EAAE,IAAI,CAAC,aAAa,CAAC,MAAM,CAAC,gBAAgB,EAAE,CAAC,CAAC,CAAC;SACxJ;QACD,WAAW,GAAG,QAAQ,CAAC,eAAe,CAAC,WAAW,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC;QAC3D,MAAM,SAAS,GAAG,IAAI,CAAC,aAAa,CAAC;QACrC,MAAM,GAAG,GAAG,IAAI,CAAC,kBAAkB,CAAC,CAAC,CAAC,CAAC;QACvC,MAAM,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;QACzB,MAAM,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;QACzB,IAAI,EAAE,EAAE,EAAE,CAAC;QACX,MAAM,MAAM,GAAG,YAAY,CAAC,gBAAgB,CAAC,gBAAgB,EAAE,OAAO,CAAC,CAAC;QACxE,MAAM,YAAY,GAAG,IAAI,CAAC,EAAE,GAAG,GAAG,GAAG,WAAW,CAAC;QACjD,MAAM,SAAS,GAAG,MAAM,CAAC,SAAS,CAAC,CAAK,yBAAyB;QACjE,MAAM,WAAW,GAAG,MAAM,CAAC,iBAAiB,CAAC,CAAC,yBAAyB;QACvE,MAAM,QAAQ,GAAG,MAAM,CAAC,cAAc,CAAC,CAAC,wBAAwB;QAChE,MAAM,cAAc,GAAG,MAAM,CAAC,oBAAoB,CAAC;QACnD,MAAM,IAAI,GAAG,QAAQ,CAAC,MAAM,EAAE,CAAC;QAC/B,MAAM,IAAI,GAAG,QAAQ,CAAC,MAAM,EAAE,CAAC;QAC/B,MAAM,MAAM,GAAG,QAAQ,CAAC,MAAM,EAAE,CAAC;QACjC,IAAI,OAAO,GAAG,CAAC,CAAC;QAChB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,WAAW,EAAE,CAAC,EAAE,EAAE;YACrC,IAAI,CAAC,GAAG,CAAC,IAAI,WAAW;gBACtB,OAAO,GAAG,CAAC,GAAG,YAAY,CAAC;;gBAE3B,OAAO,GAAG,CAAC,CAAC,GAAG,WAAW,CAAC,GAAG,YAAY,CAAC;YAC7C,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,OAAO,CAAC,CAAC;YACvB,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,OAAO,CAAC,CAAC;YACvB,MAAM,GAAG,GAAG,SAAS,CAAC,WAAW,CAAC,EAAE,GAAG,EAAE,EAAE,EAAE,GAAG,EAAE,EAAE,EAAE,CAAC,CAAC;YACxD,MAAM,CAAC,QAAQ,CAAC,GAAG,CAAC,CAAC;YAErB,IAAI,SAAS;gBACX,SAAS,CAAC,IAAI,CAAC,CAAC,GAAG,WAAW,CAAC,CAAC;YAElC,IAAI,WAAW,EAAE;gBACf,SAAS,CAAC,MAAM,CAAC,WAAW,CAAC,CAAC,EAAE,GAAG,EAAE,EAAE,EAAE,GAAG,EAAE,EAAE,GAAG,EAAE,IAAI,CAAC,CAAC;gBAC3D,WAAW,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC;aACxB;YACD,IAAI,QAAQ,EAAE;gBACZ,QAAQ,CAAC,MAAM,CAAC,CAAC,GAAG,WAAW,EAAE,CAAC,CAAC,CAAC;aACrC;YACD,IAAI,cAAc,EAAE;gBAClB,SAAS,CAAC,MAAM,CAAC,WAAW,CAAC,CAAC,EAAE,EAAE,EAAE,EAAE,CAAC,EAAE,IAAI,CAAC,CAAC;gBAC/C,SAAS,CAAC,MAAM,CAAC,WAAW,CAAC,CAAC,EAAE,GAAG,EAAE,EAAE,CAAC,EAAE,GAAG,EAAE,EAAE,EAAE,EAAE,IAAI,CAAC,CAAC;gBAC3D,IAAI,CAAC,gBAAgB,CAAC,IAAI,EAAE,MAAM,CAAC,CAAC;gBACpC,cAAc,CAAC,IAAI,CAAC,MAAM,CAAC,CAAC;aAC7B;SACF;QACD,OAAO,MAAM,CAAC;IAChB,CAAC;IAED,yEAAyE;IAClE,yBAAyB,CAAC,OAAwB;QACvD,OAAO,OAAO,CAAC,YAAY,CAAC,IAAI,CAAC,CAAC;IACpC,CAAC;IACD;;;OAGG;IACI,gBAAgB,CAAC,SAAiB;QACvC,MAAM,GAAG,GAAG,IAAI,CAAC,cAAc,CAAC,iBAAiB,CAAC,SAAS,CAAC,CAAC;QAC7D,MAAM,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;QACzB,MAAM,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;QACzB,MAAM,SAAS,GAAG,IAAI,CAAC,aAAa,CAAC;QACrC,MAAM,MAAM,GAAG,SAAS,CAAC,WAAW,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC;QAC/C,MAAM,OAAO,GAAG,SAAS,CAAC,MAAM,CAAC,WAAW,CAAC,EAAE,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;QACvD,MAAM,QAAQ,GAAG,SAAS,CAAC,MAAM,CAAC,WAAW,CAAC,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC,CAAC;QACxD,OAAO,IAAI,CAAC,MAAM,CAAC,KAAK,CAAC,MAAM,CAAC,MAAM,EAAE,OAAO,EAAE,QAAQ,CAAC,CAAC,CAAC;IAC9D,CAAC;IACD,6CAA6C;IACtC,WAAW,CAAC,KAAc,EAAE,SAAqB;QACtD,IAAI,SAAS,GAAG,IAAI,CAAC,aAAa,CAAC;QACnC,IAAI,SAAS,EAAE;YACb,SAAS,GAAG,SAAS,CAAC,0BAA0B,CAAC,SAAS,CAAC,CAAC;SAC7D;QAED,KAAK,CAAC,oBAAoB,CAAC,SAAS,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;QAClD,KAAK,CAAC,oBAAoB,CAAC,SAAS,EAAE,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;QACjD,KAAK,CAAC,oBAAoB,CAAC,SAAS,EAAE,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;QACjD,KAAK,CAAC,oBAAoB,CAAC,SAAS,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;QAEhD,KAAK,CAAC,oBAAoB,CAAC,SAAS,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;QACjD,KAAK,CAAC,oBAAoB,CAAC,SAAS,EAAE,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;QAChD,KAAK,CAAC,oBAAoB,CAAC,SAAS,EAAE,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;QAChD,KAAK,CAAC,oBAAoB,CAAC,SAAS,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;IAEjD,CAAC;IACD;;;OAGG;IACI,iBAAiB,CAAC,SAAiB,EAAE,SAAiB,EAAE,MAAgB;QAC7E,6BAA6B;QAC7B,MAAM,YAAY,GAAG,IAAI,CAAC,kBAAkB,CAAC,SAAS,CAAC,CAAC;QACxD,MAAM,UAAU,GAAG,IAAI,CAAC,kBAAkB,CAAC,SAAS,CAAC,CAAC;QACtD,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,OAAO,IAAI,CAAC,aAAa,CAAC,WAAW,CAAC,QAAQ,GAAG,MAAM,EAAE,QAAQ,GAAG,MAAM,EAAE,MAAM,EAAE,MAAM,CAAC,CAAC;IAC9F,CAAC;IACD;;;OAGG;IACI,4BAA4B,CAAC,SAAiB,EAAE,SAAiB,EAAE,MAAkC;QAC1G,MAAM,YAAY,GAAG,IAAI,CAAC,kBAAkB,CAAC,SAAS,CAAC,CAAC;QACxD,MAAM,UAAU,GAAG,IAAI,CAAC,kBAAkB,CAAC,SAAS,CAAC,CAAC;QACtD,MAAM,MAAM,GAAG,IAAI,CAAC,EAAE,GAAG,GAAG,CAAC;QAC7B,MAAM,IAAI,GAAG,IAAI,CAAC,cAAc,CAAC,YAAY,CAAC;QAC9C,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,OAAO,yBAAyB,CAAC,sBAAsB,CACrD,IAAI,CAAC,aAAa,CAAC,WAAW,CAAC,QAAQ,GAAG,MAAM,EAAE,QAAQ,GAAG,MAAM,EAAE,MAAM,CAAC,EAC5E,IAAI,CAAC,aAAa,CAAC,MAAM,CAAC,WAAW,CAAC,CAAC,MAAM,GAAG,QAAQ,EAAE,MAAM,GAAG,QAAQ,EAAE,CAAC,CAAC,EAAI,sFAAsF;QACzK,IAAI,CAAC,aAAa,CAAC,MAAM,CAAC,WAAW,CAAC,CAAC,IAAI,GAAG,QAAQ,GAAG,MAAM,EAAE,CAAC,IAAI,GAAG,QAAQ,GAAG,MAAM,EAAE,IAAI,GAAG,MAAM,CAAC,EAC1G,MAAM,CAAC,CAAC;IACZ,CAAC;IACD;;;;;OAKG;IACH,IAAW,cAAc;QACvB,OAAO,IAAI,CAAC,MAAM,IAAI,IAAI,CAAC,cAAc,CAAC,mBAAmB,CAAC;IAChE,CAAC;IACD;;;;OAIG;IACI,wBAAwB;QAC7B,wEAAwE;QACxE,MAAM,EAAE,GAAG,IAAI,CAAC,aAAa,CAAC,MAAM,CAAC,gBAAgB,EAAE,CAAC;QACxD,MAAM,EAAE,GAAG,IAAI,CAAC,aAAa,CAAC,MAAM,CAAC,gBAAgB,EAAE,CAAC;QACxD,MAAM,EAAE,GAAG,IAAI,CAAC,aAAa,CAAC,MAAM,CAAC,gBAAgB,EAAE,CAAC;QACxD,MAAM,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC;QAC/B,IAAI,KAAK,GAAG,IAAI,CAAC,EAAE,GAAG,GAAG,GAAG,KAAK,CAAC;QAClC,IAAI,CAAC,IAAI,CAAC,cAAc,CAAC,gBAAgB,CAAC,GAAG,CAAC;YAC5C,KAAK,IAAI,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,cAAc,CAAC,YAAY,CAAC,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,cAAc,CAAC,UAAU,CAAC,CAAC,CAAC,CAAC;QAC9H,MAAM,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,EAAE,EAAE,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,cAAc,CAAC,YAAY,CAAC,CAAC;QAE/E,OAAO,QAAQ,CAAC,MAAM,CAAC,KAAK,EAAE,KAAK,CAAC,CAAC;IACvC,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 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 { LineString3d } from \"../curve/LineString3d\";\r\nimport { Loop } from \"../curve/Loop\";\r\nimport { StrokeOptions } from \"../curve/StrokeOptions\";\r\nimport { Geometry } from \"../Geometry\";\r\nimport { AngleSweep } from \"../geometry3d/AngleSweep\";\r\nimport { GeometryHandler, UVSurface } from \"../geometry3d/GeometryHandler\";\r\nimport { Matrix3d } from \"../geometry3d/Matrix3d\";\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 Sphere is\r\n *\r\n * * A unit sphere (but read on ....)\r\n * * mapped by an arbitrary (possibly skewed, non-uniform scaled) transform\r\n * * hence possibly the final geometry is ellipsoidal\r\n * @public\r\n */\r\nexport class Sphere extends SolidPrimitive implements UVSurface {\r\n  /** String name for schema properties */\r\n  public readonly solidPrimitiveType = \"sphere\";\r\n\r\n  private _localToWorld: Transform;  // unit sphere maps to world through the transform0 part of this map.\r\n  private _latitudeSweep: AngleSweep;\r\n  /** Return the latitude (in radians) all fractional v. */\r\n  public vFractionToRadians(v: number): number {\r\n    return this._latitudeSweep.fractionToRadians(v);\r\n  }\r\n  /** Return the longitude (in radians) all fractional u. */\r\n  public uFractionToRadians(u: number): number {\r\n    return u * Math.PI * 2.0;\r\n  }\r\n\r\n  private constructor(localToWorld: Transform, latitudeSweep: AngleSweep, capped: boolean) {\r\n    super(capped);\r\n    this._localToWorld = localToWorld;\r\n    this._latitudeSweep = latitudeSweep ? latitudeSweep : AngleSweep.createFullLatitude();\r\n  }\r\n  /** return a deep clone */\r\n  public clone(): Sphere {\r\n    return new Sphere(this._localToWorld.clone(), this._latitudeSweep.clone(), this.capped);\r\n  }\r\n  /** Transform the sphere in place.\r\n   * * Fails if the transform is singular.\r\n   */\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  /** Return a transformed clone. */\r\n  public cloneTransformed(transform: Transform): Sphere | undefined {\r\n    const sphere1 = this.clone();\r\n    transform.multiplyTransformTransform(sphere1._localToWorld, sphere1._localToWorld);\r\n    if (transform.matrix.determinant() < 0.0) {\r\n      if (sphere1._latitudeSweep !== undefined) {\r\n        sphere1._latitudeSweep.reverseInPlace();\r\n      }\r\n    }\r\n    return sphere1;\r\n  }\r\n  /** Return a coordinate frame (right handed, unit axes)\r\n   * * origin at sphere center\r\n   * * equator 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 latitude sweep as fraction of south pole to north pole. */\r\n  public get latitudeSweepFraction(): number { return this._latitudeSweep.sweepRadians / Math.PI; }\r\n  /** Create from center and radius, with optional restricted latitudes. */\r\n  public static createCenterRadius(center: Point3d, radius: number, latitudeSweep?: AngleSweep): Sphere {\r\n    const localToWorld = Transform.createOriginAndMatrix(center, Matrix3d.createUniformScale(radius));\r\n    return new Sphere(localToWorld,\r\n      latitudeSweep ? latitudeSweep : AngleSweep.createFullLatitude(), false);\r\n  }\r\n  /** Create an ellipsoid which is a unit sphere mapped to position by an (arbitrary, possibly skewed and scaled) transform. */\r\n  public static createEllipsoid(localToWorld: Transform, latitudeSweep: AngleSweep, capped: boolean): Sphere | undefined {\r\n    return new Sphere(localToWorld, latitudeSweep, capped);\r\n  }\r\n\r\n  /** Create a sphere from the typical parameters of the Dgn file */\r\n  public static createDgnSphere(center: Point3d, vectorX: Vector3d, vectorZ: Vector3d, radiusXY: number, radiusZ: number,\r\n    latitudeSweep: AngleSweep,\r\n    capped: boolean): Sphere | undefined {\r\n    const vectorY = vectorX.rotate90Around(vectorZ);\r\n    if (vectorY && !vectorX.isParallelTo(vectorZ)) {\r\n      const matrix = Matrix3d.createColumns(vectorX, vectorY, vectorZ);\r\n      matrix.scaleColumns(radiusXY, radiusXY, radiusZ, matrix);\r\n      const frame = Transform.createOriginAndMatrix(center, matrix);\r\n      return new Sphere(frame, latitudeSweep.clone(), capped);\r\n    }\r\n    return undefined;\r\n  }\r\n\r\n  /** Create a sphere from the typical parameters of the Dgn file */\r\n  public static createFromAxesAndScales(center: Point3d, axes: undefined | Matrix3d, radiusX: number, radiusY: number, radiusZ: number,\r\n    latitudeSweep: AngleSweep | undefined,\r\n    capped: boolean): Sphere | undefined {\r\n    const localToWorld = Transform.createOriginAndMatrix(center, axes);\r\n    localToWorld.matrix.scaleColumnsInPlace(radiusX, radiusY, radiusZ);\r\n    return new Sphere(localToWorld, latitudeSweep ? latitudeSweep.clone() : AngleSweep.createFullLatitude(), capped);\r\n  }\r\n\r\n  /** return (copy of) sphere center */\r\n  public cloneCenter(): Point3d { return this._localToWorld.getOrigin(); }\r\n  /** return the (full length, i.e. scaled by radius) X vector from the sphere transform */\r\n  public cloneVectorX(): Vector3d { return this._localToWorld.matrix.columnX(); }\r\n  /** return the (full length, i.e. scaled by radius) Y vector from the sphere transform */\r\n  public cloneVectorY(): Vector3d { return this._localToWorld.matrix.columnY(); }\r\n  /** return the (full length, i.e. scaled by radius) Z vector from the sphere transform */\r\n  public cloneVectorZ(): Vector3d { return this._localToWorld.matrix.columnZ(); }\r\n  /** return (a copy of) the sphere's angle sweep. */\r\n  public cloneLatitudeSweep(): AngleSweep { return this._latitudeSweep.clone(); }\r\n  /** Test if the geometry is a true sphere taking the transform (which might have nonuniform scaling) is applied. */\r\n  public trueSphereRadius(): number | undefined {\r\n    const factors = this._localToWorld.matrix.factorRigidWithSignedScale();\r\n    if (!factors) return undefined;\r\n    if (factors && factors.scale > 0)\r\n      return factors.scale;\r\n    return undefined;\r\n  }\r\n  /**\r\n   * Return the larger of the primary xyz axis radii\r\n   */\r\n  public maxAxisRadius(): number {\r\n    const matrix = this._localToWorld.matrix;\r\n    return Geometry.maxXYZ(matrix.columnXMagnitude(), matrix.columnYMagnitude(), matrix.columnYMagnitude());\r\n  }\r\n  /**\r\n   * Return a (clone of) the sphere's local to world transformation.\r\n   */\r\n  public cloneLocalToWorld(): Transform { return this._localToWorld.clone(); }\r\n  /** Test if `other` is a `Sphere` */\r\n  public isSameGeometryClass(other: any): boolean { return other instanceof Sphere; }\r\n  /** Test for same geometry in `other` */\r\n  public override isAlmostEqual(other: GeometryQuery): boolean {\r\n    if (other instanceof Sphere) {\r\n      if (this.capped !== other.capped) return false;\r\n      if (!this._localToWorld.isAlmostEqual(other._localToWorld)) return false;\r\n      return true;\r\n    }\r\n    return false;\r\n  }\r\n  /**\r\n   *  return strokes for a cross-section (elliptic arc) at specified fraction v along the axis.\r\n   * * if strokeOptions is supplied, it is applied to the equator radii.\r\n   * @param v fractional position along the cone axis\r\n   * @param strokes stroke count or options.\r\n   */\r\n  public strokeConstantVSection(v: number, fixedStrokeCount: number | undefined,\r\n    options?: StrokeOptions): LineString3d {\r\n    let strokeCount = 16;\r\n    if (fixedStrokeCount !== undefined && Number.isFinite(fixedStrokeCount)) {\r\n      strokeCount = fixedStrokeCount;\r\n    } else if (options instanceof StrokeOptions) {\r\n      strokeCount = options.applyTolerancesToArc(Geometry.maxXY(this._localToWorld.matrix.columnXMagnitude(), this._localToWorld.matrix.columnYMagnitude()));\r\n    }\r\n    strokeCount = Geometry.clampToStartEnd(strokeCount, 4, 64);\r\n    const transform = this._localToWorld;\r\n    const phi = this.vFractionToRadians(v);\r\n    const c1 = Math.cos(phi);\r\n    const s1 = Math.sin(phi);\r\n    let c0, s0;\r\n    const result = LineString3d.createForStrokes(fixedStrokeCount, options);\r\n    const deltaRadians = Math.PI * 2.0 / strokeCount;\r\n    const fractions = result.fractions;     // possibly undefined !!!\r\n    const derivatives = result.packedDerivatives; // possibly undefined !!!\r\n    const uvParams = result.packedUVParams; // possibly undefined !!\r\n    const surfaceNormals = result.packedSurfaceNormals;\r\n    const dXdu = Vector3d.create();\r\n    const dXdv = Vector3d.create();\r\n    const normal = Vector3d.create();\r\n    let radians = 0;\r\n    for (let i = 0; i <= strokeCount; i++) {\r\n      if (i * 2 <= strokeCount)\r\n        radians = i * deltaRadians;\r\n      else\r\n        radians = (i - strokeCount) * deltaRadians;\r\n      c0 = Math.cos(radians);\r\n      s0 = Math.sin(radians);\r\n      const xyz = transform.multiplyXYZ(c1 * c0, c1 * s0, s1);\r\n      result.addPoint(xyz);\r\n\r\n      if (fractions)\r\n        fractions.push(i / strokeCount);\r\n\r\n      if (derivatives) {\r\n        transform.matrix.multiplyXYZ(-c1 * s0, c1 * c0, 0.0, dXdu);\r\n        derivatives.push(dXdu);\r\n      }\r\n      if (uvParams) {\r\n        uvParams.pushXY(i / strokeCount, v);\r\n      }\r\n      if (surfaceNormals) {\r\n        transform.matrix.multiplyXYZ(-s0, c0, 0, dXdu);\r\n        transform.matrix.multiplyXYZ(-s1 * c0, -s1 * s0, c1, dXdv);\r\n        dXdu.unitCrossProduct(dXdv, normal);\r\n        surfaceNormals.push(normal);\r\n      }\r\n    }\r\n    return result;\r\n  }\r\n\r\n  /** Second step of double dispatch:  call `handler.handleSphere(this)` */\r\n  public dispatchToGeometryHandler(handler: GeometryHandler): any {\r\n    return handler.handleSphere(this);\r\n  }\r\n  /**\r\n   * Return the Arc3d section at vFraction.  For the sphere, this is a latitude circle.\r\n   * @param vFraction fractional position along the sweep direction\r\n   */\r\n  public constantVSection(vFraction: number): CurveCollection | undefined {\r\n    const phi = this._latitudeSweep.fractionToRadians(vFraction);\r\n    const s1 = Math.sin(phi);\r\n    const c1 = Math.cos(phi);\r\n    const transform = this._localToWorld;\r\n    const center = transform.multiplyXYZ(0, 0, s1);\r\n    const vector0 = transform.matrix.multiplyXYZ(c1, 0, 0);\r\n    const vector90 = transform.matrix.multiplyXYZ(0, c1, 0);\r\n    return Loop.create(Arc3d.create(center, vector0, vector90));\r\n  }\r\n  /** Extend a range to contain this sphere. */\r\n  public extendRange(range: Range3d, transform?: Transform): void {\r\n    let placement = this._localToWorld;\r\n    if (transform) {\r\n      placement = transform.multiplyTransformTransform(placement);\r\n    }\r\n\r\n    range.extendTransformedXYZ(placement, -1, -1, -1);\r\n    range.extendTransformedXYZ(placement, 1, -1, -1);\r\n    range.extendTransformedXYZ(placement, -1, 1, -1);\r\n    range.extendTransformedXYZ(placement, 1, 1, -1);\r\n\r\n    range.extendTransformedXYZ(placement, -1, -1, 1);\r\n    range.extendTransformedXYZ(placement, 1, -1, 1);\r\n    range.extendTransformedXYZ(placement, -1, 1, 1);\r\n    range.extendTransformedXYZ(placement, 1, 1, 1);\r\n\r\n  }\r\n  /** Evaluate as a uv surface\r\n   * @param uFraction fractional position in minor (phi)\r\n   * @param vFraction fractional position on major (theta) arc\r\n   */\r\n  public uvFractionToPoint(uFraction: number, vFraction: number, result?: Point3d): Point3d {\r\n    // sphere with radius 1 . . .\r\n    const thetaRadians = this.uFractionToRadians(uFraction);\r\n    const phiRadians = this.vFractionToRadians(vFraction);\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    return this._localToWorld.multiplyXYZ(cosTheta * cosPhi, sinTheta * cosPhi, sinPhi, 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(uFraction: number, vFraction: number, result?: Plane3dByOriginAndVectors): Plane3dByOriginAndVectors {\r\n    const thetaRadians = this.uFractionToRadians(uFraction);\r\n    const phiRadians = this.vFractionToRadians(vFraction);\r\n    const fTheta = Math.PI * 2.0;\r\n    const fPhi = this._latitudeSweep.sweepRadians;\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    return Plane3dByOriginAndVectors.createOriginAndVectors(\r\n      this._localToWorld.multiplyXYZ(cosTheta * cosPhi, sinTheta * cosPhi, sinPhi),\r\n      this._localToWorld.matrix.multiplyXYZ(-fTheta * sinTheta, fTheta * cosTheta, 0),   // !!! note cosTheta term is omitted -- scale is wrong, but remains non-zero at poles.\r\n      this._localToWorld.matrix.multiplyXYZ(-fPhi * cosTheta * sinPhi, -fPhi * sinTheta * sinPhi, fPhi * cosPhi),\r\n      result);\r\n  }\r\n  /**\r\n   * * A sphere is can be closed two ways:\r\n   *   * full sphere (no caps needed for closure)\r\n   *   * incomplete but with caps\r\n   * @return true if this is a closed volume.\r\n   */\r\n  public get isClosedVolume(): boolean {\r\n    return this.capped || this._latitudeSweep.isFullLatitudeSweep;\r\n  }\r\n  /**\r\n   * Directional distance query\r\n   * * u direction is around longitude circle at maximum distance from axis.\r\n   * * v direction is on a line of longitude between the latitude limits.\r\n   */\r\n  public maxIsoParametricDistance(): Vector2d {\r\n    // approximate radius at equator .. if elliptic, this is not exact . . .\r\n    const rX = this._localToWorld.matrix.columnXMagnitude();\r\n    const rY = this._localToWorld.matrix.columnYMagnitude();\r\n    const rZ = this._localToWorld.matrix.columnZMagnitude();\r\n    const rMaxU = Math.max(rX, rY);\r\n    let dMaxU = Math.PI * 2.0 * rMaxU;\r\n    if (!this._latitudeSweep.isRadiansInSweep(0.0))\r\n      dMaxU *= Math.max(Math.cos(Math.abs(this._latitudeSweep.startRadians)), Math.cos(Math.abs(this._latitudeSweep.endRadians)));\r\n    const dMaxV = Math.max(rMaxU, rZ) * Math.abs(this._latitudeSweep.sweepRadians);\r\n\r\n    return Vector2d.create(dMaxU, dMaxV);\r\n  }\r\n}\r\n"]}

package/lib/esm/solid/SweepContour.d.ts

@@ -0,0 +1,82 @@
+import { UnionOfConvexClipPlaneSets } from "../clipping/UnionOfConvexClipPlaneSets";
+import { AnyCurve } from "../curve/CurveChain";
+import { CurveCollection } from "../curve/CurveCollection";
+import { StrokeOptions } from "../curve/StrokeOptions";
+import { Vector3d } from "../geometry3d/Point3dVector3d";
+import { Ray3d } from "../geometry3d/Ray3d";
+import { Transform } from "../geometry3d/Transform";
+import { IndexedPolyface } from "../polyface/Polyface";
+import { PolyfaceBuilder } from "../polyface/PolyfaceBuilder";
+import { MultiLineStringDataVariant } from "../topology/Triangulation";
+/**
+ * Sweepable contour with Transform for local to world interaction.
+ * * The surface/solid classes `LinearSweep`, `RotationalSweep`, `RuledSweep` use this for their swept contours.
+ * @public
+ */
+export declare class SweepContour {
+    /** The underlying curve collection, in its world coordinates position. */
+    curves: CurveCollection;
+    /** coordinate frame that in which the curves are all in the xy plane. */
+    localToWorld: Transform;
+    /** Axis used only in rotational case. */
+    axis: Ray3d | undefined;
+    private constructor();
+    /** Create for linear sweep.
+     * * The optional default normal may be useful for guiding coordinate frame setup.
+     * * the contour is CAPTURED.
+     */
+    static createForLinearSweep(contour: CurveCollection, defaultNormal?: Vector3d): SweepContour | undefined;
+    /** Create for linear sweep.
+     * * The optional default normal may be useful for guiding coordinate frame setup.
+     * * the points are captured into linestrings and Loops as needed.
+     */
+    static createForPolygon(points: MultiLineStringDataVariant, defaultNormal?: Vector3d): SweepContour | undefined;
+    /** Create for rotational sweep.
+     * * The axis ray is retained.
+     * * the contour is CAPTURED.
+     */
+    static createForRotation(contour: CurveCollection, axis: Ray3d): SweepContour | undefined;
+    /** Return (Reference to) the curves */
+    getCurves(): CurveCollection;
+    /** Apply `transform` to the curves, axis.
+     * * The local to world frame is reconstructed for the transformed curves.
+     */
+    tryTransformInPlace(transform: Transform): boolean;
+    /** Return a deep clone. */
+    clone(): SweepContour;
+    /** Return a transformed clone. */
+    cloneTransformed(transform: Transform): SweepContour | undefined;
+    /** Test for near equality of cures and local frame. */
+    isAlmostEqual(other: any): boolean;
+    private _xyStrokes?;
+    private _facets?;
+    get xyStrokes(): AnyCurve | undefined;
+    /**
+     * build the (cached) internal facets.
+     * @param options options for stroking the curves.
+     */
+    buildFacets(options: StrokeOptions | undefined): void;
+    /** delete existing facets.
+     * * This protects against PolyfaceBuilder reusing facets constructed with different options settings.
+     */
+    purgeFacets(): void;
+    /** Emit facets to a builder.
+     * This method may cache and reuse facets over multiple calls.
+     */
+    emitFacets(builder: PolyfaceBuilder, reverse: boolean, transform?: Transform): void;
+    /** Emit facets to a function
+     * This method may cache and reuse facets over multiple calls.
+     */
+    announceFacets(announce: (facets: IndexedPolyface) => void, options: StrokeOptions | undefined): void;
+    /**
+     * Triangulate the region.
+     * Create a UnionOfConvexClipPlaneSets that clips to the swept region.
+     * * If sweepVector is not given, the sweep direction is perpendicular to the plane of the contour.
+     * * If sweepVector is given, it is the sweep direction and does not have to be perpendicular to the contour.
+     * * cap0 and cap1 indicate construction of clip planes parallel to the contour plane.
+     * * If cap1 is true, the cap plane is at `anyPointOnPlane + sweepVector`.  That is, the sweep vector indicates both direction and distance.
+     * * caps are NOT created of sweepVector is not given.
+     */
+    sweepToUnionOfConvexClipPlaneSets(sweepVector?: Vector3d, cap0?: boolean, cap1?: boolean): UnionOfConvexClipPlaneSets | undefined;
+}
+//# sourceMappingURL=SweepContour.d.ts.map

package/lib/esm/solid/SweepContour.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"SweepContour.d.ts","sourceRoot":"","sources":["../../../src/solid/SweepContour.ts"],"names":[],"mappings":"AAWA,OAAO,EAAE,0BAA0B,EAAE,MAAM,wCAAwC,CAAC;AACpF,OAAO,EAAE,QAAQ,EAAE,MAAM,qBAAqB,CAAC;AAC/C,OAAO,EAAE,eAAe,EAAE,MAAM,0BAA0B,CAAC;AAK3D,OAAO,EAAE,aAAa,EAAE,MAAM,wBAAwB,CAAC;AAEvD,OAAO,EAAW,QAAQ,EAAE,MAAM,+BAA+B,CAAC;AAElE,OAAO,EAAE,KAAK,EAAE,MAAM,qBAAqB,CAAC;AAC5C,OAAO,EAAE,SAAS,EAAE,MAAM,yBAAyB,CAAC;AACpD,OAAO,EAAE,eAAe,EAAE,MAAM,sBAAsB,CAAC;AACvD,OAAO,EAAE,eAAe,EAAE,MAAM,6BAA6B,CAAC;AAE9D,OAAO,EAAE,0BAA0B,EAAgB,MAAM,2BAA2B,CAAC;AAErF;;;;GAIG;AACH,qBAAa,YAAY;IACvB,0EAA0E;IACnE,MAAM,EAAE,eAAe,CAAC;IAC/B,yEAAyE;IAClE,YAAY,EAAE,SAAS,CAAC;IAC/B,yCAAyC;IAClC,IAAI,EAAE,KAAK,GAAG,SAAS,CAAC;IAE/B,OAAO;IAKP;;;OAGG;WACW,oBAAoB,CAAC,OAAO,EAAE,eAAe,EAAE,aAAa,CAAC,EAAE,QAAQ,GAAG,YAAY,GAAG,SAAS;IAQhH;;;OAGG;WACW,gBAAgB,CAAC,MAAM,EAAE,0BAA0B,EAAE,aAAa,CAAC,EAAE,QAAQ,GAAG,YAAY,GAAG,SAAS;IAsBtH;;;OAGG;WACW,iBAAiB,CAAC,OAAO,EAAE,eAAe,EAAE,IAAI,EAAE,KAAK,GAAG,YAAY,GAAG,SAAS;IAQhG,uCAAuC;IAChC,SAAS,IAAI,eAAe;IACnC;;OAEG;IACI,mBAAmB,CAAC,SAAS,EAAE,SAAS,GAAG,OAAO;IAezD,2BAA2B;IACpB,KAAK,IAAI,YAAY;IAG5B,kCAAkC;IAC3B,gBAAgB,CAAC,SAAS,EAAE,SAAS,GAAG,YAAY,GAAG,SAAS;IAMvE,uDAAuD;IAChD,aAAa,CAAC,KAAK,EAAE,GAAG,GAAG,OAAO;IAOzC,OAAO,CAAC,UAAU,CAAC,CAAW;IAC9B,OAAO,CAAC,OAAO,CAAC,CAAkB;IAClC,IAAW,SAAS,IAAI,QAAQ,GAAG,SAAS,CAA2B;IACvE;;;OAGG;IACI,WAAW,CAAC,OAAO,EAAE,aAAa,GAAG,SAAS,GAAG,IAAI;IAsD5D;;OAEG;IACI,WAAW;IAIlB;;OAEG;IACI,UAAU,CAAC,OAAO,EAAE,eAAe,EAAE,OAAO,EAAE,OAAO,EAAE,SAAS,CAAC,EAAE,SAAS;IAMnF;;OAEG;IACI,cAAc,CAAC,QAAQ,EAAE,CAAC,MAAM,EAAE,eAAe,KAAK,IAAI,EAAE,OAAO,EAAE,aAAa,GAAG,SAAS;IAMrG;;;;;;;;OAQG;IACI,iCAAiC,CAAC,WAAW,CAAC,EAAE,QAAQ,EAAE,IAAI,GAAE,OAAe,EAAE,IAAI,GAAE,OAAe,GAAG,0BAA0B,GAAG,SAAS;CAyCvJ"}