@math.gl/geometry 5.0.0-alpha.2 → 5.0.0-alpha.20

package/dist/culling/lib/shapes/implicit-shape.d.ts

@@ -0,0 +1,56 @@
+import { Matrix4, Vector3 } from '@math.gl/core';
+import type { CullingResult } from "../../constants.js";
+import { AxisAlignedBoundingBox } from "../bounding-volumes/axis-aligned-bounding-box.js";
+import { BoundingSphere } from "../bounding-volumes/bounding-sphere.js";
+import type { BoundingVolume } from "../bounding-volumes/bounding-volume.js";
+import type { Plane } from "../plane.js";
+import type { Ray } from "../ray.js";
+export type RayIntersection = {
+    distance: number;
+    point: Vector3;
+    normal: Vector3;
+    entering: boolean;
+};
+export interface ImplicitShape extends BoundingVolume {
+    readonly type: 'box' | 'capsule' | 'cylinder' | 'plane' | 'sphere';
+    readonly matrix: Matrix4;
+    clone(): ImplicitShape;
+    equals(shape: ImplicitShape): boolean;
+    containsPoint(point: readonly number[]): boolean;
+    intersectRay(ray: Ray): RayIntersection | undefined;
+    getAxisAlignedBoundingBox(): AxisAlignedBoundingBox | undefined;
+    getBoundingSphere(): BoundingSphere | undefined;
+}
+export type LocalRayIntersection = {
+    distance: number;
+    normal: readonly number[];
+};
+/** Shared affine transform, support mapping and query implementation for finite convex shapes. */
+export declare abstract class ImplicitShapeBase implements ImplicitShape {
+    abstract readonly type: ImplicitShape['type'];
+    readonly matrix: Matrix4;
+    protected inverseMatrix: Matrix4;
+    constructor(matrix?: readonly number[]);
+    abstract clone(): ImplicitShapeBase;
+    abstract containsLocalPoint(point: readonly number[]): boolean;
+    abstract distanceToLocalPoint(point: readonly number[]): number;
+    abstract supportLocal(direction: readonly number[]): readonly number[];
+    abstract intersectLocalRay(origin: readonly number[], direction: readonly number[]): LocalRayIntersection | undefined;
+    equals(shape: ImplicitShape): boolean;
+    transform(transform: readonly number[]): this;
+    containsPoint(point: readonly number[]): boolean;
+    distanceSquaredTo(point: readonly number[]): number;
+    /** Exact for rigid/uniform transforms; conservative estimate for non-uniform scale or shear. */
+    distanceTo(point: readonly number[]): number;
+    intersectPlane(plane: Plane): CullingResult;
+    intersectRay(ray: Ray): RayIntersection | undefined;
+    getAxisAlignedBoundingBox(): AxisAlignedBoundingBox | undefined;
+    getBoundingSphere(): BoundingSphere | undefined;
+    protected getSupportPoint(worldDirection: readonly number[]): Vector3;
+}
+export declare function validateAffineMatrix(matrix: readonly number[]): void;
+export declare function solveQuadratic(a: number, b: number, c: number): number[];
+export declare function closestDistanceToSegment2D(px: number, py: number, ax: number, ay: number, bx: number, by: number): number;
+/** Distance in the radial/Y half-plane to a spherical profile arc selected by normal Y range. */
+export declare function closestDistanceToCircleArc2D(px: number, py: number, centerY: number, radius: number, minimumNormalY: number, maximumNormalY: number): number;
+//# sourceMappingURL=implicit-shape.d.ts.map

package/dist/culling/lib/shapes/implicit-shape.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"implicit-shape.d.ts","sourceRoot":"","sources":["../../../../src/culling/lib/shapes/implicit-shape.ts"],"names":[],"mappings":"AAIA,OAAO,EAAC,OAAO,EAAE,OAAO,EAAC,MAAM,eAAe,CAAC;AAC/C,OAAO,KAAK,EAAC,aAAa,EAAC,2BAAwB;AACnD,OAAO,EAAC,sBAAsB,EAAC,yDAAsD;AACrF,OAAO,EAAC,cAAc,EAAC,+CAA4C;AACnE,OAAO,KAAK,EAAC,cAAc,EAAC,+CAA4C;AACxE,OAAO,KAAK,EAAC,KAAK,EAAC,oBAAiB;AACpC,OAAO,KAAK,EAAC,GAAG,EAAC,kBAAe;AAEhC,MAAM,MAAM,eAAe,GAAG;IAC5B,QAAQ,EAAE,MAAM,CAAC;IACjB,KAAK,EAAE,OAAO,CAAC;IACf,MAAM,EAAE,OAAO,CAAC;IAChB,QAAQ,EAAE,OAAO,CAAC;CACnB,CAAC;AAEF,MAAM,WAAW,aAAc,SAAQ,cAAc;IACnD,QAAQ,CAAC,IAAI,EAAE,KAAK,GAAG,SAAS,GAAG,UAAU,GAAG,OAAO,GAAG,QAAQ,CAAC;IACnE,QAAQ,CAAC,MAAM,EAAE,OAAO,CAAC;IACzB,KAAK,IAAI,aAAa,CAAC;IACvB,MAAM,CAAC,KAAK,EAAE,aAAa,GAAG,OAAO,CAAC;IACtC,aAAa,CAAC,KAAK,EAAE,SAAS,MAAM,EAAE,GAAG,OAAO,CAAC;IACjD,YAAY,CAAC,GAAG,EAAE,GAAG,GAAG,eAAe,GAAG,SAAS,CAAC;IACpD,yBAAyB,IAAI,sBAAsB,GAAG,SAAS,CAAC;IAChE,iBAAiB,IAAI,cAAc,GAAG,SAAS,CAAC;CACjD;AAED,MAAM,MAAM,oBAAoB,GAAG;IAAC,QAAQ,EAAE,MAAM,CAAC;IAAC,MAAM,EAAE,SAAS,MAAM,EAAE,CAAA;CAAC,CAAC;AAEjF,kGAAkG;AAClG,8BAAsB,iBAAkB,YAAW,aAAa;IAC9D,QAAQ,CAAC,QAAQ,CAAC,IAAI,EAAE,aAAa,CAAC,MAAM,CAAC,CAAC;IAC9C,QAAQ,CAAC,MAAM,EAAE,OAAO,CAAC;IACzB,SAAS,CAAC,aAAa,EAAE,OAAO,CAAC;gBAErB,MAAM,CAAC,EAAE,SAAS,MAAM,EAAE;IAOtC,QAAQ,CAAC,KAAK,IAAI,iBAAiB;IACnC,QAAQ,CAAC,kBAAkB,CAAC,KAAK,EAAE,SAAS,MAAM,EAAE,GAAG,OAAO;IAC9D,QAAQ,CAAC,oBAAoB,CAAC,KAAK,EAAE,SAAS,MAAM,EAAE,GAAG,MAAM;IAC/D,QAAQ,CAAC,YAAY,CAAC,SAAS,EAAE,SAAS,MAAM,EAAE,GAAG,SAAS,MAAM,EAAE;IACtE,QAAQ,CAAC,iBAAiB,CACxB,MAAM,EAAE,SAAS,MAAM,EAAE,EACzB,SAAS,EAAE,SAAS,MAAM,EAAE,GAC3B,oBAAoB,GAAG,SAAS;IAEnC,MAAM,CAAC,KAAK,EAAE,aAAa,GAAG,OAAO;IAMrC,SAAS,CAAC,SAAS,EAAE,SAAS,MAAM,EAAE,GAAG,IAAI;IAS7C,aAAa,CAAC,KAAK,EAAE,SAAS,MAAM,EAAE,GAAG,OAAO;IAIhD,iBAAiB,CAAC,KAAK,EAAE,SAAS,MAAM,EAAE,GAAG,MAAM;IAKnD,gGAAgG;IAChG,UAAU,CAAC,KAAK,EAAE,SAAS,MAAM,EAAE,GAAG,MAAM;IAO5C,cAAc,CAAC,KAAK,EAAE,KAAK,GAAG,aAAa;IAU3C,YAAY,CAAC,GAAG,EAAE,GAAG,GAAG,eAAe,GAAG,SAAS;IAkBnD,yBAAyB,IAAI,sBAAsB,GAAG,SAAS;IAa/D,iBAAiB,IAAI,cAAc,GAAG,SAAS;IAM/C,SAAS,CAAC,eAAe,CAAC,cAAc,EAAE,SAAS,MAAM,EAAE,GAAG,OAAO;CAStE;AAED,wBAAgB,oBAAoB,CAAC,MAAM,EAAE,SAAS,MAAM,EAAE,GAAG,IAAI,CAcpE;AAED,wBAAgB,cAAc,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,GAAG,MAAM,EAAE,CAMxE;AAED,wBAAgB,0BAA0B,CACxC,EAAE,EAAE,MAAM,EACV,EAAE,EAAE,MAAM,EACV,EAAE,EAAE,MAAM,EACV,EAAE,EAAE,MAAM,EACV,EAAE,EAAE,MAAM,EACV,EAAE,EAAE,MAAM,GACT,MAAM,CAQR;AAED,iGAAiG;AACjG,wBAAgB,4BAA4B,CAC1C,EAAE,EAAE,MAAM,EACV,EAAE,EAAE,MAAM,EACV,OAAO,EAAE,MAAM,EACf,MAAM,EAAE,MAAM,EACd,cAAc,EAAE,MAAM,EACtB,cAAc,EAAE,MAAM,GACrB,MAAM,CAeR"}

package/dist/culling/lib/shapes/implicit-shape.js

@@ -0,0 +1,142 @@
+// math.gl
+// SPDX-License-Identifier: MIT
+// SPDX-FileCopyrightText: Copyright (c) vis.gl contributors
+import { Matrix4, Vector3 } from '@math.gl/core';
+import { AxisAlignedBoundingBox } from "../bounding-volumes/axis-aligned-bounding-box.js";
+import { BoundingSphere } from "../bounding-volumes/bounding-sphere.js";
+/** Shared affine transform, support mapping and query implementation for finite convex shapes. */
+export class ImplicitShapeBase {
+    constructor(matrix) {
+        this.matrix = new Matrix4();
+        if (matrix)
+            this.matrix.copy(matrix);
+        validateAffineMatrix(this.matrix);
+        this.inverseMatrix = this.matrix.clone().invert();
+    }
+    equals(shape) {
+        if (this.type !== shape.type)
+            return false;
+        for (let i = 0; i < 16; i++)
+            if (this.matrix[i] !== shape.matrix[i])
+                return false;
+        return true;
+    }
+    transform(transform) {
+        const next = new Matrix4().copy(transform);
+        validateAffineMatrix(next);
+        this.matrix.multiplyLeft(next);
+        validateAffineMatrix(this.matrix);
+        this.inverseMatrix.copy(this.matrix).invert();
+        return this;
+    }
+    containsPoint(point) {
+        return this.containsLocalPoint(this.inverseMatrix.transformAsPoint(point));
+    }
+    distanceSquaredTo(point) {
+        const distance = this.distanceTo(point);
+        return distance * distance;
+    }
+    /** Exact for rigid/uniform transforms; conservative estimate for non-uniform scale or shear. */
+    distanceTo(point) {
+        const localPoint = this.inverseMatrix.transformAsPoint(point);
+        const localDistance = this.distanceToLocalPoint(localPoint);
+        const scales = this.matrix.getScale();
+        return localDistance * Math.min(scales[0], scales[1], scales[2]);
+    }
+    intersectPlane(plane) {
+        const maximum = this.getSupportPoint(plane.normal);
+        const minimum = this.getSupportPoint([-plane.normal[0], -plane.normal[1], -plane.normal[2]]);
+        const maxDistance = plane.getPointDistance(maximum);
+        const minDistance = plane.getPointDistance(minimum);
+        if (minDistance > 0)
+            return 'inside';
+        if (maxDistance < 0)
+            return 'outside';
+        return 'intersecting';
+    }
+    intersectRay(ray) {
+        const localOrigin = this.inverseMatrix.transformAsPoint(ray.origin);
+        const localDirection = this.inverseMatrix.transformAsVector(ray.direction);
+        const hit = this.intersectLocalRay(localOrigin, localDirection);
+        if (!hit || hit.distance < 0)
+            return undefined;
+        const point = new Vector3(ray.direction).scale(hit.distance).add(ray.origin);
+        const inverse = this.inverseMatrix;
+        const normal = new Vector3(inverse[0] * hit.normal[0] + inverse[1] * hit.normal[1] + inverse[2] * hit.normal[2], inverse[4] * hit.normal[0] + inverse[5] * hit.normal[1] + inverse[6] * hit.normal[2], inverse[8] * hit.normal[0] + inverse[9] * hit.normal[1] + inverse[10] * hit.normal[2]).normalize();
+        return { distance: hit.distance, point, normal, entering: ray.direction.dot(normal) < 0 };
+    }
+    getAxisAlignedBoundingBox() {
+        const minimum = new Vector3();
+        const maximum = new Vector3();
+        for (let axis = 0; axis < 3; axis++) {
+            const direction = [0, 0, 0];
+            direction[axis] = 1;
+            maximum[axis] = this.getSupportPoint(direction)[axis];
+            direction[axis] = -1;
+            minimum[axis] = this.getSupportPoint(direction)[axis];
+        }
+        return new AxisAlignedBoundingBox(minimum, maximum);
+    }
+    getBoundingSphere() {
+        const box = this.getAxisAlignedBoundingBox();
+        if (!box)
+            return undefined;
+        return new BoundingSphere().fromCornerPoints(box.minimum, box.maximum);
+    }
+    getSupportPoint(worldDirection) {
+        const matrix = this.matrix;
+        const localDirection = [
+            matrix[0] * worldDirection[0] + matrix[1] * worldDirection[1] + matrix[2] * worldDirection[2],
+            matrix[4] * worldDirection[0] + matrix[5] * worldDirection[1] + matrix[6] * worldDirection[2],
+            matrix[8] * worldDirection[0] + matrix[9] * worldDirection[1] + matrix[10] * worldDirection[2]
+        ];
+        return new Vector3(this.supportLocal(localDirection)).transform(this.matrix);
+    }
+}
+export function validateAffineMatrix(matrix) {
+    if (matrix.length !== 16 || matrix.some(value => !Number.isFinite(value))) {
+        throw new Error('Shape matrix must contain 16 finite values');
+    }
+    if (Math.abs(matrix[3]) > 1e-12 ||
+        Math.abs(matrix[7]) > 1e-12 ||
+        Math.abs(matrix[11]) > 1e-12 ||
+        Math.abs(matrix[15] - 1) > 1e-12) {
+        throw new Error('Shape matrix must be affine');
+    }
+    if (Math.abs(new Matrix4().copy(matrix).determinant()) < 1e-15)
+        throw new Error('Shape matrix must be invertible');
+}
+export function solveQuadratic(a, b, c) {
+    if (Math.abs(a) < 1e-15)
+        return Math.abs(b) < 1e-15 ? [] : [-c / b];
+    const discriminant = b * b - 4 * a * c;
+    if (discriminant < 0)
+        return [];
+    const root = Math.sqrt(Math.max(0, discriminant));
+    return [(-b - root) / (2 * a), (-b + root) / (2 * a)].sort((left, right) => left - right);
+}
+export function closestDistanceToSegment2D(px, py, ax, ay, bx, by) {
+    const dx = bx - ax;
+    const dy = by - ay;
+    const denominator = dx * dx + dy * dy;
+    const t = denominator
+        ? Math.max(0, Math.min(1, ((px - ax) * dx + (py - ay) * dy) / denominator))
+        : 0;
+    return Math.hypot(px - (ax + t * dx), py - (ay + t * dy));
+}
+/** Distance in the radial/Y half-plane to a spherical profile arc selected by normal Y range. */
+export function closestDistanceToCircleArc2D(px, py, centerY, radius, minimumNormalY, maximumNormalY) {
+    if (radius === 0)
+        return Math.hypot(px, py - centerY);
+    const dx = px;
+    const dy = py - centerY;
+    const length = Math.hypot(dx, dy);
+    const normalY = length ? dy / length : 0;
+    if (normalY >= minimumNormalY && normalY <= maximumNormalY) {
+        return Math.abs(length - radius);
+    }
+    const minimumRadius = radius * Math.sqrt(Math.max(0, 1 - minimumNormalY ** 2));
+    const maximumRadius = radius * Math.sqrt(Math.max(0, 1 - maximumNormalY ** 2));
+    return Math.min(Math.hypot(px - minimumRadius, py - (centerY + radius * minimumNormalY)), Math.hypot(px - maximumRadius, py - (centerY + radius * maximumNormalY)));
+}
+//# sourceMappingURL=implicit-shape.js.map

package/dist/culling/lib/shapes/implicit-shape.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"implicit-shape.js","sourceRoot":"","sources":["../../../../src/culling/lib/shapes/implicit-shape.ts"],"names":[],"mappings":"AAAA,UAAU;AACV,+BAA+B;AAC/B,4DAA4D;AAE5D,OAAO,EAAC,OAAO,EAAE,OAAO,EAAC,MAAM,eAAe,CAAC;AAE/C,OAAO,EAAC,sBAAsB,EAAC,yDAAsD;AACrF,OAAO,EAAC,cAAc,EAAC,+CAA4C;AAyBnE,kGAAkG;AAClG,MAAM,OAAgB,iBAAiB;IAKrC,YAAY,MAA0B;QACpC,IAAI,CAAC,MAAM,GAAG,IAAI,OAAO,EAAE,CAAC;QAC5B,IAAI,MAAM;YAAE,IAAI,CAAC,MAAM,CAAC,IAAI,CAAC,MAAM,CAAC,CAAC;QACrC,oBAAoB,CAAC,IAAI,CAAC,MAAM,CAAC,CAAC;QAClC,IAAI,CAAC,aAAa,GAAG,IAAI,CAAC,MAAM,CAAC,KAAK,EAAE,CAAC,MAAM,EAAE,CAAC;IACpD,CAAC;IAWD,MAAM,CAAC,KAAoB;QACzB,IAAI,IAAI,CAAC,IAAI,KAAK,KAAK,CAAC,IAAI;YAAE,OAAO,KAAK,CAAC;QAC3C,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,EAAE,EAAE,CAAC,EAAE;YAAE,IAAI,IAAI,CAAC,MAAM,CAAC,CAAC,CAAC,KAAK,KAAK,CAAC,MAAM,CAAC,CAAC,CAAC;gBAAE,OAAO,KAAK,CAAC;QAClF,OAAO,IAAI,CAAC;IACd,CAAC;IAED,SAAS,CAAC,SAA4B;QACpC,MAAM,IAAI,GAAG,IAAI,OAAO,EAAE,CAAC,IAAI,CAAC,SAAS,CAAC,CAAC;QAC3C,oBAAoB,CAAC,IAAI,CAAC,CAAC;QAC3B,IAAI,CAAC,MAAM,CAAC,YAAY,CAAC,IAAI,CAAC,CAAC;QAC/B,oBAAoB,CAAC,IAAI,CAAC,MAAM,CAAC,CAAC;QAClC,IAAI,CAAC,aAAa,CAAC,IAAI,CAAC,IAAI,CAAC,MAAM,CAAC,CAAC,MAAM,EAAE,CAAC;QAC9C,OAAO,IAAI,CAAC;IACd,CAAC;IAED,aAAa,CAAC,KAAwB;QACpC,OAAO,IAAI,CAAC,kBAAkB,CAAC,IAAI,CAAC,aAAa,CAAC,gBAAgB,CAAC,KAAK,CAAsB,CAAC,CAAC;IAClG,CAAC;IAED,iBAAiB,CAAC,KAAwB;QACxC,MAAM,QAAQ,GAAG,IAAI,CAAC,UAAU,CAAC,KAAK,CAAC,CAAC;QACxC,OAAO,QAAQ,GAAG,QAAQ,CAAC;IAC7B,CAAC;IAED,gGAAgG;IAChG,UAAU,CAAC,KAAwB;QACjC,MAAM,UAAU,GAAG,IAAI,CAAC,aAAa,CAAC,gBAAgB,CAAC,KAAK,CAAC,CAAC;QAC9D,MAAM,aAAa,GAAG,IAAI,CAAC,oBAAoB,CAAC,UAA+B,CAAC,CAAC;QACjF,MAAM,MAAM,GAAG,IAAI,CAAC,MAAM,CAAC,QAAQ,EAAE,CAAC;QACtC,OAAO,aAAa,GAAG,IAAI,CAAC,GAAG,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC;IACnE,CAAC;IAED,cAAc,CAAC,KAAY;QACzB,MAAM,OAAO,GAAG,IAAI,CAAC,eAAe,CAAC,KAAK,CAAC,MAAM,CAAC,CAAC;QACnD,MAAM,OAAO,GAAG,IAAI,CAAC,eAAe,CAAC,CAAC,CAAC,KAAK,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,CAAC,KAAK,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,CAAC,KAAK,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QAC7F,MAAM,WAAW,GAAG,KAAK,CAAC,gBAAgB,CAAC,OAAO,CAAC,CAAC;QACpD,MAAM,WAAW,GAAG,KAAK,CAAC,gBAAgB,CAAC,OAAO,CAAC,CAAC;QACpD,IAAI,WAAW,GAAG,CAAC;YAAE,OAAO,QAAQ,CAAC;QACrC,IAAI,WAAW,GAAG,CAAC;YAAE,OAAO,SAAS,CAAC;QACtC,OAAO,cAAc,CAAC;IACxB,CAAC;IAED,YAAY,CAAC,GAAQ;QACnB,MAAM,WAAW,GAAG,IAAI,CAAC,aAAa,CAAC,gBAAgB,CAAC,GAAG,CAAC,MAAM,CAAC,CAAC;QACpE,MAAM,cAAc,GAAG,IAAI,CAAC,aAAa,CAAC,iBAAiB,CAAC,GAAG,CAAC,SAAS,CAAC,CAAC;QAC3E,MAAM,GAAG,GAAG,IAAI,CAAC,iBAAiB,CAChC,WAAgC,EAChC,cAAmC,CACpC,CAAC;QACF,IAAI,CAAC,GAAG,IAAI,GAAG,CAAC,QAAQ,GAAG,CAAC;YAAE,OAAO,SAAS,CAAC;QAC/C,MAAM,KAAK,GAAG,IAAI,OAAO,CAAC,GAAG,CAAC,SAAS,CAAC,CAAC,KAAK,CAAC,GAAG,CAAC,QAAQ,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,MAAM,CAAC,CAAC;QAC7E,MAAM,OAAO,GAAG,IAAI,CAAC,aAAa,CAAC;QACnC,MAAM,MAAM,GAAG,IAAI,OAAO,CACxB,OAAO,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,MAAM,CAAC,CAAC,CAAC,GAAG,OAAO,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,MAAM,CAAC,CAAC,CAAC,GAAG,OAAO,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,MAAM,CAAC,CAAC,CAAC,EACpF,OAAO,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,MAAM,CAAC,CAAC,CAAC,GAAG,OAAO,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,MAAM,CAAC,CAAC,CAAC,GAAG,OAAO,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,MAAM,CAAC,CAAC,CAAC,EACpF,OAAO,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,MAAM,CAAC,CAAC,CAAC,GAAG,OAAO,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,MAAM,CAAC,CAAC,CAAC,GAAG,OAAO,CAAC,EAAE,CAAC,GAAG,GAAG,CAAC,MAAM,CAAC,CAAC,CAAC,CACtF,CAAC,SAAS,EAAE,CAAC;QACd,OAAO,EAAC,QAAQ,EAAE,GAAG,CAAC,QAAQ,EAAE,KAAK,EAAE,MAAM,EAAE,QAAQ,EAAE,GAAG,CAAC,SAAS,CAAC,GAAG,CAAC,MAAM,CAAC,GAAG,CAAC,EAAC,CAAC;IAC1F,CAAC;IAED,yBAAyB;QACvB,MAAM,OAAO,GAAG,IAAI,OAAO,EAAE,CAAC;QAC9B,MAAM,OAAO,GAAG,IAAI,OAAO,EAAE,CAAC;QAC9B,KAAK,IAAI,IAAI,GAAG,CAAC,EAAE,IAAI,GAAG,CAAC,EAAE,IAAI,EAAE,EAAE,CAAC;YACpC,MAAM,SAAS,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;YAC5B,SAAS,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC;YACpB,OAAO,CAAC,IAAI,CAAC,GAAG,IAAI,CAAC,eAAe,CAAC,SAAS,CAAC,CAAC,IAAI,CAAC,CAAC;YACtD,SAAS,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC;YACrB,OAAO,CAAC,IAAI,CAAC,GAAG,IAAI,CAAC,eAAe,CAAC,SAAS,CAAC,CAAC,IAAI,CAAC,CAAC;QACxD,CAAC;QACD,OAAO,IAAI,sBAAsB,CAAC,OAAO,EAAE,OAAO,CAAC,CAAC;IACtD,CAAC;IAED,iBAAiB;QACf,MAAM,GAAG,GAAG,IAAI,CAAC,yBAAyB,EAAE,CAAC;QAC7C,IAAI,CAAC,GAAG;YAAE,OAAO,SAAS,CAAC;QAC3B,OAAO,IAAI,cAAc,EAAE,CAAC,gBAAgB,CAAC,GAAG,CAAC,OAAO,EAAE,GAAG,CAAC,OAAO,CAAC,CAAC;IACzE,CAAC;IAES,eAAe,CAAC,cAAiC;QACzD,MAAM,MAAM,GAAG,IAAI,CAAC,MAAM,CAAC;QAC3B,MAAM,cAAc,GAAG;YACrB,MAAM,CAAC,CAAC,CAAC,GAAG,cAAc,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,GAAG,cAAc,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,GAAG,cAAc,CAAC,CAAC,CAAC;YAC7F,MAAM,CAAC,CAAC,CAAC,GAAG,cAAc,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,GAAG,cAAc,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,GAAG,cAAc,CAAC,CAAC,CAAC;YAC7F,MAAM,CAAC,CAAC,CAAC,GAAG,cAAc,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,GAAG,cAAc,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC,EAAE,CAAC,GAAG,cAAc,CAAC,CAAC,CAAC;SAC/F,CAAC;QACF,OAAO,IAAI,OAAO,CAAC,IAAI,CAAC,YAAY,CAAC,cAAc,CAAC,CAAC,CAAC,SAAS,CAAC,IAAI,CAAC,MAAM,CAAC,CAAC;IAC/E,CAAC;CACF;AAED,MAAM,UAAU,oBAAoB,CAAC,MAAyB;IAC5D,IAAI,MAAM,CAAC,MAAM,KAAK,EAAE,IAAI,MAAM,CAAC,IAAI,CAAC,KAAK,CAAC,EAAE,CAAC,CAAC,MAAM,CAAC,QAAQ,CAAC,KAAK,CAAC,CAAC,EAAE,CAAC;QAC1E,MAAM,IAAI,KAAK,CAAC,4CAA4C,CAAC,CAAC;IAChE,CAAC;IACD,IACE,IAAI,CAAC,GAAG,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK;QAC3B,IAAI,CAAC,GAAG,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK;QAC3B,IAAI,CAAC,GAAG,CAAC,MAAM,CAAC,EAAE,CAAC,CAAC,GAAG,KAAK;QAC5B,IAAI,CAAC,GAAG,CAAC,MAAM,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,GAAG,KAAK,EAChC,CAAC;QACD,MAAM,IAAI,KAAK,CAAC,6BAA6B,CAAC,CAAC;IACjD,CAAC;IACD,IAAI,IAAI,CAAC,GAAG,CAAC,IAAI,OAAO,EAAE,CAAC,IAAI,CAAC,MAAM,CAAC,CAAC,WAAW,EAAE,CAAC,GAAG,KAAK;QAC5D,MAAM,IAAI,KAAK,CAAC,iCAAiC,CAAC,CAAC;AACvD,CAAC;AAED,MAAM,UAAU,cAAc,CAAC,CAAS,EAAE,CAAS,EAAE,CAAS;IAC5D,IAAI,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,KAAK;QAAE,OAAO,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;IACpE,MAAM,YAAY,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;IACvC,IAAI,YAAY,GAAG,CAAC;QAAE,OAAO,EAAE,CAAC;IAChC,MAAM,IAAI,GAAG,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,YAAY,CAAC,CAAC,CAAC;IAClD,OAAO,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,IAAI,EAAE,KAAK,EAAE,EAAE,CAAC,IAAI,GAAG,KAAK,CAAC,CAAC;AAC5F,CAAC;AAED,MAAM,UAAU,0BAA0B,CACxC,EAAU,EACV,EAAU,EACV,EAAU,EACV,EAAU,EACV,EAAU,EACV,EAAU;IAEV,MAAM,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC;IACnB,MAAM,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC;IACnB,MAAM,WAAW,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC;IACtC,MAAM,CAAC,GAAG,WAAW;QACnB,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,GAAG,EAAE,CAAC,GAAG,EAAE,GAAG,CAAC,EAAE,GAAG,EAAE,CAAC,GAAG,EAAE,CAAC,GAAG,WAAW,CAAC,CAAC;QAC3E,CAAC,CAAC,CAAC,CAAC;IACN,OAAO,IAAI,CAAC,KAAK,CAAC,EAAE,GAAG,CAAC,EAAE,GAAG,CAAC,GAAG,EAAE,CAAC,EAAE,EAAE,GAAG,CAAC,EAAE,GAAG,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC;AAC5D,CAAC;AAED,iGAAiG;AACjG,MAAM,UAAU,4BAA4B,CAC1C,EAAU,EACV,EAAU,EACV,OAAe,EACf,MAAc,EACd,cAAsB,EACtB,cAAsB;IAEtB,IAAI,MAAM,KAAK,CAAC;QAAE,OAAO,IAAI,CAAC,KAAK,CAAC,EAAE,EAAE,EAAE,GAAG,OAAO,CAAC,CAAC;IACtD,MAAM,EAAE,GAAG,EAAE,CAAC;IACd,MAAM,EAAE,GAAG,EAAE,GAAG,OAAO,CAAC;IACxB,MAAM,MAAM,GAAG,IAAI,CAAC,KAAK,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC;IAClC,MAAM,OAAO,GAAG,MAAM,CAAC,CAAC,CAAC,EAAE,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC;IACzC,IAAI,OAAO,IAAI,cAAc,IAAI,OAAO,IAAI,cAAc,EAAE,CAAC;QAC3D,OAAO,IAAI,CAAC,GAAG,CAAC,MAAM,GAAG,MAAM,CAAC,CAAC;IACnC,CAAC;IACD,MAAM,aAAa,GAAG,MAAM,GAAG,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,GAAG,cAAc,IAAI,CAAC,CAAC,CAAC,CAAC;IAC/E,MAAM,aAAa,GAAG,MAAM,GAAG,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,GAAG,cAAc,IAAI,CAAC,CAAC,CAAC,CAAC;IAC/E,OAAO,IAAI,CAAC,GAAG,CACb,IAAI,CAAC,KAAK,CAAC,EAAE,GAAG,aAAa,EAAE,EAAE,GAAG,CAAC,OAAO,GAAG,MAAM,GAAG,cAAc,CAAC,CAAC,EACxE,IAAI,CAAC,KAAK,CAAC,EAAE,GAAG,aAAa,EAAE,EAAE,GAAG,CAAC,OAAO,GAAG,MAAM,GAAG,cAAc,CAAC,CAAC,CACzE,CAAC;AACJ,CAAC"}

package/dist/culling/lib/shapes/plane-shape.d.ts

@@ -0,0 +1,27 @@
+import type { CullingResult } from "../../constants.js";
+import { ImplicitShapeBase, type ImplicitShape, type LocalRayIntersection } from "./implicit-shape.js";
+import type { AxisAlignedBoundingBox } from "../bounding-volumes/axis-aligned-bounding-box.js";
+import type { BoundingSphere } from "../bounding-volumes/bounding-sphere.js";
+import type { Plane } from "../plane.js";
+export type PlaneShapeProps = {
+    sizeX?: number;
+    sizeZ?: number;
+    matrix?: readonly number[];
+};
+/** Analytic glTF plane. The surface faces +Y and the negative-Y half-space is inside. */
+export declare class PlaneShape extends ImplicitShapeBase {
+    readonly type: "plane";
+    readonly sizeX?: number;
+    readonly sizeZ?: number;
+    constructor(props?: PlaneShapeProps);
+    clone(): PlaneShape;
+    equals(shape: ImplicitShape): boolean;
+    containsLocalPoint(point: readonly number[]): boolean;
+    distanceToLocalPoint(point: readonly number[]): number;
+    supportLocal(): readonly number[];
+    intersectLocalRay(origin: readonly number[], direction: readonly number[]): LocalRayIntersection | undefined;
+    intersectPlane(plane: Plane): CullingResult;
+    getAxisAlignedBoundingBox(): AxisAlignedBoundingBox | undefined;
+    getBoundingSphere(): BoundingSphere | undefined;
+}
+//# sourceMappingURL=plane-shape.d.ts.map

package/dist/culling/lib/shapes/plane-shape.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"plane-shape.d.ts","sourceRoot":"","sources":["../../../../src/culling/lib/shapes/plane-shape.ts"],"names":[],"mappings":"AAKA,OAAO,KAAK,EAAC,aAAa,EAAC,2BAAwB;AACnD,OAAO,EAAC,iBAAiB,EAAE,KAAK,aAAa,EAAE,KAAK,oBAAoB,EAAC,4BAAyB;AAClG,OAAO,KAAK,EAAC,sBAAsB,EAAC,yDAAsD;AAC1F,OAAO,KAAK,EAAC,cAAc,EAAC,+CAA4C;AACxE,OAAO,KAAK,EAAC,KAAK,EAAC,oBAAiB;AAEpC,MAAM,MAAM,eAAe,GAAG;IAAC,KAAK,CAAC,EAAE,MAAM,CAAC;IAAC,KAAK,CAAC,EAAE,MAAM,CAAC;IAAC,MAAM,CAAC,EAAE,SAAS,MAAM,EAAE,CAAA;CAAC,CAAC;AAE3F,yFAAyF;AACzF,qBAAa,UAAW,SAAQ,iBAAiB;IAC/C,QAAQ,CAAC,IAAI,EAAG,OAAO,CAAU;IACjC,QAAQ,CAAC,KAAK,CAAC,EAAE,MAAM,CAAC;IACxB,QAAQ,CAAC,KAAK,CAAC,EAAE,MAAM,CAAC;gBACZ,KAAK,GAAE,eAAoB;IASvC,KAAK,IAAI,UAAU;IAGV,MAAM,CAAC,KAAK,EAAE,aAAa,GAAG,OAAO;IAQ9C,kBAAkB,CAAC,KAAK,EAAE,SAAS,MAAM,EAAE,GAAG,OAAO;IAGrD,oBAAoB,CAAC,KAAK,EAAE,SAAS,MAAM,EAAE,GAAG,MAAM;IAGtD,YAAY,IAAI,SAAS,MAAM,EAAE;IAGjC,iBAAiB,CACf,MAAM,EAAE,SAAS,MAAM,EAAE,EACzB,SAAS,EAAE,SAAS,MAAM,EAAE,GAC3B,oBAAoB,GAAG,SAAS;IAU1B,cAAc,CAAC,KAAK,EAAE,KAAK,GAAG,aAAa;IAW3C,yBAAyB,IAAI,sBAAsB,GAAG,SAAS;IAG/D,iBAAiB,IAAI,cAAc,GAAG,SAAS;CAGzD"}

package/dist/culling/lib/shapes/plane-shape.js

@@ -0,0 +1,71 @@
+// math.gl
+// SPDX-License-Identifier: MIT
+// SPDX-FileCopyrightText: Copyright (c) vis.gl contributors
+import { Vector3 } from '@math.gl/core';
+import { ImplicitShapeBase } from "./implicit-shape.js";
+/** Analytic glTF plane. The surface faces +Y and the negative-Y half-space is inside. */
+export class PlaneShape extends ImplicitShapeBase {
+    constructor(props = {}) {
+        super(props.matrix);
+        this.type = 'plane';
+        this.sizeX = props.sizeX;
+        this.sizeZ = props.sizeZ;
+        if (this.sizeX !== undefined && (!Number.isFinite(this.sizeX) || this.sizeX <= 0))
+            throw new Error('Plane sizeX must be greater than zero');
+        if (this.sizeZ !== undefined && (!Number.isFinite(this.sizeZ) || this.sizeZ <= 0))
+            throw new Error('Plane sizeZ must be greater than zero');
+    }
+    clone() {
+        return new PlaneShape({ sizeX: this.sizeX, sizeZ: this.sizeZ, matrix: this.matrix });
+    }
+    equals(shape) {
+        return (shape instanceof PlaneShape &&
+            super.equals(shape) &&
+            this.sizeX === shape.sizeX &&
+            this.sizeZ === shape.sizeZ);
+    }
+    containsLocalPoint(point) {
+        return point[1] <= 1e-12;
+    }
+    distanceToLocalPoint(point) {
+        return Math.max(0, point[1]);
+    }
+    supportLocal() {
+        throw new Error('An unbounded plane half-space has no finite support point');
+    }
+    intersectLocalRay(origin, direction) {
+        if (Math.abs(direction[1]) < 1e-15)
+            return undefined;
+        const distance = -origin[1] / direction[1];
+        if (distance < 0)
+            return undefined;
+        const x = origin[0] + distance * direction[0];
+        const z = origin[2] + distance * direction[2];
+        if (this.sizeX !== undefined && Math.abs(x) > this.sizeX / 2 + 1e-12)
+            return undefined;
+        if (this.sizeZ !== undefined && Math.abs(z) > this.sizeZ / 2 + 1e-12)
+            return undefined;
+        return { distance, normal: [0, 1, 0] };
+    }
+    intersectPlane(plane) {
+        const boundaryPoint = new Vector3(0, 0, 0).transform(this.matrix);
+        const inverse = this.inverseMatrix;
+        const boundaryNormal = new Vector3(inverse[1], inverse[5], inverse[9]).normalize();
+        const alignment = boundaryNormal.dot(plane.normal);
+        if (Math.abs(Math.abs(alignment) - 1) > 1e-10)
+            return 'intersecting';
+        const boundaryDistance = plane.getPointDistance(boundaryPoint);
+        if (alignment < 0 && boundaryDistance > 0)
+            return 'inside';
+        if (alignment > 0 && boundaryDistance < 0)
+            return 'outside';
+        return 'intersecting';
+    }
+    getAxisAlignedBoundingBox() {
+        return undefined;
+    }
+    getBoundingSphere() {
+        return undefined;
+    }
+}
+//# sourceMappingURL=plane-shape.js.map

package/dist/culling/lib/shapes/plane-shape.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"plane-shape.js","sourceRoot":"","sources":["../../../../src/culling/lib/shapes/plane-shape.ts"],"names":[],"mappings":"AAAA,UAAU;AACV,+BAA+B;AAC/B,4DAA4D;AAE5D,OAAO,EAAC,OAAO,EAAC,MAAM,eAAe,CAAC;AAEtC,OAAO,EAAC,iBAAiB,EAAgD,4BAAyB;AAOlG,yFAAyF;AACzF,MAAM,OAAO,UAAW,SAAQ,iBAAiB;IAI/C,YAAY,QAAyB,EAAE;QACrC,KAAK,CAAC,KAAK,CAAC,MAAM,CAAC,CAAC;QAJb,SAAI,GAAG,OAAgB,CAAC;QAK/B,IAAI,CAAC,KAAK,GAAG,KAAK,CAAC,KAAK,CAAC;QACzB,IAAI,CAAC,KAAK,GAAG,KAAK,CAAC,KAAK,CAAC;QACzB,IAAI,IAAI,CAAC,KAAK,KAAK,SAAS,IAAI,CAAC,CAAC,MAAM,CAAC,QAAQ,CAAC,IAAI,CAAC,KAAK,CAAC,IAAI,IAAI,CAAC,KAAK,IAAI,CAAC,CAAC;YAC/E,MAAM,IAAI,KAAK,CAAC,uCAAuC,CAAC,CAAC;QAC3D,IAAI,IAAI,CAAC,KAAK,KAAK,SAAS,IAAI,CAAC,CAAC,MAAM,CAAC,QAAQ,CAAC,IAAI,CAAC,KAAK,CAAC,IAAI,IAAI,CAAC,KAAK,IAAI,CAAC,CAAC;YAC/E,MAAM,IAAI,KAAK,CAAC,uCAAuC,CAAC,CAAC;IAC7D,CAAC;IACD,KAAK;QACH,OAAO,IAAI,UAAU,CAAC,EAAC,KAAK,EAAE,IAAI,CAAC,KAAK,EAAE,KAAK,EAAE,IAAI,CAAC,KAAK,EAAE,MAAM,EAAE,IAAI,CAAC,MAAM,EAAC,CAAC,CAAC;IACrF,CAAC;IACQ,MAAM,CAAC,KAAoB;QAClC,OAAO,CACL,KAAK,YAAY,UAAU;YAC3B,KAAK,CAAC,MAAM,CAAC,KAAK,CAAC;YACnB,IAAI,CAAC,KAAK,KAAK,KAAK,CAAC,KAAK;YAC1B,IAAI,CAAC,KAAK,KAAK,KAAK,CAAC,KAAK,CAC3B,CAAC;IACJ,CAAC;IACD,kBAAkB,CAAC,KAAwB;QACzC,OAAO,KAAK,CAAC,CAAC,CAAC,IAAI,KAAK,CAAC;IAC3B,CAAC;IACD,oBAAoB,CAAC,KAAwB;QAC3C,OAAO,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC;IAC/B,CAAC;IACD,YAAY;QACV,MAAM,IAAI,KAAK,CAAC,2DAA2D,CAAC,CAAC;IAC/E,CAAC;IACD,iBAAiB,CACf,MAAyB,EACzB,SAA4B;QAE5B,IAAI,IAAI,CAAC,GAAG,CAAC,SAAS,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK;YAAE,OAAO,SAAS,CAAC;QACrD,MAAM,QAAQ,GAAG,CAAC,MAAM,CAAC,CAAC,CAAC,GAAG,SAAS,CAAC,CAAC,CAAC,CAAC;QAC3C,IAAI,QAAQ,GAAG,CAAC;YAAE,OAAO,SAAS,CAAC;QACnC,MAAM,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,GAAG,QAAQ,GAAG,SAAS,CAAC,CAAC,CAAC,CAAC;QAC9C,MAAM,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,GAAG,QAAQ,GAAG,SAAS,CAAC,CAAC,CAAC,CAAC;QAC9C,IAAI,IAAI,CAAC,KAAK,KAAK,SAAS,IAAI,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,KAAK,GAAG,CAAC,GAAG,KAAK;YAAE,OAAO,SAAS,CAAC;QACvF,IAAI,IAAI,CAAC,KAAK,KAAK,SAAS,IAAI,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,KAAK,GAAG,CAAC,GAAG,KAAK;YAAE,OAAO,SAAS,CAAC;QACvF,OAAO,EAAC,QAAQ,EAAE,MAAM,EAAE,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,EAAC,CAAC;IACvC,CAAC;IACQ,cAAc,CAAC,KAAY;QAClC,MAAM,aAAa,GAAG,IAAI,OAAO,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,SAAS,CAAC,IAAI,CAAC,MAAM,CAAC,CAAC;QAClE,MAAM,OAAO,GAAG,IAAI,CAAC,aAAa,CAAC;QACnC,MAAM,cAAc,GAAG,IAAI,OAAO,CAAC,OAAO,CAAC,CAAC,CAAC,EAAE,OAAO,CAAC,CAAC,CAAC,EAAE,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC,SAAS,EAAE,CAAC;QACnF,MAAM,SAAS,GAAG,cAAc,CAAC,GAAG,CAAC,KAAK,CAAC,MAAM,CAAC,CAAC;QACnD,IAAI,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,SAAS,CAAC,GAAG,CAAC,CAAC,GAAG,KAAK;YAAE,OAAO,cAAc,CAAC;QACrE,MAAM,gBAAgB,GAAG,KAAK,CAAC,gBAAgB,CAAC,aAAa,CAAC,CAAC;QAC/D,IAAI,SAAS,GAAG,CAAC,IAAI,gBAAgB,GAAG,CAAC;YAAE,OAAO,QAAQ,CAAC;QAC3D,IAAI,SAAS,GAAG,CAAC,IAAI,gBAAgB,GAAG,CAAC;YAAE,OAAO,SAAS,CAAC;QAC5D,OAAO,cAAc,CAAC;IACxB,CAAC;IACQ,yBAAyB;QAChC,OAAO,SAAS,CAAC;IACnB,CAAC;IACQ,iBAAiB;QACxB,OAAO,SAAS,CAAC;IACnB,CAAC;CACF"}

package/dist/culling/lib/shapes/sphere-shape.d.ts

@@ -0,0 +1,18 @@
+import { ImplicitShapeBase, type ImplicitShape, type LocalRayIntersection } from "./implicit-shape.js";
+export type SphereShapeProps = {
+    radius?: number;
+    matrix?: readonly number[];
+};
+/** Analytic glTF sphere shape. */
+export declare class SphereShape extends ImplicitShapeBase {
+    readonly type: "sphere";
+    readonly radius: number;
+    constructor(props?: SphereShapeProps);
+    clone(): SphereShape;
+    equals(shape: ImplicitShape): boolean;
+    containsLocalPoint(point: readonly number[]): boolean;
+    distanceToLocalPoint(point: readonly number[]): number;
+    supportLocal(direction: readonly number[]): readonly number[];
+    intersectLocalRay(origin: readonly number[], direction: readonly number[]): LocalRayIntersection | undefined;
+}
+//# sourceMappingURL=sphere-shape.d.ts.map

package/dist/culling/lib/shapes/sphere-shape.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"sphere-shape.d.ts","sourceRoot":"","sources":["../../../../src/culling/lib/shapes/sphere-shape.ts"],"names":[],"mappings":"AAIA,OAAO,EACL,iBAAiB,EAEjB,KAAK,aAAa,EAClB,KAAK,oBAAoB,EAC1B,4BAAyB;AAI1B,MAAM,MAAM,gBAAgB,GAAG;IAAC,MAAM,CAAC,EAAE,MAAM,CAAC;IAAC,MAAM,CAAC,EAAE,SAAS,MAAM,EAAE,CAAA;CAAC,CAAC;AAE7E,kCAAkC;AAClC,qBAAa,WAAY,SAAQ,iBAAiB;IAChD,QAAQ,CAAC,IAAI,EAAG,QAAQ,CAAU;IAClC,QAAQ,CAAC,MAAM,EAAE,MAAM,CAAC;gBACZ,KAAK,GAAE,gBAAqB;IAMxC,KAAK,IAAI,WAAW;IAGX,MAAM,CAAC,KAAK,EAAE,aAAa,GAAG,OAAO;IAG9C,kBAAkB,CAAC,KAAK,EAAE,SAAS,MAAM,EAAE,GAAG,OAAO;IAGrD,oBAAoB,CAAC,KAAK,EAAE,SAAS,MAAM,EAAE,GAAG,MAAM;IAGtD,YAAY,CAAC,SAAS,EAAE,SAAS,MAAM,EAAE,GAAG,SAAS,MAAM,EAAE;IAI7D,iBAAiB,CACf,MAAM,EAAE,SAAS,MAAM,EAAE,EACzB,SAAS,EAAE,SAAS,MAAM,EAAE,GAC3B,oBAAoB,GAAG,SAAS;CAWpC"}

package/dist/culling/lib/shapes/sphere-shape.js

@@ -0,0 +1,43 @@
+// math.gl
+// SPDX-License-Identifier: MIT
+// SPDX-FileCopyrightText: Copyright (c) vis.gl contributors
+import { ImplicitShapeBase, solveQuadratic } from "./implicit-shape.js";
+const EPSILON = 1e-6;
+/** Analytic glTF sphere shape. */
+export class SphereShape extends ImplicitShapeBase {
+    constructor(props = {}) {
+        super(props.matrix);
+        this.type = 'sphere';
+        this.radius = props.radius ?? 0.5;
+        if (!Number.isFinite(this.radius) || this.radius <= 0)
+            throw new Error('Sphere radius must be greater than zero');
+    }
+    clone() {
+        return new SphereShape({ radius: this.radius, matrix: this.matrix });
+    }
+    equals(shape) {
+        return shape instanceof SphereShape && super.equals(shape) && this.radius === shape.radius;
+    }
+    containsLocalPoint(point) {
+        return Math.hypot(...point) <= this.radius + EPSILON;
+    }
+    distanceToLocalPoint(point) {
+        return Math.max(0, Math.hypot(...point) - this.radius);
+    }
+    supportLocal(direction) {
+        const length = Math.hypot(...direction);
+        return length ? direction.map(value => (value / length) * this.radius) : [this.radius, 0, 0];
+    }
+    intersectLocalRay(origin, direction) {
+        const roots = solveQuadratic(dot(direction, direction), 2 * dot(origin, direction), dot(origin, origin) - this.radius * this.radius);
+        const distance = roots.find(value => value >= 0);
+        if (distance === undefined)
+            return undefined;
+        const point = origin.map((value, i) => value + distance * direction[i]);
+        return { distance, normal: point.map(value => value / this.radius) };
+    }
+}
+function dot(a, b) {
+    return a[0] * b[0] + a[1] * b[1] + a[2] * b[2];
+}
+//# sourceMappingURL=sphere-shape.js.map

package/dist/culling/lib/shapes/sphere-shape.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"sphere-shape.js","sourceRoot":"","sources":["../../../../src/culling/lib/shapes/sphere-shape.ts"],"names":[],"mappings":"AAAA,UAAU;AACV,+BAA+B;AAC/B,4DAA4D;AAE5D,OAAO,EACL,iBAAiB,EACjB,cAAc,EAGf,4BAAyB;AAE1B,MAAM,OAAO,GAAG,IAAI,CAAC;AAIrB,kCAAkC;AAClC,MAAM,OAAO,WAAY,SAAQ,iBAAiB;IAGhD,YAAY,QAA0B,EAAE;QACtC,KAAK,CAAC,KAAK,CAAC,MAAM,CAAC,CAAC;QAHb,SAAI,GAAG,QAAiB,CAAC;QAIhC,IAAI,CAAC,MAAM,GAAG,KAAK,CAAC,MAAM,IAAI,GAAG,CAAC;QAClC,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,IAAI,CAAC,MAAM,CAAC,IAAI,IAAI,CAAC,MAAM,IAAI,CAAC;YACnD,MAAM,IAAI,KAAK,CAAC,yCAAyC,CAAC,CAAC;IAC/D,CAAC;IACD,KAAK;QACH,OAAO,IAAI,WAAW,CAAC,EAAC,MAAM,EAAE,IAAI,CAAC,MAAM,EAAE,MAAM,EAAE,IAAI,CAAC,MAAM,EAAC,CAAC,CAAC;IACrE,CAAC;IACQ,MAAM,CAAC,KAAoB;QAClC,OAAO,KAAK,YAAY,WAAW,IAAI,KAAK,CAAC,MAAM,CAAC,KAAK,CAAC,IAAI,IAAI,CAAC,MAAM,KAAK,KAAK,CAAC,MAAM,CAAC;IAC7F,CAAC;IACD,kBAAkB,CAAC,KAAwB;QACzC,OAAO,IAAI,CAAC,KAAK,CAAC,GAAG,KAAK,CAAC,IAAI,IAAI,CAAC,MAAM,GAAG,OAAO,CAAC;IACvD,CAAC;IACD,oBAAoB,CAAC,KAAwB;QAC3C,OAAO,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,IAAI,CAAC,KAAK,CAAC,GAAG,KAAK,CAAC,GAAG,IAAI,CAAC,MAAM,CAAC,CAAC;IACzD,CAAC;IACD,YAAY,CAAC,SAA4B;QACvC,MAAM,MAAM,GAAG,IAAI,CAAC,KAAK,CAAC,GAAG,SAAS,CAAC,CAAC;QACxC,OAAO,MAAM,CAAC,CAAC,CAAC,SAAS,CAAC,GAAG,CAAC,KAAK,CAAC,EAAE,CAAC,CAAC,KAAK,GAAG,MAAM,CAAC,GAAG,IAAI,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,MAAM,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;IAC/F,CAAC;IACD,iBAAiB,CACf,MAAyB,EACzB,SAA4B;QAE5B,MAAM,KAAK,GAAG,cAAc,CAC1B,GAAG,CAAC,SAAS,EAAE,SAAS,CAAC,EACzB,CAAC,GAAG,GAAG,CAAC,MAAM,EAAE,SAAS,CAAC,EAC1B,GAAG,CAAC,MAAM,EAAE,MAAM,CAAC,GAAG,IAAI,CAAC,MAAM,GAAG,IAAI,CAAC,MAAM,CAChD,CAAC;QACF,MAAM,QAAQ,GAAG,KAAK,CAAC,IAAI,CAAC,KAAK,CAAC,EAAE,CAAC,KAAK,IAAI,CAAC,CAAC,CAAC;QACjD,IAAI,QAAQ,KAAK,SAAS;YAAE,OAAO,SAAS,CAAC;QAC7C,MAAM,KAAK,GAAG,MAAM,CAAC,GAAG,CAAC,CAAC,KAAK,EAAE,CAAC,EAAE,EAAE,CAAC,KAAK,GAAG,QAAQ,GAAG,SAAS,CAAC,CAAC,CAAC,CAAC,CAAC;QACxE,OAAO,EAAC,QAAQ,EAAE,MAAM,EAAE,KAAK,CAAC,GAAG,CAAC,KAAK,CAAC,EAAE,CAAC,KAAK,GAAG,IAAI,CAAC,MAAM,CAAC,EAAC,CAAC;IACrE,CAAC;CACF;AAED,SAAS,GAAG,CAAC,CAAoB,EAAE,CAAoB;IACrD,OAAO,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;AACjD,CAAC"}

package/dist/culling/oriented-box-2d.cjs

@@ -0,0 +1,73 @@
+var __defProp = Object.defineProperty;
+var __getOwnPropDesc = Object.getOwnPropertyDescriptor;
+var __getOwnPropNames = Object.getOwnPropertyNames;
+var __hasOwnProp = Object.prototype.hasOwnProperty;
+var __export = (target, all) => {
+  for (var name in all)
+    __defProp(target, name, { get: all[name], enumerable: true });
+};
+var __copyProps = (to, from, except, desc) => {
+  if (from && typeof from === "object" || typeof from === "function") {
+    for (let key of __getOwnPropNames(from))
+      if (!__hasOwnProp.call(to, key) && key !== except)
+        __defProp(to, key, { get: () => from[key], enumerable: !(desc = __getOwnPropDesc(from, key)) || desc.enumerable });
+  }
+  return to;
+};
+var __toCommonJS = (mod) => __copyProps(__defProp({}, "__esModule", { value: true }), mod);
+
+// dist/culling/oriented-box-2d.js
+var oriented_box_2d_exports = {};
+__export(oriented_box_2d_exports, {
+  intersectOrientedBoxes2D: () => intersectOrientedBoxes2D
+});
+module.exports = __toCommonJS(oriented_box_2d_exports);
+function intersectOrientedBoxes2D(first, second, epsilon = 0) {
+  validateBox(first);
+  validateBox(second);
+  if (!Number.isFinite(epsilon) || epsilon < 0) {
+    throw new RangeError("Box separation allowance must be finite and nonnegative");
+  }
+  const firstDirectionScale = Math.max(Math.abs(first.direction[0]), Math.abs(first.direction[1]));
+  const secondDirectionScale = Math.max(Math.abs(second.direction[0]), Math.abs(second.direction[1]));
+  const firstWidthX = first.direction[0] / firstDirectionScale;
+  const firstWidthY = first.direction[1] / firstDirectionScale;
+  const secondWidthX = second.direction[0] / secondDirectionScale;
+  const secondWidthY = second.direction[1] / secondDirectionScale;
+  const firstLength = Math.hypot(firstWidthX, firstWidthY);
+  const secondLength = Math.hypot(secondWidthX, secondWidthY);
+  const firstX = firstWidthX / firstLength;
+  const firstY = firstWidthY / firstLength;
+  const secondX = secondWidthX / secondLength;
+  const secondY = secondWidthY / secondLength;
+  const cosine = Math.min(1, Math.abs(firstX * secondX + firstY * secondY));
+  const sine = Math.min(1, Math.abs(firstX * secondY - firstY * secondX));
+  return overlapsOnAxis(first, second, firstX, firstY, first.halfSize[0], second.halfSize[0] * cosine, second.halfSize[1] * sine, epsilon) && overlapsOnAxis(first, second, -firstY, firstX, first.halfSize[1], second.halfSize[0] * sine, second.halfSize[1] * cosine, epsilon) && overlapsOnAxis(first, second, secondX, secondY, second.halfSize[0], first.halfSize[0] * cosine, first.halfSize[1] * sine, epsilon) && overlapsOnAxis(first, second, -secondY, secondX, second.halfSize[1], first.halfSize[0] * sine, first.halfSize[1] * cosine, epsilon);
+}
+function overlapsOnAxis(first, second, axisX, axisY, ownRadius, otherWidthRadius, otherHeightRadius, epsilon) {
+  const displacementX = second.center[0] - first.center[0];
+  const displacementY = second.center[1] - first.center[1];
+  const projectionX = Number.isFinite(displacementX) ? displacementX * axisX : second.center[0] * axisX - first.center[0] * axisX;
+  const projectionY = Number.isFinite(displacementY) ? displacementY * axisY : second.center[1] * axisY - first.center[1] * axisY;
+  const scale = Math.max(ownRadius, otherWidthRadius, otherHeightRadius, epsilon, Number.isFinite(projectionX) ? Math.abs(projectionX) : Math.max(Math.abs(first.center[0] * axisX), Math.abs(second.center[0] * axisX)), Number.isFinite(projectionY) ? Math.abs(projectionY) : Math.max(Math.abs(first.center[1] * axisY), Math.abs(second.center[1] * axisY))) || 1;
+  const normalizedX = Number.isFinite(projectionX) ? projectionX / scale : second.center[0] * axisX / scale - first.center[0] * axisX / scale;
+  const normalizedY = Number.isFinite(projectionY) ? projectionY / scale : second.center[1] * axisY / scale - first.center[1] * axisY / scale;
+  return Math.abs(normalizedX + normalizedY) <= ownRadius / scale + otherWidthRadius / scale + otherHeightRadius / scale + epsilon / scale;
+}
+function validateBox(box) {
+  validateVector(box.center);
+  validateVector(box.halfSize);
+  validateVector(box.direction);
+  if (box.halfSize[0] < 0 || box.halfSize[1] < 0) {
+    throw new RangeError("Box half-sizes must be nonnegative");
+  }
+  if (box.direction[0] === 0 && box.direction[1] === 0) {
+    throw new RangeError("Box direction must be nonzero");
+  }
+}
+function validateVector(vector) {
+  if (vector.length !== 2 || !Number.isFinite(vector[0]) || !Number.isFinite(vector[1])) {
+    throw new RangeError("Box vectors must have exactly two finite components");
+  }
+}
+//# sourceMappingURL=oriented-box-2d.cjs.map

package/dist/culling/oriented-box-2d.cjs.map

@@ -0,0 +1,6 @@
+{
+  "version": 3,
+  "sources": ["../../src/culling/oriented-box-2d.ts"],
+  "mappings": ";;;;;;;;;;;;;;;;;;;AAAA;;;;;AAyBM,SAAU,yBACd,OACA,QACA,UAAU,GAAC;AAEX,cAAY,KAAK;AACjB,cAAY,MAAM;AAClB,MAAI,CAAC,OAAO,SAAS,OAAO,KAAK,UAAU,GAAG;AAC5C,UAAM,IAAI,WAAW,yDAAyD;EAChF;AAGA,QAAM,sBAAsB,KAAK,IAAI,KAAK,IAAI,MAAM,UAAU,CAAC,CAAC,GAAG,KAAK,IAAI,MAAM,UAAU,CAAC,CAAC,CAAC;AAC/F,QAAM,uBAAuB,KAAK,IAChC,KAAK,IAAI,OAAO,UAAU,CAAC,CAAC,GAC5B,KAAK,IAAI,OAAO,UAAU,CAAC,CAAC,CAAC;AAE/B,QAAM,cAAc,MAAM,UAAU,CAAC,IAAI;AACzC,QAAM,cAAc,MAAM,UAAU,CAAC,IAAI;AACzC,QAAM,eAAe,OAAO,UAAU,CAAC,IAAI;AAC3C,QAAM,eAAe,OAAO,UAAU,CAAC,IAAI;AAC3C,QAAM,cAAc,KAAK,MAAM,aAAa,WAAW;AACvD,QAAM,eAAe,KAAK,MAAM,cAAc,YAAY;AAC1D,QAAM,SAAS,cAAc;AAC7B,QAAM,SAAS,cAAc;AAC7B,QAAM,UAAU,eAAe;AAC/B,QAAM,UAAU,eAAe;AAG/B,QAAM,SAAS,KAAK,IAAI,GAAG,KAAK,IAAI,SAAS,UAAU,SAAS,OAAO,CAAC;AACxE,QAAM,OAAO,KAAK,IAAI,GAAG,KAAK,IAAI,SAAS,UAAU,SAAS,OAAO,CAAC;AACtE,SACE,eACE,OACA,QACA,QACA,QACA,MAAM,SAAS,CAAC,GAChB,OAAO,SAAS,CAAC,IAAI,QACrB,OAAO,SAAS,CAAC,IAAI,MACrB,OAAO,KAET,eACE,OACA,QACA,CAAC,QACD,QACA,MAAM,SAAS,CAAC,GAChB,OAAO,SAAS,CAAC,IAAI,MACrB,OAAO,SAAS,CAAC,IAAI,QACrB,OAAO,KAET,eACE,OACA,QACA,SACA,SACA,OAAO,SAAS,CAAC,GACjB,MAAM,SAAS,CAAC,IAAI,QACpB,MAAM,SAAS,CAAC,IAAI,MACpB,OAAO,KAET,eACE,OACA,QACA,CAAC,SACD,SACA,OAAO,SAAS,CAAC,GACjB,MAAM,SAAS,CAAC,IAAI,MACpB,MAAM,SAAS,CAAC,IAAI,QACpB,OAAO;AAGb;AAGA,SAAS,eACP,OACA,QACA,OACA,OACA,WACA,kBACA,mBACA,SAAe;AAEf,QAAM,gBAAgB,OAAO,OAAO,CAAC,IAAI,MAAM,OAAO,CAAC;AACvD,QAAM,gBAAgB,OAAO,OAAO,CAAC,IAAI,MAAM,OAAO,CAAC;AAGvD,QAAM,cAAc,OAAO,SAAS,aAAa,IAC7C,gBAAgB,QAChB,OAAO,OAAO,CAAC,IAAI,QAAQ,MAAM,OAAO,CAAC,IAAI;AACjD,QAAM,cAAc,OAAO,SAAS,aAAa,IAC7C,gBAAgB,QAChB,OAAO,OAAO,CAAC,IAAI,QAAQ,MAAM,OAAO,CAAC,IAAI;AACjD,QAAM,QACJ,KAAK,IACH,WACA,kBACA,mBACA,SACA,OAAO,SAAS,WAAW,IACvB,KAAK,IAAI,WAAW,IACpB,KAAK,IAAI,KAAK,IAAI,MAAM,OAAO,CAAC,IAAI,KAAK,GAAG,KAAK,IAAI,OAAO,OAAO,CAAC,IAAI,KAAK,CAAC,GAClF,OAAO,SAAS,WAAW,IACvB,KAAK,IAAI,WAAW,IACpB,KAAK,IAAI,KAAK,IAAI,MAAM,OAAO,CAAC,IAAI,KAAK,GAAG,KAAK,IAAI,OAAO,OAAO,CAAC,IAAI,KAAK,CAAC,CAAC,KAChF;AACP,QAAM,cAAc,OAAO,SAAS,WAAW,IAC3C,cAAc,QACb,OAAO,OAAO,CAAC,IAAI,QAAS,QAAS,MAAM,OAAO,CAAC,IAAI,QAAS;AACrE,QAAM,cAAc,OAAO,SAAS,WAAW,IAC3C,cAAc,QACb,OAAO,OAAO,CAAC,IAAI,QAAS,QAAS,MAAM,OAAO,CAAC,IAAI,QAAS;AACrE,SACE,KAAK,IAAI,cAAc,WAAW,KAClC,YAAY,QAAQ,mBAAmB,QAAQ,oBAAoB,QAAQ,UAAU;AAEzF;AAGA,SAAS,YAAY,KAAkB;AACrC,iBAAe,IAAI,MAAM;AACzB,iBAAe,IAAI,QAAQ;AAC3B,iBAAe,IAAI,SAAS;AAC5B,MAAI,IAAI,SAAS,CAAC,IAAI,KAAK,IAAI,SAAS,CAAC,IAAI,GAAG;AAC9C,UAAM,IAAI,WAAW,oCAAoC;EAC3D;AACA,MAAI,IAAI,UAAU,CAAC,MAAM,KAAK,IAAI,UAAU,CAAC,MAAM,GAAG;AACpD,UAAM,IAAI,WAAW,+BAA+B;EACtD;AACF;AAGA,SAAS,eAAe,QAAmC;AACzD,MAAI,OAAO,WAAW,KAAK,CAAC,OAAO,SAAS,OAAO,CAAC,CAAC,KAAK,CAAC,OAAO,SAAS,OAAO,CAAC,CAAC,GAAG;AACrF,UAAM,IAAI,WAAW,qDAAqD;EAC5E;AACF;",
+  "names": []
+}

package/dist/culling/oriented-box-2d.d.ts

@@ -0,0 +1,22 @@
+/** A rectangle in a shared two-dimensional Euclidean coordinate space. */
+export type OrientedBox2D = {
+    /** Two finite center coordinates; arrays, typed arrays and core vectors are accepted. */
+    readonly center: Readonly<ArrayLike<number>>;
+    /** Nonnegative half-width and half-height, in the center's coordinate units. */
+    readonly halfSize: Readonly<ArrayLike<number>>;
+    /** Nonzero finite width-axis direction. Its magnitude does not change the size. */
+    readonly direction: Readonly<ArrayLike<number>>;
+};
+/** Test closed rectangles using their four separating axes, without allocating query objects.
+ * Directions are normalized locally; the height axis is the perpendicular of the width axis.
+ * Points and line segments (zero half-sizes) are supported. Inputs are never modified.
+ * Subtracting centers before projection avoids projecting large translated corner coordinates.
+ * @param first First rectangle in the shared coordinate space.
+ * @param second Second rectangle in the shared coordinate space.
+ * @param epsilon Nonnegative finite separation allowance in coordinate units; default zero.
+ * Touching edges/vertices intersect. A gap no larger than epsilon on every tested axis also
+ * intersects; this is a SAT allowance, not a Euclidean-distance query.
+ * @throws RangeError For invalid vectors, sizes, directions or epsilon.
+ */
+export declare function intersectOrientedBoxes2D(first: OrientedBox2D, second: OrientedBox2D, epsilon?: number): boolean;
+//# sourceMappingURL=oriented-box-2d.d.ts.map

package/dist/culling/oriented-box-2d.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"oriented-box-2d.d.ts","sourceRoot":"","sources":["../../src/culling/oriented-box-2d.ts"],"names":[],"mappings":"AAIA,0EAA0E;AAC1E,MAAM,MAAM,aAAa,GAAG;IAC1B,yFAAyF;IACzF,QAAQ,CAAC,MAAM,EAAE,QAAQ,CAAC,SAAS,CAAC,MAAM,CAAC,CAAC,CAAC;IAC7C,gFAAgF;IAChF,QAAQ,CAAC,QAAQ,EAAE,QAAQ,CAAC,SAAS,CAAC,MAAM,CAAC,CAAC,CAAC;IAC/C,mFAAmF;IACnF,QAAQ,CAAC,SAAS,EAAE,QAAQ,CAAC,SAAS,CAAC,MAAM,CAAC,CAAC,CAAC;CACjD,CAAC;AAEF;;;;;;;;;;GAUG;AACH,wBAAgB,wBAAwB,CACtC,KAAK,EAAE,aAAa,EACpB,MAAM,EAAE,aAAa,EACrB,OAAO,SAAI,GACV,OAAO,CAqET"}

package/dist/culling/oriented-box-2d.js

@@ -0,0 +1,86 @@
+// math.gl
+// SPDX-License-Identifier: MIT
+// SPDX-FileCopyrightText: Copyright (c) vis.gl contributors
+/** Test closed rectangles using their four separating axes, without allocating query objects.
+ * Directions are normalized locally; the height axis is the perpendicular of the width axis.
+ * Points and line segments (zero half-sizes) are supported. Inputs are never modified.
+ * Subtracting centers before projection avoids projecting large translated corner coordinates.
+ * @param first First rectangle in the shared coordinate space.
+ * @param second Second rectangle in the shared coordinate space.
+ * @param epsilon Nonnegative finite separation allowance in coordinate units; default zero.
+ * Touching edges/vertices intersect. A gap no larger than epsilon on every tested axis also
+ * intersects; this is a SAT allowance, not a Euclidean-distance query.
+ * @throws RangeError For invalid vectors, sizes, directions or epsilon.
+ */
+export function intersectOrientedBoxes2D(first, second, epsilon = 0) {
+    validateBox(first);
+    validateBox(second);
+    if (!Number.isFinite(epsilon) || epsilon < 0) {
+        throw new RangeError('Box separation allowance must be finite and nonnegative');
+    }
+    // Scale before normalizing so both huge and subnormal directions remain usable.
+    const firstDirectionScale = Math.max(Math.abs(first.direction[0]), Math.abs(first.direction[1]));
+    const secondDirectionScale = Math.max(Math.abs(second.direction[0]), Math.abs(second.direction[1]));
+    const firstWidthX = first.direction[0] / firstDirectionScale;
+    const firstWidthY = first.direction[1] / firstDirectionScale;
+    const secondWidthX = second.direction[0] / secondDirectionScale;
+    const secondWidthY = second.direction[1] / secondDirectionScale;
+    const firstLength = Math.hypot(firstWidthX, firstWidthY);
+    const secondLength = Math.hypot(secondWidthX, secondWidthY);
+    const firstX = firstWidthX / firstLength;
+    const firstY = firstWidthY / firstLength;
+    const secondX = secondWidthX / secondLength;
+    const secondY = secondWidthY / secondLength;
+    // Clamp unit-axis dot products against rounding slightly above one.
+    const cosine = Math.min(1, Math.abs(firstX * secondX + firstY * secondY));
+    const sine = Math.min(1, Math.abs(firstX * secondY - firstY * secondX));
+    return (overlapsOnAxis(first, second, firstX, firstY, first.halfSize[0], second.halfSize[0] * cosine, second.halfSize[1] * sine, epsilon) &&
+        overlapsOnAxis(first, second, -firstY, firstX, first.halfSize[1], second.halfSize[0] * sine, second.halfSize[1] * cosine, epsilon) &&
+        overlapsOnAxis(first, second, secondX, secondY, second.halfSize[0], first.halfSize[0] * cosine, first.halfSize[1] * sine, epsilon) &&
+        overlapsOnAxis(first, second, -secondY, secondX, second.halfSize[1], first.halfSize[0] * sine, first.halfSize[1] * cosine, epsilon));
+}
+/** Scale each projection independently so an irrelevant huge extent cannot erase a small gap. */
+function overlapsOnAxis(first, second, axisX, axisY, ownRadius, otherWidthRadius, otherHeightRadius, epsilon) {
+    const displacementX = second.center[0] - first.center[0];
+    const displacementY = second.center[1] - first.center[1];
+    // Keep a representable center difference before multiplying. If subtraction overflows,
+    // project the finite centers first; a zero axis component then stays exactly zero.
+    const projectionX = Number.isFinite(displacementX)
+        ? displacementX * axisX
+        : second.center[0] * axisX - first.center[0] * axisX;
+    const projectionY = Number.isFinite(displacementY)
+        ? displacementY * axisY
+        : second.center[1] * axisY - first.center[1] * axisY;
+    const scale = Math.max(ownRadius, otherWidthRadius, otherHeightRadius, epsilon, Number.isFinite(projectionX)
+        ? Math.abs(projectionX)
+        : Math.max(Math.abs(first.center[0] * axisX), Math.abs(second.center[0] * axisX)), Number.isFinite(projectionY)
+        ? Math.abs(projectionY)
+        : Math.max(Math.abs(first.center[1] * axisY), Math.abs(second.center[1] * axisY))) || 1;
+    const normalizedX = Number.isFinite(projectionX)
+        ? projectionX / scale
+        : (second.center[0] * axisX) / scale - (first.center[0] * axisX) / scale;
+    const normalizedY = Number.isFinite(projectionY)
+        ? projectionY / scale
+        : (second.center[1] * axisY) / scale - (first.center[1] * axisY) / scale;
+    return (Math.abs(normalizedX + normalizedY) <=
+        ownRadius / scale + otherWidthRadius / scale + otherHeightRadius / scale + epsilon / scale);
+}
+/** Validate structural boxes before entering the scalar query kernel. */
+function validateBox(box) {
+    validateVector(box.center);
+    validateVector(box.halfSize);
+    validateVector(box.direction);
+    if (box.halfSize[0] < 0 || box.halfSize[1] < 0) {
+        throw new RangeError('Box half-sizes must be nonnegative');
+    }
+    if (box.direction[0] === 0 && box.direction[1] === 0) {
+        throw new RangeError('Box direction must be nonzero');
+    }
+}
+/** Require exactly two finite components, without copying the vector. */
+function validateVector(vector) {
+    if (vector.length !== 2 || !Number.isFinite(vector[0]) || !Number.isFinite(vector[1])) {
+        throw new RangeError('Box vectors must have exactly two finite components');
+    }
+}
+//# sourceMappingURL=oriented-box-2d.js.map

package/dist/culling/oriented-box-2d.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"oriented-box-2d.js","sourceRoot":"","sources":["../../src/culling/oriented-box-2d.ts"],"names":[],"mappings":"AAAA,UAAU;AACV,+BAA+B;AAC/B,4DAA4D;AAY5D;;;;;;;;;;GAUG;AACH,MAAM,UAAU,wBAAwB,CACtC,KAAoB,EACpB,MAAqB,EACrB,OAAO,GAAG,CAAC;IAEX,WAAW,CAAC,KAAK,CAAC,CAAC;IACnB,WAAW,CAAC,MAAM,CAAC,CAAC;IACpB,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,OAAO,CAAC,IAAI,OAAO,GAAG,CAAC,EAAE,CAAC;QAC7C,MAAM,IAAI,UAAU,CAAC,yDAAyD,CAAC,CAAC;IAClF,CAAC;IAED,gFAAgF;IAChF,MAAM,mBAAmB,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,SAAS,CAAC,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,SAAS,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IACjG,MAAM,oBAAoB,GAAG,IAAI,CAAC,GAAG,CACnC,IAAI,CAAC,GAAG,CAAC,MAAM,CAAC,SAAS,CAAC,CAAC,CAAC,CAAC,EAC7B,IAAI,CAAC,GAAG,CAAC,MAAM,CAAC,SAAS,CAAC,CAAC,CAAC,CAAC,CAC9B,CAAC;IACF,MAAM,WAAW,GAAG,KAAK,CAAC,SAAS,CAAC,CAAC,CAAC,GAAG,mBAAmB,CAAC;IAC7D,MAAM,WAAW,GAAG,KAAK,CAAC,SAAS,CAAC,CAAC,CAAC,GAAG,mBAAmB,CAAC;IAC7D,MAAM,YAAY,GAAG,MAAM,CAAC,SAAS,CAAC,CAAC,CAAC,GAAG,oBAAoB,CAAC;IAChE,MAAM,YAAY,GAAG,MAAM,CAAC,SAAS,CAAC,CAAC,CAAC,GAAG,oBAAoB,CAAC;IAChE,MAAM,WAAW,GAAG,IAAI,CAAC,KAAK,CAAC,WAAW,EAAE,WAAW,CAAC,CAAC;IACzD,MAAM,YAAY,GAAG,IAAI,CAAC,KAAK,CAAC,YAAY,EAAE,YAAY,CAAC,CAAC;IAC5D,MAAM,MAAM,GAAG,WAAW,GAAG,WAAW,CAAC;IACzC,MAAM,MAAM,GAAG,WAAW,GAAG,WAAW,CAAC;IACzC,MAAM,OAAO,GAAG,YAAY,GAAG,YAAY,CAAC;IAC5C,MAAM,OAAO,GAAG,YAAY,GAAG,YAAY,CAAC;IAE5C,oEAAoE;IACpE,MAAM,MAAM,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,MAAM,GAAG,OAAO,GAAG,MAAM,GAAG,OAAO,CAAC,CAAC,CAAC;IAC1E,MAAM,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,MAAM,GAAG,OAAO,GAAG,MAAM,GAAG,OAAO,CAAC,CAAC,CAAC;IACxE,OAAO,CACL,cAAc,CACZ,KAAK,EACL,MAAM,EACN,MAAM,EACN,MAAM,EACN,KAAK,CAAC,QAAQ,CAAC,CAAC,CAAC,EACjB,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,GAAG,MAAM,EAC3B,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,GAAG,IAAI,EACzB,OAAO,CACR;QACD,cAAc,CACZ,KAAK,EACL,MAAM,EACN,CAAC,MAAM,EACP,MAAM,EACN,KAAK,CAAC,QAAQ,CAAC,CAAC,CAAC,EACjB,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,GAAG,IAAI,EACzB,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,GAAG,MAAM,EAC3B,OAAO,CACR;QACD,cAAc,CACZ,KAAK,EACL,MAAM,EACN,OAAO,EACP,OAAO,EACP,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,EAClB,KAAK,CAAC,QAAQ,CAAC,CAAC,CAAC,GAAG,MAAM,EAC1B,KAAK,CAAC,QAAQ,CAAC,CAAC,CAAC,GAAG,IAAI,EACxB,OAAO,CACR;QACD,cAAc,CACZ,KAAK,EACL,MAAM,EACN,CAAC,OAAO,EACR,OAAO,EACP,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,EAClB,KAAK,CAAC,QAAQ,CAAC,CAAC,CAAC,GAAG,IAAI,EACxB,KAAK,CAAC,QAAQ,CAAC,CAAC,CAAC,GAAG,MAAM,EAC1B,OAAO,CACR,CACF,CAAC;AACJ,CAAC;AAED,iGAAiG;AACjG,SAAS,cAAc,CACrB,KAAoB,EACpB,MAAqB,EACrB,KAAa,EACb,KAAa,EACb,SAAiB,EACjB,gBAAwB,EACxB,iBAAyB,EACzB,OAAe;IAEf,MAAM,aAAa,GAAG,MAAM,CAAC,MAAM,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC;IACzD,MAAM,aAAa,GAAG,MAAM,CAAC,MAAM,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC;IACzD,uFAAuF;IACvF,mFAAmF;IACnF,MAAM,WAAW,GAAG,MAAM,CAAC,QAAQ,CAAC,aAAa,CAAC;QAChD,CAAC,CAAC,aAAa,GAAG,KAAK;QACvB,CAAC,CAAC,MAAM,CAAC,MAAM,CAAC,CAAC,CAAC,GAAG,KAAK,GAAG,KAAK,CAAC,MAAM,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC;IACvD,MAAM,WAAW,GAAG,MAAM,CAAC,QAAQ,CAAC,aAAa,CAAC;QAChD,CAAC,CAAC,aAAa,GAAG,KAAK;QACvB,CAAC,CAAC,MAAM,CAAC,MAAM,CAAC,CAAC,CAAC,GAAG,KAAK,GAAG,KAAK,CAAC,MAAM,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC;IACvD,MAAM,KAAK,GACT,IAAI,CAAC,GAAG,CACN,SAAS,EACT,gBAAgB,EAChB,iBAAiB,EACjB,OAAO,EACP,MAAM,CAAC,QAAQ,CAAC,WAAW,CAAC;QAC1B,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,WAAW,CAAC;QACvB,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,MAAM,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,MAAM,CAAC,MAAM,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,EACnF,MAAM,CAAC,QAAQ,CAAC,WAAW,CAAC;QAC1B,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,WAAW,CAAC;QACvB,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,MAAM,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,MAAM,CAAC,MAAM,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CACpF,IAAI,CAAC,CAAC;IACT,MAAM,WAAW,GAAG,MAAM,CAAC,QAAQ,CAAC,WAAW,CAAC;QAC9C,CAAC,CAAC,WAAW,GAAG,KAAK;QACrB,CAAC,CAAC,CAAC,MAAM,CAAC,MAAM,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,GAAG,KAAK,GAAG,CAAC,KAAK,CAAC,MAAM,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,GAAG,KAAK,CAAC;IAC3E,MAAM,WAAW,GAAG,MAAM,CAAC,QAAQ,CAAC,WAAW,CAAC;QAC9C,CAAC,CAAC,WAAW,GAAG,KAAK;QACrB,CAAC,CAAC,CAAC,MAAM,CAAC,MAAM,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,GAAG,KAAK,GAAG,CAAC,KAAK,CAAC,MAAM,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,GAAG,KAAK,CAAC;IAC3E,OAAO,CACL,IAAI,CAAC,GAAG,CAAC,WAAW,GAAG,WAAW,CAAC;QACnC,SAAS,GAAG,KAAK,GAAG,gBAAgB,GAAG,KAAK,GAAG,iBAAiB,GAAG,KAAK,GAAG,OAAO,GAAG,KAAK,CAC3F,CAAC;AACJ,CAAC;AAED,yEAAyE;AACzE,SAAS,WAAW,CAAC,GAAkB;IACrC,cAAc,CAAC,GAAG,CAAC,MAAM,CAAC,CAAC;IAC3B,cAAc,CAAC,GAAG,CAAC,QAAQ,CAAC,CAAC;IAC7B,cAAc,CAAC,GAAG,CAAC,SAAS,CAAC,CAAC;IAC9B,IAAI,GAAG,CAAC,QAAQ,CAAC,CAAC,CAAC,GAAG,CAAC,IAAI,GAAG,CAAC,QAAQ,CAAC,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC;QAC/C,MAAM,IAAI,UAAU,CAAC,oCAAoC,CAAC,CAAC;IAC7D,CAAC;IACD,IAAI,GAAG,CAAC,SAAS,CAAC,CAAC,CAAC,KAAK,CAAC,IAAI,GAAG,CAAC,SAAS,CAAC,CAAC,CAAC,KAAK,CAAC,EAAE,CAAC;QACrD,MAAM,IAAI,UAAU,CAAC,+BAA+B,CAAC,CAAC;IACxD,CAAC;AACH,CAAC;AAED,yEAAyE;AACzE,SAAS,cAAc,CAAC,MAAmC;IACzD,IAAI,MAAM,CAAC,MAAM,KAAK,CAAC,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC;QACtF,MAAM,IAAI,UAAU,CAAC,qDAAqD,CAAC,CAAC;IAC9E,CAAC;AACH,CAAC"}