@itwin/core-geometry 5.2.0-dev.7 → 5.2.0

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

@@ -1 +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;IACd,MAAM,CAAU;IAChB,SAAS,CAAS;IAClB,SAAS,CAAS;IAElB,KAAK,CAAiC,CAAC,qDAAqD;IAEpG,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,CAAC;YAC9B,MAAM,OAAO,GAA4B,EAAE,CAAC;YAC5C,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,EAAE,CAAC,EAAE,EAAE,CAAC;gBAC9B,OAAO,CAAC,IAAI,CAAC,SAAS,CAAC,CAAC;YAC1B,CAAC;YACD,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,OAAO,CAAC,CAAC;QAC3B,CAAC;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,CAAC;YACZ,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;QACxC,CAAC;aAAM,CAAC;YACN,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;QACjE,CAAC;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,CAAC;YACZ,MAAM,GAAG,EAAE,CAAC;YACZ,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC;QAC5B,CAAC;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;IACjB,OAAO,CAAuB;IAC9B,QAAQ,CAAsB;IACtC,gEAAgE;IAChE,IAAW,SAAS,KAA0B,OAAO,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC;IACrE,YAAoB,MAA4B,EAAE,OAA4B;QAC5E,IAAI,CAAC,OAAO,GAAG,MAAM,CAAC;QACtB,IAAI,CAAC,QAAQ,GAAG,OAAO,CAAC;IAC1B,CAAC;IACD,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,CAAC;YAC3B,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;QAC3C,CAAC;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,CAAC;YAC9B,KAAK,IAAI,CAAC,GAAG,EAAE,EAAE,CAAC,IAAI,EAAE,EAAE,CAAC,EAAE,EAAE,CAAC;gBAC9B,MAAM,UAAU,GAAG,IAAI,CAAC,QAAQ,CAAC,cAAc,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;gBACtD,IAAI,UAAU,KAAK,SAAS,EAAE,CAAC;oBAC7B,KAAK,MAAM,CAAC,IAAI,UAAU,EAAE,CAAC;wBAC3B,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC;4BACV,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;wBACb,CAAC;oBACH,CAAC;gBACH,CAAC;YACH,CAAC;QACH,CAAC;IACH,CAAC;CACF","sourcesContent":["/*---------------------------------------------------------------------------------------------\n* Copyright (c) Bentley Systems, Incorporated. All rights reserved.\n* See LICENSE.md in the project root for license terms and full copyright notice.\n*--------------------------------------------------------------------------------------------*/\n/** @packageDocumentation\n * @module RangeSearch\n */\n\nimport { Geometry } from \"../../Geometry\";\nimport { IndexedXYZCollection } from \"../../geometry3d/IndexedXYZCollection\";\nimport { Point3d } from \"../../geometry3d/Point3dVector3d\";\nimport { Range2d, Range3d } from \"../../geometry3d/Range\";\nimport { LowAndHighXY } from \"../../geometry3d/XYZProps\";\n\n/** Type for a value which may be either (a) undefined or (b) an array of type []. */\nexport type OptionalArray<T> = T[] | undefined;\n/**\n * Arrays of type T values distributed by xy position when entered.\n * @internal\n */\nexport class XYIndexGrid<T> {\n  private _range: Range2d;\n  private _numXEdge: number;\n  private _numYEdge: number;\n\n  private _data: Array<Array<OptionalArray<T>>>; // _indices[j][i][k] = index k value in row j block i\n\n  protected constructor(range: Range2d, numX: number, numY: number) {\n    this._range = range;\n    this._numXEdge = Math.max(numX, 1);\n    this._numYEdge = Math.max(numY, 1);\n    this._data = [];\n    for (let j = 0; j < numY; j++) {\n      const thisRow: Array<OptionalArray<T>> = [];\n      for (let i = 0; i < numX; i++) {\n        thisRow.push(undefined);\n      }\n      this._data.push(thisRow);\n    }\n  }\n  /** Return the number of x edges in the grid */\n  public get numXEdge(): number { return this._numXEdge; }\n  /** Return the number of y edges in the grid */\n  public get numYEdge(): number { return this._numYEdge; }\n\n  /** Return the `i` index of cells containing x coordinate */\n  public xIndex(x: number): number {\n    const fraction = (x - this._range.low.x) / (this._range.high.x - this._range.low.x);\n    const q = Math.floor(fraction * this._numXEdge);\n    if (q < 0)\n      return 0;\n    if (q > this._numXEdge - 1)\n      return this._numXEdge - 1;\n    return q;\n  }\n  /** Return the `j` index of cells containing x coordinate */\n  public yIndex(y: number): number {\n    const fraction = (y - this._range.low.y) / (this._range.high.y - this._range.low.y);\n    const q = Math.floor(fraction * this._numYEdge);\n    if (q < 0)\n      return 0;\n    if (q > this._numYEdge - 1)\n      return this._numYEdge - 1;\n    return q;\n  }\n  /**\n   * Construct an array with cells mapped to a range, with counts determined by estimated total count and target number of entries per cell.\n   * @param range\n   * @param totalEntries\n   * @param targetEntriesPerCell\n   */\n  public static createWithEstimatedCounts<T>(range: LowAndHighXY, totalEntries: number, targetEntriesPerCell: number): XYIndexGrid<T> | undefined {\n    if (range.low.x >= range.high.x || range.low.y >= range.high.y)\n      return undefined;\n    const range2d = Range2d.createXYXY(range.low.x, range.low.y, range.high.x, range.high.y);\n\n    const dx = range2d.xLength();\n    const dy = range2d.yLength();\n    // numX / numY = dX / dY  (numX,numY will be integers  . . .)\n    // numX = dX * numY / dY\n    let numX: number;\n    let numY: number;\n    if (dy > dx) {\n      numY = Math.ceil(Math.sqrt(dy * totalEntries / (targetEntriesPerCell * dx)));\n      numX = Math.ceil(totalEntries / numY);\n    } else {\n      numX = Math.ceil(Math.sqrt(dx * totalEntries / (targetEntriesPerCell * dy)));\n      numY = Math.ceil(totalEntries / (numX * targetEntriesPerCell));\n    }\n    return new XYIndexGrid(range2d, numX, numY);\n  }\n  /**\n   * Add (save) a new data value to the grid cell containing x,y\n   * @param x\n   * @param y\n   * @param value\n   */\n  public addDataAtXY(x: number, y: number, value: T) {\n    const i = this.xIndex(x);\n    const j = this.yIndex(y);\n    let dataJI = this._data[j][i];\n    if (!dataJI) {\n      dataJI = [];\n      this._data[j][i] = dataJI;\n    }\n    dataJI.push(value);\n    return;\n  }\n  /**\n   * Get the (reference to the possibly null array of) data values for the cell indicated by xy.\n   * @param x\n   * @param y\n   */\n  public getDataAtXY(x: number, y: number): OptionalArray<T> {\n    const i = this.xIndex(x);\n    const j = this.yIndex(y);\n    return this._data[j][i];\n  }\n\n  /**\n   * Get the (reference to the possibly null array of) data values for the cell indicated by indices in the x and y direction\n   * @param xIndex\n   * @param yIndex\n   */\n  public getDataAtIndex(xIndex: number, yIndex: number): OptionalArray<T> {\n    if (xIndex < 0)\n      return undefined;\n    if (xIndex >= this._numXEdge)\n      return undefined;\n    if (yIndex < 0)\n      return undefined;\n    if (yIndex >= this._numYEdge)\n      return undefined;\n\n    return this._data[yIndex][xIndex];\n  }\n  /** Return true if (xIndex, yIndex) is a valid cell index. */\n  public isValidIndex(xIndex: number, yIndex: number): boolean {\n    if (xIndex < 0)\n      return false;\n    if (xIndex >= this._numXEdge)\n      return false;\n    if (yIndex < 0)\n      return false;\n    if (yIndex >= this._numYEdge)\n      return false;\n    return true;\n  }\n}\n/** Manage buckets of points for fast search.\n * @internal\n */\nexport class XYPointBuckets {\n  private _points: IndexedXYZCollection;\n  private _buckets: XYIndexGrid<number>;\n  /** Return the underlying grid with indices recorded by block */\n  public get indexGrid(): XYIndexGrid<number> { return this._buckets; }\n  private constructor(points: IndexedXYZCollection, buckets: XYIndexGrid<number>) {\n    this._points = points;\n    this._buckets = buckets;\n  }\n  /** Create an XYIndex grid with all indices of all `points` entered */\n  public static create(points: IndexedXYZCollection, targetPointsPerCell: number): XYPointBuckets | undefined {\n    const n = points.length;\n    if (points.length < 1)\n      return undefined;\n    const range = points.getRange();\n    range.expandInPlace(Geometry.smallMetricDistance * 1000.0);\n    const buckets = XYIndexGrid.createWithEstimatedCounts<number>(range, points.length, targetPointsPerCell);\n    if (buckets === undefined)\n      return undefined;\n    const result = new XYPointBuckets(points, buckets);\n    const point = Point3d.create();\n    for (let i = 0; i < n; i++) {\n      points.getPoint3dAtUncheckedPointIndex(i, point);\n      buckets.addDataAtXY(point.x, point.y, i);\n    }\n    return result;\n  }\n  /** call the `announce` function with the index and coordinates of all points in given range.\n   * * continue the search if `announce` returns true.\n   * * terminate the search if `announce` returns false;\n   */\n  public announcePointsInRange(range: Range2d | Range3d, announce: (index: number, x: number, y: number, z: number) => boolean) {\n    const i0 = this._buckets.xIndex(range.low.x);\n    const i1 = this._buckets.xIndex(range.high.x);\n    const j0 = this._buckets.yIndex(range.low.y);\n    const j1 = this._buckets.yIndex(range.high.y);\n    const n = this._points.length;\n    for (let i = i0; i <= i1; i++) {\n      for (let j = j0; j <= j1; j++) {\n        const candidates = this._buckets.getDataAtIndex(i, j);\n        if (candidates !== undefined) {\n          for (const k of candidates) {\n            if (k < n) {\n              const x = this._points.getXAtUncheckedPointIndex(k);\n              const y = this._points.getYAtUncheckedPointIndex(k);\n              const z = this._points.getZAtUncheckedPointIndex(k);\n              if (range.containsXY(x, y))\n                if (!announce(k, x, y, z))\n                  return;\n            }\n          }\n        }\n      }\n    }\n  }\n}\n"]}
+{"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;IACd,MAAM,CAAU;IAChB,SAAS,CAAS;IAClB,SAAS,CAAS;IAElB,KAAK,CAAiC,CAAC,qDAAqD;IAEpG,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,CAAC;YAC9B,MAAM,OAAO,GAA4B,EAAE,CAAC;YAC5C,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,EAAE,CAAC,EAAE,EAAE,CAAC;gBAC9B,OAAO,CAAC,IAAI,CAAC,SAAS,CAAC,CAAC;YAC1B,CAAC;YACD,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,OAAO,CAAC,CAAC;QAC3B,CAAC;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,CAAC;YACZ,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;QACxC,CAAC;aAAM,CAAC;YACN,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;QACjE,CAAC;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,CAAC;YACZ,MAAM,GAAG,EAAE,CAAC;YACZ,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC;QAC5B,CAAC;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;IACjB,OAAO,CAAuB;IAC9B,QAAQ,CAAsB;IACtC,gEAAgE;IAChE,IAAW,SAAS,KAA0B,OAAO,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC;IACrE,YAAoB,MAA4B,EAAE,OAA4B;QAC5E,IAAI,CAAC,OAAO,GAAG,MAAM,CAAC;QACtB,IAAI,CAAC,QAAQ,GAAG,OAAO,CAAC;IAC1B,CAAC;IACD,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,CAAC;YAC3B,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;QAC3C,CAAC;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,CAAC;YAC9B,KAAK,IAAI,CAAC,GAAG,EAAE,EAAE,CAAC,IAAI,EAAE,EAAE,CAAC,EAAE,EAAE,CAAC;gBAC9B,MAAM,UAAU,GAAG,IAAI,CAAC,QAAQ,CAAC,cAAc,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;gBACtD,IAAI,UAAU,KAAK,SAAS,EAAE,CAAC;oBAC7B,KAAK,MAAM,CAAC,IAAI,UAAU,EAAE,CAAC;wBAC3B,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC;4BACV,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;wBACb,CAAC;oBACH,CAAC;gBACH,CAAC;YACH,CAAC;QACH,CAAC;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"]}