@itwin/core-geometry 5.0.0-dev.3 → 5.0.0-dev.34

package/lib/cjs/polyface/PolyfaceBuilder.js

@@ -5,11 +5,14 @@
 *--------------------------------------------------------------------------------------------*/
 Object.defineProperty(exports, "__esModule", { value: true });
 exports.PolyfaceBuilder = void 0;
+/** @packageDocumentation
+ * @module Polyface
+ */
+const Arc3d_1 = require("../curve/Arc3d");
 const ConstructCurveBetweenCurves_1 = require("../curve/ConstructCurveBetweenCurves");
 const CurveCollection_1 = require("../curve/CurveCollection");
 const CurveFactory_1 = require("../curve/CurveFactory");
 const CurvePrimitive_1 = require("../curve/CurvePrimitive");
-const GeometryQuery_1 = require("../curve/GeometryQuery");
 const LineString3d_1 = require("../curve/LineString3d");
 const ParityRegion_1 = require("../curve/ParityRegion");
 const CylindricalRange_1 = require("../curve/Query/CylindricalRange");
@@ -23,7 +26,6 @@ const GeometryHandler_1 = require("../geometry3d/GeometryHandler");
 const GrowableFloat64Array_1 = require("../geometry3d/GrowableFloat64Array");
 const GrowableXYArray_1 = require("../geometry3d/GrowableXYArray");
 const GrowableXYZArray_1 = require("../geometry3d/GrowableXYZArray");
-const IndexedXYZCollection_1 = require("../geometry3d/IndexedXYZCollection");
 const Matrix3d_1 = require("../geometry3d/Matrix3d");
 const Plane3dByOriginAndVectors_1 = require("../geometry3d/Plane3dByOriginAndVectors");
 const Point2dVector2d_1 = require("../geometry3d/Point2dVector2d");
@@ -108,6 +110,7 @@ class FacetSector {
         if (Math.abs(uvw.z) < tol)
             uvw.z = 0.0;
     }
+    static { this._edgeVector = Point3dVector3d_1.Vector3d.create(); }
     /**
      * Given two sectors with xyz and sectionDerivative (u derivative), use the edge from A to B as v direction
      * in-surface derivative. Compute cross products (and normalize).
@@ -129,7 +132,6 @@ class FacetSector {
         }
     }
 }
-FacetSector._edgeVector = Point3dVector3d_1.Vector3d.create();
 /**
  * * Simple construction for strongly typed GeometryQuery objects:
  *
@@ -340,6 +342,9 @@ class PolyfaceBuilder extends GeometryHandler_1.NullGeometryHandler {
     findOrAddParamXY(x, y) {
         return this.addParamXY(x, y);
     }
+    static { this._workPointFindOrAddA = Point3dVector3d_1.Point3d.create(); }
+    static { this._workVectorFindOrAdd = Point3dVector3d_1.Vector3d.create(); }
+    static { this._workUVFindOrAdd = Point2dVector2d_1.Point2d.create(); }
     /**
      * Add a point to the polyface. The implementation is free to either create a new point or return the index of a
      * prior point with the same coordinates.
@@ -1080,17 +1085,28 @@ class PolyfaceBuilder extends GeometryHandler_1.NullGeometryHandler {
     }
     /** Construct facets for a rotational sweep. */
     addRotationalSweep(surface) {
-        const contour = surface.getCurves();
-        const section0 = StrokeCountChain_1.StrokeCountSection.createForParityRegionOrChain(contour, this._options);
+        const contour = surface.getSweepContourRef();
+        const section0 = StrokeCountChain_1.StrokeCountSection.createForParityRegionOrChain(contour.getCurves(), this._options);
         const baseStrokes = section0.getStrokes();
+        // ensure sweep is positive for buildRotationalNormalsInLineStrings
         const axis = surface.cloneAxisRay();
+        const sweepAngle = surface.getSweep();
+        if (sweepAngle.radians < 0.0) {
+            axis.direction.scaleInPlace(-1);
+            sweepAngle.setRadians(-sweepAngle.radians);
+        }
+        // swingVector points in the direction of positive sweep
         const perpendicularVector = CylindricalRange_1.CylindricalRangeQuery.computeMaxVectorFromRay(axis, baseStrokes);
         const swingVector = axis.direction.crossProduct(perpendicularVector);
+        // ensure contour computed normal is aligned with swingVector for buildRotationalNormalsInLineStrings
+        const contourNormalAgreesWithSwingDir = contour.localToWorld.matrix.dotColumnZ(swingVector) > 0;
+        if (!contourNormalAgreesWithSwingDir)
+            baseStrokes.reverseInPlace();
         if (this._options.needNormals)
             CylindricalRange_1.CylindricalRangeQuery.buildRotationalNormalsInLineStrings(baseStrokes, axis, swingVector);
         const maxDistance = perpendicularVector.magnitude();
-        const maxPath = Math.abs(maxDistance * surface.getSweep().radians);
-        let numStep = StrokeOptions_1.StrokeOptions.applyAngleTol(this._options, 1, surface.getSweep().radians, undefined);
+        const maxPath = Math.abs(maxDistance * sweepAngle.radians);
+        let numStep = StrokeOptions_1.StrokeOptions.applyAngleTol(this._options, 1, sweepAngle.radians, undefined);
         numStep = StrokeOptions_1.StrokeOptions.applyMaxEdgeLength(this._options, numStep, maxPath);
         for (let i = 1; i <= numStep; i++) {
             const transformA = surface.getFractionalRotationTransform((i - 1) / numStep);
@@ -1098,11 +1114,10 @@ class PolyfaceBuilder extends GeometryHandler_1.NullGeometryHandler {
             this.addBetweenRotatedStrokeSets(baseStrokes, transformA, i - 1, transformB, i);
         }
         if (surface.capped) {
-            const capContour = surface.getSweepContourRef();
-            capContour.purgeFacets();
-            capContour.emitFacets(this, true, undefined);
-            // final loop pass left transformA at end
-            capContour.emitFacets(this, false, surface.getFractionalRotationTransform(1.0));
+            contour.purgeFacets();
+            const reverseNearCap = contourNormalAgreesWithSwingDir;
+            contour.emitFacets(this, reverseNearCap, undefined);
+            contour.emitFacets(this, !reverseNearCap, surface.getFractionalRotationTransform(1.0));
         }
     }
     /** Construct facets for any planar region. */
@@ -1166,9 +1181,9 @@ class PolyfaceBuilder extends GeometryHandler_1.NullGeometryHandler {
             const normalIndices = ls.ensureEmptyNormalIndices();
             const normalIndex0 = this.findOrAddNormalInLineString(ls, 0, transform);
             normalIndices.push(normalIndex0);
-            let normalIndexA = normalIndex0;
-            let normalIndexB;
             if (n > 1) {
+                let normalIndexA = normalIndex0;
+                let normalIndexB;
                 for (let i = 1; i + 1 < n; i++) {
                     normalIndexB = this.findOrAddNormalInLineString(ls, i, transform, normalIndexA);
                     normalIndices.push(normalIndexB);
@@ -1183,9 +1198,9 @@ class PolyfaceBuilder extends GeometryHandler_1.NullGeometryHandler {
             const uvIndices = ls.ensureEmptyUVIndices();
             const uvIndex0 = this.findOrAddParamInLineString(ls, 0, vParam);
             uvIndices.push(uvIndex0);
-            let uvIndexA = uvIndex0;
-            let uvIndexB;
             if (n > 1) {
+                let uvIndexA = uvIndex0;
+                let uvIndexB;
                 for (let i = 1; i + 1 < n; i++) {
                     uvIndexB = this.findOrAddParamInLineString(ls, i, vParam, uvIndexA);
                     uvIndices.push(uvIndexB);
@@ -1820,56 +1835,39 @@ class PolyfaceBuilder extends GeometryHandler_1.NullGeometryHandler {
         const context = GreedyTriangulationBetweenLineStrings_1.GreedyTriangulationBetweenLineStrings.createContext();
         context.emitTriangles(resolveToIndexedXYZCollectionOrCarrier(pointsA), resolveToIndexedXYZCollectionOrCarrier(pointsB), (triangle) => { this.addTriangleFacet(triangle.points); });
     }
-    addMiteredPipesFromPoints(centerline, sectionData, numFacetAround = 12) {
-        const sections = CurveFactory_1.CurveFactory.createMiteredPipeSections(centerline, sectionData);
-        const pointA0 = Point3dVector3d_1.Point3d.create();
-        const pointA1 = Point3dVector3d_1.Point3d.create();
-        const pointB0 = Point3dVector3d_1.Point3d.create();
-        const pointB1 = Point3dVector3d_1.Point3d.create();
-        if (numFacetAround < 3)
-            numFacetAround = 3;
-        const df = 1.0 / numFacetAround;
-        for (let i = 1; i < sections.length; i++) {
-            const arcA = sections[i - 1];
-            const arcB = sections[i];
-            arcA.fractionToPoint(0.0, pointA0);
-            arcB.fractionToPoint(0.0, pointB0);
-            for (let k = 1; k <= numFacetAround; k++, pointA0.setFromPoint3d(pointA1), pointB0.setFromPoint3d(pointB1)) {
-                const f = k * df;
-                arcA.fractionToPoint(f, pointA1);
-                arcB.fractionToPoint(f, pointB1);
-                this.addQuadFacet([pointA0, pointB0, pointB1, pointA1]);
-            }
-        }
-    }
     /**
      * Add quad facets along a mitered pipe that follows a centerline curve.
+     * * At the end of each pipe segment, the pipe is mitered by the plane that bisects the angle between successive
+     * centerline segments.
      * * Circular or elliptical pipe cross sections can be specified by supplying either a radius, a pair of semi-axis
-     * lengths, or a full Arc3d:
-     *    * For semi-axis length input, x corresponds to an ellipse local axis nominally situated parallel to the xy-plane.
-     *    * For Arc3d input, the center is translated to the centerline start point to act as initial cross section.
-     * @param centerline centerline of pipe. If curved, it will be stroked using the builder's StrokeOptions.
-     * @param sectionData circle radius, ellipse semi-axis lengths, or full Arc3d.
+     * lengths, or an Arc3d:
+     *    * For semi-axis length input, x and y correspond to ellipse local axes perpendicular to each other and to the
+     * start tangent.
+     *    * For Arc3d input, the center is translated to the centerline start point. For best results, ensure this arc
+     * is perpendicular to the centerline start tangent.
+     * * This function internally calls [[CurveFactory.createMiteredSweepSections]], passing in the `startTangent`
+     * option to preserve the rotation of Arc3d-type `sectionData`.
+     * @param centerline centerline of pipe. If curved, it will be stroked using the builder's StrokeOptions, otherwise
+     * for best results, ensure no successive duplicate points with e.g., [[GrowableXYZArray.createCompressed]].
+     * @param sectionData circle radius, ellipse semi-axis lengths, or Arc3d.
      * @param numFacetAround how many equal parameter-space chords around each section.
+     * @param capped if `true`, add a cap at each end of the pipe; defaults to `false`.
      */
-    addMiteredPipes(centerline, sectionData, numFacetAround = 12) {
-        if (Array.isArray(centerline)) {
-            this.addMiteredPipesFromPoints(new Point3dArrayCarrier_1.Point3dArrayCarrier(centerline), sectionData, numFacetAround);
-        }
-        else if (centerline instanceof GrowableXYZArray_1.GrowableXYZArray) {
-            this.addMiteredPipesFromPoints(centerline, sectionData, numFacetAround);
-        }
-        else if (centerline instanceof IndexedXYZCollection_1.IndexedXYZCollection) {
-            this.addMiteredPipesFromPoints(centerline, sectionData, numFacetAround);
-        }
-        else if (centerline instanceof LineString3d_1.LineString3d) {
-            this.addMiteredPipesFromPoints(centerline.packedPoints, sectionData, numFacetAround);
-        }
-        else if (centerline instanceof GeometryQuery_1.GeometryQuery) {
-            const linestring = LineString3d_1.LineString3d.create();
-            centerline.emitStrokes(linestring, this._options);
-            this.addMiteredPipesFromPoints(linestring.packedPoints, sectionData, numFacetAround);
-        }
+    addMiteredPipes(centerline, sectionData, numFacetAround = 12, capped = false) {
+        const arc = CurveFactory_1.CurveFactory.createArcFromSectionData(centerline, sectionData);
+        if (!arc)
+            return;
+        if (numFacetAround < 3)
+            numFacetAround = 3;
+        const section = LineString3d_1.LineString3d.create();
+        section.appendFractionalStrokePoints(arc, numFacetAround, 0.0, 1.0, true);
+        const strokeOptions = this._options.clone();
+        const outputSelect = CurveFactory_1.MiteredSweepOutputSelect.AlsoMesh;
+        const startTangent = sectionData instanceof Arc3d_1.Arc3d ? arc.binormalVector() : undefined; // preserve arc orientation
+        const options = { strokeOptions, capped, outputSelect, startTangent };
+        const sections = CurveFactory_1.CurveFactory.createMiteredSweepSections(centerline, section, options);
+        if (sections && sections.mesh)
+            this.addIndexedPolyface(sections.mesh);
     }
     /** Return the polyface index array indices corresponding to the given edge, or `undefined` if error. */
     getEdgeIndices(edge) {
@@ -1980,9 +1978,6 @@ class PolyfaceBuilder extends GeometryHandler_1.NullGeometryHandler {
     }
 }
 exports.PolyfaceBuilder = PolyfaceBuilder;
-PolyfaceBuilder._workPointFindOrAddA = Point3dVector3d_1.Point3d.create();
-PolyfaceBuilder._workVectorFindOrAdd = Point3dVector3d_1.Vector3d.create();
-PolyfaceBuilder._workUVFindOrAdd = Point2dVector2d_1.Point2d.create();
 function resolveToIndexedXYZCollectionOrCarrier(points) {
     if (Array.isArray(points))
         return new Point3dArrayCarrier_1.Point3dArrayCarrier(points);