@itwin/core-geometry 4.3.0-dev.8 → 4.3.0

package/lib/cjs/geometry3d/PointHelpers.d.ts

@@ -18,7 +18,7 @@ export declare class NumberArray {
     /** Return true if arrays have identical counts and equal entries (using `!==` comparison) */
     static isExactEqual(dataA: any[] | Float64Array | undefined, dataB: any[] | Float64Array | undefined): boolean;
     /** Return true if arrays have identical counts and entries equal within tolerance */
-    static isAlmostEqual(dataA: number[] | Float64Array | undefined, dataB: number[] | Float64Array | undefined, tolerance: number): boolean;
+    static isAlmostEqual(dataA: number[] | Float64Array | undefined, dataB: number[] | Float64Array | undefined, tolerance?: number): boolean;
     /** return the sum of numbers in an array.  Note that "PreciseSum" may be more accurate. */
     static sum(data: number[] | Float64Array): number;
     /** test if coordinate x appears (to tolerance by `Geometry.isSameCoordinate`) in this array of numbers */
@@ -42,8 +42,23 @@ export declare class NumberArray {
      * @param step max permitted step
      */
     static createArrayWithMaxStepSize(low: number, high: number, step: number): number[];
-    /** copy numbers from variant sources to number[]. */
+    /** Copy numbers from variant sources to number[]. */
     static create(source: number[] | Float64Array): number[];
+    /** Copy number[][]. */
+    static copy2d(source: number[][]): number[][];
+    /** Copy number[][][]. */
+    static copy3d(source: number[][][]): number[][][];
+    /** Copy numbers from Float64Array to number[][].
+     * @param numPerBlock block size
+     */
+    static unpack2d(source: Float64Array, numPerBlock: number): number[][] | undefined;
+    /** Copy numbers from Float64Array to number[][][].
+     * @param numPerRow row size
+     * @param numPerBlock block size
+     */
+    static unpack3d(source: Float64Array, numPerRow: number, numPerBlock: number): number[][][] | undefined;
+    /** Copy numbers from 1d/2d/3d array to Float64Array. */
+    static pack(source: number[] | number[][] | number[][][]): Float64Array;
     /** Return a copy of the knots array, with multiplicity of first and last knots raised or lowered to expectedMultiplicity. */
     static cloneWithStartAndEndMultiplicity(knots: number[] | undefined, target0: number, target1: number): number[];
     /** Compute the linear combination s of the numbers and scales.
@@ -96,13 +111,30 @@ export declare class Vector3dArray {
  * @public
  */
 export declare class Point4dArray {
-    /** pack each point and its corresponding weight into a buffer of xyzw xyzw ... */
+    /**
+     * Copy each weighted point and its corresponding weight into a packed buffer.
+     * @param data array of weighted xyz
+     * @param weights scalar weight array
+     * @param result optional destination array. If insufficiently sized, a new array is returned.
+     * @return packed weighted point array
+     */
     static packPointsAndWeightsToFloat64Array(data: Point3d[] | Float64Array | number[], weights: number[] | Float64Array, result?: Float64Array): Float64Array | undefined;
-    /** pack x,y,z,w in Float64Array. */
+    /**
+     * Copy 4d points into a packed buffer.
+     * @param data array of xyzw
+     * @param result optional destination array. If insufficiently sized, a new array is returned.
+     * @return packed point array
+     */
     static packToFloat64Array(data: Point4d[], result?: Float64Array): Float64Array;
     /** unpack from  ... to array of Point4d */
     static unpackToPoint4dArray(data: Float64Array): Point4d[];
-    /** unpack from xyzw xyzw... array to array of Point3d and array of weight.
+    /**
+     * Unpack packed 4D data to a Point3d array and an array of weights.
+     * * `WeightStyle` of `data` is not assumed. If input data is of form [a,b,c,d], default output arrays will have form [a,b,c] and [d].
+     * @param data input 4D points (packed)
+     * @param points output 3D data
+     * @param weights output weights (w portion of input)
+     * @param pointFormatter optional xyz formatter. By default, returns a Point3d created from the xyz portion of the input.
      */
     static unpackFloat64ArrayToPointsAndWeights(data: Float64Array, points: Point3d[], weights: number[], pointFormatter?: (x: number, y: number, z: number) => any): void;
     private static _workPoint4d;
@@ -112,7 +144,7 @@ export declare class Point4dArray {
      * @param xyzw array of x,y,z,w points.
      */
     static multiplyInPlace(transform: Transform, xyzw: Float64Array): void;
-    /** test for near equality of all corresponding numeric values, treated as coordinates. */
+    /** Test arrays for near equality of all corresponding numeric values, treated as coordinates. */
     static isAlmostEqual(dataA: Point4d[] | Float64Array | undefined, dataB: Point4d[] | Float64Array | undefined): boolean;
     /** return true iff all xyzw points' altitudes are within tolerance of the plane.*/
     static isCloseToPlane(data: Point4d[] | Float64Array, plane: Plane3dByOriginAndUnitNormal, tolerance?: number): boolean;
@@ -122,8 +154,13 @@ export declare class Point4dArray {
  * @public
  */
 export declare class Point3dArray {
-    /** pack x,y,z to `Float64Array` */
-    static packToFloat64Array(data: Point3d[]): Float64Array;
+    /**
+     * Copy 3d points into a packed buffer.
+     * @param data array of xyz
+     * @param result optional destination array. If insufficiently sized, a new array is returned.
+     * @return packed point array
+     */
+    static packToFloat64Array(data: Point3d[], result?: Float64Array): Float64Array;
     /**
      * Compute the 8 weights of trilinear mapping
      * By appropriate choice of weights, this can be used for both point and derivative mappings.
@@ -187,14 +224,15 @@ export declare class Point3dArray {
      */
     static unpackNumbersToNestedArrays(data: Float64Array, numPerBlock: number): any[];
     /**
-     * return an 3-dimensional array containing all the values of `data` in arrays numPerRow blocks of numPerBlock
+     * Return a 3-dimensional array containing all the values of `data` in rows of numPerRow blocks of size numPerBlock.
      * @param data simple array of numbers
-     * @param numPerBlock number of values in each block at first level down
+     * @param numPerBlock number of values in each block
+     * @param numPerRow number of blocks per row
      */
     static unpackNumbersToNestedArraysIJK(data: Float64Array, numPerBlock: number, numPerRow: number): any[];
     /**  multiply a transform times each x,y,z triple and replace the x,y,z in the packed array */
     static multiplyInPlace(transform: Transform, xyz: Float64Array): void;
-    /** Apply Geometry.isAlmostEqual to corresponding coordinates */
+    /** Test arrays for near equality of all corresponding numeric values, treated as coordinates. */
     static isAlmostEqual(dataA: Point3d[] | Float64Array | undefined, dataB: Point3d[] | Float64Array | undefined): boolean;
     /** return simple average of all coordinates.   (000 if empty array) */
     static centroid(points: IndexedXYZCollection | Point3d[], result?: Point3d): Point3d;
@@ -255,18 +293,18 @@ export declare class Point3dArray {
      */
     static cloneDeepXYZPoint3dArrays(data: MultiLineStringDataVariant): any[];
     /**
-     * return perpendicular distance from points[indexB] to the segment points[indexA] to points[indexC].
-     * * extrapolation option when projection is outside of fraction range 0..1 are:
-     *   * false ==> measure distance to closest endpoint
-     *   * true ==> measure distance to extended line segment.
-     * (no index checking!)
+     * Return perpendicular distance from points[indexB] to the segment from points[indexA] to points[indexC].
+     * * Extrapolation options when the projection is outside of the fraction range [0,1] are:
+     *   * false ==> return distance to closest endpoint
+     *   * true ==> return distance to extended line segment
+     * * There is no index checking!
      */
     static distanceIndexedPointBToSegmentAC(points: Point3d[], indexA: number, indexB: number, indexC: number, extrapolate: boolean): number;
     /** Computes the hull of the XY projection of points.
-     * * Returns the hull as an array of Point3d
-     * * Optionally returns non-hull points in `insidePoints[]`
-     * * If both arrays empty if less than 3 points.
-     * *
+     * @param points input points, z-coordinates ignored
+     * @param hullPoints (output) points on the convex hull (cloned from input points)
+     * @param insidePoints (output) points not on the convex hull (cloned from input points)
+     * @param addClosurePoint whether to append the first hull point to `hullPoints`
      */
     static computeConvexHullXY(points: Point3d[], hullPoints: Point3d[], insidePoints: Point3d[], addClosurePoint?: boolean): void;
     /**

package/lib/cjs/geometry3d/PointHelpers.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"PointHelpers.d.ts","sourceRoot":"","sources":["../../../src/geometry3d/PointHelpers.ts"],"names":[],"mappings":"AAWA,OAAO,EAAE,OAAO,EAAE,MAAM,uBAAuB,CAAC;AAChD,OAAO,EAAE,0BAA0B,EAAE,MAAM,2BAA2B,CAAC;AACvE,OAAO,EAAE,oBAAoB,EAAE,MAAM,wBAAwB,CAAC;AAC9D,OAAO,EAAE,4BAA4B,EAAE,MAAM,gCAAgC,CAAC;AAC9E,OAAO,EAAE,OAAO,EAAE,MAAM,mBAAmB,CAAC;AAE5C,OAAO,EAAE,OAAO,EAAE,QAAQ,EAAE,GAAG,EAAE,MAAM,mBAAmB,CAAC;AAE3D,OAAO,EAAE,SAAS,EAAE,MAAM,aAAa,CAAC;AACxC,OAAO,EAAE,KAAK,EAAE,MAAM,EAAE,QAAQ,EAAE,MAAM,YAAY,CAAC;AAiBrD;;;GAGG;AACH,qBAAa,WAAW;IACtB;;OAEG;WACW,UAAU,CAAC,IAAI,EAAE,MAAM,EAAE,GAAG,MAAM;IAgBhD,6FAA6F;WAC/E,YAAY,CAAC,KAAK,EAAE,GAAG,EAAE,GAAG,YAAY,GAAG,SAAS,EAAE,KAAK,EAAE,GAAG,EAAE,GAAG,YAAY,GAAG,SAAS,GAAG,OAAO;IAWrH,qFAAqF;WACvE,aAAa,CACzB,KAAK,EAAE,MAAM,EAAE,GAAG,YAAY,GAAG,SAAS,EAC1C,KAAK,EAAE,MAAM,EAAE,GAAG,YAAY,GAAG,SAAS,EAC1C,SAAS,EAAE,MAAM,GAAG,OAAO;IAW7B,2FAA2F;WAC7E,GAAG,CAAC,IAAI,EAAE,MAAM,EAAE,GAAG,YAAY,GAAG,MAAM;IAKxD,0GAA0G;WAC5F,mBAAmB,CAAC,CAAC,EAAE,MAAM,EAAE,IAAI,EAAE,MAAM,EAAE,GAAG,SAAS,GAAG,OAAO;IAMjF,2DAA2D;WAC7C,WAAW,CAAC,MAAM,EAAE,MAAM,EAAE,GAAG,MAAM;IAcnD,yDAAyD;WAC3C,SAAS,CAAC,EAAE,EAAE,MAAM,EAAE,EAAE,EAAE,MAAM,GAAG,MAAM;IAKvD;;OAEG;WACW,UAAU,CAAC,KAAK,EAAE,MAAM,EAAE,GAAG,YAAY,EAAE,KAAK,EAAE,MAAM,EAAE,GAAG,YAAY,GAAG,MAAM;IAOhG;;OAEG;WACW,iBAAiB,CAAC,KAAK,EAAE,YAAY,EAAE,KAAK,EAAE,YAAY,GAAG,MAAM;IAMjF;;;;;OAKG;WACW,0BAA0B,CAAC,GAAG,EAAE,MAAM,EAAE,IAAI,EAAE,MAAM,EAAE,IAAI,EAAE,MAAM,GAAG,MAAM,EAAE;IAc3F,qDAAqD;WACvC,MAAM,CAAC,MAAM,EAAE,MAAM,EAAE,GAAG,YAAY,GAAG,MAAM,EAAE;IAO/D,6HAA6H;WAC/G,gCAAgC,CAAC,KAAK,EAAE,MAAM,EAAE,GAAG,SAAS,EAAE,OAAO,EAAE,MAAM,EAAE,OAAO,EAAE,MAAM,GAAG,MAAM,EAAE;IAqBvH;;;;OAIG;WACW,iBAAiB,CAAC,IAAI,EAAE,MAAM,EAAE,EAAE,MAAM,EAAE,MAAM,EAAE,GAAG,MAAM;IAQzE;;;;;OAKG;WACW,yBAAyB,CAAC,MAAM,EAAE,MAAM,EAAE,EAAE,MAAM,EAAE,MAAM,EAAE,GAAG,MAAM;CAcpF;AAED;;;GAGG;AACH,qBAAa,YAAY;IACvB,uEAAuE;WACzD,aAAa,CAAC,KAAK,EAAE,SAAS,GAAG,OAAO,EAAE,EAAE,KAAK,EAAE,SAAS,GAAG,OAAO,EAAE,GAAG,OAAO;IAYhG;;;OAGG;WACW,iBAAiB,CAAC,IAAI,EAAE,OAAO,EAAE,GAAG,OAAO,EAAE;IAG3D;;;OAGG;WACW,qCAAqC,CAAC,IAAI,EAAE,KAAK,EAAE,GAAG,MAAM;CAc3E;AAED;;;GAGG;AACH,qBAAa,aAAa;IACxB,uEAAuE;WACzD,aAAa,CAAC,KAAK,EAAE,SAAS,GAAG,QAAQ,EAAE,EAAE,KAAK,EAAE,SAAS,GAAG,QAAQ,EAAE,GAAG,OAAO;IAWlG;;;OAGG;WACW,kBAAkB,CAAC,IAAI,EAAE,MAAM,EAAE,GAAG,QAAQ,EAAE;CAG7D;AAED;;;GAGG;AACH,qBAAa,YAAY;IACvB,kFAAkF;WACpE,kCAAkC,CAAC,IAAI,EAAE,OAAO,EAAE,GAAG,YAAY,GAAG,MAAM,EAAE,EAAE,OAAO,EAAE,MAAM,EAAE,GAAG,YAAY,EAC1H,MAAM,CAAC,EAAE,YAAY,GAAG,YAAY,GAAG,SAAS;IAkClD,oCAAoC;WACtB,kBAAkB,CAAC,IAAI,EAAE,OAAO,EAAE,EAAE,MAAM,CAAC,EAAE,YAAY,GAAG,YAAY;IAWtF,2CAA2C;WAC7B,oBAAoB,CAAC,IAAI,EAAE,YAAY,GAAG,OAAO,EAAE;IAOjE;OACG;WACW,oCAAoC,CAAC,IAAI,EAAE,YAAY,EAAE,MAAM,EAAE,OAAO,EAAE,EAAE,OAAO,EAAE,MAAM,EAAE,EACzG,cAAc,GAAE,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,KAAK,GAA0C;IAQjG,OAAO,CAAC,MAAM,CAAC,YAAY,CAAoB;IAC/C;;;;OAIG;WACW,eAAe,CAAC,SAAS,EAAE,SAAS,EAAE,IAAI,EAAE,YAAY,GAAG,IAAI;IAW7E,0FAA0F;WAC5E,aAAa,CAAC,KAAK,EAAE,OAAO,EAAE,GAAG,YAAY,GAAG,SAAS,EAAE,KAAK,EAAE,OAAO,EAAE,GAAG,YAAY,GAAG,SAAS,GAAG,OAAO;IAkB9H,mFAAmF;WACrE,cAAc,CAAC,IAAI,EAAE,OAAO,EAAE,GAAG,YAAY,EAAE,KAAK,EAAE,4BAA4B,EAAE,SAAS,GAAE,MAAqC,GAAG,OAAO;CAgB7J;AACD;;;GAGG;AAEH,qBAAa,YAAY;IACvB,mCAAmC;WACrB,kBAAkB,CAAC,IAAI,EAAE,OAAO,EAAE,GAAG,YAAY;IAU/D;;;;;;;;;;OAUG;WACW,wBAAwB,CAAC,OAAO,EAAE,YAAY,EAAE,EAAE,EAAE,MAAM,EAAE,EAAE,EAAE,MAAM,EAAE,EAAE,EAAE,MAAM,EAAE,EAAE,EAAE,MAAM,EAAE,EAAE,EAAE,MAAM,EAAE,EAAE,EAAE,MAAM;IAWpI;;;;;;OAMG;WACW,YAAY,CAAC,OAAO,EAAE,YAAY,EAAE,MAAM,EAAE,OAAO,EAAE,GAAG,MAAM;IAQ5E;;;;;;OAMG;WACW,YAAY,CAAC,OAAO,EAAE,YAAY,EAAE,MAAM,EAAE,OAAO,EAAE,GAAG,MAAM;IAQ5E;;;;;;OAMG;WACW,YAAY,CAAC,OAAO,EAAE,YAAY,EAAE,MAAM,EAAE,OAAO,EAAE,GAAG,MAAM;IAQ5E,OAAO,CAAC,MAAM,CAAC,UAAU,CAAuB;IAChD,OAAO,CAAC,MAAM,CAAC,SAAS,CAAuB;IAC/C,OAAO,CAAC,MAAM,CAAC,SAAS,CAAuB;IAC/C,OAAO,CAAC,MAAM,CAAC,SAAS,CAAuB;IAC/C;;;;OAIG;WACW,sBAAsB,CAAC,MAAM,EAAE,OAAO,EAAE,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,MAAM,CAAC,EAAE,OAAO,GAAG,OAAO;IAYnH;;;;;;OAMG;WACW,oCAAoC,CAAC,MAAM,EAAE,OAAO,EAAE,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,MAAM,CAAC,EAAE,SAAS,GAAG,SAAS;IAWrI,0EAA0E;WAC5D,2BAA2B,CAAC,IAAI,EAAE,YAAY,GAAG,MAAM,EAAE,GAAG,OAAO,EAAE;IAQnF;;;;OAIG;WACW,2BAA2B,CAAC,IAAI,EAAE,YAAY,EAAE,WAAW,EAAE,MAAM,GAAG,GAAG,EAAE;IAmBzF;;;;OAIG;WACW,8BAA8B,CAAC,IAAI,EAAE,YAAY,EAAE,WAAW,EAAE,MAAM,EAAE,SAAS,EAAE,MAAM,GAAG,GAAG,EAAE;IAuB/G,8FAA8F;WAChF,eAAe,CAAC,SAAS,EAAE,SAAS,EAAE,GAAG,EAAE,YAAY,GAAG,IAAI;IAU5E,gEAAgE;WAClD,aAAa,CAAC,KAAK,EAAE,OAAO,EAAE,GAAG,YAAY,GAAG,SAAS,EAAE,KAAK,EAAE,OAAO,EAAE,GAAG,YAAY,GAAG,SAAS,GAAG,OAAO;IAmB9H,uEAAuE;WACzD,QAAQ,CAAC,MAAM,EAAE,oBAAoB,GAAG,OAAO,EAAE,EAAE,MAAM,CAAC,EAAE,OAAO,GAAG,OAAO;IAiB3F,iEAAiE;WACnD,uBAAuB,CAAC,MAAM,EAAE,OAAO,EAAE,EAAE,UAAU,EAAE,GAAG,EAAE,SAAS,EAAE,QAAQ,GAAG,MAAM,GAAG,SAAS;IAgBlH,gIAAgI;WAClH,wCAAwC,CAAC,MAAM,EAAE,OAAO,EAAE,EAAE,UAAU,EAAE,OAAO,EAAE,MAAM,EAAE,QAAQ,EAAE,SAAS,EAAE,QAAQ,GAAG,MAAM,GAAG,SAAS;IAmBzJ,oEAAoE;WACtD,iBAAiB,CAAC,IAAI,EAAE,MAAM,EAAE,EAAE,UAAU,EAAE,MAAM,GAAG,MAAM;IAgB3E,8EAA8E;WAChE,cAAc,CAAC,IAAI,EAAE,OAAO,EAAE,GAAG,YAAY,EAAE,KAAK,EAAE,4BAA4B,EAAE,SAAS,GAAE,MAAqC,GAAG,OAAO;IAiB5J;;;OAGG;WACW,cAAc,CAAC,IAAI,EAAE,OAAO,EAAE,GAAG,YAAY,EAAE,cAAc,GAAE,OAAe,EAAE,cAAc,CAAC,EAAE,MAAM,GAAG,MAAM;IA0B9H;;;;;OAKG;WACY,kBAAkB,CAAC,MAAM,EAAE,OAAO,EAAE,EAAE,SAAS,GAAE,MAAqC,GAAG,MAAM;IAe9G;;;OAGG;WACW,iBAAiB,CAAC,IAAI,EAAE,QAAQ,EAAE,GAAG,YAAY,GAAG,OAAO,EAAE;IAiB3E;;;OAGG;WACW,iBAAiB,CAAC,IAAI,EAAE,MAAM,EAAE,GAAG,OAAO,EAAE;IAG1D;;;;OAIG;WACW,sBAAsB,CAAC,MAAM,EAAE,OAAO,EAAE,EAAE,aAAa,EAAE,MAAM,GAAG,OAAO,EAAE;IAczF,mDAAmD;IACnD,OAAO,CAAC,MAAM,CAAC,UAAU;IAEzB;;;OAGG;WACW,yBAAyB,CAAC,IAAI,EAAE,0BAA0B,GAAG,MAAM,EAAE,EAAE;IAKrF;;OAEG;WACW,0BAA0B,CAAC,IAAI,EAAE,QAAQ,EAAE,GAAG,MAAM,EAAE,EAAE;IAYtE;;OAEG;WACW,2BAA2B,CAAC,IAAI,EAAE,QAAQ,EAAE,GAAG,YAAY;IAiBzE;;;OAGG;WACW,yBAAyB,CAAC,IAAI,EAAE,0BAA0B,GAAG,GAAG,EAAE;IAMhF;;;;;;OAMG;WACW,gCAAgC,CAAC,MAAM,EAAE,OAAO,EAAE,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,WAAW,EAAE,OAAO,GAAG,MAAM;IAqB/I;;;;;OAKG;WACW,mBAAmB,CAAC,MAAM,EAAE,OAAO,EAAE,EAAE,UAAU,EAAE,OAAO,EAAE,EAAE,YAAY,EAAE,OAAO,EAAE,EAAE,eAAe,GAAE,OAAe;IAiDrI;;;OAGG;WACW,YAAY,CAAC,IAAI,EAAE,OAAO,EAAE,GAAG;QAAE,SAAS,EAAE,OAAO,CAAC;QAAC,SAAS,EAAE,OAAO,CAAC;QAAC,SAAS,EAAE,OAAO,CAAC;QAAC,SAAS,EAAE,OAAO,CAAA;KAAE,GAAG,SAAS;CAc5I"}
+{"version":3,"file":"PointHelpers.d.ts","sourceRoot":"","sources":["../../../src/geometry3d/PointHelpers.ts"],"names":[],"mappings":"AAWA,OAAO,EAAE,OAAO,EAAE,MAAM,uBAAuB,CAAC;AAChD,OAAO,EAAE,0BAA0B,EAAE,MAAM,2BAA2B,CAAC;AACvE,OAAO,EAAE,oBAAoB,EAAE,MAAM,wBAAwB,CAAC;AAC9D,OAAO,EAAE,4BAA4B,EAAE,MAAM,gCAAgC,CAAC;AAC9E,OAAO,EAAE,OAAO,EAAE,MAAM,mBAAmB,CAAC;AAE5C,OAAO,EAAE,OAAO,EAAE,QAAQ,EAAE,GAAG,EAAE,MAAM,mBAAmB,CAAC;AAE3D,OAAO,EAAE,SAAS,EAAE,MAAM,aAAa,CAAC;AACxC,OAAO,EAAE,KAAK,EAAE,MAAM,EAAE,QAAQ,EAAE,MAAM,YAAY,CAAC;AAiBrD;;;GAGG;AACH,qBAAa,WAAW;IACtB;;OAEG;WACW,UAAU,CAAC,IAAI,EAAE,MAAM,EAAE,GAAG,MAAM;IAgBhD,6FAA6F;WAC/E,YAAY,CAAC,KAAK,EAAE,GAAG,EAAE,GAAG,YAAY,GAAG,SAAS,EAAE,KAAK,EAAE,GAAG,EAAE,GAAG,YAAY,GAAG,SAAS,GAAG,OAAO;IAWrH,qFAAqF;WACvE,aAAa,CACzB,KAAK,EAAE,MAAM,EAAE,GAAG,YAAY,GAAG,SAAS,EAC1C,KAAK,EAAE,MAAM,EAAE,GAAG,YAAY,GAAG,SAAS,EAC1C,SAAS,GAAE,MAAqC,GAAG,OAAO;IAY5D,2FAA2F;WAC7E,GAAG,CAAC,IAAI,EAAE,MAAM,EAAE,GAAG,YAAY,GAAG,MAAM;IAKxD,0GAA0G;WAC5F,mBAAmB,CAAC,CAAC,EAAE,MAAM,EAAE,IAAI,EAAE,MAAM,EAAE,GAAG,SAAS,GAAG,OAAO;IAMjF,2DAA2D;WAC7C,WAAW,CAAC,MAAM,EAAE,MAAM,EAAE,GAAG,MAAM;IAcnD,yDAAyD;WAC3C,SAAS,CAAC,EAAE,EAAE,MAAM,EAAE,EAAE,EAAE,MAAM,GAAG,MAAM;IAKvD;;OAEG;WACW,UAAU,CAAC,KAAK,EAAE,MAAM,EAAE,GAAG,YAAY,EAAE,KAAK,EAAE,MAAM,EAAE,GAAG,YAAY,GAAG,MAAM;IAOhG;;OAEG;WACW,iBAAiB,CAAC,KAAK,EAAE,YAAY,EAAE,KAAK,EAAE,YAAY,GAAG,MAAM;IAMjF;;;;;OAKG;WACW,0BAA0B,CAAC,GAAG,EAAE,MAAM,EAAE,IAAI,EAAE,MAAM,EAAE,IAAI,EAAE,MAAM,GAAG,MAAM,EAAE;IAc3F,qDAAqD;WACvC,MAAM,CAAC,MAAM,EAAE,MAAM,EAAE,GAAG,YAAY,GAAG,MAAM,EAAE;IAO/D,uBAAuB;WACT,MAAM,CAAC,MAAM,EAAE,MAAM,EAAE,EAAE,GAAG,MAAM,EAAE,EAAE;IAWpD,yBAAyB;WACX,MAAM,CAAC,MAAM,EAAE,MAAM,EAAE,EAAE,EAAE,GAAG,MAAM,EAAE,EAAE,EAAE;IAexD;;OAEG;WACW,QAAQ,CAAC,MAAM,EAAE,YAAY,EAAE,WAAW,EAAE,MAAM,GAAG,MAAM,EAAE,EAAE,GAAG,SAAS;IAMzF;;;OAGG;WACW,QAAQ,CAAC,MAAM,EAAE,YAAY,EAAE,SAAS,EAAE,MAAM,EAAE,WAAW,EAAE,MAAM,GAAG,MAAM,EAAE,EAAE,EAAE,GAAG,SAAS;IAM9G,wDAAwD;WAC1C,IAAI,CAAC,MAAM,EAAE,MAAM,EAAE,GAAG,MAAM,EAAE,EAAE,GAAG,MAAM,EAAE,EAAE,EAAE,GAAG,YAAY;IAuC9E,6HAA6H;WAC/G,gCAAgC,CAAC,KAAK,EAAE,MAAM,EAAE,GAAG,SAAS,EAAE,OAAO,EAAE,MAAM,EAAE,OAAO,EAAE,MAAM,GAAG,MAAM,EAAE;IAqBvH;;;;OAIG;WACW,iBAAiB,CAAC,IAAI,EAAE,MAAM,EAAE,EAAE,MAAM,EAAE,MAAM,EAAE,GAAG,MAAM;IAQzE;;;;;OAKG;WACW,yBAAyB,CAAC,MAAM,EAAE,MAAM,EAAE,EAAE,MAAM,EAAE,MAAM,EAAE,GAAG,MAAM;CAcpF;AAED;;;GAGG;AACH,qBAAa,YAAY;IACvB,uEAAuE;WACzD,aAAa,CAAC,KAAK,EAAE,SAAS,GAAG,OAAO,EAAE,EAAE,KAAK,EAAE,SAAS,GAAG,OAAO,EAAE,GAAG,OAAO;IAYhG;;;OAGG;WACW,iBAAiB,CAAC,IAAI,EAAE,OAAO,EAAE,GAAG,OAAO,EAAE;IAG3D;;;OAGG;WACW,qCAAqC,CAAC,IAAI,EAAE,KAAK,EAAE,GAAG,MAAM;CAc3E;AAED;;;GAGG;AACH,qBAAa,aAAa;IACxB,uEAAuE;WACzD,aAAa,CAAC,KAAK,EAAE,SAAS,GAAG,QAAQ,EAAE,EAAE,KAAK,EAAE,SAAS,GAAG,QAAQ,EAAE,GAAG,OAAO;IAWlG;;;OAGG;WACW,kBAAkB,CAAC,IAAI,EAAE,MAAM,EAAE,GAAG,QAAQ,EAAE;CAG7D;AAED;;;GAGG;AACH,qBAAa,YAAY;IACvB;;;;;;OAMG;WACW,kCAAkC,CAAC,IAAI,EAAE,OAAO,EAAE,GAAG,YAAY,GAAG,MAAM,EAAE,EAAE,OAAO,EAAE,MAAM,EAAE,GAAG,YAAY,EAAE,MAAM,CAAC,EAAE,YAAY,GAAG,YAAY,GAAG,SAAS;IAkC9K;;;;;OAKG;WACW,kBAAkB,CAAC,IAAI,EAAE,OAAO,EAAE,EAAE,MAAM,CAAC,EAAE,YAAY,GAAG,YAAY;IAatF,2CAA2C;WAC7B,oBAAoB,CAAC,IAAI,EAAE,YAAY,GAAG,OAAO,EAAE;IAOjE;;;;;;;OAOG;WACW,oCAAoC,CAAC,IAAI,EAAE,YAAY,EAAE,MAAM,EAAE,OAAO,EAAE,EAAE,OAAO,EAAE,MAAM,EAAE,EACzG,cAAc,GAAE,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,KAAK,GAA0C;IAQjG,OAAO,CAAC,MAAM,CAAC,YAAY,CAAoB;IAC/C;;;;OAIG;WACW,eAAe,CAAC,SAAS,EAAE,SAAS,EAAE,IAAI,EAAE,YAAY,GAAG,IAAI;IAW7E,iGAAiG;WACnF,aAAa,CAAC,KAAK,EAAE,OAAO,EAAE,GAAG,YAAY,GAAG,SAAS,EAAE,KAAK,EAAE,OAAO,EAAE,GAAG,YAAY,GAAG,SAAS,GAAG,OAAO;IAgC9H,mFAAmF;WACrE,cAAc,CAAC,IAAI,EAAE,OAAO,EAAE,GAAG,YAAY,EAAE,KAAK,EAAE,4BAA4B,EAAE,SAAS,GAAE,MAAqC,GAAG,OAAO;CAgB7J;AACD;;;GAGG;AAEH,qBAAa,YAAY;IACvB;;;;;OAKG;WACW,kBAAkB,CAAC,IAAI,EAAE,OAAO,EAAE,EAAE,MAAM,CAAC,EAAE,YAAY,GAAG,YAAY;IAYtF;;;;;;;;;;OAUG;WACW,wBAAwB,CAAC,OAAO,EAAE,YAAY,EAAE,EAAE,EAAE,MAAM,EAAE,EAAE,EAAE,MAAM,EAAE,EAAE,EAAE,MAAM,EAAE,EAAE,EAAE,MAAM,EAAE,EAAE,EAAE,MAAM,EAAE,EAAE,EAAE,MAAM;IAWpI;;;;;;OAMG;WACW,YAAY,CAAC,OAAO,EAAE,YAAY,EAAE,MAAM,EAAE,OAAO,EAAE,GAAG,MAAM;IAQ5E;;;;;;OAMG;WACW,YAAY,CAAC,OAAO,EAAE,YAAY,EAAE,MAAM,EAAE,OAAO,EAAE,GAAG,MAAM;IAQ5E;;;;;;OAMG;WACW,YAAY,CAAC,OAAO,EAAE,YAAY,EAAE,MAAM,EAAE,OAAO,EAAE,GAAG,MAAM;IAQ5E,OAAO,CAAC,MAAM,CAAC,UAAU,CAAuB;IAChD,OAAO,CAAC,MAAM,CAAC,SAAS,CAAuB;IAC/C,OAAO,CAAC,MAAM,CAAC,SAAS,CAAuB;IAC/C,OAAO,CAAC,MAAM,CAAC,SAAS,CAAuB;IAC/C;;;;OAIG;WACW,sBAAsB,CAAC,MAAM,EAAE,OAAO,EAAE,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,MAAM,CAAC,EAAE,OAAO,GAAG,OAAO;IAYnH;;;;;;OAMG;WACW,oCAAoC,CAAC,MAAM,EAAE,OAAO,EAAE,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,MAAM,CAAC,EAAE,SAAS,GAAG,SAAS;IAWrI,0EAA0E;WAC5D,2BAA2B,CAAC,IAAI,EAAE,YAAY,GAAG,MAAM,EAAE,GAAG,OAAO,EAAE;IAQnF;;;;OAIG;WACW,2BAA2B,CAAC,IAAI,EAAE,YAAY,EAAE,WAAW,EAAE,MAAM,GAAG,GAAG,EAAE;IAmBzF;;;;;OAKG;WACW,8BAA8B,CAAC,IAAI,EAAE,YAAY,EAAE,WAAW,EAAE,MAAM,EAAE,SAAS,EAAE,MAAM,GAAG,GAAG,EAAE;IAuB/G,8FAA8F;WAChF,eAAe,CAAC,SAAS,EAAE,SAAS,EAAE,GAAG,EAAE,YAAY,GAAG,IAAI;IAU5E,iGAAiG;WACnF,aAAa,CAAC,KAAK,EAAE,OAAO,EAAE,GAAG,YAAY,GAAG,SAAS,EAAE,KAAK,EAAE,OAAO,EAAE,GAAG,YAAY,GAAG,SAAS,GAAG,OAAO;IAgC9H,uEAAuE;WACzD,QAAQ,CAAC,MAAM,EAAE,oBAAoB,GAAG,OAAO,EAAE,EAAE,MAAM,CAAC,EAAE,OAAO,GAAG,OAAO;IAiB3F,iEAAiE;WACnD,uBAAuB,CAAC,MAAM,EAAE,OAAO,EAAE,EAAE,UAAU,EAAE,GAAG,EAAE,SAAS,EAAE,QAAQ,GAAG,MAAM,GAAG,SAAS;IAgBlH,gIAAgI;WAClH,wCAAwC,CAAC,MAAM,EAAE,OAAO,EAAE,EAAE,UAAU,EAAE,OAAO,EAAE,MAAM,EAAE,QAAQ,EAAE,SAAS,EAAE,QAAQ,GAAG,MAAM,GAAG,SAAS;IAmBzJ,oEAAoE;WACtD,iBAAiB,CAAC,IAAI,EAAE,MAAM,EAAE,EAAE,UAAU,EAAE,MAAM,GAAG,MAAM;IAgB3E,8EAA8E;WAChE,cAAc,CAAC,IAAI,EAAE,OAAO,EAAE,GAAG,YAAY,EAAE,KAAK,EAAE,4BAA4B,EAAE,SAAS,GAAE,MAAqC,GAAG,OAAO;IAiB5J;;;OAGG;WACW,cAAc,CAAC,IAAI,EAAE,OAAO,EAAE,GAAG,YAAY,EAAE,cAAc,GAAE,OAAe,EAAE,cAAc,CAAC,EAAE,MAAM,GAAG,MAAM;IA0B9H;;;;;OAKG;WACW,kBAAkB,CAAC,MAAM,EAAE,OAAO,EAAE,EAAE,SAAS,GAAE,MAAqC,GAAG,MAAM;IAe7G;;;OAGG;WACW,iBAAiB,CAAC,IAAI,EAAE,QAAQ,EAAE,GAAG,YAAY,GAAG,OAAO,EAAE;IAiB3E;;;OAGG;WACW,iBAAiB,CAAC,IAAI,EAAE,MAAM,EAAE,GAAG,OAAO,EAAE;IAG1D;;;;OAIG;WACW,sBAAsB,CAAC,MAAM,EAAE,OAAO,EAAE,EAAE,aAAa,EAAE,MAAM,GAAG,OAAO,EAAE;IAczF,mDAAmD;IACnD,OAAO,CAAC,MAAM,CAAC,UAAU;IAEzB;;;OAGG;WACW,yBAAyB,CAAC,IAAI,EAAE,0BAA0B,GAAG,MAAM,EAAE,EAAE;IAKrF;;OAEG;WACW,0BAA0B,CAAC,IAAI,EAAE,QAAQ,EAAE,GAAG,MAAM,EAAE,EAAE;IAoBtE;;OAEG;WACW,2BAA2B,CAAC,IAAI,EAAE,QAAQ,EAAE,GAAG,YAAY;IAqBzE;;;OAGG;WACW,yBAAyB,CAAC,IAAI,EAAE,0BAA0B,GAAG,GAAG,EAAE;IAMhF;;;;;;OAMG;WACW,gCAAgC,CAAC,MAAM,EAAE,OAAO,EAAE,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,WAAW,EAAE,OAAO,GAAG,MAAM;IAoB/I;;;;;OAKG;WACW,mBAAmB,CAAC,MAAM,EAAE,OAAO,EAAE,EAAE,UAAU,EAAE,OAAO,EAAE,EAAE,YAAY,EAAE,OAAO,EAAE,EAAE,eAAe,GAAE,OAAe;IAiDrI;;;OAGG;WACW,YAAY,CAAC,IAAI,EAAE,OAAO,EAAE,GAAG;QAAE,SAAS,EAAE,OAAO,CAAC;QAAC,SAAS,EAAE,OAAO,CAAC;QAAC,SAAS,EAAE,OAAO,CAAC;QAAC,SAAS,EAAE,OAAO,CAAA;KAAE,GAAG,SAAS;CAc5I"}

package/lib/cjs/geometry3d/PointHelpers.js

@@ -68,12 +68,12 @@ class NumberArray {
         return (dataA === undefined && dataB === undefined);
     }
     /** Return true if arrays have identical counts and entries equal within tolerance */
-    static isAlmostEqual(dataA, dataB, tolerance) {
+    static isAlmostEqual(dataA, dataB, tolerance = Geometry_1.Geometry.smallMetricDistance) {
         if (dataA && dataB) {
             if (dataA.length !== dataB.length)
                 return false;
             for (let i = 0; i < dataA.length; i++)
-                if (Math.abs(dataA[i] - dataB[i]) >= tolerance)
+                if (Math.abs(dataA[i] - dataB[i]) > tolerance)
                     return false;
             return true;
         }
@@ -159,13 +159,97 @@ class NumberArray {
         result.push(high);
         return result;
     }
-    /** copy numbers from variant sources to number[]. */
+    /** Copy numbers from variant sources to number[]. */
     static create(source) {
         const result = [];
         for (const q of source)
             result.push(q);
         return result;
     }
+    /** Copy number[][]. */
+    static copy2d(source) {
+        const result = [];
+        for (const row of source) {
+            const newRow = [];
+            for (const entry of row)
+                newRow.push(entry);
+            result.push(newRow);
+        }
+        return result;
+    }
+    /** Copy number[][][]. */
+    static copy3d(source) {
+        const result = [];
+        for (const row of source) {
+            const newRow = [];
+            for (const block of row) {
+                const newBlock = [];
+                for (const entry of block)
+                    newBlock.push(entry);
+                newRow.push(newBlock);
+            }
+            result.push(newRow);
+        }
+        return result;
+    }
+    /** Copy numbers from Float64Array to number[][].
+     * @param numPerBlock block size
+     */
+    static unpack2d(source, numPerBlock) {
+        if (numPerBlock < 1)
+            return undefined;
+        return Point3dArray.unpackNumbersToNestedArrays(source, numPerBlock);
+    }
+    /** Copy numbers from Float64Array to number[][][].
+     * @param numPerRow row size
+     * @param numPerBlock block size
+     */
+    static unpack3d(source, numPerRow, numPerBlock) {
+        if (numPerBlock < 1 || numPerRow < 1)
+            return undefined;
+        return Point3dArray.unpackNumbersToNestedArraysIJK(source, numPerBlock, numPerRow);
+    }
+    /** Copy numbers from 1d/2d/3d array to Float64Array. */
+    static pack(source) {
+        const numRows = source.length;
+        let numPerRow = 0;
+        let numPerBlock = 0;
+        let numCoords = 0;
+        if (numRows > 0) {
+            numCoords = numRows;
+            if (Array.isArray(source[0])) {
+                numPerRow = source[0].length;
+                if (numPerRow > 0) {
+                    numCoords *= numPerRow;
+                    if (Array.isArray(source[0][0])) {
+                        numPerBlock = source[0][0].length;
+                        if (numPerBlock > 0)
+                            numCoords *= numPerBlock;
+                    }
+                }
+            }
+        }
+        const result = new Float64Array(numCoords);
+        if (numPerBlock > 0) {
+            const src3d = source;
+            for (let i = 0, c = 0; i < numRows; ++i)
+                for (let j = 0; j < numPerRow; ++j)
+                    for (let k = 0; k < numPerBlock; ++k)
+                        result[c++] = src3d[i][j][k];
+        }
+        else if (numPerRow > 0) {
+            const src2d = source;
+            for (let i = 0, c = 0; i < numRows; ++i)
+                for (let j = 0; j < numPerRow; ++j)
+                    result[c++] = src2d[i][j];
+        }
+        else if (numRows > 0) {
+            const src1d = source;
+            for (let i = 0, c = 0; i < numRows; ++i)
+                result[c++] = src1d[i];
+        }
+        return result;
+    }
     /** Return a copy of the knots array, with multiplicity of first and last knots raised or lowered to expectedMultiplicity. */
     static cloneWithStartAndEndMultiplicity(knots, target0, target1) {
         const result = [];
@@ -299,16 +383,23 @@ exports.Vector3dArray = Vector3dArray;
  * @public
  */
 class Point4dArray {
-    /** pack each point and its corresponding weight into a buffer of xyzw xyzw ... */
+    /**
+     * Copy each weighted point and its corresponding weight into a packed buffer.
+     * @param data array of weighted xyz
+     * @param weights scalar weight array
+     * @param result optional destination array. If insufficiently sized, a new array is returned.
+     * @return packed weighted point array
+     */
     static packPointsAndWeightsToFloat64Array(data, weights, result) {
+        let points;
         if (Array.isArray(data) && data[0] instanceof Point3dVector3d_1.Point3d) {
-            const points = data;
+            points = data;
             if (points.length !== weights.length)
                 return undefined;
-            result = result ? result : new Float64Array(4 * points.length);
-            let i = 0;
-            let k = 0;
-            for (k = 0; k < points.length; k++) {
+            const numValues = 4 * points.length;
+            if (!result || result.length < numValues)
+                result = new Float64Array(numValues);
+            for (let i = 0, k = 0; k < points.length; k++) {
                 result[i++] = points[k].x;
                 result[i++] = points[k].y;
                 result[i++] = points[k].z;
@@ -316,28 +407,32 @@ class Point4dArray {
             }
             return result;
         }
-        else {
-            const points = data;
-            const numPoints = weights.length;
-            if (points.length !== 3 * numPoints)
-                return undefined;
-            let i = 0;
-            let k;
-            result = result ? result : new Float64Array(4 * numPoints);
-            for (k = 0; k < numPoints; k++) {
-                const k0 = 3 * k;
-                result[i++] = points[k0];
-                result[i++] = points[k0 + 1];
-                result[i++] = points[k0 + 2];
-                result[i++] = weights[k];
-            }
-            return result;
+        points = data;
+        const numPoints = weights.length;
+        if (points.length !== 3 * numPoints)
+            return undefined;
+        const numValues1 = 4 * numPoints;
+        if (!result || result.length < numValues1)
+            result = new Float64Array(numValues1);
+        for (let i = 0, k = 0; k < numPoints; k++) {
+            const k0 = 3 * k;
+            result[i++] = points[k0];
+            result[i++] = points[k0 + 1];
+            result[i++] = points[k0 + 2];
+            result[i++] = weights[k];
         }
-        return undefined;
+        return result;
     }
-    /** pack x,y,z,w in Float64Array. */
+    /**
+     * Copy 4d points into a packed buffer.
+     * @param data array of xyzw
+     * @param result optional destination array. If insufficiently sized, a new array is returned.
+     * @return packed point array
+     */
     static packToFloat64Array(data, result) {
-        result = result ? result : new Float64Array(4 * data.length);
+        const numValues = 4 * data.length;
+        if (!result || result.length < numValues)
+            result = new Float64Array(numValues);
         let i = 0;
         for (const p of data) {
             result[i++] = p.x;
@@ -355,7 +450,13 @@ class Point4dArray {
         }
         return result;
     }
-    /** unpack from xyzw xyzw... array to array of Point3d and array of weight.
+    /**
+     * Unpack packed 4D data to a Point3d array and an array of weights.
+     * * `WeightStyle` of `data` is not assumed. If input data is of form [a,b,c,d], default output arrays will have form [a,b,c] and [d].
+     * @param data input 4D points (packed)
+     * @param points output 3D data
+     * @param weights output weights (w portion of input)
+     * @param pointFormatter optional xyz formatter. By default, returns a Point3d created from the xyz portion of the input.
      */
     static unpackFloat64ArrayToPointsAndWeights(data, points, weights, pointFormatter = (x, y, z) => Point3dVector3d_1.Point3d.create(x, y, z)) {
         points.length = 0;
@@ -381,21 +482,36 @@ class Point4dArray {
             xyzw[i + 3] = xyzw1.w;
         }
     }
-    /** test for near equality of all corresponding numeric values, treated as coordinates. */
+    /** Test arrays for near equality of all corresponding numeric values, treated as coordinates. */
     static isAlmostEqual(dataA, dataB) {
         if (dataA && dataB) {
-            if (dataA.length !== dataB.length)
-                return false;
             if (dataA instanceof Float64Array && dataB instanceof Float64Array) {
+                if (dataA.length !== dataB.length)
+                    return false;
                 for (let i = 0; i < dataA.length; i++)
                     if (!Geometry_1.Geometry.isSameCoordinate(dataA[i], dataB[i]))
                         return false;
             }
             else if (Array.isArray(dataA) && Array.isArray(dataB)) {
+                if (dataA.length !== dataB.length)
+                    return false;
                 for (let i = 0; i < dataA.length; i++)
                     if (!dataA[i].isAlmostEqual(dataB[i]))
                         return false;
             }
+            else { // different types
+                const points = dataA instanceof Float64Array ? dataB : dataA;
+                const numbers = dataA instanceof Float64Array ? dataA : dataB;
+                if (numbers.length !== points.length * 4)
+                    return false;
+                for (let iPoint = 0; iPoint < points.length; ++iPoint) {
+                    if (!Geometry_1.Geometry.isSameCoordinate(points[iPoint].x, numbers[4 * iPoint]) ||
+                        !Geometry_1.Geometry.isSameCoordinate(points[iPoint].y, numbers[4 * iPoint + 1]) ||
+                        !Geometry_1.Geometry.isSameCoordinate(points[iPoint].z, numbers[4 * iPoint + 2]) ||
+                        !Geometry_1.Geometry.isSameCoordinate(points[iPoint].w, numbers[4 * iPoint + 3]))
+                        return false;
+                }
+            }
             return true;
         }
         // if both are null it is equal, otherwise unequal
@@ -426,9 +542,16 @@ exports.Point4dArray = Point4dArray;
  * @public
  */
 class Point3dArray {
-    /** pack x,y,z to `Float64Array` */
-    static packToFloat64Array(data) {
-        const result = new Float64Array(3 * data.length);
+    /**
+     * Copy 3d points into a packed buffer.
+     * @param data array of xyz
+     * @param result optional destination array. If insufficiently sized, a new array is returned.
+     * @return packed point array
+     */
+    static packToFloat64Array(data, result) {
+        const numValues = 3 * data.length;
+        if (!result || result.length < numValues)
+            result = new Float64Array(numValues);
         let i = 0;
         for (const p of data) {
             result[i++] = p.x;
@@ -564,9 +687,10 @@ class Point3dArray {
         return result;
     }
     /**
-     * return an 3-dimensional array containing all the values of `data` in arrays numPerRow blocks of numPerBlock
+     * Return a 3-dimensional array containing all the values of `data` in rows of numPerRow blocks of size numPerBlock.
      * @param data simple array of numbers
-     * @param numPerBlock number of values in each block at first level down
+     * @param numPerBlock number of values in each block
+     * @param numPerRow number of blocks per row
      */
     static unpackNumbersToNestedArraysIJK(data, numPerBlock, numPerRow) {
         const result = [];
@@ -602,21 +726,35 @@ class Point3dArray {
             xyz[i + 2] = xyz1.z;
         }
     }
-    /** Apply Geometry.isAlmostEqual to corresponding coordinates */
+    /** Test arrays for near equality of all corresponding numeric values, treated as coordinates. */
     static isAlmostEqual(dataA, dataB) {
         if (dataA && dataB) {
-            if (dataA.length !== dataB.length)
-                return false;
             if (dataA instanceof Float64Array && dataB instanceof Float64Array) {
+                if (dataA.length !== dataB.length)
+                    return false;
                 for (let i = 0; i < dataA.length; i++)
                     if (!Geometry_1.Geometry.isSameCoordinate(dataA[i], dataB[i]))
                         return false;
             }
             else if (Array.isArray(dataA) && Array.isArray(dataB)) {
+                if (dataA.length !== dataB.length)
+                    return false;
                 for (let i = 0; i < dataA.length; i++)
                     if (!dataA[i].isAlmostEqual(dataB[i]))
                         return false;
             }
+            else { // different types
+                const points = dataA instanceof Float64Array ? dataB : dataA;
+                const numbers = dataA instanceof Float64Array ? dataA : dataB;
+                if (numbers.length !== points.length * 3)
+                    return false;
+                for (let iPoint = 0; iPoint < points.length; ++iPoint) {
+                    if (!Geometry_1.Geometry.isSameCoordinate(points[iPoint].x, numbers[3 * iPoint]) ||
+                        !Geometry_1.Geometry.isSameCoordinate(points[iPoint].y, numbers[3 * iPoint + 1]) ||
+                        !Geometry_1.Geometry.isSameCoordinate(points[iPoint].z, numbers[3 * iPoint + 2]))
+                        return false;
+                }
+            }
             return true;
         }
         // if both are null it is equal, otherwise unequal
@@ -825,7 +963,16 @@ class Point3dArray {
                 result.push([p.x, p.y, p.z]);
             }
             else if (Array.isArray(p)) {
-                result.push([p[0], p[1], p.length > 2 ? p[2] : 0.0]);
+                const x = p.length > 0 ? p[0] : 0.0;
+                const y = p.length > 1 ? p[1] : 0.0;
+                const z = p.length > 2 ? p[2] : 0.0;
+                result.push([x, y, z]);
+            }
+            else {
+                const x = p.x !== undefined ? p.x : 0.0;
+                const y = p.y !== undefined ? p.y : 0.0;
+                const z = p.z !== undefined ? p.z : 0.0;
+                result.push([x, y, z]);
             }
         }
         return result;
@@ -843,9 +990,14 @@ class Point3dArray {
                 result[i++] = p.z;
             }
             else if (Array.isArray(p)) {
-                result[i++] = p[0];
-                result[i++] = p[1];
-                result[i++] = p.length > 2 ? p[2] : 0.0; // allow missing z
+                result[i++] = p.length > 0 ? p[0] : 0.0;
+                result[i++] = p.length > 1 ? p[1] : 0.0;
+                result[i++] = p.length > 2 ? p[2] : 0.0;
+            }
+            else {
+                result[i++] = p.x !== undefined ? p.x : 0.0;
+                result[i++] = p.y !== undefined ? p.y : 0.0;
+                result[i++] = p.z !== undefined ? p.z : 0.0;
             }
         }
         return result;
@@ -860,37 +1012,36 @@ class Point3dArray {
         return collector.claimResult();
     }
     /**
-     * return perpendicular distance from points[indexB] to the segment points[indexA] to points[indexC].
-     * * extrapolation option when projection is outside of fraction range 0..1 are:
-     *   * false ==> measure distance to closest endpoint
-     *   * true ==> measure distance to extended line segment.
-     * (no index checking!)
+     * Return perpendicular distance from points[indexB] to the segment from points[indexA] to points[indexC].
+     * * Extrapolation options when the projection is outside of the fraction range [0,1] are:
+     *   * false ==> return distance to closest endpoint
+     *   * true ==> return distance to extended line segment
+     * * There is no index checking!
      */
     static distanceIndexedPointBToSegmentAC(points, indexA, indexB, indexC, extrapolate) {
         const vectorU = Point3dVector3d_1.Vector3d.createStartEnd(points[indexA], points[indexC]);
         const vectorV = Point3dVector3d_1.Vector3d.createStartEnd(points[indexA], points[indexB]);
         const uDotU = vectorU.dotProduct(vectorU);
         const uDotV = vectorU.dotProduct(vectorV);
-        let fraction = Geometry_1.Geometry.conditionalDivideFraction(uDotV, uDotU);
+        const fraction = Geometry_1.Geometry.conditionalDivideFraction(uDotV, uDotU);
         if (fraction === undefined)
-            fraction = 0.0;
+            return vectorV.magnitude(); // AC is degenerate; return ||B-A||
         if (!extrapolate) {
             if (fraction > 1.0)
-                fraction = 1.0;
+                return points[indexB].distance(points[indexC]); // return ||B-C||
             if (fraction < 0.0)
-                fraction = 0.0;
+                return vectorV.magnitude(); // return ||B-A||
         }
-        let h2 = vectorV.magnitudeSquared() - fraction * fraction * uDotU;
-        // h2 should never be negative except for quirky tolerance ..
-        if (h2 < 0.0)
-            h2 = 0.0;
-        return Math.sqrt(h2);
+        // return distance to projection on (extended) segment
+        const h2 = vectorV.magnitudeSquared() - fraction * fraction * uDotU;
+        // h2 should never be negative except for quirky tolerance...
+        return h2 <= 0.0 ? 0.0 : Math.sqrt(h2);
     }
     /** Computes the hull of the XY projection of points.
-     * * Returns the hull as an array of Point3d
-     * * Optionally returns non-hull points in `insidePoints[]`
-     * * If both arrays empty if less than 3 points.
-     * *
+     * @param points input points, z-coordinates ignored
+     * @param hullPoints (output) points on the convex hull (cloned from input points)
+     * @param insidePoints (output) points not on the convex hull (cloned from input points)
+     * @param addClosurePoint whether to append the first hull point to `hullPoints`
      */
     static computeConvexHullXY(points, hullPoints, insidePoints, addClosurePoint = false) {
         hullPoints.length = 0;