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

package/lib/esm/polyface/multiclip/XYPointBuckets.d.ts

@@ -0,0 +1,73 @@
+/** @packageDocumentation
+ * @module RangeSearch
+ */
+import { IndexedXYZCollection } from "../../geometry3d/IndexedXYZCollection";
+import { Range2d, Range3d } from "../../geometry3d/Range";
+import { LowAndHighXY } from "../../geometry3d/XYZProps";
+/** Type for a value which may be either (a) undefined or (b) an array of type []. */
+export declare type OptionalArray<T> = T[] | undefined;
+/**
+ * Arrays of type T values distributed by xy position when entered.
+ * @internal
+ */
+export declare class XYIndexGrid<T> {
+    private _range;
+    private _numXEdge;
+    private _numYEdge;
+    private _data;
+    protected constructor(range: Range2d, numX: number, numY: number);
+    /** Return the number of x edges in the grid */
+    get numXEdge(): number;
+    /** Return the number of y edges in the grid */
+    get numYEdge(): number;
+    /** Return the `i` index of cells containing x coordinate */
+    xIndex(x: number): number;
+    /** Return the `j` index of cells containing x coordinate */
+    yIndex(y: number): number;
+    /**
+     * Construct an array with cells mapped to a range, with counts determined by estimated total count and target number of entries per cell.
+     * @param range
+     * @param totalEntries
+     * @param targetEntriesPerCell
+     */
+    static createWithEstimatedCounts<T>(range: LowAndHighXY, totalEntries: number, targetEntriesPerCell: number): XYIndexGrid<T> | undefined;
+    /**
+     * Add (save) a new data value to the grid cell containing x,y
+     * @param x
+     * @param y
+     * @param value
+     */
+    addDataAtXY(x: number, y: number, value: T): void;
+    /**
+     * Get the (reference to the possibly null array of) data values for the cell indicated by xy.
+     * @param x
+     * @param y
+     */
+    getDataAtXY(x: number, y: number): OptionalArray<T>;
+    /**
+     * Get the (reference to the possibly null array of) data values for the cell indicated by indices in the x and y direction
+     * @param xIndex
+     * @param yIndex
+     */
+    getDataAtIndex(xIndex: number, yIndex: number): OptionalArray<T>;
+    /** Return true if (xIndex, yIndex) is a valid cell index. */
+    isValidIndex(xIndex: number, yIndex: number): boolean;
+}
+/** Manage buckets of points for fast search.
+ * @internal
+ */
+export declare class XYPointBuckets {
+    private _points;
+    private _buckets;
+    /** Return the underlying grid with indices recorded by block */
+    get indexGrid(): XYIndexGrid<number>;
+    private constructor();
+    /** Create an XYIndex grid with all indices of all `points` entered */
+    static create(points: IndexedXYZCollection, targetPointsPerCell: number): XYPointBuckets | undefined;
+    /** call the `announce` function with the index and coordinates of all points in given range.
+     * * continue the search if `announce` returns true.
+     * * terminate the search if `announce` returns false;
+     */
+    announcePointsInRange(range: Range2d | Range3d, announce: (index: number, x: number, y: number, z: number) => boolean): void;
+}
+//# sourceMappingURL=XYPointBuckets.d.ts.map

package/lib/esm/polyface/multiclip/XYPointBuckets.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"XYPointBuckets.d.ts","sourceRoot":"","sources":["../../../../src/polyface/multiclip/XYPointBuckets.ts"],"names":[],"mappings":"AAIA;;GAEG;AAGH,OAAO,EAAE,oBAAoB,EAAE,MAAM,uCAAuC,CAAC;AAE7E,OAAO,EAAE,OAAO,EAAE,OAAO,EAAE,MAAM,wBAAwB,CAAC;AAC1D,OAAO,EAAE,YAAY,EAAE,MAAM,2BAA2B,CAAC;AAEzD,qFAAqF;AACrF,oBAAY,aAAa,CAAC,CAAC,IAAI,CAAC,EAAE,GAAG,SAAS,CAAC;AAC/C;;;GAGG;AACH,qBAAa,WAAW,CAAC,CAAC;IACxB,OAAO,CAAC,MAAM,CAAU;IACxB,OAAO,CAAC,SAAS,CAAS;IAC1B,OAAO,CAAC,SAAS,CAAS;IAE1B,OAAO,CAAC,KAAK,CAAiC;IAE9C,SAAS,aAAa,KAAK,EAAE,OAAO,EAAE,IAAI,EAAE,MAAM,EAAE,IAAI,EAAE,MAAM;IAahE,+CAA+C;IAC/C,IAAW,QAAQ,IAAI,MAAM,CAA2B;IACxD,+CAA+C;IAC/C,IAAW,QAAQ,IAAI,MAAM,CAA2B;IAExD,4DAA4D;IACrD,MAAM,CAAC,CAAC,EAAE,MAAM,GAAG,MAAM;IAShC,4DAA4D;IACrD,MAAM,CAAC,CAAC,EAAE,MAAM,GAAG,MAAM;IAShC;;;;;OAKG;WACW,yBAAyB,CAAC,CAAC,EAAE,KAAK,EAAE,YAAY,EAAE,YAAY,EAAE,MAAM,EAAE,oBAAoB,EAAE,MAAM,GAAG,WAAW,CAAC,CAAC,CAAC,GAAG,SAAS;IAoB/I;;;;;OAKG;IACI,WAAW,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,KAAK,EAAE,CAAC;IAWjD;;;;OAIG;IACI,WAAW,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,GAAG,aAAa,CAAC,CAAC,CAAC;IAM1D;;;;OAIG;IACI,cAAc,CAAC,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,GAAG,aAAa,CAAC,CAAC,CAAC;IAYvE,6DAA6D;IACtD,YAAY,CAAC,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,GAAG,OAAO;CAW7D;AACD;;GAEG;AACH,qBAAa,cAAc;IACzB,OAAO,CAAC,OAAO,CAAuB;IACtC,OAAO,CAAC,QAAQ,CAAsB;IACtC,gEAAgE;IAChE,IAAW,SAAS,IAAI,WAAW,CAAC,MAAM,CAAC,CAA0B;IACrE,OAAO;IAIP,sEAAsE;WACxD,MAAM,CAAC,MAAM,EAAE,oBAAoB,EAAE,mBAAmB,EAAE,MAAM,GAAG,cAAc,GAAG,SAAS;IAiB3G;;;OAGG;IACI,qBAAqB,CAAC,KAAK,EAAE,OAAO,GAAG,OAAO,EAAE,QAAQ,EAAE,CAAC,KAAK,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,KAAK,OAAO;CAwB7H"}

package/lib/esm/polyface/multiclip/XYPointBuckets.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 RangeSearch
+ */
+import { Geometry } from "../../Geometry";
+import { Point3d } from "../../geometry3d/Point3dVector3d";
+import { Range2d } from "../../geometry3d/Range";
+/**
+ * Arrays of type T values distributed by xy position when entered.
+ * @internal
+ */
+export class XYIndexGrid {
+    constructor(range, numX, numY) {
+        this._range = range;
+        this._numXEdge = Math.max(numX, 1);
+        this._numYEdge = Math.max(numY, 1);
+        this._data = [];
+        for (let j = 0; j < numY; j++) {
+            const thisRow = [];
+            for (let i = 0; i < numX; i++) {
+                thisRow.push(undefined);
+            }
+            this._data.push(thisRow);
+        }
+    }
+    /** Return the number of x edges in the grid */
+    get numXEdge() { return this._numXEdge; }
+    /** Return the number of y edges in the grid */
+    get numYEdge() { return this._numYEdge; }
+    /** Return the `i` index of cells containing x coordinate */
+    xIndex(x) {
+        const fraction = (x - this._range.low.x) / (this._range.high.x - this._range.low.x);
+        const q = Math.floor(fraction * this._numXEdge);
+        if (q < 0)
+            return 0;
+        if (q > this._numXEdge - 1)
+            return this._numXEdge - 1;
+        return q;
+    }
+    /** Return the `j` index of cells containing x coordinate */
+    yIndex(y) {
+        const fraction = (y - this._range.low.y) / (this._range.high.y - this._range.low.y);
+        const q = Math.floor(fraction * this._numYEdge);
+        if (q < 0)
+            return 0;
+        if (q > this._numYEdge - 1)
+            return this._numYEdge - 1;
+        return q;
+    }
+    /**
+     * Construct an array with cells mapped to a range, with counts determined by estimated total count and target number of entries per cell.
+     * @param range
+     * @param totalEntries
+     * @param targetEntriesPerCell
+     */
+    static createWithEstimatedCounts(range, totalEntries, targetEntriesPerCell) {
+        if (range.low.x >= range.high.x || range.low.y >= range.high.y)
+            return undefined;
+        const range2d = Range2d.createXYXY(range.low.x, range.low.y, range.high.x, range.high.y);
+        const dx = range2d.xLength();
+        const dy = range2d.yLength();
+        // numX / numY = dX / dY  (numX,numY will be integers  . . .)
+        // numX = dX * numY / dY
+        let numX;
+        let numY;
+        if (dy > dx) {
+            numY = Math.ceil(Math.sqrt(dy * totalEntries / (targetEntriesPerCell * dx)));
+            numX = Math.ceil(totalEntries / numY);
+        }
+        else {
+            numX = Math.ceil(Math.sqrt(dx * totalEntries / (targetEntriesPerCell * dy)));
+            numY = Math.ceil(totalEntries / (numX * targetEntriesPerCell));
+        }
+        return new XYIndexGrid(range2d, numX, numY);
+    }
+    /**
+     * Add (save) a new data value to the grid cell containing x,y
+     * @param x
+     * @param y
+     * @param value
+     */
+    addDataAtXY(x, y, value) {
+        const i = this.xIndex(x);
+        const j = this.yIndex(y);
+        let dataJI = this._data[j][i];
+        if (!dataJI) {
+            dataJI = [];
+            this._data[j][i] = dataJI;
+        }
+        dataJI.push(value);
+        return;
+    }
+    /**
+     * Get the (reference to the possibly null array of) data values for the cell indicated by xy.
+     * @param x
+     * @param y
+     */
+    getDataAtXY(x, y) {
+        const i = this.xIndex(x);
+        const j = this.yIndex(y);
+        return this._data[j][i];
+    }
+    /**
+     * Get the (reference to the possibly null array of) data values for the cell indicated by indices in the x and y direction
+     * @param xIndex
+     * @param yIndex
+     */
+    getDataAtIndex(xIndex, yIndex) {
+        if (xIndex < 0)
+            return undefined;
+        if (xIndex >= this._numXEdge)
+            return undefined;
+        if (yIndex < 0)
+            return undefined;
+        if (yIndex >= this._numYEdge)
+            return undefined;
+        return this._data[yIndex][xIndex];
+    }
+    /** Return true if (xIndex, yIndex) is a valid cell index. */
+    isValidIndex(xIndex, yIndex) {
+        if (xIndex < 0)
+            return false;
+        if (xIndex >= this._numXEdge)
+            return false;
+        if (yIndex < 0)
+            return false;
+        if (yIndex >= this._numYEdge)
+            return false;
+        return true;
+    }
+}
+/** Manage buckets of points for fast search.
+ * @internal
+ */
+export class XYPointBuckets {
+    constructor(points, buckets) {
+        this._points = points;
+        this._buckets = buckets;
+    }
+    /** Return the underlying grid with indices recorded by block */
+    get indexGrid() { return this._buckets; }
+    /** Create an XYIndex grid with all indices of all `points` entered */
+    static create(points, targetPointsPerCell) {
+        const n = points.length;
+        if (points.length < 1)
+            return undefined;
+        const range = points.getRange();
+        range.expandInPlace(Geometry.smallMetricDistance * 1000.0);
+        const buckets = XYIndexGrid.createWithEstimatedCounts(range, points.length, targetPointsPerCell);
+        if (buckets === undefined)
+            return undefined;
+        const result = new XYPointBuckets(points, buckets);
+        const point = Point3d.create();
+        for (let i = 0; i < n; i++) {
+            points.getPoint3dAtUncheckedPointIndex(i, point);
+            buckets.addDataAtXY(point.x, point.y, i);
+        }
+        return result;
+    }
+    /** call the `announce` function with the index and coordinates of all points in given range.
+     * * continue the search if `announce` returns true.
+     * * terminate the search if `announce` returns false;
+     */
+    announcePointsInRange(range, announce) {
+        const i0 = this._buckets.xIndex(range.low.x);
+        const i1 = this._buckets.xIndex(range.high.x);
+        const j0 = this._buckets.yIndex(range.low.y);
+        const j1 = this._buckets.yIndex(range.high.y);
+        const n = this._points.length;
+        for (let i = i0; i <= i1; i++) {
+            for (let j = j0; j <= j1; j++) {
+                const candidates = this._buckets.getDataAtIndex(i, j);
+                if (candidates !== undefined) {
+                    for (const k of candidates) {
+                        if (k < n) {
+                            const x = this._points.getXAtUncheckedPointIndex(k);
+                            const y = this._points.getYAtUncheckedPointIndex(k);
+                            const z = this._points.getZAtUncheckedPointIndex(k);
+                            if (range.containsXY(x, y))
+                                if (!announce(k, x, y, z))
+                                    return;
+                        }
+                    }
+                }
+            }
+        }
+    }
+}
+//# sourceMappingURL=XYPointBuckets.js.map

package/lib/esm/polyface/multiclip/XYPointBuckets.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"XYPointBuckets.js","sourceRoot":"","sources":["../../../../src/polyface/multiclip/XYPointBuckets.ts"],"names":[],"mappings":"AAAA;;;+FAG+F;AAC/F;;GAEG;AAEH,OAAO,EAAE,QAAQ,EAAE,MAAM,gBAAgB,CAAC;AAE1C,OAAO,EAAE,OAAO,EAAE,MAAM,kCAAkC,CAAC;AAC3D,OAAO,EAAE,OAAO,EAAW,MAAM,wBAAwB,CAAC;AAK1D;;;GAGG;AACH,MAAM,OAAO,WAAW;IAOtB,YAAsB,KAAc,EAAE,IAAY,EAAE,IAAY;QAC9D,IAAI,CAAC,MAAM,GAAG,KAAK,CAAC;QACpB,IAAI,CAAC,SAAS,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,EAAE,CAAC,CAAC,CAAC;QACnC,IAAI,CAAC,SAAS,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,EAAE,CAAC,CAAC,CAAC;QACnC,IAAI,CAAC,KAAK,GAAG,EAAE,CAAC;QAChB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,EAAE,CAAC,EAAE,EAAE;YAC7B,MAAM,OAAO,GAA4B,EAAE,CAAC;YAC5C,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,EAAE,CAAC,EAAE,EAAE;gBAC7B,OAAO,CAAC,IAAI,CAAC,SAAS,CAAC,CAAC;aACzB;YACD,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,OAAO,CAAC,CAAC;SAC1B;IACH,CAAC;IACD,+CAA+C;IAC/C,IAAW,QAAQ,KAAa,OAAO,IAAI,CAAC,SAAS,CAAC,CAAC,CAAC;IACxD,+CAA+C;IAC/C,IAAW,QAAQ,KAAa,OAAO,IAAI,CAAC,SAAS,CAAC,CAAC,CAAC;IAExD,4DAA4D;IACrD,MAAM,CAAC,CAAS;QACrB,MAAM,QAAQ,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,IAAI,CAAC,MAAM,CAAC,IAAI,CAAC,CAAC,GAAG,IAAI,CAAC,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;QACpF,MAAM,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,QAAQ,GAAG,IAAI,CAAC,SAAS,CAAC,CAAC;QAChD,IAAI,CAAC,GAAG,CAAC;YACP,OAAO,CAAC,CAAC;QACX,IAAI,CAAC,GAAG,IAAI,CAAC,SAAS,GAAG,CAAC;YACxB,OAAO,IAAI,CAAC,SAAS,GAAG,CAAC,CAAC;QAC5B,OAAO,CAAC,CAAC;IACX,CAAC;IACD,4DAA4D;IACrD,MAAM,CAAC,CAAS;QACrB,MAAM,QAAQ,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,IAAI,CAAC,MAAM,CAAC,IAAI,CAAC,CAAC,GAAG,IAAI,CAAC,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;QACpF,MAAM,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,QAAQ,GAAG,IAAI,CAAC,SAAS,CAAC,CAAC;QAChD,IAAI,CAAC,GAAG,CAAC;YACP,OAAO,CAAC,CAAC;QACX,IAAI,CAAC,GAAG,IAAI,CAAC,SAAS,GAAG,CAAC;YACxB,OAAO,IAAI,CAAC,SAAS,GAAG,CAAC,CAAC;QAC5B,OAAO,CAAC,CAAC;IACX,CAAC;IACD;;;;;OAKG;IACI,MAAM,CAAC,yBAAyB,CAAI,KAAmB,EAAE,YAAoB,EAAE,oBAA4B;QAChH,IAAI,KAAK,CAAC,GAAG,CAAC,CAAC,IAAI,KAAK,CAAC,IAAI,CAAC,CAAC,IAAI,KAAK,CAAC,GAAG,CAAC,CAAC,IAAI,KAAK,CAAC,IAAI,CAAC,CAAC;YAC5D,OAAO,SAAS,CAAC;QACnB,MAAM,OAAO,GAAG,OAAO,CAAC,UAAU,CAAC,KAAK,CAAC,GAAG,CAAC,CAAC,EAAE,KAAK,CAAC,GAAG,CAAC,CAAC,EAAE,KAAK,CAAC,IAAI,CAAC,CAAC,EAAE,KAAK,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;QAEzF,MAAM,EAAE,GAAG,OAAO,CAAC,OAAO,EAAE,CAAC;QAC7B,MAAM,EAAE,GAAG,OAAO,CAAC,OAAO,EAAE,CAAC;QAC7B,6DAA6D;QAC7D,wBAAwB;QACxB,IAAI,IAAY,CAAC;QACjB,IAAI,IAAY,CAAC;QACjB,IAAI,EAAE,GAAG,EAAE,EAAE;YACX,IAAI,GAAG,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,EAAE,GAAG,YAAY,GAAG,CAAC,oBAAoB,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC;YAC7E,IAAI,GAAG,IAAI,CAAC,IAAI,CAAC,YAAY,GAAG,IAAI,CAAC,CAAC;SACvC;aAAM;YACL,IAAI,GAAG,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,EAAE,GAAG,YAAY,GAAG,CAAC,oBAAoB,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC;YAC7E,IAAI,GAAG,IAAI,CAAC,IAAI,CAAC,YAAY,GAAG,CAAC,IAAI,GAAG,oBAAoB,CAAC,CAAC,CAAC;SAChE;QACD,OAAO,IAAI,WAAW,CAAC,OAAO,EAAE,IAAI,EAAE,IAAI,CAAC,CAAC;IAC9C,CAAC;IACD;;;;;OAKG;IACI,WAAW,CAAC,CAAS,EAAE,CAAS,EAAE,KAAQ;QAC/C,MAAM,CAAC,GAAG,IAAI,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC;QACzB,MAAM,CAAC,GAAG,IAAI,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC;QACzB,IAAI,MAAM,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QAC9B,IAAI,CAAC,MAAM,EAAE;YACX,MAAM,GAAG,EAAE,CAAC;YACZ,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC;SAC3B;QACD,MAAM,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC;QACnB,OAAO;IACT,CAAC;IACD;;;;OAIG;IACI,WAAW,CAAC,CAAS,EAAE,CAAS;QACrC,MAAM,CAAC,GAAG,IAAI,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC;QACzB,MAAM,CAAC,GAAG,IAAI,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC;QACzB,OAAO,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IAC1B,CAAC;IAED;;;;OAIG;IACI,cAAc,CAAC,MAAc,EAAE,MAAc;QAClD,IAAI,MAAM,GAAG,CAAC;YACZ,OAAO,SAAS,CAAC;QACnB,IAAI,MAAM,IAAI,IAAI,CAAC,SAAS;YAC1B,OAAO,SAAS,CAAC;QACnB,IAAI,MAAM,GAAG,CAAC;YACZ,OAAO,SAAS,CAAC;QACnB,IAAI,MAAM,IAAI,IAAI,CAAC,SAAS;YAC1B,OAAO,SAAS,CAAC;QAEnB,OAAO,IAAI,CAAC,KAAK,CAAC,MAAM,CAAC,CAAC,MAAM,CAAC,CAAC;IACpC,CAAC;IACD,6DAA6D;IACtD,YAAY,CAAC,MAAc,EAAE,MAAc;QAChD,IAAI,MAAM,GAAG,CAAC;YACZ,OAAO,KAAK,CAAC;QACf,IAAI,MAAM,IAAI,IAAI,CAAC,SAAS;YAC1B,OAAO,KAAK,CAAC;QACf,IAAI,MAAM,GAAG,CAAC;YACZ,OAAO,KAAK,CAAC;QACf,IAAI,MAAM,IAAI,IAAI,CAAC,SAAS;YAC1B,OAAO,KAAK,CAAC;QACf,OAAO,IAAI,CAAC;IACd,CAAC;CACF;AACD;;GAEG;AACH,MAAM,OAAO,cAAc;IAKzB,YAAoB,MAA4B,EAAE,OAA4B;QAC5E,IAAI,CAAC,OAAO,GAAG,MAAM,CAAC;QACtB,IAAI,CAAC,QAAQ,GAAG,OAAO,CAAC;IAC1B,CAAC;IALD,gEAAgE;IAChE,IAAW,SAAS,KAA0B,OAAO,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC;IAKrE,sEAAsE;IAC/D,MAAM,CAAC,MAAM,CAAC,MAA4B,EAAE,mBAA2B;QAC5E,MAAM,CAAC,GAAG,MAAM,CAAC,MAAM,CAAC;QACxB,IAAI,MAAM,CAAC,MAAM,GAAG,CAAC;YACnB,OAAO,SAAS,CAAC;QACnB,MAAM,KAAK,GAAG,MAAM,CAAC,QAAQ,EAAE,CAAC;QAChC,KAAK,CAAC,aAAa,CAAC,QAAQ,CAAC,mBAAmB,GAAG,MAAM,CAAC,CAAC;QAC3D,MAAM,OAAO,GAAG,WAAW,CAAC,yBAAyB,CAAS,KAAK,EAAE,MAAM,CAAC,MAAM,EAAE,mBAAmB,CAAC,CAAC;QACzG,IAAI,OAAO,KAAK,SAAS;YACvB,OAAO,SAAS,CAAC;QACnB,MAAM,MAAM,GAAG,IAAI,cAAc,CAAC,MAAM,EAAE,OAAO,CAAC,CAAC;QACnD,MAAM,KAAK,GAAG,OAAO,CAAC,MAAM,EAAE,CAAC;QAC/B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE;YAC1B,MAAM,CAAC,+BAA+B,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC;YACjD,OAAO,CAAC,WAAW,CAAC,KAAK,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;SAC1C;QACD,OAAO,MAAM,CAAC;IAChB,CAAC;IACD;;;OAGG;IACI,qBAAqB,CAAC,KAAwB,EAAE,QAAqE;QAC1H,MAAM,EAAE,GAAG,IAAI,CAAC,QAAQ,CAAC,MAAM,CAAC,KAAK,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;QAC7C,MAAM,EAAE,GAAG,IAAI,CAAC,QAAQ,CAAC,MAAM,CAAC,KAAK,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;QAC9C,MAAM,EAAE,GAAG,IAAI,CAAC,QAAQ,CAAC,MAAM,CAAC,KAAK,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;QAC7C,MAAM,EAAE,GAAG,IAAI,CAAC,QAAQ,CAAC,MAAM,CAAC,KAAK,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;QAC9C,MAAM,CAAC,GAAG,IAAI,CAAC,OAAO,CAAC,MAAM,CAAC;QAC9B,KAAK,IAAI,CAAC,GAAG,EAAE,EAAE,CAAC,IAAI,EAAE,EAAE,CAAC,EAAE,EAAE;YAC7B,KAAK,IAAI,CAAC,GAAG,EAAE,EAAE,CAAC,IAAI,EAAE,EAAE,CAAC,EAAE,EAAE;gBAC7B,MAAM,UAAU,GAAG,IAAI,CAAC,QAAQ,CAAC,cAAc,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;gBACtD,IAAI,UAAU,KAAK,SAAS,EAAE;oBAC5B,KAAK,MAAM,CAAC,IAAI,UAAU,EAAE;wBAC1B,IAAI,CAAC,GAAG,CAAC,EAAE;4BACT,MAAM,CAAC,GAAG,IAAI,CAAC,OAAO,CAAC,yBAAyB,CAAC,CAAC,CAAC,CAAC;4BACpD,MAAM,CAAC,GAAG,IAAI,CAAC,OAAO,CAAC,yBAAyB,CAAC,CAAC,CAAC,CAAC;4BACpD,MAAM,CAAC,GAAG,IAAI,CAAC,OAAO,CAAC,yBAAyB,CAAC,CAAC,CAAC,CAAC;4BACpD,IAAI,KAAK,CAAC,UAAU,CAAC,CAAC,EAAE,CAAC,CAAC;gCACxB,IAAI,CAAC,QAAQ,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;oCACvB,OAAO;yBACZ;qBACF;iBACF;aACF;SACF;IACH,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 RangeSearch\r\n */\r\n\r\nimport { Geometry } from \"../../Geometry\";\r\nimport { IndexedXYZCollection } from \"../../geometry3d/IndexedXYZCollection\";\r\nimport { Point3d } from \"../../geometry3d/Point3dVector3d\";\r\nimport { Range2d, Range3d } from \"../../geometry3d/Range\";\r\nimport { LowAndHighXY } from \"../../geometry3d/XYZProps\";\r\n\r\n/** Type for a value which may be either (a) undefined or (b) an array of type []. */\r\nexport type OptionalArray<T> = T[] | undefined;\r\n/**\r\n * Arrays of type T values distributed by xy position when entered.\r\n * @internal\r\n */\r\nexport class XYIndexGrid<T> {\r\n  private _range: Range2d;\r\n  private _numXEdge: number;\r\n  private _numYEdge: number;\r\n\r\n  private _data: Array<Array<OptionalArray<T>>>; // _indices[j][i][k] = index k value in row j block i\r\n\r\n  protected constructor(range: Range2d, numX: number, numY: number) {\r\n    this._range = range;\r\n    this._numXEdge = Math.max(numX, 1);\r\n    this._numYEdge = Math.max(numY, 1);\r\n    this._data = [];\r\n    for (let j = 0; j < numY; j++) {\r\n      const thisRow: Array<OptionalArray<T>> = [];\r\n      for (let i = 0; i < numX; i++) {\r\n        thisRow.push(undefined);\r\n      }\r\n      this._data.push(thisRow);\r\n    }\r\n  }\r\n  /** Return the number of x edges in the grid */\r\n  public get numXEdge(): number { return this._numXEdge; }\r\n  /** Return the number of y edges in the grid */\r\n  public get numYEdge(): number { return this._numYEdge; }\r\n\r\n  /** Return the `i` index of cells containing x coordinate */\r\n  public xIndex(x: number): number {\r\n    const fraction = (x - this._range.low.x) / (this._range.high.x - this._range.low.x);\r\n    const q = Math.floor(fraction * this._numXEdge);\r\n    if (q < 0)\r\n      return 0;\r\n    if (q > this._numXEdge - 1)\r\n      return this._numXEdge - 1;\r\n    return q;\r\n  }\r\n  /** Return the `j` index of cells containing x coordinate */\r\n  public yIndex(y: number): number {\r\n    const fraction = (y - this._range.low.y) / (this._range.high.y - this._range.low.y);\r\n    const q = Math.floor(fraction * this._numYEdge);\r\n    if (q < 0)\r\n      return 0;\r\n    if (q > this._numYEdge - 1)\r\n      return this._numYEdge - 1;\r\n    return q;\r\n  }\r\n  /**\r\n   * Construct an array with cells mapped to a range, with counts determined by estimated total count and target number of entries per cell.\r\n   * @param range\r\n   * @param totalEntries\r\n   * @param targetEntriesPerCell\r\n   */\r\n  public static createWithEstimatedCounts<T>(range: LowAndHighXY, totalEntries: number, targetEntriesPerCell: number): XYIndexGrid<T> | undefined {\r\n    if (range.low.x >= range.high.x || range.low.y >= range.high.y)\r\n      return undefined;\r\n    const range2d = Range2d.createXYXY(range.low.x, range.low.y, range.high.x, range.high.y);\r\n\r\n    const dx = range2d.xLength();\r\n    const dy = range2d.yLength();\r\n    // numX / numY = dX / dY  (numX,numY will be integers  . . .)\r\n    // numX = dX * numY / dY\r\n    let numX: number;\r\n    let numY: number;\r\n    if (dy > dx) {\r\n      numY = Math.ceil(Math.sqrt(dy * totalEntries / (targetEntriesPerCell * dx)));\r\n      numX = Math.ceil(totalEntries / numY);\r\n    } else {\r\n      numX = Math.ceil(Math.sqrt(dx * totalEntries / (targetEntriesPerCell * dy)));\r\n      numY = Math.ceil(totalEntries / (numX * targetEntriesPerCell));\r\n    }\r\n    return new XYIndexGrid(range2d, numX, numY);\r\n  }\r\n  /**\r\n   * Add (save) a new data value to the grid cell containing x,y\r\n   * @param x\r\n   * @param y\r\n   * @param value\r\n   */\r\n  public addDataAtXY(x: number, y: number, value: T) {\r\n    const i = this.xIndex(x);\r\n    const j = this.yIndex(y);\r\n    let dataJI = this._data[j][i];\r\n    if (!dataJI) {\r\n      dataJI = [];\r\n      this._data[j][i] = dataJI;\r\n    }\r\n    dataJI.push(value);\r\n    return;\r\n  }\r\n  /**\r\n   * Get the (reference to the possibly null array of) data values for the cell indicated by xy.\r\n   * @param x\r\n   * @param y\r\n   */\r\n  public getDataAtXY(x: number, y: number): OptionalArray<T> {\r\n    const i = this.xIndex(x);\r\n    const j = this.yIndex(y);\r\n    return this._data[j][i];\r\n  }\r\n\r\n  /**\r\n   * Get the (reference to the possibly null array of) data values for the cell indicated by indices in the x and y direction\r\n   * @param xIndex\r\n   * @param yIndex\r\n   */\r\n  public getDataAtIndex(xIndex: number, yIndex: number): OptionalArray<T> {\r\n    if (xIndex < 0)\r\n      return undefined;\r\n    if (xIndex >= this._numXEdge)\r\n      return undefined;\r\n    if (yIndex < 0)\r\n      return undefined;\r\n    if (yIndex >= this._numYEdge)\r\n      return undefined;\r\n\r\n    return this._data[yIndex][xIndex];\r\n  }\r\n  /** Return true if (xIndex, yIndex) is a valid cell index. */\r\n  public isValidIndex(xIndex: number, yIndex: number): boolean {\r\n    if (xIndex < 0)\r\n      return false;\r\n    if (xIndex >= this._numXEdge)\r\n      return false;\r\n    if (yIndex < 0)\r\n      return false;\r\n    if (yIndex >= this._numYEdge)\r\n      return false;\r\n    return true;\r\n  }\r\n}\r\n/** Manage buckets of points for fast search.\r\n * @internal\r\n */\r\nexport class XYPointBuckets {\r\n  private _points: IndexedXYZCollection;\r\n  private _buckets: XYIndexGrid<number>;\r\n  /** Return the underlying grid with indices recorded by block */\r\n  public get indexGrid(): XYIndexGrid<number> { return this._buckets; }\r\n  private constructor(points: IndexedXYZCollection, buckets: XYIndexGrid<number>) {\r\n    this._points = points;\r\n    this._buckets = buckets;\r\n  }\r\n  /** Create an XYIndex grid with all indices of all `points` entered */\r\n  public static create(points: IndexedXYZCollection, targetPointsPerCell: number): XYPointBuckets | undefined {\r\n    const n = points.length;\r\n    if (points.length < 1)\r\n      return undefined;\r\n    const range = points.getRange();\r\n    range.expandInPlace(Geometry.smallMetricDistance * 1000.0);\r\n    const buckets = XYIndexGrid.createWithEstimatedCounts<number>(range, points.length, targetPointsPerCell);\r\n    if (buckets === undefined)\r\n      return undefined;\r\n    const result = new XYPointBuckets(points, buckets);\r\n    const point = Point3d.create();\r\n    for (let i = 0; i < n; i++) {\r\n      points.getPoint3dAtUncheckedPointIndex(i, point);\r\n      buckets.addDataAtXY(point.x, point.y, i);\r\n    }\r\n    return result;\r\n  }\r\n  /** call the `announce` function with the index and coordinates of all points in given range.\r\n   * * continue the search if `announce` returns true.\r\n   * * terminate the search if `announce` returns false;\r\n   */\r\n  public announcePointsInRange(range: Range2d | Range3d, announce: (index: number, x: number, y: number, z: number) => boolean) {\r\n    const i0 = this._buckets.xIndex(range.low.x);\r\n    const i1 = this._buckets.xIndex(range.high.x);\r\n    const j0 = this._buckets.yIndex(range.low.y);\r\n    const j1 = this._buckets.yIndex(range.high.y);\r\n    const n = this._points.length;\r\n    for (let i = i0; i <= i1; i++) {\r\n      for (let j = j0; j <= j1; j++) {\r\n        const candidates = this._buckets.getDataAtIndex(i, j);\r\n        if (candidates !== undefined) {\r\n          for (const k of candidates) {\r\n            if (k < n) {\r\n              const x = this._points.getXAtUncheckedPointIndex(k);\r\n              const y = this._points.getYAtUncheckedPointIndex(k);\r\n              const z = this._points.getZAtUncheckedPointIndex(k);\r\n              if (range.containsXY(x, y))\r\n                if (!announce(k, x, y, z))\r\n                  return;\r\n            }\r\n          }\r\n        }\r\n      }\r\n    }\r\n  }\r\n}\r\n"]}