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

package/lib/esm/curve/Query/StrokeCountMap.js

@@ -0,0 +1,138 @@
+/*---------------------------------------------------------------------------------------------
+* Copyright (c) Bentley Systems, Incorporated. All rights reserved.
+* See LICENSE.md in the project root for license terms and full copyright notice.
+*--------------------------------------------------------------------------------------------*/
+/** @packageDocumentation
+ * @module Curve
+ */
+import { Geometry } from "../../Geometry";
+/**
+ * data carrier interface for per-primitive stroke counts and distances used by PolyfaceBuilder.
+ * * For a simple primitive (Line segment or arc) that is stroked with uniform fraction, the members are:
+ *   * `numStroke` = total number of strokes
+ *   * `curveLength` = length of this curve
+ *   * `a0` = external mapped coordinate for fraction 0 on this primitive or component
+ *   * `a1` = external mapped coordinate for fraction 1 on this primitive or component
+ *
+ * * For linestring and bspline curve, those numbers are totals for the overall curve, and breakdown within
+ *     the components (line segments or bezier spans) is recorded on the optional array `componentData[]`
+ *   * Members of the array are annotated with componentIndex within the linestring or bspline curve
+ * @public
+ */
+export class StrokeCountMap {
+    /**
+     * Constructor.  Initialize all fields from arguments.
+     * * Callers that expect to announce numStroke and curveLength for multiple components send an empty componentData array.
+     * * Callers that do not have multiple components send undefined component data.
+     * @param numStroke
+     * @param curveLength
+     * @param a0
+     * @param a1
+     * @param componentData
+     */
+    constructor(numStroke, curveLength, a0, a1, componentData) {
+        this.numStroke = numStroke;
+        this.curveLength = curveLength;
+        this.a0 = a0;
+        this.a1 = a1;
+        this.componentData = componentData;
+    }
+    /**
+     * Create a `StrokeCountMap` with curve primitive and optional componentData array.
+     * @param primitive
+     * @param numStroke
+     * @param curveLength
+     * @param a0
+     * @param a1
+     * @param componentData
+     */
+    static createWithCurvePrimitive(primitive, numStroke, curveLength, a0, a1, componentData) {
+        const result = new StrokeCountMap(numStroke, curveLength, a0, a1, componentData);
+        result.primitive = primitive;
+        return result;
+    }
+    /**
+     * Create a `StrokeCountMap` with `componentIndex` (but no primitive or componentData array)
+     * @param index
+     * @param numStroke
+     * @param curveLength
+     * @param a0
+     * @param a1
+     */
+    static createWithComponentIndex(componentIndex = 0, numStroke = 0, curveLength = 0, a0 = 0, a1 = 0) {
+        const result = new StrokeCountMap(numStroke, curveLength, a0, a1);
+        result.componentIndex = componentIndex;
+        return result;
+    }
+    /**
+     * create a StrokeCountMap, optionally
+     * * (a) use parent a1 as new a0
+     * * (b) attach a (usually empty) array for component counts.
+     * @param parentMap optional map whose a1 becomes a0 in the new map.
+     * @param componentData optional array of component StrokeCountMaps.
+     */
+    static createWithCurvePrimitiveAndOptionalParent(curvePrimitive, parentMap, componentData) {
+        const a0 = parentMap ? parentMap.a1 : 0.0;
+        const result = new StrokeCountMap(0, 0, a0, a0, componentData);
+        result.primitive = curvePrimitive;
+        return result;
+    }
+    /**
+     * Apply stroke count and curve length from a component to a parent map.
+     * If componentData is present, install the new count and length with distance limits
+     * @param parentMap map to be updated.
+     * @param numStroke number of strokes on new child curve
+     * @param curveLength curve length for new child curve.
+     */
+    addToCountAndLength(numStroke, curveLength) {
+        const a2 = this.a1 + curveLength;
+        if (this.componentData) {
+            this.componentData.push(new StrokeCountMap(numStroke, curveLength, this.a1, a2));
+        }
+        this.numStroke += numStroke;
+        this.curveLength += curveLength;
+        this.a1 = a2;
+    }
+    /** return true if `other` has the same component structure as `this`
+     * * testing recurses through corresponding members of componentData arrays.
+     */
+    isCompatibleComponentStructure(other, enforceCounts) {
+        if (enforceCounts && this.numStroke !== other.numStroke)
+            return false;
+        if (this.componentData === undefined && other.componentData === undefined)
+            return true;
+        if (this.componentData && other.componentData) {
+            // both have components. Recurse . . ..
+            if (this.componentData.length !== other.componentData.length)
+                return false;
+            const n = this.componentData.length;
+            for (let i = 0; i < n; i++)
+                if (!this.componentData[i].isCompatibleComponentStructure(other.componentData[i], enforceCounts))
+                    return false;
+            return true;
+        }
+        // one has componentData, the other not.
+        return false;
+    }
+    /**
+     * * clone all data from root.
+     * * clone componentData arrays recursively.
+     */
+    clone() {
+        const a = new StrokeCountMap(this.numStroke, this.curveLength, this.a0, this.a1);
+        if (this.componentData) {
+            a.componentData = [];
+            for (const child of this.componentData)
+                a.componentData.push(child.clone());
+        }
+        return a;
+    }
+    /**
+     * interpolate in the a0,a1 mapping.
+     * @param fraction fractional position between a0 and a1
+     */
+    fractionToA(fraction) {
+        return Geometry.interpolate(this.a0, fraction, this.a1);
+    }
+}
+//# sourceMappingURL=StrokeCountMap.js.map

package/lib/esm/curve/Query/StrokeCountMap.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"StrokeCountMap.js","sourceRoot":"","sources":["../../../../src/curve/Query/StrokeCountMap.ts"],"names":[],"mappings":"AAAA;;;+FAG+F;AAC/F;;GAEG;AAEH,OAAO,EAAE,QAAQ,EAAE,MAAM,gBAAgB,CAAC;AAG1C;;;;;;;;;;;;GAYG;AACH,MAAM,OAAO,cAAc;IAezB;;;;;;;;;OASG;IACH,YAAoB,SAAiB,EAAE,WAAmB,EAAE,EAAU,EAAE,EAAU,EAAE,aAAgC;QAClH,IAAI,CAAC,SAAS,GAAG,SAAS,CAAC;QAC3B,IAAI,CAAC,WAAW,GAAG,WAAW,CAAC;QAC/B,IAAI,CAAC,EAAE,GAAG,EAAE,CAAC;QACb,IAAI,CAAC,EAAE,GAAG,EAAE,CAAC;QACb,IAAI,CAAC,aAAa,GAAG,aAAa,CAAC;IACrC,CAAC;IACD;;;;;;;;OAQG;IACI,MAAM,CAAC,wBAAwB,CAAC,SAAyB,EAAE,SAAiB,EAAE,WAAmB,EAAE,EAAU,EAAE,EAAU,EAAE,aAAgC;QAChK,MAAM,MAAM,GAAG,IAAI,cAAc,CAAC,SAAS,EAAE,WAAW,EAAE,EAAE,EAAE,EAAE,EAAE,aAAa,CAAC,CAAC;QACjF,MAAM,CAAC,SAAS,GAAG,SAAS,CAAC;QAC7B,OAAO,MAAM,CAAC;IAChB,CAAC;IACD;;;;;;;OAOG;IACI,MAAM,CAAC,wBAAwB,CAAC,iBAAyB,CAAC,EAAE,YAAoB,CAAC,EAAE,cAAsB,CAAC,EAAE,KAAa,CAAC,EAAE,KAAa,CAAC;QAC/I,MAAM,MAAM,GAAG,IAAI,cAAc,CAAC,SAAS,EAAE,WAAW,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC;QAClE,MAAM,CAAC,cAAc,GAAG,cAAc,CAAC;QACvC,OAAO,MAAM,CAAC;IAChB,CAAC;IAED;;;;;;OAMG;IACI,MAAM,CAAC,yCAAyC,CAAC,cAA8B,EAAE,SAA0B,EAAE,aAAgC;QAClJ,MAAM,EAAE,GAAG,SAAS,CAAC,CAAC,CAAC,SAAS,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC;QAC1C,MAAM,MAAM,GAAG,IAAI,cAAc,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,aAAa,CAAC,CAAC;QAC/D,MAAM,CAAC,SAAS,GAAG,cAAc,CAAC;QAClC,OAAO,MAAM,CAAC;IAEhB,CAAC;IACD;;;;;;OAMG;IACI,mBAAmB,CAAC,SAAiB,EAAE,WAAmB;QAC/D,MAAM,EAAE,GAAG,IAAI,CAAC,EAAE,GAAG,WAAW,CAAC;QACjC,IAAI,IAAI,CAAC,aAAa,EAAE;YACtB,IAAI,CAAC,aAAa,CAAC,IAAI,CACrB,IAAI,cAAc,CAAC,SAAS,EAAE,WAAW,EAAE,IAAI,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC,CAAC;SAC5D;QACD,IAAI,CAAC,SAAS,IAAI,SAAS,CAAC;QAC5B,IAAI,CAAC,WAAW,IAAI,WAAW,CAAC;QAChC,IAAI,CAAC,EAAE,GAAG,EAAE,CAAC;IACf,CAAC;IACD;;OAEG;IACI,8BAA8B,CAAC,KAAqB,EAAE,aAAsB;QACjF,IAAI,aAAa,IAAI,IAAI,CAAC,SAAS,KAAK,KAAK,CAAC,SAAS;YACrD,OAAO,KAAK,CAAC;QACf,IAAI,IAAI,CAAC,aAAa,KAAK,SAAS,IAAI,KAAK,CAAC,aAAa,KAAK,SAAS;YACvE,OAAO,IAAI,CAAC;QACd,IAAI,IAAI,CAAC,aAAa,IAAI,KAAK,CAAC,aAAa,EAAE;YAC7C,uCAAuC;YACvC,IAAI,IAAI,CAAC,aAAa,CAAC,MAAM,KAAK,KAAK,CAAC,aAAa,CAAC,MAAM;gBAC1D,OAAO,KAAK,CAAC;YACf,MAAM,CAAC,GAAG,IAAI,CAAC,aAAa,CAAC,MAAM,CAAC;YACpC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE;gBACxB,IAAI,CAAC,IAAI,CAAC,aAAa,CAAC,CAAC,CAAC,CAAC,8BAA8B,CAAC,KAAK,CAAC,aAAa,CAAC,CAAC,CAAC,EAAE,aAAa,CAAC;oBAC9F,OAAO,KAAK,CAAC;YACjB,OAAO,IAAI,CAAC;SACb;QACD,wCAAwC;QACxC,OAAO,KAAK,CAAC;IACf,CAAC;IACD;;;OAGG;IACI,KAAK;QACV,MAAM,CAAC,GAAG,IAAI,cAAc,CAAC,IAAI,CAAC,SAAS,EAAE,IAAI,CAAC,WAAW,EAAE,IAAI,CAAC,EAAE,EAAE,IAAI,CAAC,EAAE,CAAC,CAAC;QACjF,IAAI,IAAI,CAAC,aAAa,EAAE;YACtB,CAAC,CAAC,aAAa,GAAG,EAAE,CAAC;YACrB,KAAK,MAAM,KAAK,IAAI,IAAI,CAAC,aAAa;gBACpC,CAAC,CAAC,aAAa,CAAC,IAAI,CAAC,KAAK,CAAC,KAAK,EAAE,CAAC,CAAC;SACvC;QACD,OAAO,CAAC,CAAC;IACX,CAAC;IACD;;;OAGG;IACI,WAAW,CAAC,QAAgB;QACjC,OAAO,QAAQ,CAAC,WAAW,CAAC,IAAI,CAAC,EAAE,EAAE,QAAQ,EAAE,IAAI,CAAC,EAAE,CAAC,CAAC;IAC1D,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 Curve\r\n */\r\n\r\nimport { Geometry } from \"../../Geometry\";\r\nimport { CurvePrimitive } from \"../CurvePrimitive\";\r\n\r\n/**\r\n * data carrier interface for per-primitive stroke counts and distances used by PolyfaceBuilder.\r\n * * For a simple primitive (Line segment or arc) that is stroked with uniform fraction, the members are:\r\n *   * `numStroke` = total number of strokes\r\n *   * `curveLength` = length of this curve\r\n *   * `a0` = external mapped coordinate for fraction 0 on this primitive or component\r\n *   * `a1` = external mapped coordinate for fraction 1 on this primitive or component\r\n *\r\n * * For linestring and bspline curve, those numbers are totals for the overall curve, and breakdown within\r\n *     the components (line segments or bezier spans) is recorded on the optional array `componentData[]`\r\n *   * Members of the array are annotated with componentIndex within the linestring or bspline curve\r\n * @public\r\n */\r\nexport class StrokeCountMap {\r\n  /** number of strokes expected in this interval. */\r\n  public numStroke: number;\r\n  /** Length of the curve interval. */\r\n  public curveLength: number;\r\n  /** start coordinate (in user-defined space) for fraction 0 on this primitive or component */\r\n  public a0: number;\r\n  /** end coordinate (in user-defined space) for fraction 0 on this primitive or component */\r\n  public a1: number;\r\n  /** further StrokeCountMap's within this interval (e.g. for individual segments of a linestring.) */\r\n  public componentData?: StrokeCountMap[];\r\n  /** The curve that this map represents */\r\n  public primitive?: CurvePrimitive;\r\n  /** this curve's index within its parent. */\r\n  public componentIndex?: number;\r\n  /**\r\n   * Constructor.  Initialize all fields from arguments.\r\n   * * Callers that expect to announce numStroke and curveLength for multiple components send an empty componentData array.\r\n   * * Callers that do not have multiple components send undefined component data.\r\n   * @param numStroke\r\n   * @param curveLength\r\n   * @param a0\r\n   * @param a1\r\n   * @param componentData\r\n   */\r\n  private constructor(numStroke: number, curveLength: number, a0: number, a1: number, componentData?: StrokeCountMap[]) {\r\n    this.numStroke = numStroke;\r\n    this.curveLength = curveLength;\r\n    this.a0 = a0;\r\n    this.a1 = a1;\r\n    this.componentData = componentData;\r\n  }\r\n  /**\r\n   * Create a `StrokeCountMap` with curve primitive and optional componentData array.\r\n   * @param primitive\r\n   * @param numStroke\r\n   * @param curveLength\r\n   * @param a0\r\n   * @param a1\r\n   * @param componentData\r\n   */\r\n  public static createWithCurvePrimitive(primitive: CurvePrimitive, numStroke: number, curveLength: number, a0: number, a1: number, componentData?: StrokeCountMap[]) {\r\n    const result = new StrokeCountMap(numStroke, curveLength, a0, a1, componentData);\r\n    result.primitive = primitive;\r\n    return result;\r\n  }\r\n  /**\r\n   * Create a `StrokeCountMap` with `componentIndex` (but no primitive or componentData array)\r\n   * @param index\r\n   * @param numStroke\r\n   * @param curveLength\r\n   * @param a0\r\n   * @param a1\r\n   */\r\n  public static createWithComponentIndex(componentIndex: number = 0, numStroke: number = 0, curveLength: number = 0, a0: number = 0, a1: number = 0) {\r\n    const result = new StrokeCountMap(numStroke, curveLength, a0, a1);\r\n    result.componentIndex = componentIndex;\r\n    return result;\r\n  }\r\n\r\n  /**\r\n   * create a StrokeCountMap, optionally\r\n   * * (a) use parent a1 as new a0\r\n   * * (b) attach a (usually empty) array for component counts.\r\n   * @param parentMap optional map whose a1 becomes a0 in the new map.\r\n   * @param componentData optional array of component StrokeCountMaps.\r\n   */\r\n  public static createWithCurvePrimitiveAndOptionalParent(curvePrimitive: CurvePrimitive, parentMap?: StrokeCountMap, componentData?: StrokeCountMap[]): StrokeCountMap {\r\n    const a0 = parentMap ? parentMap.a1 : 0.0;\r\n    const result = new StrokeCountMap(0, 0, a0, a0, componentData);\r\n    result.primitive = curvePrimitive;\r\n    return result;\r\n\r\n  }\r\n  /**\r\n   * Apply stroke count and curve length from a component to a parent map.\r\n   * If componentData is present, install the new count and length with distance limits\r\n   * @param parentMap map to be updated.\r\n   * @param numStroke number of strokes on new child curve\r\n   * @param curveLength curve length for new child curve.\r\n   */\r\n  public addToCountAndLength(numStroke: number, curveLength: number) {\r\n    const a2 = this.a1 + curveLength;\r\n    if (this.componentData) {\r\n      this.componentData.push(\r\n        new StrokeCountMap(numStroke, curveLength, this.a1, a2));\r\n    }\r\n    this.numStroke += numStroke;\r\n    this.curveLength += curveLength;\r\n    this.a1 = a2;\r\n  }\r\n  /** return true if `other` has the same component structure as `this`\r\n   * * testing recurses through corresponding members of componentData arrays.\r\n   */\r\n  public isCompatibleComponentStructure(other: StrokeCountMap, enforceCounts: boolean): boolean {\r\n    if (enforceCounts && this.numStroke !== other.numStroke)\r\n      return false;\r\n    if (this.componentData === undefined && other.componentData === undefined)\r\n      return true;\r\n    if (this.componentData && other.componentData) {\r\n      // both have components. Recurse . . ..\r\n      if (this.componentData.length !== other.componentData.length)\r\n        return false;\r\n      const n = this.componentData.length;\r\n      for (let i = 0; i < n; i++)\r\n        if (!this.componentData[i].isCompatibleComponentStructure(other.componentData[i], enforceCounts))\r\n          return false;\r\n      return true;\r\n    }\r\n    // one has componentData, the other not.\r\n    return false;\r\n  }\r\n  /**\r\n   * * clone all data from root.\r\n   * * clone componentData arrays recursively.\r\n   */\r\n  public clone(): StrokeCountMap {\r\n    const a = new StrokeCountMap(this.numStroke, this.curveLength, this.a0, this.a1);\r\n    if (this.componentData) {\r\n      a.componentData = [];\r\n      for (const child of this.componentData)\r\n        a.componentData.push(child.clone());\r\n    }\r\n    return a;\r\n  }\r\n  /**\r\n   * interpolate in the a0,a1 mapping.\r\n   * @param fraction fractional position between a0 and a1\r\n   */\r\n  public fractionToA(fraction: number) {\r\n    return Geometry.interpolate(this.a0, fraction, this.a1);\r\n  }\r\n}\r\n"]}

package/lib/esm/curve/RegionMomentsXY.d.ts

@@ -0,0 +1,61 @@
+/** @packageDocumentation
+ * @module Curve
+ */
+import { BSplineCurve3d } from "../bspline/BSplineCurve";
+import { BSplineCurve3dH } from "../bspline/BSplineCurve3dH";
+import { NullGeometryHandler } from "../geometry3d/GeometryHandler";
+import { MomentData } from "../geometry4d/MomentData";
+import { Arc3d } from "./Arc3d";
+import { CurvePrimitive } from "./CurvePrimitive";
+import { LineSegment3d } from "./LineSegment3d";
+import { LineString3d } from "./LineString3d";
+import { Loop } from "./Loop";
+import { ParityRegion } from "./ParityRegion";
+import { TransitionSpiral3d } from "./spiral/TransitionSpiral3d";
+import { UnionRegion } from "./UnionRegion";
+/**
+ * Implementation class for computing XY area moments.
+ * @internal
+ */
+export declare class RegionMomentsXY extends NullGeometryHandler {
+    private _activeMomentData?;
+    private _point0;
+    private _point1;
+    /** Accumulate (independent) integrations over
+     * * origin to chord of the arc.
+     * * origin to the "cap" between the chord and arc.
+     */
+    handleArc3d(arc: Arc3d): void;
+    /** Accumulate integrals over the (triangular) areas from the origin to each line segment */
+    handleLineString3d(ls: LineString3d): void;
+    /** Accumulate integrals over the (triangular) area from the origin to this line segment */
+    handleLineSegment3d(segment: LineSegment3d): void;
+    /** Accumulate integrals from origin to all primitives in the chain. */
+    handleLoop(loop: Loop): MomentData | undefined;
+    /**
+     * ASSUMPTIONS FOR ORIENTATION AND CONTAINMENT ISSUES
+     * * Largest area is outer
+     * * All others are interior (and not overlapping)
+     * Hence
+     * * Outer area sign must be positive -- negate all integrations as needed
+     * * Outer area signs must be positive -- negate all integrations as needed
+     * @param region
+     */
+    handleParityRegion(region: ParityRegion): MomentData | undefined;
+    /** Accumulate (as simple addition) products over each component of the union region. */
+    handleUnionRegion(region: UnionRegion): MomentData | undefined;
+    private _strokeOptions?;
+    private getStrokeOptions;
+    /** Single curve primitive (not loop . . .).
+     * * stroke the curve
+     * * accumulate stroke array.
+     */
+    handleCurvePrimitive(cp: CurvePrimitive): void;
+    /** handle strongly typed  BSplineCurve3d  as generic curve primitive */
+    handleBSplineCurve3d(g: BSplineCurve3d): void;
+    /** handle strongly typed  BSplineCurve3dH  as generic curve primitive */
+    handleBSplineCurve3dH(g: BSplineCurve3dH): void;
+    /** handle strongly typed  TransitionSpiral as generic curve primitive  */
+    handleTransitionSpiral(g: TransitionSpiral3d): void;
+}
+//# sourceMappingURL=RegionMomentsXY.d.ts.map

package/lib/esm/curve/RegionMomentsXY.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"RegionMomentsXY.d.ts","sourceRoot":"","sources":["../../../src/curve/RegionMomentsXY.ts"],"names":[],"mappings":"AAKA;;GAEG;AAIH,OAAO,EAAE,cAAc,EAAE,MAAM,yBAAyB,CAAC;AACzD,OAAO,EAAE,eAAe,EAAE,MAAM,4BAA4B,CAAC;AAE7D,OAAO,EAAE,mBAAmB,EAAE,MAAM,+BAA+B,CAAC;AAEpE,OAAO,EAAE,UAAU,EAAE,MAAM,0BAA0B,CAAC;AACtD,OAAO,EAAE,KAAK,EAAE,MAAM,SAAS,CAAC;AAChC,OAAO,EAAE,cAAc,EAAE,MAAM,kBAAkB,CAAC;AAClD,OAAO,EAAE,aAAa,EAAE,MAAM,iBAAiB,CAAC;AAChD,OAAO,EAAE,YAAY,EAAE,MAAM,gBAAgB,CAAC;AAC9C,OAAO,EAAE,IAAI,EAAE,MAAM,QAAQ,CAAC;AAC9B,OAAO,EAAE,YAAY,EAAE,MAAM,gBAAgB,CAAC;AAE9C,OAAO,EAAE,kBAAkB,EAAE,MAAM,6BAA6B,CAAC;AACjE,OAAO,EAAE,WAAW,EAAE,MAAM,eAAe,CAAC;AAE5C;;;GAGG;AACH,qBAAa,eAAgB,SAAQ,mBAAmB;IACtD,OAAO,CAAC,iBAAiB,CAAC,CAAa;IACvC,OAAO,CAAC,OAAO,CAAoB;IACnC,OAAO,CAAC,OAAO,CAAoB;IAEnC;;;OAGG;IACa,WAAW,CAAC,GAAG,EAAE,KAAK,GAAG,IAAI;IAgC7C,4FAA4F;IAC5E,kBAAkB,CAAC,EAAE,EAAE,YAAY,GAAG,IAAI;IAI1D,2FAA2F;IAC3E,mBAAmB,CAAC,OAAO,EAAE,aAAa,GAAG,IAAI;IAMjE,uEAAuE;IACvD,UAAU,CAAC,IAAI,EAAE,IAAI,GAAG,UAAU,GAAG,SAAS;IAQ9D;;;;;;;;OAQG;IACa,kBAAkB,CAAC,MAAM,EAAE,YAAY,GAAG,UAAU,GAAG,SAAS;IAgChF,wFAAwF;IACxE,iBAAiB,CAAC,MAAM,EAAE,WAAW,GAAG,UAAU,GAAG,SAAS;IAY9E,OAAO,CAAC,cAAc,CAAC,CAAgB;IACvC,OAAO,CAAC,gBAAgB;IASxB;;;OAGG;IACI,oBAAoB,CAAC,EAAE,EAAE,cAAc;IAM9C,wEAAwE;IACxD,oBAAoB,CAAC,CAAC,EAAE,cAAc;IACtD,yEAAyE;IACzD,qBAAqB,CAAC,CAAC,EAAE,eAAe;IACxD,0EAA0E;IAC1D,sBAAsB,CAAC,CAAC,EAAE,kBAAkB;CAE7D"}

package/lib/esm/curve/RegionMomentsXY.js

@@ -0,0 +1,156 @@
+/*---------------------------------------------------------------------------------------------
+* Copyright (c) Bentley Systems, Incorporated. All rights reserved.
+* See LICENSE.md in the project root for license terms and full copyright notice.
+*--------------------------------------------------------------------------------------------*/
+import { Angle } from "../geometry3d/Angle";
+import { NullGeometryHandler } from "../geometry3d/GeometryHandler";
+import { Point3d } from "../geometry3d/Point3dVector3d";
+import { MomentData } from "../geometry4d/MomentData";
+import { LineString3d } from "./LineString3d";
+import { Loop } from "./Loop";
+import { StrokeOptions } from "./StrokeOptions";
+/**
+ * Implementation class for computing XY area moments.
+ * @internal
+ */
+export class RegionMomentsXY extends NullGeometryHandler {
+    constructor() {
+        super(...arguments);
+        this._point0 = Point3d.create();
+        this._point1 = Point3d.create();
+    }
+    /** Accumulate (independent) integrations over
+     * * origin to chord of the arc.
+     * * origin to the "cap" between the chord and arc.
+     */
+    handleArc3d(arc) {
+        const momentData = this._activeMomentData;
+        const sweepRadians = arc.sweep.sweepRadians;
+        const alphaRadians = sweepRadians * 0.5;
+        // from https://apps.dtic.mil/dtic/tr/fulltext/u2/274936.pdf page 71  for radius = 1
+        let s = Math.sin(alphaRadians);
+        let c = Math.cos(alphaRadians);
+        let s1 = Math.sin(sweepRadians);
+        if (Angle.isFullCircleRadians(sweepRadians)) {
+            s = 0.0;
+            c = -1.0;
+            s1 = 0.0;
+        }
+        const q = 2 * s * s * s * c / (alphaRadians - s * c);
+        const s3 = s * s * s;
+        const s6 = s3 * s3;
+        const area = 0.5 * (sweepRadians - s1);
+        const inertiaXX = 0.25 * area * (1.0 - q / 3.0);
+        const inertiaYY1 = 0.25 * area * (1.0 + q);
+        const inertiaYY = inertiaYY1 - 4.0 * s6 / (9.0 * area);
+        const productXX = inertiaYY;
+        const productYY = inertiaXX;
+        const centerToCentroid = 4.0 * s * s * s / (3.0 * (sweepRadians - s1));
+        const midRadians = arc.sweep.fractionToRadians(0.5);
+        const centralPlane = arc.radiansToRotatedBasis(midRadians);
+        const centroid = centralPlane.origin.plusScaled(centralPlane.vectorU, centerToCentroid);
+        momentData.accumulateXYProductsInCentroidalFrame(productXX, 0.0, productYY, area, centroid, centralPlane.vectorU, centralPlane.vectorV);
+        const pointB = arc.fractionToPoint(0.0);
+        const pointC = arc.fractionToPoint(1.0);
+        momentData.accumulateTriangleMomentsXY(undefined, pointB, pointC);
+    }
+    /** Accumulate integrals over the (triangular) areas from the origin to each line segment */
+    handleLineString3d(ls) {
+        const momentData = this._activeMomentData;
+        momentData.accumulateTriangleToLineStringMomentsXY(undefined, ls.packedPoints);
+    }
+    /** Accumulate integrals over the (triangular) area from the origin to this line segment */
+    handleLineSegment3d(segment) {
+        const momentData = this._activeMomentData;
+        segment.startPoint(this._point0);
+        segment.endPoint(this._point1);
+        momentData.accumulateTriangleMomentsXY(undefined, this._point0, this._point1);
+    }
+    /** Accumulate integrals from origin to all primitives in the chain. */
+    handleLoop(loop) {
+        const momentData = this._activeMomentData = MomentData.create();
+        momentData.needOrigin = false;
+        for (const child of loop.children)
+            child.dispatchToGeometryHandler(this);
+        this._activeMomentData = undefined;
+        return momentData;
+    }
+    /**
+     * ASSUMPTIONS FOR ORIENTATION AND CONTAINMENT ISSUES
+     * * Largest area is outer
+     * * All others are interior (and not overlapping)
+     * Hence
+     * * Outer area sign must be positive -- negate all integrations as needed
+     * * Outer area signs must be positive -- negate all integrations as needed
+     * @param region
+     */
+    handleParityRegion(region) {
+        const allChildMoments = [];
+        let maxAbsArea = 0.0;
+        let largestChildMoments;
+        for (const child of region.children) {
+            if (child instanceof Loop) {
+                const childMoments = this.handleLoop(child);
+                if (childMoments) {
+                    allChildMoments.push(childMoments);
+                    const q = Math.abs(childMoments.quantitySum);
+                    if (q > maxAbsArea) {
+                        maxAbsArea = q;
+                        largestChildMoments = childMoments;
+                    }
+                }
+            }
+        }
+        if (largestChildMoments) {
+            const summedMoments = MomentData.create();
+            const sign0 = largestChildMoments.signFactor(1.0);
+            summedMoments.accumulateProducts(largestChildMoments, sign0);
+            for (const childMoments of allChildMoments) {
+                if (childMoments !== largestChildMoments) {
+                    const sign1 = childMoments.signFactor(-1.0);
+                    summedMoments.accumulateProducts(childMoments, sign1);
+                }
+            }
+            return summedMoments;
+        }
+        return undefined;
+    }
+    /** Accumulate (as simple addition) products over each component of the union region. */
+    handleUnionRegion(region) {
+        const summedMoments = MomentData.create();
+        for (const child of region.children) {
+            const childMoments = child.dispatchToGeometryHandler(this);
+            if (childMoments) {
+                const sign0 = childMoments.signFactor(1.0);
+                summedMoments.accumulateProducts(childMoments, sign0);
+            }
+        }
+        return summedMoments;
+    }
+    getStrokeOptions() {
+        if (this._strokeOptions)
+            return this._strokeOptions;
+        const options = StrokeOptions.createForCurves();
+        // this is unusually fine for stroking, but appropriate for sum.
+        options.angleTol = Angle.createDegrees(5.0);
+        this._strokeOptions = options;
+        return options;
+    }
+    /** Single curve primitive (not loop . . .).
+     * * stroke the curve
+     * * accumulate stroke array.
+     */
+    handleCurvePrimitive(cp) {
+        const strokes = LineString3d.create();
+        const options = this.getStrokeOptions();
+        cp.emitStrokes(strokes, options);
+        this.handleLineString3d(strokes);
+    }
+    /** handle strongly typed  BSplineCurve3d  as generic curve primitive */
+    handleBSplineCurve3d(g) { return this.handleCurvePrimitive(g); }
+    /** handle strongly typed  BSplineCurve3dH  as generic curve primitive */
+    handleBSplineCurve3dH(g) { return this.handleCurvePrimitive(g); }
+    /** handle strongly typed  TransitionSpiral as generic curve primitive  */
+    handleTransitionSpiral(g) { return this.handleCurvePrimitive(g); }
+}
+//# sourceMappingURL=RegionMomentsXY.js.map

package/lib/esm/curve/RegionMomentsXY.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"RegionMomentsXY.js","sourceRoot":"","sources":["../../../src/curve/RegionMomentsXY.ts"],"names":[],"mappings":"AAAA;;;+FAG+F;AAU/F,OAAO,EAAE,KAAK,EAAE,MAAM,qBAAqB,CAAC;AAC5C,OAAO,EAAE,mBAAmB,EAAE,MAAM,+BAA+B,CAAC;AACpE,OAAO,EAAE,OAAO,EAAE,MAAM,+BAA+B,CAAC;AACxD,OAAO,EAAE,UAAU,EAAE,MAAM,0BAA0B,CAAC;AAItD,OAAO,EAAE,YAAY,EAAE,MAAM,gBAAgB,CAAC;AAC9C,OAAO,EAAE,IAAI,EAAE,MAAM,QAAQ,CAAC;AAE9B,OAAO,EAAE,aAAa,EAAE,MAAM,iBAAiB,CAAC;AAIhD;;;GAGG;AACH,MAAM,OAAO,eAAgB,SAAQ,mBAAmB;IAAxD;;QAEU,YAAO,GAAG,OAAO,CAAC,MAAM,EAAE,CAAC;QAC3B,YAAO,GAAG,OAAO,CAAC,MAAM,EAAE,CAAC;IA4IrC,CAAC;IA1IC;;;OAGG;IACa,WAAW,CAAC,GAAU;QACpC,MAAM,UAAU,GAAG,IAAI,CAAC,iBAAkB,CAAC;QAC3C,MAAM,YAAY,GAAG,GAAG,CAAC,KAAK,CAAC,YAAY,CAAC;QAC5C,MAAM,YAAY,GAAG,YAAY,GAAG,GAAG,CAAC;QACxC,oFAAoF;QACpF,IAAI,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,YAAY,CAAC,CAAC;QAC/B,IAAI,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,YAAY,CAAC,CAAC;QAC/B,IAAI,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,YAAY,CAAC,CAAC;QAChC,IAAI,KAAK,CAAC,mBAAmB,CAAC,YAAY,CAAC,EAAE;YAC3C,CAAC,GAAG,GAAG,CAAC;YACR,CAAC,GAAG,CAAC,GAAG,CAAC;YACT,EAAE,GAAG,GAAG,CAAC;SACV;QACD,MAAM,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,YAAY,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC;QACrD,MAAM,EAAE,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;QACrB,MAAM,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC;QACnB,MAAM,IAAI,GAAG,GAAG,GAAG,CAAC,YAAY,GAAG,EAAE,CAAC,CAAC;QACvC,MAAM,SAAS,GAAG,IAAI,GAAG,IAAI,GAAG,CAAC,GAAG,GAAG,CAAC,GAAG,GAAG,CAAC,CAAC;QAChD,MAAM,UAAU,GAAG,IAAI,GAAG,IAAI,GAAG,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC;QAC3C,MAAM,SAAS,GAAG,UAAU,GAAG,GAAG,GAAG,EAAE,GAAG,CAAC,GAAG,GAAG,IAAI,CAAC,CAAC;QACvD,MAAM,SAAS,GAAG,SAAS,CAAC;QAC5B,MAAM,SAAS,GAAG,SAAS,CAAC;QAC5B,MAAM,gBAAgB,GAAG,GAAG,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC,YAAY,GAAG,EAAE,CAAC,CAAC,CAAC;QACvE,MAAM,UAAU,GAAG,GAAG,CAAC,KAAK,CAAC,iBAAiB,CAAC,GAAG,CAAC,CAAC;QACpD,MAAM,YAAY,GAAG,GAAG,CAAC,qBAAqB,CAAC,UAAU,CAAC,CAAC;QAC3D,MAAM,QAAQ,GAAG,YAAY,CAAC,MAAM,CAAC,UAAU,CAAC,YAAY,CAAC,OAAO,EAAE,gBAAgB,CAAC,CAAC;QACxF,UAAU,CAAC,qCAAqC,CAAC,SAAS,EAAE,GAAG,EAAE,SAAS,EAAE,IAAI,EAAE,QAAQ,EAAE,YAAY,CAAC,OAAO,EAAE,YAAY,CAAC,OAAO,CAAC,CAAC;QACxI,MAAM,MAAM,GAAG,GAAG,CAAC,eAAe,CAAC,GAAG,CAAC,CAAC;QACxC,MAAM,MAAM,GAAG,GAAG,CAAC,eAAe,CAAC,GAAG,CAAC,CAAC;QACxC,UAAU,CAAC,2BAA2B,CAAC,SAAS,EAAE,MAAM,EAAE,MAAM,CAAC,CAAC;IACpE,CAAC;IAED,4FAA4F;IAC5E,kBAAkB,CAAC,EAAgB;QACjD,MAAM,UAAU,GAAG,IAAI,CAAC,iBAAkB,CAAC;QAC3C,UAAU,CAAC,uCAAuC,CAAC,SAAS,EAAE,EAAE,CAAC,YAAY,CAAC,CAAC;IACjF,CAAC;IACD,2FAA2F;IAC3E,mBAAmB,CAAC,OAAsB;QACxD,MAAM,UAAU,GAAG,IAAI,CAAC,iBAAkB,CAAC;QAC3C,OAAO,CAAC,UAAU,CAAC,IAAI,CAAC,OAAO,CAAC,CAAC;QACjC,OAAO,CAAC,QAAQ,CAAC,IAAI,CAAC,OAAO,CAAC,CAAC;QAC/B,UAAU,CAAC,2BAA2B,CAAC,SAAS,EAAE,IAAI,CAAC,OAAO,EAAE,IAAI,CAAC,OAAO,CAAC,CAAC;IAChF,CAAC;IACD,uEAAuE;IACvD,UAAU,CAAC,IAAU;QACnC,MAAM,UAAU,GAAG,IAAI,CAAC,iBAAiB,GAAG,UAAU,CAAC,MAAM,EAAE,CAAC;QAChE,UAAU,CAAC,UAAU,GAAG,KAAK,CAAC;QAC9B,KAAK,MAAM,KAAK,IAAI,IAAI,CAAC,QAAQ;YAC/B,KAAK,CAAC,yBAAyB,CAAC,IAAI,CAAC,CAAC;QACxC,IAAI,CAAC,iBAAiB,GAAG,SAAS,CAAC;QACnC,OAAO,UAAU,CAAC;IACpB,CAAC;IACD;;;;;;;;OAQG;IACa,kBAAkB,CAAC,MAAoB;QACrD,MAAM,eAAe,GAAiB,EAAE,CAAC;QACzC,IAAI,UAAU,GAAG,GAAG,CAAC;QACrB,IAAI,mBAA2C,CAAC;QAChD,KAAK,MAAM,KAAK,IAAI,MAAM,CAAC,QAAQ,EAAE;YACnC,IAAI,KAAK,YAAY,IAAI,EAAE;gBACzB,MAAM,YAAY,GAAG,IAAI,CAAC,UAAU,CAAC,KAAK,CAAC,CAAC;gBAC5C,IAAI,YAAY,EAAE;oBAChB,eAAe,CAAC,IAAI,CAAC,YAAY,CAAC,CAAC;oBACnC,MAAM,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,YAAY,CAAC,WAAW,CAAC,CAAC;oBAC7C,IAAI,CAAC,GAAG,UAAU,EAAE;wBAClB,UAAU,GAAG,CAAC,CAAC;wBACf,mBAAmB,GAAG,YAAY,CAAC;qBACpC;iBAEF;aACF;SACF;QACD,IAAI,mBAAmB,EAAE;YACvB,MAAM,aAAa,GAAG,UAAU,CAAC,MAAM,EAAE,CAAC;YAC1C,MAAM,KAAK,GAAG,mBAAmB,CAAC,UAAU,CAAC,GAAG,CAAC,CAAC;YAClD,aAAa,CAAC,kBAAkB,CAAC,mBAAmB,EAAE,KAAK,CAAC,CAAC;YAC7D,KAAK,MAAM,YAAY,IAAI,eAAe,EAAE;gBAC1C,IAAI,YAAY,KAAK,mBAAmB,EAAE;oBACxC,MAAM,KAAK,GAAG,YAAY,CAAC,UAAU,CAAC,CAAC,GAAG,CAAC,CAAC;oBAC5C,aAAa,CAAC,kBAAkB,CAAC,YAAY,EAAE,KAAK,CAAC,CAAC;iBACvD;aACF;YACD,OAAO,aAAa,CAAC;SACtB;QACD,OAAO,SAAS,CAAC;IACnB,CAAC;IACD,wFAAwF;IACxE,iBAAiB,CAAC,MAAmB;QACnD,MAAM,aAAa,GAAG,UAAU,CAAC,MAAM,EAAE,CAAC;QAC1C,KAAK,MAAM,KAAK,IAAI,MAAM,CAAC,QAAQ,EAAE;YACnC,MAAM,YAAY,GAAG,KAAK,CAAC,yBAAyB,CAAC,IAAI,CAAC,CAAC;YAC3D,IAAI,YAAY,EAAE;gBAChB,MAAM,KAAK,GAAG,YAAY,CAAC,UAAU,CAAC,GAAG,CAAC,CAAC;gBAC3C,aAAa,CAAC,kBAAkB,CAAC,YAAY,EAAE,KAAK,CAAC,CAAC;aACvD;SACF;QACD,OAAO,aAAa,CAAC;IACvB,CAAC;IAGO,gBAAgB;QACtB,IAAI,IAAI,CAAC,cAAc;YACrB,OAAO,IAAI,CAAC,cAAc,CAAC;QAC7B,MAAM,OAAO,GAAG,aAAa,CAAC,eAAe,EAAE,CAAC;QAChD,gEAAgE;QAChE,OAAO,CAAC,QAAQ,GAAG,KAAK,CAAC,aAAa,CAAC,GAAG,CAAC,CAAC;QAC5C,IAAI,CAAC,cAAc,GAAG,OAAO,CAAC;QAC9B,OAAO,OAAO,CAAC;IACjB,CAAC;IACD;;;OAGG;IACI,oBAAoB,CAAC,EAAkB;QAC5C,MAAM,OAAO,GAAG,YAAY,CAAC,MAAM,EAAE,CAAC;QACtC,MAAM,OAAO,GAAG,IAAI,CAAC,gBAAgB,EAAE,CAAC;QACxC,EAAE,CAAC,WAAW,CAAC,OAAO,EAAE,OAAO,CAAC,CAAC;QACjC,IAAI,CAAC,kBAAkB,CAAC,OAAO,CAAC,CAAC;IACnC,CAAC;IACD,wEAAwE;IACxD,oBAAoB,CAAC,CAAiB,IAAI,OAAO,IAAI,CAAC,oBAAoB,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IAChG,yEAAyE;IACzD,qBAAqB,CAAC,CAAkB,IAAI,OAAO,IAAI,CAAC,oBAAoB,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IAClG,0EAA0E;IAC1D,sBAAsB,CAAC,CAAqB,IAAI,OAAO,IAAI,CAAC,oBAAoB,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;CAEvG","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 Curve\r\n */\r\n\r\n// import { Geometry, Angle, AngleSweep } from \"../Geometry\";\r\n\r\nimport { BSplineCurve3d } from \"../bspline/BSplineCurve\";\r\nimport { BSplineCurve3dH } from \"../bspline/BSplineCurve3dH\";\r\nimport { Angle } from \"../geometry3d/Angle\";\r\nimport { NullGeometryHandler } from \"../geometry3d/GeometryHandler\";\r\nimport { Point3d } from \"../geometry3d/Point3dVector3d\";\r\nimport { MomentData } from \"../geometry4d/MomentData\";\r\nimport { Arc3d } from \"./Arc3d\";\r\nimport { CurvePrimitive } from \"./CurvePrimitive\";\r\nimport { LineSegment3d } from \"./LineSegment3d\";\r\nimport { LineString3d } from \"./LineString3d\";\r\nimport { Loop } from \"./Loop\";\r\nimport { ParityRegion } from \"./ParityRegion\";\r\nimport { StrokeOptions } from \"./StrokeOptions\";\r\nimport { TransitionSpiral3d } from \"./spiral/TransitionSpiral3d\";\r\nimport { UnionRegion } from \"./UnionRegion\";\r\n\r\n/**\r\n * Implementation class for computing XY area moments.\r\n * @internal\r\n */\r\nexport class RegionMomentsXY extends NullGeometryHandler {\r\n  private _activeMomentData?: MomentData;\r\n  private _point0 = Point3d.create();\r\n  private _point1 = Point3d.create();\r\n\r\n  /** Accumulate (independent) integrations over\r\n   * * origin to chord of the arc.\r\n   * * origin to the \"cap\" between the chord and arc.\r\n   */\r\n  public override handleArc3d(arc: Arc3d): void {\r\n    const momentData = this._activeMomentData!;\r\n    const sweepRadians = arc.sweep.sweepRadians;\r\n    const alphaRadians = sweepRadians * 0.5;\r\n    // from https://apps.dtic.mil/dtic/tr/fulltext/u2/274936.pdf page 71  for radius = 1\r\n    let s = Math.sin(alphaRadians);\r\n    let c = Math.cos(alphaRadians);\r\n    let s1 = Math.sin(sweepRadians);\r\n    if (Angle.isFullCircleRadians(sweepRadians)) {\r\n      s = 0.0;\r\n      c = -1.0;\r\n      s1 = 0.0;\r\n    }\r\n    const q = 2 * s * s * s * c / (alphaRadians - s * c);\r\n    const s3 = s * s * s;\r\n    const s6 = s3 * s3;\r\n    const area = 0.5 * (sweepRadians - s1);\r\n    const inertiaXX = 0.25 * area * (1.0 - q / 3.0);\r\n    const inertiaYY1 = 0.25 * area * (1.0 + q);\r\n    const inertiaYY = inertiaYY1 - 4.0 * s6 / (9.0 * area);\r\n    const productXX = inertiaYY;\r\n    const productYY = inertiaXX;\r\n    const centerToCentroid = 4.0 * s * s * s / (3.0 * (sweepRadians - s1));\r\n    const midRadians = arc.sweep.fractionToRadians(0.5);\r\n    const centralPlane = arc.radiansToRotatedBasis(midRadians);\r\n    const centroid = centralPlane.origin.plusScaled(centralPlane.vectorU, centerToCentroid);\r\n    momentData.accumulateXYProductsInCentroidalFrame(productXX, 0.0, productYY, area, centroid, centralPlane.vectorU, centralPlane.vectorV);\r\n    const pointB = arc.fractionToPoint(0.0);\r\n    const pointC = arc.fractionToPoint(1.0);\r\n    momentData.accumulateTriangleMomentsXY(undefined, pointB, pointC);\r\n  }\r\n\r\n  /** Accumulate integrals over the (triangular) areas from the origin to each line segment */\r\n  public override handleLineString3d(ls: LineString3d): void {\r\n    const momentData = this._activeMomentData!;\r\n    momentData.accumulateTriangleToLineStringMomentsXY(undefined, ls.packedPoints);\r\n  }\r\n  /** Accumulate integrals over the (triangular) area from the origin to this line segment */\r\n  public override handleLineSegment3d(segment: LineSegment3d): void {\r\n    const momentData = this._activeMomentData!;\r\n    segment.startPoint(this._point0);\r\n    segment.endPoint(this._point1);\r\n    momentData.accumulateTriangleMomentsXY(undefined, this._point0, this._point1);\r\n  }\r\n  /** Accumulate integrals from origin to all primitives in the chain. */\r\n  public override handleLoop(loop: Loop): MomentData | undefined {\r\n    const momentData = this._activeMomentData = MomentData.create();\r\n    momentData.needOrigin = false;\r\n    for (const child of loop.children)\r\n      child.dispatchToGeometryHandler(this);\r\n    this._activeMomentData = undefined;\r\n    return momentData;\r\n  }\r\n  /**\r\n   * ASSUMPTIONS FOR ORIENTATION AND CONTAINMENT ISSUES\r\n   * * Largest area is outer\r\n   * * All others are interior (and not overlapping)\r\n   * Hence\r\n   * * Outer area sign must be positive -- negate all integrations as needed\r\n   * * Outer area signs must be positive -- negate all integrations as needed\r\n   * @param region\r\n   */\r\n  public override handleParityRegion(region: ParityRegion): MomentData | undefined {\r\n    const allChildMoments: MomentData[] = [];\r\n    let maxAbsArea = 0.0;\r\n    let largestChildMoments: MomentData | undefined;\r\n    for (const child of region.children) {\r\n      if (child instanceof Loop) {\r\n        const childMoments = this.handleLoop(child);\r\n        if (childMoments) {\r\n          allChildMoments.push(childMoments);\r\n          const q = Math.abs(childMoments.quantitySum);\r\n          if (q > maxAbsArea) {\r\n            maxAbsArea = q;\r\n            largestChildMoments = childMoments;\r\n          }\r\n\r\n        }\r\n      }\r\n    }\r\n    if (largestChildMoments) {\r\n      const summedMoments = MomentData.create();\r\n      const sign0 = largestChildMoments.signFactor(1.0);\r\n      summedMoments.accumulateProducts(largestChildMoments, sign0);\r\n      for (const childMoments of allChildMoments) {\r\n        if (childMoments !== largestChildMoments) {\r\n          const sign1 = childMoments.signFactor(-1.0);\r\n          summedMoments.accumulateProducts(childMoments, sign1);\r\n        }\r\n      }\r\n      return summedMoments;\r\n    }\r\n    return undefined;\r\n  }\r\n  /** Accumulate (as simple addition) products over each component of the union region. */\r\n  public override handleUnionRegion(region: UnionRegion): MomentData | undefined {\r\n    const summedMoments = MomentData.create();\r\n    for (const child of region.children) {\r\n      const childMoments = child.dispatchToGeometryHandler(this);\r\n      if (childMoments) {\r\n        const sign0 = childMoments.signFactor(1.0);\r\n        summedMoments.accumulateProducts(childMoments, sign0);\r\n      }\r\n    }\r\n    return summedMoments;\r\n  }\r\n\r\n  private _strokeOptions?: StrokeOptions;\r\n  private getStrokeOptions(): StrokeOptions {\r\n    if (this._strokeOptions)\r\n      return this._strokeOptions;\r\n    const options = StrokeOptions.createForCurves();\r\n    // this is unusually fine for stroking, but appropriate for sum.\r\n    options.angleTol = Angle.createDegrees(5.0);\r\n    this._strokeOptions = options;\r\n    return options;\r\n  }\r\n  /** Single curve primitive (not loop . . .).\r\n   * * stroke the curve\r\n   * * accumulate stroke array.\r\n   */\r\n  public handleCurvePrimitive(cp: CurvePrimitive) {\r\n    const strokes = LineString3d.create();\r\n    const options = this.getStrokeOptions();\r\n    cp.emitStrokes(strokes, options);\r\n    this.handleLineString3d(strokes);\r\n  }\r\n  /** handle strongly typed  BSplineCurve3d  as generic curve primitive */\r\n  public override handleBSplineCurve3d(g: BSplineCurve3d) { return this.handleCurvePrimitive(g); }\r\n  /** handle strongly typed  BSplineCurve3dH  as generic curve primitive */\r\n  public override handleBSplineCurve3dH(g: BSplineCurve3dH) { return this.handleCurvePrimitive(g); }\r\n  /** handle strongly typed  TransitionSpiral as generic curve primitive  */\r\n  public override handleTransitionSpiral(g: TransitionSpiral3d) { return this.handleCurvePrimitive(g); }\r\n\r\n}\r\n"]}