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

package/lib/esm/polyface/BoxTopology.js

@@ -0,0 +1,99 @@
+/*---------------------------------------------------------------------------------------------
+* Copyright (c) Bentley Systems, Incorporated. All rights reserved.
+* See LICENSE.md in the project root for license terms and full copyright notice.
+*--------------------------------------------------------------------------------------------*/
+/** @packageDocumentation
+ * @module Polyface
+ */
+import { Point3d } from "../geometry3d/Point3dVector3d";
+//
+//      2------------------3
+//      | \     F4       / |
+//      |   6----------7   |
+//      |   |          |   |   (BOTTOM = F0)
+//      |F5 |   F1     |F3 |
+//      |   |          |   |
+//      |   4----------5   |
+//      | /     F2       \ |
+//      0------------------1
+//
+/**
+ * methods for gathering points and indices of a box (8 corners, 6 faces, 12 edges)
+ * @internal
+ */
+export class BoxTopology {
+    static pointsClone() {
+        const clones = [];
+        for (const p of this.points)
+            clones.push(p.clone());
+        return clones;
+    }
+}
+/**
+ * static readonly array with the coordinates of the 8 unit cube corners in standard order, which is:
+ * x varies fastest
+ * * The point indices for the x edges are (0 to 1), (2 to 3), (4 to 5), (6 to 7)
+ * * The point indices for the y edges are (0 to 2), (1 to 3), (4 to 6), (5 to 7)
+ * * The point indices for the z edges are (0 to 4), (1 to 5), (2 to 6), (3 to 7)
+ * * These indices are tabulated in the `axisEdgeVertex[axis][edge][vertex]` array
+ */
+BoxTopology.points = [
+    Point3d.create(0, 0, 0),
+    Point3d.create(1, 0, 0),
+    Point3d.create(0, 1, 0),
+    Point3d.create(1, 1, 0),
+    Point3d.create(0, 0, 1),
+    Point3d.create(1, 0, 1),
+    Point3d.create(0, 1, 1),
+    Point3d.create(1, 1, 1),
+];
+/** IN faceId pair, the first component for bottom and top caps is `primaryCapId` */
+BoxTopology.primaryCapId = -1;
+/** Indices of vertices around faces, in CCW from the outside. */
+BoxTopology.cornerIndexCCW = [
+    [1, 0, 2, 3],
+    [4, 5, 7, 6],
+    [0, 1, 5, 4],
+    [1, 3, 7, 5],
+    [3, 2, 6, 7],
+    [2, 0, 4, 6]
+];
+/**  // [partnerFace[faceIndex][k] = index of k'th adjacent face  */
+BoxTopology.partnerFace = [
+    [5, 4, 3, 2],
+    [2, 3, 4, 5],
+    [0, 3, 1, 5],
+    [0, 4, 1, 2],
+    [0, 5, 1, 3],
+    [0, 2, 1, 4],
+];
+/** face id as used in SolidPrimitive methods */
+BoxTopology.faceId = [
+    [BoxTopology.primaryCapId, 0],
+    [BoxTopology.primaryCapId, 1],
+    [0, 0],
+    [0, 1],
+    [0, 2],
+    [0, 3]
+];
+/**
+ * Table to look up axis indices of edges and normals in box faces.
+ * faceDirections[faceIndex] =[[edge0AxisIndex, edge1AxisIndex, normalAxisIndex],[direction sign for along the axis]
+ */
+BoxTopology.faceDirections = [
+    [[0, 1, 2], [-1, 1, -1]],
+    [[0, 1, 2], [1, 1, 1]],
+    [[0, 2, 1], [1, -1, 1]],
+    [[1, 2, 0], [1, 1, 1]],
+    [[0, 2, 1], [-1, 1, 1]],
+    [[1, 2, 0], [-1, 1, -1]]
+];
+/** There are 4 edges in each axis direction.
+ *  * axisEdgeVertex[axisIndex][edgeIndex 0..3][*] = vertex index at end of edge in axisIndex direction.
+ */
+BoxTopology.axisEdgeVertex = [
+    [[0, 1], [2, 3], [4, 5], [6, 7]],
+    [[0, 2], [1, 3], [4, 6], [5, 7]],
+    [[0, 4], [1, 5], [2, 6], [3, 7]]
+];
+//# sourceMappingURL=BoxTopology.js.map

package/lib/esm/polyface/BoxTopology.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"BoxTopology.js","sourceRoot":"","sources":["../../../src/polyface/BoxTopology.ts"],"names":[],"mappings":"AAAA;;;+FAG+F;AAE/F;;GAEG;AAEH,OAAO,EAAE,OAAO,EAAE,MAAM,+BAA+B,CAAC;AAExD,EAAE;AACF,4BAA4B;AAC5B,4BAA4B;AAC5B,4BAA4B;AAC5B,4CAA4C;AAC5C,4BAA4B;AAC5B,4BAA4B;AAC5B,4BAA4B;AAC5B,4BAA4B;AAC5B,4BAA4B;AAC5B,EAAE;AACF;;;GAGG;AACH,MAAM,OAAO,WAAW;IAmBf,MAAM,CAAC,WAAW;QACvB,MAAM,MAAM,GAAG,EAAE,CAAC;QAClB,KAAK,MAAM,CAAC,IAAI,IAAI,CAAC,MAAM;YACzB,MAAM,CAAC,IAAI,CAAC,CAAC,CAAC,KAAK,EAAE,CAAC,CAAC;QACzB,OAAO,MAAM,CAAC;IAChB,CAAC;;AAvBD;;;;;;;GAOG;AACoB,kBAAM,GAAG;IAC9B,OAAO,CAAC,MAAM,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;IACvB,OAAO,CAAC,MAAM,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;IACvB,OAAO,CAAC,MAAM,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;IACvB,OAAO,CAAC,MAAM,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;IACvB,OAAO,CAAC,MAAM,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;IACvB,OAAO,CAAC,MAAM,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;IACvB,OAAO,CAAC,MAAM,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;IACvB,OAAO,CAAC,MAAM,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;CACxB,CAAC;AAOF,oFAAoF;AAC7D,wBAAY,GAAG,CAAC,CAAC,CAAC;AACzC,iEAAiE;AAC1C,0BAAc,GACrC;IACE,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;IACZ,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;IACZ,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;IACZ,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;IACZ,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;IACZ,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;CAAC,CAAC;AAChB,oEAAoE;AAC7C,uBAAW,GAClC;IACE,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;IACZ,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;IACZ,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;IACZ,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;IACZ,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;IACZ,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;CACb,CAAC;AACF,gDAAgD;AACzB,kBAAM,GAAG;IAC9B,CAAC,WAAW,CAAC,YAAY,EAAE,CAAC,CAAC;IAC7B,CAAC,WAAW,CAAC,YAAY,EAAE,CAAC,CAAC;IAC7B,CAAC,CAAC,EAAE,CAAC,CAAC;IACN,CAAC,CAAC,EAAE,CAAC,CAAC;IACN,CAAC,CAAC,EAAE,CAAC,CAAC;IACN,CAAC,CAAC,EAAE,CAAC,CAAC;CAAC,CAAC;AACV;;;GAGG;AACoB,0BAAc,GACrC;IACE,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;IACxB,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;IACtB,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;IACvB,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;IACtB,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;IACvB,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;CAAC,CAAC;AAC5B;;GAEG;AACoB,0BAAc,GACrC;IACE,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;IAChC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;IAChC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;CAAC,CAAC","sourcesContent":["/*---------------------------------------------------------------------------------------------\r\n* Copyright (c) Bentley Systems, Incorporated. All rights reserved.\r\n* See LICENSE.md in the project root for license terms and full copyright notice.\r\n*--------------------------------------------------------------------------------------------*/\r\n\r\n/** @packageDocumentation\r\n * @module Polyface\r\n */\r\n\r\nimport { Point3d } from \"../geometry3d/Point3dVector3d\";\r\n\r\n//\r\n//      2------------------3\r\n//      | \\     F4       / |\r\n//      |   6----------7   |\r\n//      |   |          |   |   (BOTTOM = F0)\r\n//      |F5 |   F1     |F3 |\r\n//      |   |          |   |\r\n//      |   4----------5   |\r\n//      | /     F2       \\ |\r\n//      0------------------1\r\n//\r\n/**\r\n * methods for gathering points and indices of a box (8 corners, 6 faces, 12 edges)\r\n * @internal\r\n */\r\nexport class BoxTopology {\r\n  /**\r\n   * static readonly array with the coordinates of the 8 unit cube corners in standard order, which is:\r\n   * x varies fastest\r\n   * * The point indices for the x edges are (0 to 1), (2 to 3), (4 to 5), (6 to 7)\r\n   * * The point indices for the y edges are (0 to 2), (1 to 3), (4 to 6), (5 to 7)\r\n   * * The point indices for the z edges are (0 to 4), (1 to 5), (2 to 6), (3 to 7)\r\n   * * These indices are tabulated in the `axisEdgeVertex[axis][edge][vertex]` array\r\n   */\r\n  public static readonly points = [\r\n    Point3d.create(0, 0, 0),\r\n    Point3d.create(1, 0, 0),\r\n    Point3d.create(0, 1, 0),\r\n    Point3d.create(1, 1, 0),\r\n    Point3d.create(0, 0, 1),\r\n    Point3d.create(1, 0, 1),\r\n    Point3d.create(0, 1, 1),\r\n    Point3d.create(1, 1, 1),\r\n  ];\r\n  public static pointsClone(): Point3d[] {\r\n    const clones = [];\r\n    for (const p of this.points)\r\n      clones.push(p.clone());\r\n    return clones;\r\n  }\r\n  /** IN faceId pair, the first component for bottom and top caps is `primaryCapId` */\r\n  public static readonly primaryCapId = -1;\r\n  /** Indices of vertices around faces, in CCW from the outside. */\r\n  public static readonly cornerIndexCCW =\r\n  [\r\n    [1, 0, 2, 3],\r\n    [4, 5, 7, 6],\r\n    [0, 1, 5, 4],\r\n    [1, 3, 7, 5],\r\n    [3, 2, 6, 7],\r\n    [2, 0, 4, 6]];\r\n  /**  // [partnerFace[faceIndex][k] = index of k'th adjacent face  */\r\n  public static readonly partnerFace =\r\n  [\r\n    [5, 4, 3, 2],\r\n    [2, 3, 4, 5],\r\n    [0, 3, 1, 5],\r\n    [0, 4, 1, 2],\r\n    [0, 5, 1, 3],\r\n    [0, 2, 1, 4],\r\n  ];\r\n  /** face id as used in SolidPrimitive methods */\r\n  public static readonly faceId = [\r\n    [BoxTopology.primaryCapId, 0],\r\n    [BoxTopology.primaryCapId, 1],\r\n    [0, 0],\r\n    [0, 1],\r\n    [0, 2],\r\n    [0, 3]];\r\n  /**\r\n   * Table to look up axis indices of edges and normals in box faces.\r\n   * faceDirections[faceIndex] =[[edge0AxisIndex, edge1AxisIndex, normalAxisIndex],[direction sign for along the axis]\r\n   */\r\n  public static readonly faceDirections =\r\n  [\r\n    [[0, 1, 2], [-1, 1, -1]],\r\n    [[0, 1, 2], [1, 1, 1]],\r\n    [[0, 2, 1], [1, -1, 1]],\r\n    [[1, 2, 0], [1, 1, 1]],\r\n    [[0, 2, 1], [-1, 1, 1]],\r\n    [[1, 2, 0], [-1, 1, -1]]];\r\n  /** There are 4 edges in each axis direction.\r\n   *  * axisEdgeVertex[axisIndex][edgeIndex 0..3][*] = vertex index at end of edge in axisIndex direction.\r\n   */\r\n  public static readonly axisEdgeVertex =\r\n  [\r\n    [[0, 1], [2, 3], [4, 5], [6, 7]],\r\n    [[0, 2], [1, 3], [4, 6], [5, 7]],\r\n    [[0, 4], [1, 5], [2, 6], [3, 7]]];\r\n}\r\n"]}

package/lib/esm/polyface/FacetFaceData.d.ts

@@ -0,0 +1,44 @@
+/** @packageDocumentation
+ * @module Polyface
+ */
+import { Point2d } from "../geometry3d/Point2dVector2d";
+import { Range2d } from "../geometry3d/Range";
+import { IndexedPolyface } from "./Polyface";
+/**
+ * Data for a face in a polyface containing facets.
+ * This is built up cooperatively by the PolyfaceBuilder and its
+ * callers, and stored as a FaceData array in PolyfaceData.
+ * @public
+ */
+export declare class FacetFaceData {
+    private _paramDistanceRange;
+    private _paramRange;
+    /** (property accessor) Return a reference to the distance-scaled parameter range. */
+    get paramDistanceRange(): Range2d;
+    /** (property accessor) Return a reference to the parameter range. */
+    get paramRange(): Range2d;
+    private constructor();
+    /** Create a FacetFaceData with null ranges. */
+    static createNull(): FacetFaceData;
+    /** Create a deep copy of this FacetFaceData object. */
+    clone(result?: FacetFaceData): FacetFaceData;
+    /** Restore this FacetFaceData to its null constructor state. */
+    setNull(): void;
+    /** Return distance-based parameter from stored parameter value. */
+    convertParamXYToDistance(x: number, y: number, result?: Point2d): Point2d;
+    /** Return normalized (0-1) parameter from stored parameter value. */
+    convertParamXYToNormalized(x: number, y: number, result?: Point2d): Point2d;
+    /** Return distance-based parameter from stored parameter value. */
+    convertParamToDistance(param: Point2d, result?: Point2d): Point2d;
+    /** Return normalized (0-1) parameter from stored parameter value. */
+    convertParamToNormalized(param: Point2d, result?: Point2d): Point2d;
+    /** Scale distance parameters. */
+    scaleDistances(distanceScale: number): void;
+    /**
+     * Sets the param and paramDistance range of this FacetFaceData based on the newly terminated facets that make it up.
+     * Takes the polyface itself, the first and last indexes of the facets to be included in the face.
+     * Returns true on success, false otherwise.
+     */
+    setParamDistanceRangeFromNewFaceData(polyface: IndexedPolyface, facetStart: number, facetEnd: number): boolean;
+}
+//# sourceMappingURL=FacetFaceData.d.ts.map

package/lib/esm/polyface/FacetFaceData.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"FacetFaceData.d.ts","sourceRoot":"","sources":["../../../src/polyface/FacetFaceData.ts"],"names":[],"mappings":"AAIA;;GAEG;AAEH,OAAO,EAAE,OAAO,EAAY,MAAM,+BAA+B,CAAC;AAElE,OAAO,EAAE,OAAO,EAAE,MAAM,qBAAqB,CAAC;AAC9C,OAAO,EAAE,eAAe,EAAE,MAAM,YAAY,CAAC;AAG7C;;;;;GAKG;AACH,qBAAa,aAAa;IACxB,OAAO,CAAC,mBAAmB,CAAU;IACrC,OAAO,CAAC,WAAW,CAAU;IAC7B,qFAAqF;IACrF,IAAW,kBAAkB,IAAI,OAAO,CAAqC;IAC7E,qEAAqE;IACrE,IAAW,UAAU,IAAI,OAAO,CAA6B;IAC7D,OAAO;IAIP,+CAA+C;WACjC,UAAU,IAAI,aAAa;IAGzC,uDAAuD;IAChD,KAAK,CAAC,MAAM,CAAC,EAAE,aAAa,GAAG,aAAa;IAQnD,gEAAgE;IACzD,OAAO;IAId,mEAAmE;IAC5D,wBAAwB,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,MAAM,CAAC,EAAE,OAAO,GAAG,OAAO;IAShF,qEAAqE;IAC9D,0BAA0B,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,MAAM,CAAC,EAAE,OAAO,GAAG,OAAO;IAOlF,mEAAmE;IAC5D,sBAAsB,CAAC,KAAK,EAAE,OAAO,EAAE,MAAM,CAAC,EAAE,OAAO,GAAG,OAAO;IAGxE,qEAAqE;IAC9D,wBAAwB,CAAC,KAAK,EAAE,OAAO,EAAE,MAAM,CAAC,EAAE,OAAO,GAAG,OAAO;IAG1E,iCAAiC;IAC1B,cAAc,CAAC,aAAa,EAAE,MAAM;IAM3C;;;;OAIG;IACI,oCAAoC,CAAC,QAAQ,EAAE,eAAe,EAAE,UAAU,EAAE,MAAM,EAAE,QAAQ,EAAE,MAAM,GAAG,OAAO;CA0DtH"}

package/lib/esm/polyface/FacetFaceData.js

@@ -0,0 +1,142 @@
+/*---------------------------------------------------------------------------------------------
+* Copyright (c) Bentley Systems, Incorporated. All rights reserved.
+* See LICENSE.md in the project root for license terms and full copyright notice.
+*--------------------------------------------------------------------------------------------*/
+/** @packageDocumentation
+ * @module Polyface
+ */
+import { Point2d, Vector2d } from "../geometry3d/Point2dVector2d";
+import { Point3d } from "../geometry3d/Point3dVector3d";
+import { Range2d } from "../geometry3d/Range";
+import { IndexedPolyfaceVisitor } from "./IndexedPolyfaceVisitor";
+/**
+ * Data for a face in a polyface containing facets.
+ * This is built up cooperatively by the PolyfaceBuilder and its
+ * callers, and stored as a FaceData array in PolyfaceData.
+ * @public
+ */
+export class FacetFaceData {
+    constructor(distanceRange, paramRange) {
+        this._paramDistanceRange = distanceRange;
+        this._paramRange = paramRange;
+    }
+    /** (property accessor) Return a reference to the distance-scaled parameter range. */
+    get paramDistanceRange() { return this._paramDistanceRange; }
+    /** (property accessor) Return a reference to the parameter range. */
+    get paramRange() { return this._paramRange; }
+    /** Create a FacetFaceData with null ranges. */
+    static createNull() {
+        return new FacetFaceData(Range2d.createNull(), Range2d.createNull());
+    }
+    /** Create a deep copy of this FacetFaceData object. */
+    clone(result) {
+        if (result) {
+            this._paramDistanceRange.clone(result._paramDistanceRange);
+            this._paramRange.clone(result._paramRange);
+            return result;
+        }
+        return new FacetFaceData(this._paramDistanceRange.clone(), this._paramRange.clone());
+    }
+    /** Restore this FacetFaceData to its null constructor state. */
+    setNull() {
+        this._paramDistanceRange.setNull();
+        this._paramRange.setNull();
+    }
+    /** Return distance-based parameter from stored parameter value. */
+    convertParamXYToDistance(x, y, result) {
+        result = result ? result : Point2d.create();
+        const paramDelta = this._paramRange.high.minus(this._paramRange.low);
+        result.x = (0 === paramDelta.x) ? x : (this._paramDistanceRange.low.x + (x - this._paramRange.low.x)
+            * (this._paramDistanceRange.high.x - this._paramDistanceRange.low.x) / paramDelta.x);
+        result.y = (0.0 === paramDelta.y) ? y : (this.paramDistanceRange.low.y + (y - this._paramRange.low.y)
+            * (this._paramDistanceRange.high.y - this._paramDistanceRange.low.y) / paramDelta.y);
+        return result;
+    }
+    /** Return normalized (0-1) parameter from stored parameter value. */
+    convertParamXYToNormalized(x, y, result) {
+        result = result ? result : Point2d.create();
+        const paramDelta = this._paramRange.high.minus(this._paramRange.low);
+        result.x = (0.0 === paramDelta.x) ? x : ((x - this._paramRange.low.x) / paramDelta.x);
+        result.y = (0.0 === paramDelta.y) ? y : ((y - this._paramRange.low.y) / paramDelta.y);
+        return result;
+    }
+    /** Return distance-based parameter from stored parameter value. */
+    convertParamToDistance(param, result) {
+        return this.convertParamXYToDistance(param.x, param.y, result);
+    }
+    /** Return normalized (0-1) parameter from stored parameter value. */
+    convertParamToNormalized(param, result) {
+        return this.convertParamXYToNormalized(param.x, param.y, result);
+    }
+    /** Scale distance parameters. */
+    scaleDistances(distanceScale) {
+        this._paramDistanceRange.low.x *= distanceScale;
+        this._paramDistanceRange.low.y *= distanceScale;
+        this._paramDistanceRange.high.x *= distanceScale;
+        this._paramDistanceRange.high.y *= distanceScale;
+    }
+    /**
+     * Sets the param and paramDistance range of this FacetFaceData based on the newly terminated facets that make it up.
+     * Takes the polyface itself, the first and last indexes of the facets to be included in the face.
+     * Returns true on success, false otherwise.
+     */
+    setParamDistanceRangeFromNewFaceData(polyface, facetStart, facetEnd) {
+        const dSTotal = Point2d.create();
+        const dSSquaredTotal = Point2d.create();
+        this.setNull();
+        let aveTotal = 0;
+        const visitor = IndexedPolyfaceVisitor.create(polyface, 0);
+        if (!visitor.moveToReadIndex(facetStart) || facetEnd <= facetStart)
+            return false;
+        do {
+            const numPointsInFacet = visitor.numEdgesThisFacet;
+            const visitorPoints = visitor.point;
+            const trianglePointIndexes = [];
+            const visitorParams = visitor.param;
+            const triangleParamIndexes = [];
+            if (!visitorParams)
+                return false;
+            visitorParams.extendRange(this._paramRange);
+            const dUV0 = Vector2d.create();
+            const dUV1 = Vector2d.create();
+            for (let k = 0; k < numPointsInFacet; k++) {
+                trianglePointIndexes[2] = k;
+                triangleParamIndexes[2] = k;
+                if (k > 1) {
+                    visitorParams.vectorIndexIndex(triangleParamIndexes[1], triangleParamIndexes[0], dUV0);
+                    visitorParams.vectorIndexIndex(triangleParamIndexes[1], triangleParamIndexes[2], dUV1);
+                    const delta0 = visitorPoints.getPoint3dAtUncheckedPointIndex(trianglePointIndexes[0]).minus(visitorPoints.getPoint3dAtUncheckedPointIndex(trianglePointIndexes[1]));
+                    const delta1 = visitorPoints.getPoint3dAtUncheckedPointIndex(trianglePointIndexes[1]).minus(visitorPoints.getPoint3dAtUncheckedPointIndex(trianglePointIndexes[2]));
+                    const uvCross = Math.abs(dUV0.x * dUV1.y - dUV1.x * dUV0.y);
+                    if (uvCross) {
+                        const dwDu = Point3d.createFrom(delta0);
+                        dwDu.scaleInPlace(dUV1.y);
+                        dwDu.addScaledInPlace(delta1, -dUV0.y);
+                        const dwDv = Point3d.createFrom(delta1);
+                        dwDv.scaleInPlace(dUV0.x);
+                        dwDv.addScaledInPlace(delta0, -dUV1.x);
+                        const dS = Point2d.create(dwDu.magnitude() / uvCross, dwDv.magnitude() / uvCross);
+                        dSTotal.x += dS.x;
+                        dSTotal.y += dS.y;
+                        dSSquaredTotal.x += dS.x * dS.x;
+                        dSSquaredTotal.y += dS.y * dS.y;
+                        aveTotal++;
+                    }
+                }
+                triangleParamIndexes[0] = triangleParamIndexes[1];
+                triangleParamIndexes[1] = triangleParamIndexes[2];
+                trianglePointIndexes[0] = trianglePointIndexes[1];
+                trianglePointIndexes[1] = trianglePointIndexes[2];
+            }
+        } while (visitor.moveToNextFacet() && visitor.currentReadIndex() < facetEnd);
+        if (aveTotal !== 0) {
+            const dS = Point2d.create(dSTotal.x / aveTotal, dSTotal.y / aveTotal);
+            const standardDeviation = Point2d.create(Math.sqrt(Math.abs((dSSquaredTotal.x / aveTotal) - dS.x * dS.x)), Math.sqrt(Math.abs((dSSquaredTotal.y / aveTotal) - dS.y * dS.y)));
+            // TR# 268980 - Add standard deviation to match QV....
+            this._paramDistanceRange.low.set(0, 0);
+            this._paramDistanceRange.high.set((dS.x + standardDeviation.x) * (this._paramRange.high.x - this._paramRange.low.x), (dS.y + standardDeviation.y) * (this._paramRange.high.y - this._paramRange.low.y));
+        }
+        return true;
+    }
+}
+//# sourceMappingURL=FacetFaceData.js.map

package/lib/esm/polyface/FacetFaceData.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"FacetFaceData.js","sourceRoot":"","sources":["../../../src/polyface/FacetFaceData.ts"],"names":[],"mappings":"AAAA;;;+FAG+F;AAC/F;;GAEG;AAEH,OAAO,EAAE,OAAO,EAAE,QAAQ,EAAE,MAAM,+BAA+B,CAAC;AAClE,OAAO,EAAE,OAAO,EAAE,MAAM,+BAA+B,CAAC;AACxD,OAAO,EAAE,OAAO,EAAE,MAAM,qBAAqB,CAAC;AAE9C,OAAO,EAAE,sBAAsB,EAAE,MAAM,0BAA0B,CAAC;AAElE;;;;;GAKG;AACH,MAAM,OAAO,aAAa;IAOxB,YAAoB,aAAsB,EAAE,UAAmB;QAC7D,IAAI,CAAC,mBAAmB,GAAG,aAAa,CAAC;QACzC,IAAI,CAAC,WAAW,GAAG,UAAU,CAAC;IAChC,CAAC;IAPD,qFAAqF;IACrF,IAAW,kBAAkB,KAAc,OAAO,IAAI,CAAC,mBAAmB,CAAC,CAAC,CAAC;IAC7E,qEAAqE;IACrE,IAAW,UAAU,KAAc,OAAO,IAAI,CAAC,WAAW,CAAC,CAAC,CAAC;IAK7D,+CAA+C;IACxC,MAAM,CAAC,UAAU;QACtB,OAAO,IAAI,aAAa,CAAC,OAAO,CAAC,UAAU,EAAE,EAAE,OAAO,CAAC,UAAU,EAAE,CAAC,CAAC;IACvE,CAAC;IACD,uDAAuD;IAChD,KAAK,CAAC,MAAsB;QACjC,IAAI,MAAM,EAAE;YACV,IAAI,CAAC,mBAAmB,CAAC,KAAK,CAAC,MAAM,CAAC,mBAAmB,CAAC,CAAC;YAC3D,IAAI,CAAC,WAAW,CAAC,KAAK,CAAC,MAAM,CAAC,WAAW,CAAC,CAAC;YAC3C,OAAO,MAAM,CAAC;SACf;QACD,OAAO,IAAI,aAAa,CAAC,IAAI,CAAC,mBAAmB,CAAC,KAAK,EAAE,EAAE,IAAI,CAAC,WAAW,CAAC,KAAK,EAAE,CAAC,CAAC;IACvF,CAAC;IACD,gEAAgE;IACzD,OAAO;QACZ,IAAI,CAAC,mBAAmB,CAAC,OAAO,EAAE,CAAC;QACnC,IAAI,CAAC,WAAW,CAAC,OAAO,EAAE,CAAC;IAC7B,CAAC;IACD,mEAAmE;IAC5D,wBAAwB,CAAC,CAAS,EAAE,CAAS,EAAE,MAAgB;QACpE,MAAM,GAAG,MAAM,CAAC,CAAC,CAAC,MAAM,CAAC,CAAC,CAAC,OAAO,CAAC,MAAM,EAAE,CAAC;QAC5C,MAAM,UAAU,GAAG,IAAI,CAAC,WAAW,CAAC,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,WAAW,CAAC,GAAG,CAAC,CAAC;QACrE,MAAM,CAAC,CAAC,GAAG,CAAC,CAAC,KAAK,UAAU,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,mBAAmB,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,WAAW,CAAC,GAAG,CAAC,CAAC,CAAC;cAChG,CAAC,IAAI,CAAC,mBAAmB,CAAC,IAAI,CAAC,CAAC,GAAG,IAAI,CAAC,mBAAmB,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,UAAU,CAAC,CAAC,CAAC,CAAC;QACvF,MAAM,CAAC,CAAC,GAAG,CAAC,GAAG,KAAK,UAAU,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,kBAAkB,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,WAAW,CAAC,GAAG,CAAC,CAAC,CAAC;cACjG,CAAC,IAAI,CAAC,mBAAmB,CAAC,IAAI,CAAC,CAAC,GAAG,IAAI,CAAC,mBAAmB,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,UAAU,CAAC,CAAC,CAAC,CAAC;QACvF,OAAO,MAAM,CAAC;IAChB,CAAC;IACD,qEAAqE;IAC9D,0BAA0B,CAAC,CAAS,EAAE,CAAS,EAAE,MAAgB;QACtE,MAAM,GAAG,MAAM,CAAC,CAAC,CAAC,MAAM,CAAC,CAAC,CAAC,OAAO,CAAC,MAAM,EAAE,CAAC;QAC5C,MAAM,UAAU,GAAG,IAAI,CAAC,WAAW,CAAC,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,WAAW,CAAC,GAAG,CAAC,CAAC;QACrE,MAAM,CAAC,CAAC,GAAG,CAAC,GAAG,KAAK,UAAU,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,WAAW,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,UAAU,CAAC,CAAC,CAAC,CAAC;QACtF,MAAM,CAAC,CAAC,GAAG,CAAC,GAAG,KAAK,UAAU,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,WAAW,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,UAAU,CAAC,CAAC,CAAC,CAAC;QACtF,OAAO,MAAM,CAAC;IAChB,CAAC;IACD,mEAAmE;IAC5D,sBAAsB,CAAC,KAAc,EAAE,MAAgB;QAC5D,OAAO,IAAI,CAAC,wBAAwB,CAAC,KAAK,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC;IACjE,CAAC;IACD,qEAAqE;IAC9D,wBAAwB,CAAC,KAAc,EAAE,MAAgB;QAC9D,OAAO,IAAI,CAAC,0BAA0B,CAAC,KAAK,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC;IACnE,CAAC;IACD,iCAAiC;IAC1B,cAAc,CAAC,aAAqB;QACzC,IAAI,CAAC,mBAAmB,CAAC,GAAG,CAAC,CAAC,IAAI,aAAa,CAAC;QAChD,IAAI,CAAC,mBAAmB,CAAC,GAAG,CAAC,CAAC,IAAI,aAAa,CAAC;QAChD,IAAI,CAAC,mBAAmB,CAAC,IAAI,CAAC,CAAC,IAAI,aAAa,CAAC;QACjD,IAAI,CAAC,mBAAmB,CAAC,IAAI,CAAC,CAAC,IAAI,aAAa,CAAC;IACnD,CAAC;IACD;;;;OAIG;IACI,oCAAoC,CAAC,QAAyB,EAAE,UAAkB,EAAE,QAAgB;QACzG,MAAM,OAAO,GAAG,OAAO,CAAC,MAAM,EAAE,CAAC;QACjC,MAAM,cAAc,GAAG,OAAO,CAAC,MAAM,EAAE,CAAC;QACxC,IAAI,CAAC,OAAO,EAAE,CAAC;QACf,IAAI,QAAQ,GAAG,CAAC,CAAC;QACjB,MAAM,OAAO,GAAG,sBAAsB,CAAC,MAAM,CAAC,QAAQ,EAAE,CAAC,CAAC,CAAC;QAC3D,IAAI,CAAC,OAAO,CAAC,eAAe,CAAC,UAAU,CAAC,IAAI,QAAQ,IAAI,UAAU;YAChE,OAAO,KAAK,CAAC;QACf,GAAG;YACD,MAAM,gBAAgB,GAAG,OAAO,CAAC,iBAAiB,CAAC;YACnD,MAAM,aAAa,GAAG,OAAO,CAAC,KAAK,CAAC;YACpC,MAAM,oBAAoB,GAAa,EAAE,CAAC;YAC1C,MAAM,aAAa,GAAG,OAAO,CAAC,KAAK,CAAC;YACpC,MAAM,oBAAoB,GAAa,EAAE,CAAC;YAC1C,IAAI,CAAC,aAAa;gBAChB,OAAO,KAAK,CAAC;YACf,aAAa,CAAC,WAAW,CAAC,IAAI,CAAC,WAAW,CAAC,CAAC;YAC5C,MAAM,IAAI,GAAG,QAAQ,CAAC,MAAM,EAAE,CAAC;YAC/B,MAAM,IAAI,GAAG,QAAQ,CAAC,MAAM,EAAE,CAAC;YAC/B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,gBAAgB,EAAE,CAAC,EAAE,EAAE;gBACzC,oBAAoB,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;gBAC5B,oBAAoB,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;gBAC5B,IAAI,CAAC,GAAG,CAAC,EAAE;oBACT,aAAa,CAAC,gBAAgB,CAAC,oBAAoB,CAAC,CAAC,CAAC,EAAE,oBAAoB,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC;oBACvF,aAAa,CAAC,gBAAgB,CAAC,oBAAoB,CAAC,CAAC,CAAC,EAAE,oBAAoB,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC;oBACvF,MAAM,MAAM,GAAG,aAAa,CAAC,+BAA+B,CAAC,oBAAoB,CAAC,CAAC,CAAC,CAAC,CAAC,KAAK,CAAC,aAAa,CAAC,+BAA+B,CAAC,oBAAoB,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;oBACpK,MAAM,MAAM,GAAG,aAAa,CAAC,+BAA+B,CAAC,oBAAoB,CAAC,CAAC,CAAC,CAAC,CAAC,KAAK,CAAC,aAAa,CAAC,+BAA+B,CAAC,oBAAoB,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;oBACpK,MAAM,OAAO,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC;oBAC5D,IAAI,OAAO,EAAE;wBACX,MAAM,IAAI,GAAG,OAAO,CAAC,UAAU,CAAC,MAAM,CAAC,CAAC;wBACxC,IAAI,CAAC,YAAY,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;wBAC1B,IAAI,CAAC,gBAAgB,CAAC,MAAM,EAAE,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;wBACvC,MAAM,IAAI,GAAG,OAAO,CAAC,UAAU,CAAC,MAAM,CAAC,CAAC;wBACxC,IAAI,CAAC,YAAY,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;wBAC1B,IAAI,CAAC,gBAAgB,CAAC,MAAM,EAAE,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;wBACvC,MAAM,EAAE,GAAG,OAAO,CAAC,MAAM,CAAC,IAAI,CAAC,SAAS,EAAE,GAAG,OAAO,EAAE,IAAI,CAAC,SAAS,EAAE,GAAG,OAAO,CAAC,CAAC;wBAClF,OAAO,CAAC,CAAC,IAAI,EAAE,CAAC,CAAC,CAAC;wBAClB,OAAO,CAAC,CAAC,IAAI,EAAE,CAAC,CAAC,CAAC;wBAClB,cAAc,CAAC,CAAC,IAAI,EAAE,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC;wBAChC,cAAc,CAAC,CAAC,IAAI,EAAE,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC;wBAChC,QAAQ,EAAE,CAAC;qBACZ;iBACF;gBACD,oBAAoB,CAAC,CAAC,CAAC,GAAG,oBAAoB,CAAC,CAAC,CAAC,CAAC;gBAClD,oBAAoB,CAAC,CAAC,CAAC,GAAG,oBAAoB,CAAC,CAAC,CAAC,CAAC;gBAClD,oBAAoB,CAAC,CAAC,CAAC,GAAG,oBAAoB,CAAC,CAAC,CAAC,CAAC;gBAClD,oBAAoB,CAAC,CAAC,CAAC,GAAG,oBAAoB,CAAC,CAAC,CAAC,CAAC;aACnD;SACF,QAAQ,OAAO,CAAC,eAAe,EAAE,IAAI,OAAO,CAAC,gBAAgB,EAAE,GAAG,QAAQ,EAAE;QAC7E,IAAI,QAAQ,KAAK,CAAC,EAAE;YAClB,MAAM,EAAE,GAAG,OAAO,CAAC,MAAM,CAAC,OAAO,CAAC,CAAC,GAAG,QAAQ,EAAE,OAAO,CAAC,CAAC,GAAG,QAAQ,CAAC,CAAC;YACtE,MAAM,iBAAiB,GAAG,OAAO,CAAC,MAAM,CAAC,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,cAAc,CAAC,CAAC,GAAG,QAAQ,CAAC,GAAG,EAAE,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,cAAc,CAAC,CAAC,GAAG,QAAQ,CAAC,GAAG,EAAE,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;YAC7K,sDAAsD;YACtD,IAAI,CAAC,mBAAmB,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;YACvC,IAAI,CAAC,mBAAmB,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,iBAAiB,CAAC,CAAC,CAAC,GAAG,CAAC,IAAI,CAAC,WAAW,CAAC,IAAI,CAAC,CAAC,GAAG,IAAI,CAAC,WAAW,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,GAAG,iBAAiB,CAAC,CAAC,CAAC,GAAG,CAAC,IAAI,CAAC,WAAW,CAAC,IAAI,CAAC,CAAC,GAAG,IAAI,CAAC,WAAW,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;SACzM;QACD,OAAO,IAAI,CAAC;IACd,CAAC;CACF","sourcesContent":["/*---------------------------------------------------------------------------------------------\r\n* Copyright (c) Bentley Systems, Incorporated. All rights reserved.\r\n* See LICENSE.md in the project root for license terms and full copyright notice.\r\n*--------------------------------------------------------------------------------------------*/\r\n/** @packageDocumentation\r\n * @module Polyface\r\n */\r\n\r\nimport { Point2d, Vector2d } from \"../geometry3d/Point2dVector2d\";\r\nimport { Point3d } from \"../geometry3d/Point3dVector3d\";\r\nimport { Range2d } from \"../geometry3d/Range\";\r\nimport { IndexedPolyface } from \"./Polyface\";\r\nimport { IndexedPolyfaceVisitor } from \"./IndexedPolyfaceVisitor\";\r\n\r\n/**\r\n * Data for a face in a polyface containing facets.\r\n * This is built up cooperatively by the PolyfaceBuilder and its\r\n * callers, and stored as a FaceData array in PolyfaceData.\r\n * @public\r\n */\r\nexport class FacetFaceData {\r\n  private _paramDistanceRange: Range2d;\r\n  private _paramRange: Range2d;\r\n  /** (property accessor) Return a reference to the distance-scaled parameter range. */\r\n  public get paramDistanceRange(): Range2d { return this._paramDistanceRange; }\r\n  /** (property accessor) Return a reference to the parameter range. */\r\n  public get paramRange(): Range2d { return this._paramRange; }\r\n  private constructor(distanceRange: Range2d, paramRange: Range2d) {\r\n    this._paramDistanceRange = distanceRange;\r\n    this._paramRange = paramRange;\r\n  }\r\n  /** Create a FacetFaceData with null ranges. */\r\n  public static createNull(): FacetFaceData {\r\n    return new FacetFaceData(Range2d.createNull(), Range2d.createNull());\r\n  }\r\n  /** Create a deep copy of this FacetFaceData object. */\r\n  public clone(result?: FacetFaceData): FacetFaceData {\r\n    if (result) {\r\n      this._paramDistanceRange.clone(result._paramDistanceRange);\r\n      this._paramRange.clone(result._paramRange);\r\n      return result;\r\n    }\r\n    return new FacetFaceData(this._paramDistanceRange.clone(), this._paramRange.clone());\r\n  }\r\n  /** Restore this FacetFaceData to its null constructor state. */\r\n  public setNull() {\r\n    this._paramDistanceRange.setNull();\r\n    this._paramRange.setNull();\r\n  }\r\n  /** Return distance-based parameter from stored parameter value. */\r\n  public convertParamXYToDistance(x: number, y: number, result?: Point2d): Point2d {\r\n    result = result ? result : Point2d.create();\r\n    const paramDelta = this._paramRange.high.minus(this._paramRange.low);\r\n    result.x = (0 === paramDelta.x) ? x : (this._paramDistanceRange.low.x + (x - this._paramRange.low.x)\r\n      * (this._paramDistanceRange.high.x - this._paramDistanceRange.low.x) / paramDelta.x);\r\n    result.y = (0.0 === paramDelta.y) ? y : (this.paramDistanceRange.low.y + (y - this._paramRange.low.y)\r\n      * (this._paramDistanceRange.high.y - this._paramDistanceRange.low.y) / paramDelta.y);\r\n    return result;\r\n  }\r\n  /** Return normalized (0-1) parameter from stored parameter value. */\r\n  public convertParamXYToNormalized(x: number, y: number, result?: Point2d): Point2d {\r\n    result = result ? result : Point2d.create();\r\n    const paramDelta = this._paramRange.high.minus(this._paramRange.low);\r\n    result.x = (0.0 === paramDelta.x) ? x : ((x - this._paramRange.low.x) / paramDelta.x);\r\n    result.y = (0.0 === paramDelta.y) ? y : ((y - this._paramRange.low.y) / paramDelta.y);\r\n    return result;\r\n  }\r\n  /** Return distance-based parameter from stored parameter value. */\r\n  public convertParamToDistance(param: Point2d, result?: Point2d): Point2d {\r\n    return this.convertParamXYToDistance(param.x, param.y, result);\r\n  }\r\n  /** Return normalized (0-1) parameter from stored parameter value. */\r\n  public convertParamToNormalized(param: Point2d, result?: Point2d): Point2d {\r\n    return this.convertParamXYToNormalized(param.x, param.y, result);\r\n  }\r\n  /** Scale distance parameters. */\r\n  public scaleDistances(distanceScale: number) {\r\n    this._paramDistanceRange.low.x *= distanceScale;\r\n    this._paramDistanceRange.low.y *= distanceScale;\r\n    this._paramDistanceRange.high.x *= distanceScale;\r\n    this._paramDistanceRange.high.y *= distanceScale;\r\n  }\r\n  /**\r\n   * Sets the param and paramDistance range of this FacetFaceData based on the newly terminated facets that make it up.\r\n   * Takes the polyface itself, the first and last indexes of the facets to be included in the face.\r\n   * Returns true on success, false otherwise.\r\n   */\r\n  public setParamDistanceRangeFromNewFaceData(polyface: IndexedPolyface, facetStart: number, facetEnd: number): boolean {\r\n    const dSTotal = Point2d.create();\r\n    const dSSquaredTotal = Point2d.create();\r\n    this.setNull();\r\n    let aveTotal = 0;\r\n    const visitor = IndexedPolyfaceVisitor.create(polyface, 0);\r\n    if (!visitor.moveToReadIndex(facetStart) || facetEnd <= facetStart)\r\n      return false;\r\n    do {\r\n      const numPointsInFacet = visitor.numEdgesThisFacet;\r\n      const visitorPoints = visitor.point;\r\n      const trianglePointIndexes: number[] = [];\r\n      const visitorParams = visitor.param;\r\n      const triangleParamIndexes: number[] = [];\r\n      if (!visitorParams)\r\n        return false;\r\n      visitorParams.extendRange(this._paramRange);\r\n      const dUV0 = Vector2d.create();\r\n      const dUV1 = Vector2d.create();\r\n      for (let k = 0; k < numPointsInFacet; k++) {\r\n        trianglePointIndexes[2] = k;\r\n        triangleParamIndexes[2] = k;\r\n        if (k > 1) {\r\n          visitorParams.vectorIndexIndex(triangleParamIndexes[1], triangleParamIndexes[0], dUV0);\r\n          visitorParams.vectorIndexIndex(triangleParamIndexes[1], triangleParamIndexes[2], dUV1);\r\n          const delta0 = visitorPoints.getPoint3dAtUncheckedPointIndex(trianglePointIndexes[0]).minus(visitorPoints.getPoint3dAtUncheckedPointIndex(trianglePointIndexes[1]));\r\n          const delta1 = visitorPoints.getPoint3dAtUncheckedPointIndex(trianglePointIndexes[1]).minus(visitorPoints.getPoint3dAtUncheckedPointIndex(trianglePointIndexes[2]));\r\n          const uvCross = Math.abs(dUV0.x * dUV1.y - dUV1.x * dUV0.y);\r\n          if (uvCross) {\r\n            const dwDu = Point3d.createFrom(delta0);\r\n            dwDu.scaleInPlace(dUV1.y);\r\n            dwDu.addScaledInPlace(delta1, -dUV0.y);\r\n            const dwDv = Point3d.createFrom(delta1);\r\n            dwDv.scaleInPlace(dUV0.x);\r\n            dwDv.addScaledInPlace(delta0, -dUV1.x);\r\n            const dS = Point2d.create(dwDu.magnitude() / uvCross, dwDv.magnitude() / uvCross);\r\n            dSTotal.x += dS.x;\r\n            dSTotal.y += dS.y;\r\n            dSSquaredTotal.x += dS.x * dS.x;\r\n            dSSquaredTotal.y += dS.y * dS.y;\r\n            aveTotal++;\r\n          }\r\n        }\r\n        triangleParamIndexes[0] = triangleParamIndexes[1];\r\n        triangleParamIndexes[1] = triangleParamIndexes[2];\r\n        trianglePointIndexes[0] = trianglePointIndexes[1];\r\n        trianglePointIndexes[1] = trianglePointIndexes[2];\r\n      }\r\n    } while (visitor.moveToNextFacet() && visitor.currentReadIndex() < facetEnd);\r\n    if (aveTotal !== 0) {\r\n      const dS = Point2d.create(dSTotal.x / aveTotal, dSTotal.y / aveTotal);\r\n      const standardDeviation = Point2d.create(Math.sqrt(Math.abs((dSSquaredTotal.x / aveTotal) - dS.x * dS.x)), Math.sqrt(Math.abs((dSSquaredTotal.y / aveTotal) - dS.y * dS.y)));\r\n      // TR# 268980 - Add standard deviation to match QV....\r\n      this._paramDistanceRange.low.set(0, 0);\r\n      this._paramDistanceRange.high.set((dS.x + standardDeviation.x) * (this._paramRange.high.x - this._paramRange.low.x), (dS.y + standardDeviation.y) * (this._paramRange.high.y - this._paramRange.low.y));\r\n    }\r\n    return true;\r\n  }\r\n}\r\n"]}

package/lib/esm/polyface/FacetOrientation.d.ts

@@ -0,0 +1,39 @@
+import { IndexedPolyface } from "./Polyface";
+export declare class FacetOrientationFixup {
+    private _edges;
+    private _edgeToPartnerEdge;
+    private _edgeToEdgeInComponent;
+    private _facetToFirstEdgeInComponent;
+    private _facetOrientation;
+    private _components;
+    private _visitor;
+    private _mesh;
+    private constructor();
+    private edgeIdToFacetOrientation;
+    /**
+     * RETURN FALSE IF ANY EDGE HAS 3 ORE MORE FACETS
+     */
+    private setupUnoriented;
+    private recordFacetInComponent;
+    private initializeComponent;
+    private _workArray;
+    private pushFacetEdgesOnStack;
+    /**
+     * * Run flood search from every possible seed, assigning positive and negative orientations
+     * * Halt and return false if flood detects Klein bottle effects.
+     * @return true if flood
+     */
+    private doFlood;
+    private doFacetReversals;
+    static doFixup(mesh: IndexedPolyface): boolean;
+    /** swap entries at indices in a number array.
+     * * indices are not checked for validity.
+     */
+    private static swapEntries;
+    /**
+     *
+     * @param data an array of cyclically linked loops.
+     */
+    private static extractCyclicIndices;
+}
+//# sourceMappingURL=FacetOrientation.d.ts.map

package/lib/esm/polyface/FacetOrientation.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"FacetOrientation.d.ts","sourceRoot":"","sources":["../../../src/polyface/FacetOrientation.ts"],"names":[],"mappings":"AAQA,OAAO,EAAE,eAAe,EAAmB,MAAM,YAAY,CAAC;AAmB9D,qBAAa,qBAAqB;IAChC,OAAO,CAAC,MAAM,CAAqB;IACnC,OAAO,CAAC,kBAAkB,CAAW;IACrC,OAAO,CAAC,sBAAsB,CAAW;IACzC,OAAO,CAAC,4BAA4B,CAAW;IAC/C,OAAO,CAAC,iBAAiB,CAAW;IACpC,OAAO,CAAC,WAAW,CAA0B;IAC7C,OAAO,CAAC,QAAQ,CAAkB;IAClC,OAAO,CAAC,KAAK,CAAkB;IAC/B,OAAO;IAWP,OAAO,CAAC,wBAAwB;IAIhC;;OAEG;IACH,OAAO,CAAC,eAAe;IA4CvB,OAAO,CAAC,sBAAsB;IAM9B,OAAO,CAAC,mBAAmB;IAK3B,OAAO,CAAC,UAAU,CAAgB;IAElC,OAAO,CAAC,qBAAqB;IAM7B;;;;OAIG;IACH,OAAO,CAAC,OAAO;IA2Cf,OAAO,CAAC,gBAAgB;WAWV,OAAO,CAAC,IAAI,EAAE,eAAe,GAAG,OAAO;IASrD;;OAEG;IACH,OAAO,CAAC,MAAM,CAAC,WAAW;IAK1B;;;OAGG;IACH,OAAO,CAAC,MAAM,CAAC,oBAAoB;CAQpC"}

package/lib/esm/polyface/FacetOrientation.js

@@ -0,0 +1,192 @@
+/*---------------------------------------------------------------------------------------------
+* Copyright (c) Bentley Systems, Incorporated. All rights reserved.
+* See LICENSE.md in the project root for license terms and full copyright notice.
+*--------------------------------------------------------------------------------------------*/
+/** @packageDocumentation
+ * @module Polyface
+ */
+import { SortableEdge } from "./IndexedEdgeMatcher";
+import { PolyfaceQuery } from "./PolyfaceQuery";
+class OrientedComponentData {
+    constructor(firstEdgeIndex) {
+        this.numPositive = this.numNegative = 0;
+        this.firstEdgeIndex = firstEdgeIndex;
+    }
+    // announce a new facet orientation.
+    recordOrientation(orientation) {
+        if (orientation > 0)
+            this.numPositive++;
+        else if (orientation < 0)
+            this.numNegative++;
+    }
+}
+export class FacetOrientationFixup {
+    constructor(mesh) {
+        this._workArray = [];
+        this._visitor = mesh.createVisitor(1);
+        this._edges = PolyfaceQuery.createIndexedEdges(this._visitor);
+        this._edgeToPartnerEdge = [];
+        this._edgeToEdgeInComponent = [];
+        this._facetToFirstEdgeInComponent = [];
+        this._facetOrientation = [];
+        this._components = [];
+        this._mesh = mesh;
+    }
+    edgeIdToFacetOrientation(edgeId) {
+        const facetIndex = this._edges.edges[edgeId].facetIndex;
+        return this._facetOrientation[facetIndex];
+    }
+    /**
+     * RETURN FALSE IF ANY EDGE HAS 3 ORE MORE FACETS
+     */
+    setupUnoriented() {
+        this._edges.sort();
+        const edgeArray = this._edges.edges;
+        // make each edge a singleton with itself as successor . .
+        let maxFacetIndex = -1;
+        const numEdges = this._edges.edges.length;
+        for (let i = 0; i < numEdges; i++) {
+            const facetIndex = this._edges.edges[i].facetIndex;
+            if (facetIndex > maxFacetIndex)
+                maxFacetIndex = facetIndex;
+            this._edgeToEdgeInComponent.push(i);
+            this._edgeToPartnerEdge.push(i);
+        }
+        for (let i = 0; i <= maxFacetIndex; i++) {
+            this._facetToFirstEdgeInComponent.push(-1);
+            this._facetOrientation.push(0);
+        }
+        // make each edge part of a component for only its own facet
+        // (This component exists only ephemerally -- it is not referenced by true component)
+        // (swapping entries in the cyclic loops merges loops.)
+        for (let i = 0; i < numEdges; i++) {
+            const facetIndex = this._edges.edges[i].facetIndex;
+            const j = this._facetToFirstEdgeInComponent[facetIndex];
+            if (j === -1)
+                this._facetToFirstEdgeInComponent[facetIndex] = i;
+            else {
+                FacetOrientationFixup.swapEntries(this._edgeToEdgeInComponent, i, j);
+            }
+        }
+        // make contiguous partner edges cycles in edgeToPartnerEdge.
+        // edge i is initially always a singleton.
+        for (let edgeIndex0 = 0; edgeIndex0 < numEdges; edgeIndex0++) {
+            let edgeIndex1 = edgeIndex0 + 1;
+            while (edgeIndex1 < numEdges && SortableEdge.areUndirectedPartners(edgeArray[edgeIndex0], edgeArray[edgeIndex1])) {
+                // splice the loops ...
+                FacetOrientationFixup.swapEntries(this._edgeToPartnerEdge, edgeIndex0, edgeIndex1);
+                edgeIndex1++;
+            }
+            // BUT .. everything else will fail if more than 2 anywhere .....
+            if (edgeIndex1 > edgeIndex0 + 2)
+                return false;
+        }
+        return true;
+    }
+    recordFacetInComponent(facetIndex, orientation) {
+        const componentData = this._components[this._components.length - 1];
+        this._facetOrientation[facetIndex] = orientation;
+        componentData.recordOrientation(orientation);
+    }
+    initializeComponent(edgeIndex) {
+        const facetIndex = this._edges.edges[edgeIndex].facetIndex;
+        this._components.push(new OrientedComponentData(edgeIndex));
+        this.recordFacetInComponent(facetIndex, 1);
+    }
+    pushFacetEdgesOnStack(seedEdge, stack) {
+        FacetOrientationFixup.extractCyclicIndices(this._edgeToEdgeInComponent, seedEdge, this._workArray);
+        for (const edgeIndex of this._workArray) {
+            stack.push(edgeIndex);
+        }
+    }
+    /**
+     * * Run flood search from every possible seed, assigning positive and negative orientations
+     * * Halt and return false if flood detects Klein bottle effects.
+     * @return true if flood
+     */
+    doFlood() {
+        // Finally time for flood search.
+        const edgeStack = [];
+        const edgeArray = this._edges.edges;
+        const numEdges = edgeArray.length;
+        const neighborEdges = [];
+        for (let seedEdgeIndex = 0; seedEdgeIndex < numEdges; seedEdgeIndex++) {
+            if (this.edgeIdToFacetOrientation(seedEdgeIndex) === 0) {
+                edgeStack.length = 0;
+                this.initializeComponent(seedEdgeIndex);
+                this.pushFacetEdgesOnStack(seedEdgeIndex, edgeStack);
+                // An edge on the stack is part of the current component.
+                // Its partners may or may not have been visited.
+                let baseEdgeIndex;
+                while (undefined !== (baseEdgeIndex = edgeStack.pop())) {
+                    const baseFacet = edgeArray[baseEdgeIndex].facetIndex;
+                    const baseOrientation = this._facetOrientation[baseFacet];
+                    FacetOrientationFixup.extractCyclicIndices(this._edgeToPartnerEdge, baseEdgeIndex, neighborEdges);
+                    for (const neighborEdgeIndex of neighborEdges) {
+                        if (neighborEdgeIndex !== baseEdgeIndex) {
+                            const neighborFacet = edgeArray[neighborEdgeIndex].facetIndex;
+                            const neighborOrientation = this._facetOrientation[neighborFacet];
+                            if (neighborOrientation === 0) {
+                                // first visit to this facet !
+                                // orientations of baseEdge and neighborEdge tell us how to orient it.
+                                const newOrientation = SortableEdge.areDirectedPartners(edgeArray[baseEdgeIndex], edgeArray[neighborEdgeIndex]) ? baseOrientation : -baseOrientation;
+                                this.recordFacetInComponent(neighborFacet, newOrientation);
+                                this.pushFacetEdgesOnStack(neighborEdgeIndex, edgeStack);
+                                FacetOrientationFixup.swapEntries(this._edgeToEdgeInComponent, baseEdgeIndex, neighborEdgeIndex);
+                            }
+                            else {
+                                // looking across to an already-visited facet ..
+                                const edgeOrientation = SortableEdge.relativeOrientation(edgeArray[baseEdgeIndex], edgeArray[neighborEdgeIndex]);
+                                if (edgeOrientation * baseOrientation * neighborOrientation > 0)
+                                    return false;
+                            }
+                        }
+                    }
+                }
+            }
+        }
+        // fall through when floods have completed in all components, with no Klein bottle effects
+        return true;
+    }
+    doFacetReversals() {
+        let numReverse = 0;
+        for (this._visitor.reset(); this._visitor.moveToNextFacet();) {
+            const facetIndex = this._visitor.currentReadIndex();
+            if (this._facetOrientation[facetIndex] < 0) {
+                numReverse++;
+                this._mesh.reverseSingleFacet(facetIndex);
+            }
+        }
+        return numReverse;
+    }
+    static doFixup(mesh) {
+        const context = new FacetOrientationFixup(mesh);
+        if (!context.setupUnoriented())
+            return false;
+        const ok = context.doFlood();
+        if (ok)
+            context.doFacetReversals();
+        return ok;
+    }
+    /** swap entries at indices in a number array.
+     * * indices are not checked for validity.
+     */
+    static swapEntries(data, i, j) {
+        const q = data[i];
+        data[i] = data[j];
+        data[j] = q;
+    }
+    /**
+     *
+     * @param data an array of cyclically linked loops.
+     */
+    static extractCyclicIndices(data, index0, loopIndices) {
+        loopIndices.length = 0;
+        let i = index0;
+        do {
+            loopIndices.push(i);
+            i = data[i];
+        } while (i !== index0);
+    }
+}
+//# sourceMappingURL=FacetOrientation.js.map

package/lib/esm/polyface/FacetOrientation.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"FacetOrientation.js","sourceRoot":"","sources":["../../../src/polyface/FacetOrientation.ts"],"names":[],"mappings":"AAAA;;;+FAG+F;AAC/F;;GAEG;AACH,OAAO,EAAsB,YAAY,EAAE,MAAM,sBAAsB,CAAC;AAExE,OAAO,EAAE,aAAa,EAAE,MAAM,iBAAiB,CAAC;AAEhD,MAAM,qBAAqB;IAIzB,YAAmB,cAAsB;QACvC,IAAI,CAAC,WAAW,GAAG,IAAI,CAAC,WAAW,GAAG,CAAC,CAAC;QACxC,IAAI,CAAC,cAAc,GAAG,cAAc,CAAC;IACvC,CAAC;IACD,oCAAoC;IAC7B,iBAAiB,CAAC,WAAmB;QAC1C,IAAI,WAAW,GAAG,CAAC;YACjB,IAAI,CAAC,WAAW,EAAE,CAAC;aAChB,IAAI,WAAW,GAAG,CAAC;YACtB,IAAI,CAAC,WAAW,EAAE,CAAC;IACvB,CAAC;CACF;AACD,MAAM,OAAO,qBAAqB;IAShC,YAAoB,IAAqB;QAyEjC,eAAU,GAAa,EAAE,CAAC;QAxEhC,IAAI,CAAC,QAAQ,GAAG,IAAI,CAAC,aAAa,CAAC,CAAC,CAAC,CAAC;QACtC,IAAI,CAAC,MAAM,GAAG,aAAa,CAAC,kBAAkB,CAAC,IAAI,CAAC,QAAQ,CAAC,CAAC;QAC9D,IAAI,CAAC,kBAAkB,GAAG,EAAE,CAAC;QAC7B,IAAI,CAAC,sBAAsB,GAAG,EAAE,CAAC;QACjC,IAAI,CAAC,4BAA4B,GAAG,EAAE,CAAC;QACvC,IAAI,CAAC,iBAAiB,GAAG,EAAE,CAAC;QAC5B,IAAI,CAAC,WAAW,GAAG,EAAE,CAAC;QACtB,IAAI,CAAC,KAAK,GAAG,IAAI,CAAC;IACpB,CAAC;IAEO,wBAAwB,CAAC,MAAc;QAC7C,MAAM,UAAU,GAAG,IAAI,CAAC,MAAM,CAAC,KAAK,CAAC,MAAM,CAAC,CAAC,UAAU,CAAC;QACxD,OAAO,IAAI,CAAC,iBAAiB,CAAC,UAAU,CAAC,CAAC;IAC5C,CAAC;IACD;;OAEG;IACK,eAAe;QACrB,IAAI,CAAC,MAAM,CAAC,IAAI,EAAE,CAAC;QACnB,MAAM,SAAS,GAAG,IAAI,CAAC,MAAM,CAAC,KAAK,CAAC;QACpC,0DAA0D;QAC1D,IAAI,aAAa,GAAG,CAAC,CAAC,CAAC;QACvB,MAAM,QAAQ,GAAG,IAAI,CAAC,MAAM,CAAC,KAAK,CAAC,MAAM,CAAC;QAC1C,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,QAAQ,EAAE,CAAC,EAAE,EAAE;YACjC,MAAM,UAAU,GAAG,IAAI,CAAC,MAAM,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,UAAU,CAAC;YACnD,IAAI,UAAU,GAAG,aAAa;gBAC5B,aAAa,GAAG,UAAU,CAAC;YAC7B,IAAI,CAAC,sBAAsB,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;YACpC,IAAI,CAAC,kBAAkB,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;SACjC;QACD,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,aAAa,EAAE,CAAC,EAAE,EAAE;YACvC,IAAI,CAAC,4BAA4B,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC;YAC3C,IAAI,CAAC,iBAAiB,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;SAChC;QACD,4DAA4D;QAC5D,qFAAqF;QACrF,uDAAuD;QACvD,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,QAAQ,EAAE,CAAC,EAAE,EAAE;YACjC,MAAM,UAAU,GAAG,IAAI,CAAC,MAAM,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,UAAU,CAAC;YACnD,MAAM,CAAC,GAAG,IAAI,CAAC,4BAA4B,CAAC,UAAU,CAAC,CAAC;YACxD,IAAI,CAAC,KAAK,CAAC,CAAC;gBACV,IAAI,CAAC,4BAA4B,CAAC,UAAU,CAAC,GAAG,CAAC,CAAC;iBAC/C;gBACH,qBAAqB,CAAC,WAAW,CAAC,IAAI,CAAC,sBAAsB,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;aACtE;SACF;QACD,6DAA6D;QAC7D,0CAA0C;QAC1C,KAAK,IAAI,UAAU,GAAG,CAAC,EAAE,UAAU,GAAG,QAAQ,EAAE,UAAU,EAAE,EAAE;YAC5D,IAAI,UAAU,GAAG,UAAU,GAAG,CAAC,CAAC;YAChC,OAAO,UAAU,GAAG,QAAQ,IAAI,YAAY,CAAC,qBAAqB,CAAC,SAAS,CAAC,UAAU,CAAC,EAAE,SAAS,CAAC,UAAU,CAAC,CAAC,EAAE;gBAChH,uBAAuB;gBACvB,qBAAqB,CAAC,WAAW,CAAC,IAAI,CAAC,kBAAkB,EAAE,UAAU,EAAE,UAAU,CAAC,CAAC;gBACnF,UAAU,EAAE,CAAC;aACd;YACD,iEAAiE;YACjE,IAAI,UAAU,GAAG,UAAU,GAAG,CAAC;gBAC7B,OAAO,KAAK,CAAC;SAChB;QACD,OAAO,IAAI,CAAC;IACd,CAAC;IACO,sBAAsB,CAAC,UAAkB,EAAE,WAAmB;QACpE,MAAM,aAAa,GAAG,IAAI,CAAC,WAAW,CAAC,IAAI,CAAC,WAAW,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC;QACpE,IAAI,CAAC,iBAAiB,CAAC,UAAU,CAAC,GAAG,WAAW,CAAC;QACjD,aAAa,CAAC,iBAAiB,CAAC,WAAW,CAAC,CAAC;IAC/C,CAAC;IAEO,mBAAmB,CAAC,SAAiB;QAC3C,MAAM,UAAU,GAAG,IAAI,CAAC,MAAM,CAAC,KAAK,CAAC,SAAS,CAAC,CAAC,UAAU,CAAC;QAC3D,IAAI,CAAC,WAAW,CAAC,IAAI,CAAC,IAAI,qBAAqB,CAAC,SAAS,CAAC,CAAC,CAAC;QAC5D,IAAI,CAAC,sBAAsB,CAAC,UAAU,EAAE,CAAC,CAAC,CAAC;IAC7C,CAAC;IAGO,qBAAqB,CAAC,QAAgB,EAAE,KAAe;QAC7D,qBAAqB,CAAC,oBAAoB,CAAC,IAAI,CAAC,sBAAsB,EAAE,QAAQ,EAAE,IAAI,CAAC,UAAU,CAAC,CAAC;QACnG,KAAK,MAAM,SAAS,IAAI,IAAI,CAAC,UAAU,EAAE;YACvC,KAAK,CAAC,IAAI,CAAC,SAAS,CAAC,CAAC;SACvB;IACH,CAAC;IACD;;;;OAIG;IACK,OAAO;QACb,iCAAiC;QACjC,MAAM,SAAS,GAAa,EAAE,CAAC;QAC/B,MAAM,SAAS,GAAG,IAAI,CAAC,MAAM,CAAC,KAAK,CAAC;QACpC,MAAM,QAAQ,GAAG,SAAS,CAAC,MAAM,CAAC;QAClC,MAAM,aAAa,GAAa,EAAE,CAAC;QACnC,KAAK,IAAI,aAAa,GAAG,CAAC,EAAE,aAAa,GAAG,QAAQ,EAAE,aAAa,EAAE,EAAE;YACrE,IAAI,IAAI,CAAC,wBAAwB,CAAC,aAAa,CAAC,KAAK,CAAC,EAAE;gBACtD,SAAS,CAAC,MAAM,GAAG,CAAC,CAAC;gBACrB,IAAI,CAAC,mBAAmB,CAAC,aAAa,CAAC,CAAC;gBACxC,IAAI,CAAC,qBAAqB,CAAC,aAAa,EAAE,SAAS,CAAC,CAAC;gBACrD,yDAAyD;gBACzD,iDAAiD;gBACjD,IAAI,aAAiC,CAAC;gBACtC,OAAO,SAAS,KAAK,CAAC,aAAa,GAAG,SAAS,CAAC,GAAG,EAAE,CAAC,EAAE;oBACtD,MAAM,SAAS,GAAG,SAAS,CAAC,aAAa,CAAC,CAAC,UAAU,CAAC;oBACtD,MAAM,eAAe,GAAG,IAAI,CAAC,iBAAiB,CAAC,SAAS,CAAC,CAAC;oBAC1D,qBAAqB,CAAC,oBAAoB,CAAC,IAAI,CAAC,kBAAkB,EAAE,aAAa,EAAE,aAAa,CAAC,CAAC;oBAClG,KAAK,MAAM,iBAAiB,IAAI,aAAa,EAAE;wBAC7C,IAAI,iBAAiB,KAAK,aAAa,EAAE;4BACvC,MAAM,aAAa,GAAG,SAAS,CAAC,iBAAiB,CAAC,CAAC,UAAU,CAAC;4BAC9D,MAAM,mBAAmB,GAAG,IAAI,CAAC,iBAAiB,CAAC,aAAa,CAAC,CAAC;4BAClE,IAAI,mBAAmB,KAAK,CAAC,EAAE;gCAC7B,8BAA8B;gCAC9B,sEAAsE;gCACtE,MAAM,cAAc,GAAG,YAAY,CAAC,mBAAmB,CAAC,SAAS,CAAC,aAAa,CAAC,EAAE,SAAS,CAAC,iBAAiB,CAAC,CAAC,CAAC,CAAC,CAAC,eAAe,CAAC,CAAC,CAAC,CAAC,eAAe,CAAC;gCACrJ,IAAI,CAAC,sBAAsB,CAAC,aAAa,EAAE,cAAc,CAAC,CAAC;gCAC3D,IAAI,CAAC,qBAAqB,CAAC,iBAAiB,EAAE,SAAS,CAAC,CAAC;gCACzD,qBAAqB,CAAC,WAAW,CAAC,IAAI,CAAC,sBAAsB,EAAE,aAAa,EAAE,iBAAiB,CAAC,CAAC;6BAClG;iCAAM;gCACL,gDAAgD;gCAChD,MAAM,eAAe,GAAG,YAAY,CAAC,mBAAmB,CAAC,SAAS,CAAC,aAAa,CAAC,EAAE,SAAS,CAAC,iBAAiB,CAAC,CAAC,CAAC;gCACjH,IAAI,eAAe,GAAG,eAAe,GAAG,mBAAmB,GAAG,CAAC;oCAC7D,OAAO,KAAK,CAAC;6BAChB;yBACF;qBACF;iBACF;aACF;SACF;QACD,0FAA0F;QAC1F,OAAO,IAAI,CAAC;IACd,CAAC;IACO,gBAAgB;QACtB,IAAI,UAAU,GAAG,CAAC,CAAC;QACnB,KAAK,IAAI,CAAC,QAAQ,CAAC,KAAK,EAAE,EAAE,IAAI,CAAC,QAAQ,CAAC,eAAe,EAAE,GAAG;YAC5D,MAAM,UAAU,GAAG,IAAI,CAAC,QAAQ,CAAC,gBAAgB,EAAE,CAAC;YACpD,IAAI,IAAI,CAAC,iBAAiB,CAAC,UAAU,CAAC,GAAG,CAAC,EAAE;gBAC1C,UAAU,EAAE,CAAC;gBACb,IAAI,CAAC,KAAK,CAAC,kBAAkB,CAAC,UAAU,CAAC,CAAC;aAC3C;SACF;QACD,OAAO,UAAU,CAAC;IACpB,CAAC;IACM,MAAM,CAAC,OAAO,CAAC,IAAqB;QACzC,MAAM,OAAO,GAAG,IAAI,qBAAqB,CAAC,IAAI,CAAC,CAAC;QAChD,IAAI,CAAC,OAAO,CAAC,eAAe,EAAE;YAC5B,OAAO,KAAK,CAAC;QACf,MAAM,EAAE,GAAG,OAAO,CAAC,OAAO,EAAE,CAAC;QAC7B,IAAI,EAAE;YACJ,OAAO,CAAC,gBAAgB,EAAE,CAAC;QAC7B,OAAO,EAAE,CAAC;IACZ,CAAC;IACD;;OAEG;IACK,MAAM,CAAC,WAAW,CAAC,IAAc,EAAE,CAAS,EAAE,CAAS;QAC7D,MAAM,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC;QAClB,IAAI,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC;QAClB,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;IACd,CAAC;IACD;;;OAGG;IACK,MAAM,CAAC,oBAAoB,CAAC,IAAc,EAAE,MAAc,EAAE,WAAqB;QACvF,WAAW,CAAC,MAAM,GAAG,CAAC,CAAC;QACvB,IAAI,CAAC,GAAG,MAAM,CAAC;QACf,GAAG;YACD,WAAW,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;YACpB,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC;SACb,QAAQ,CAAC,KAAK,MAAM,EAAE;IACzB,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/** @packageDocumentation\r\n * @module Polyface\r\n */\r\nimport { IndexedEdgeMatcher, SortableEdge } from \"./IndexedEdgeMatcher\";\r\nimport { IndexedPolyface, PolyfaceVisitor } from \"./Polyface\";\r\nimport { PolyfaceQuery } from \"./PolyfaceQuery\";\r\n\r\nclass OrientedComponentData {\r\n  public numPositive: number;\r\n  public numNegative: number;\r\n  public firstEdgeIndex: number;\r\n  public constructor(firstEdgeIndex: number) {\r\n    this.numPositive = this.numNegative = 0;\r\n    this.firstEdgeIndex = firstEdgeIndex;\r\n  }\r\n  // announce a new facet orientation.\r\n  public recordOrientation(orientation: number) {\r\n    if (orientation > 0)\r\n      this.numPositive++;\r\n    else if (orientation < 0)\r\n      this.numNegative++;\r\n  }\r\n}\r\nexport class FacetOrientationFixup {\r\n  private _edges: IndexedEdgeMatcher;\r\n  private _edgeToPartnerEdge: number[];\r\n  private _edgeToEdgeInComponent: number[];\r\n  private _facetToFirstEdgeInComponent: number[];\r\n  private _facetOrientation: number[];\r\n  private _components: OrientedComponentData[];\r\n  private _visitor: PolyfaceVisitor;\r\n  private _mesh: IndexedPolyface;\r\n  private constructor(mesh: IndexedPolyface) {\r\n    this._visitor = mesh.createVisitor(1);\r\n    this._edges = PolyfaceQuery.createIndexedEdges(this._visitor);\r\n    this._edgeToPartnerEdge = [];\r\n    this._edgeToEdgeInComponent = [];\r\n    this._facetToFirstEdgeInComponent = [];\r\n    this._facetOrientation = [];\r\n    this._components = [];\r\n    this._mesh = mesh;\r\n  }\r\n\r\n  private edgeIdToFacetOrientation(edgeId: number): number {\r\n    const facetIndex = this._edges.edges[edgeId].facetIndex;\r\n    return this._facetOrientation[facetIndex];\r\n  }\r\n  /**\r\n   * RETURN FALSE IF ANY EDGE HAS 3 ORE MORE FACETS\r\n   */\r\n  private setupUnoriented(): boolean {\r\n    this._edges.sort();\r\n    const edgeArray = this._edges.edges;\r\n    // make each edge a singleton with itself as successor . .\r\n    let maxFacetIndex = -1;\r\n    const numEdges = this._edges.edges.length;\r\n    for (let i = 0; i < numEdges; i++) {\r\n      const facetIndex = this._edges.edges[i].facetIndex;\r\n      if (facetIndex > maxFacetIndex)\r\n        maxFacetIndex = facetIndex;\r\n      this._edgeToEdgeInComponent.push(i);\r\n      this._edgeToPartnerEdge.push(i);\r\n    }\r\n    for (let i = 0; i <= maxFacetIndex; i++) {\r\n      this._facetToFirstEdgeInComponent.push(-1);\r\n      this._facetOrientation.push(0);\r\n    }\r\n    // make each edge part of a component for only its own facet\r\n    // (This component exists only ephemerally -- it is not referenced by true component)\r\n    // (swapping entries in the cyclic loops merges loops.)\r\n    for (let i = 0; i < numEdges; i++) {\r\n      const facetIndex = this._edges.edges[i].facetIndex;\r\n      const j = this._facetToFirstEdgeInComponent[facetIndex];\r\n      if (j === -1)\r\n        this._facetToFirstEdgeInComponent[facetIndex] = i;\r\n      else {\r\n        FacetOrientationFixup.swapEntries(this._edgeToEdgeInComponent, i, j);\r\n      }\r\n    }\r\n    // make contiguous partner edges cycles in edgeToPartnerEdge.\r\n    // edge i is initially always a singleton.\r\n    for (let edgeIndex0 = 0; edgeIndex0 < numEdges; edgeIndex0++) {\r\n      let edgeIndex1 = edgeIndex0 + 1;\r\n      while (edgeIndex1 < numEdges && SortableEdge.areUndirectedPartners(edgeArray[edgeIndex0], edgeArray[edgeIndex1])) {\r\n        // splice the loops ...\r\n        FacetOrientationFixup.swapEntries(this._edgeToPartnerEdge, edgeIndex0, edgeIndex1);\r\n        edgeIndex1++;\r\n      }\r\n      // BUT .. everything else will fail if more than 2 anywhere .....\r\n      if (edgeIndex1 > edgeIndex0 + 2)\r\n        return false;\r\n    }\r\n    return true;\r\n  }\r\n  private recordFacetInComponent(facetIndex: number, orientation: number) {\r\n    const componentData = this._components[this._components.length - 1];\r\n    this._facetOrientation[facetIndex] = orientation;\r\n    componentData.recordOrientation(orientation);\r\n  }\r\n\r\n  private initializeComponent(edgeIndex: number) {\r\n    const facetIndex = this._edges.edges[edgeIndex].facetIndex;\r\n    this._components.push(new OrientedComponentData(edgeIndex));\r\n    this.recordFacetInComponent(facetIndex, 1);\r\n  }\r\n  private _workArray: number[] = [];\r\n\r\n  private pushFacetEdgesOnStack(seedEdge: number, stack: number[]) {\r\n    FacetOrientationFixup.extractCyclicIndices(this._edgeToEdgeInComponent, seedEdge, this._workArray);\r\n    for (const edgeIndex of this._workArray) {\r\n      stack.push(edgeIndex);\r\n    }\r\n  }\r\n  /**\r\n   * * Run flood search from every possible seed, assigning positive and negative orientations\r\n   * * Halt and return false if flood detects Klein bottle effects.\r\n   * @return true if flood\r\n   */\r\n  private doFlood(): boolean {\r\n    // Finally time for flood search.\r\n    const edgeStack: number[] = [];\r\n    const edgeArray = this._edges.edges;\r\n    const numEdges = edgeArray.length;\r\n    const neighborEdges: number[] = [];\r\n    for (let seedEdgeIndex = 0; seedEdgeIndex < numEdges; seedEdgeIndex++) {\r\n      if (this.edgeIdToFacetOrientation(seedEdgeIndex) === 0) {\r\n        edgeStack.length = 0;\r\n        this.initializeComponent(seedEdgeIndex);\r\n        this.pushFacetEdgesOnStack(seedEdgeIndex, edgeStack);\r\n        // An edge on the stack is part of the current component.\r\n        // Its partners may or may not have been visited.\r\n        let baseEdgeIndex: number | undefined;\r\n        while (undefined !== (baseEdgeIndex = edgeStack.pop())) {\r\n          const baseFacet = edgeArray[baseEdgeIndex].facetIndex;\r\n          const baseOrientation = this._facetOrientation[baseFacet];\r\n          FacetOrientationFixup.extractCyclicIndices(this._edgeToPartnerEdge, baseEdgeIndex, neighborEdges);\r\n          for (const neighborEdgeIndex of neighborEdges) {\r\n            if (neighborEdgeIndex !== baseEdgeIndex) {\r\n              const neighborFacet = edgeArray[neighborEdgeIndex].facetIndex;\r\n              const neighborOrientation = this._facetOrientation[neighborFacet];\r\n              if (neighborOrientation === 0) {\r\n                // first visit to this facet !\r\n                // orientations of baseEdge and neighborEdge tell us how to orient it.\r\n                const newOrientation = SortableEdge.areDirectedPartners(edgeArray[baseEdgeIndex], edgeArray[neighborEdgeIndex]) ? baseOrientation : -baseOrientation;\r\n                this.recordFacetInComponent(neighborFacet, newOrientation);\r\n                this.pushFacetEdgesOnStack(neighborEdgeIndex, edgeStack);\r\n                FacetOrientationFixup.swapEntries(this._edgeToEdgeInComponent, baseEdgeIndex, neighborEdgeIndex);\r\n              } else {\r\n                // looking across to an already-visited facet ..\r\n                const edgeOrientation = SortableEdge.relativeOrientation(edgeArray[baseEdgeIndex], edgeArray[neighborEdgeIndex]);\r\n                if (edgeOrientation * baseOrientation * neighborOrientation > 0)\r\n                  return false;\r\n              }\r\n            }\r\n          }\r\n        }\r\n      }\r\n    }\r\n    // fall through when floods have completed in all components, with no Klein bottle effects\r\n    return true;\r\n  }\r\n  private doFacetReversals(): number {\r\n    let numReverse = 0;\r\n    for (this._visitor.reset(); this._visitor.moveToNextFacet();) {\r\n      const facetIndex = this._visitor.currentReadIndex();\r\n      if (this._facetOrientation[facetIndex] < 0) {\r\n        numReverse++;\r\n        this._mesh.reverseSingleFacet(facetIndex);\r\n      }\r\n    }\r\n    return numReverse;\r\n  }\r\n  public static doFixup(mesh: IndexedPolyface): boolean {\r\n    const context = new FacetOrientationFixup(mesh);\r\n    if (!context.setupUnoriented())\r\n      return false;\r\n    const ok = context.doFlood();\r\n    if (ok)\r\n      context.doFacetReversals();\r\n    return ok;\r\n  }\r\n  /** swap entries at indices in a number array.\r\n   * * indices are not checked for validity.\r\n   */\r\n  private static swapEntries(data: number[], i: number, j: number) {\r\n    const q = data[i];\r\n    data[i] = data[j];\r\n    data[j] = q;\r\n  }\r\n  /**\r\n   *\r\n   * @param data an array of cyclically linked loops.\r\n   */\r\n  private static extractCyclicIndices(data: number[], index0: number, loopIndices: number[]) {\r\n    loopIndices.length = 0;\r\n    let i = index0;\r\n    do {\r\n      loopIndices.push(i);\r\n      i = data[i];\r\n    } while (i !== index0);\r\n  }\r\n}\r\n"]}