@math.gl/culling 4.2.0-alpha.3 → 4.2.0-alpha.5

package/dist/lib/shapes/capsule-shape.js

@@ -0,0 +1,151 @@
+// math.gl
+// SPDX-License-Identifier: MIT
+// Copyright (c) vis.gl contributors
+import { closestDistanceToCircleArc2D, closestDistanceToSegment2D, ImplicitShapeBase, solveQuadratic } from "./implicit-shape.js";
+const EPSILON = 1e-6;
+/** Analytic glTF capsule: the convex hull of two spheres on the Y axis. */
+export class CapsuleShape extends ImplicitShapeBase {
+    constructor(props = {}) {
+        super(props.matrix);
+        this.type = 'capsule';
+        this.height = props.height ?? 1;
+        this.radiusBottom = props.radiusBottom ?? 0.5;
+        this.radiusTop = props.radiusTop ?? 0.5;
+        if (!Number.isFinite(this.height) || this.height <= 0)
+            throw new Error('Capsule height must be greater than zero');
+        if (![this.radiusBottom, this.radiusTop].every(value => Number.isFinite(value) && value >= 0))
+            throw new Error('Capsule radii must be non-negative');
+        if (this.radiusBottom === 0 && this.radiusTop === 0)
+            throw new Error('At least one capsule radius must be greater than zero');
+    }
+    clone() {
+        return new CapsuleShape({
+            height: this.height,
+            radiusBottom: this.radiusBottom,
+            radiusTop: this.radiusTop,
+            matrix: this.matrix
+        });
+    }
+    equals(shape) {
+        return (shape instanceof CapsuleShape &&
+            super.equals(shape) &&
+            this.height === shape.height &&
+            this.radiusBottom === shape.radiusBottom &&
+            this.radiusTop === shape.radiusTop);
+    }
+    containsLocalPoint(point) {
+        const profile = this.getProfile();
+        const radial = Math.hypot(point[0], point[2]);
+        if (profile.sphere)
+            return (Math.hypot(radial, point[1] - profile.sphere.center) <= profile.sphere.radius + EPSILON);
+        if (point[1] < profile.bottom.y)
+            return Math.hypot(radial, point[1] + this.height / 2) <= this.radiusBottom + EPSILON;
+        if (point[1] > profile.top.y)
+            return Math.hypot(radial, point[1] - this.height / 2) <= this.radiusTop + EPSILON;
+        const t = (point[1] - profile.bottom.y) / (profile.top.y - profile.bottom.y);
+        return (radial <= profile.bottom.radius + (profile.top.radius - profile.bottom.radius) * t + EPSILON);
+    }
+    distanceToLocalPoint(point) {
+        if (this.containsLocalPoint(point))
+            return 0;
+        const profile = this.getProfile();
+        const radial = Math.hypot(point[0], point[2]);
+        if (profile.sphere)
+            return Math.max(0, Math.hypot(radial, point[1] - profile.sphere.center) - profile.sphere.radius);
+        const delta = this.radiusTop - this.radiusBottom;
+        const sideNormalY = -delta / this.height;
+        const bottomArc = closestDistanceToCircleArc2D(radial, point[1], -this.height / 2, this.radiusBottom, -1, sideNormalY);
+        const topArc = closestDistanceToCircleArc2D(radial, point[1], this.height / 2, this.radiusTop, sideNormalY, 1);
+        const side = closestDistanceToSegment2D(radial, point[1], profile.bottom.radius, profile.bottom.y, profile.top.radius, profile.top.y);
+        return Math.min(bottomArc, topArc, side);
+    }
+    supportLocal(direction) {
+        const length = Math.hypot(...direction);
+        const unit = length ? direction.map(value => value / length) : [1, 0, 0];
+        const bottom = [
+            unit[0] * this.radiusBottom,
+            -this.height / 2 + unit[1] * this.radiusBottom,
+            unit[2] * this.radiusBottom
+        ];
+        const top = [
+            unit[0] * this.radiusTop,
+            this.height / 2 + unit[1] * this.radiusTop,
+            unit[2] * this.radiusTop
+        ];
+        return dot(top, direction) >= dot(bottom, direction) ? top : bottom;
+    }
+    intersectLocalRay(origin, direction) {
+        const profile = this.getProfile();
+        if (profile.sphere)
+            return intersectSphere(origin, direction, profile.sphere.center, profile.sphere.radius);
+        const candidates = [];
+        const slope = (profile.top.radius - profile.bottom.radius) / (profile.top.y - profile.bottom.y);
+        const intercept = profile.bottom.radius - slope * profile.bottom.y;
+        const radialAtOrigin = intercept + slope * origin[1];
+        const radialAlongRay = slope * direction[1];
+        for (const distance of solveQuadratic(direction[0] ** 2 + direction[2] ** 2 - radialAlongRay ** 2, 2 * (origin[0] * direction[0] + origin[2] * direction[2] - radialAtOrigin * radialAlongRay), origin[0] ** 2 + origin[2] ** 2 - radialAtOrigin ** 2)) {
+            const y = origin[1] + distance * direction[1];
+            if (distance >= 0 && y >= profile.bottom.y - EPSILON && y <= profile.top.y + EPSILON) {
+                const x = origin[0] + distance * direction[0];
+                const z = origin[2] + distance * direction[2];
+                const radial = Math.hypot(x, z) || 1;
+                const normal = [x / radial, -slope, z / radial];
+                const length = Math.hypot(...normal);
+                candidates.push({ distance, normal: normal.map(value => value / length) });
+            }
+        }
+        for (const bottomHit of intersectSphereRoots(origin, direction, -this.height / 2, this.radiusBottom)) {
+            const y = origin[1] + bottomHit.distance * direction[1];
+            if (y <= profile.bottom.y + EPSILON)
+                candidates.push(bottomHit);
+        }
+        for (const topHit of intersectSphereRoots(origin, direction, this.height / 2, this.radiusTop)) {
+            const y = origin[1] + topHit.distance * direction[1];
+            if (y >= profile.top.y - EPSILON)
+                candidates.push(topHit);
+        }
+        return candidates.sort((a, b) => a.distance - b.distance)[0];
+    }
+    getProfile() {
+        const delta = this.radiusTop - this.radiusBottom;
+        if (Math.abs(delta) >= this.height) {
+            const topWins = this.radiusTop >= this.radiusBottom;
+            return {
+                sphere: {
+                    center: topWins ? this.height / 2 : -this.height / 2,
+                    radius: topWins ? this.radiusTop : this.radiusBottom
+                },
+                bottom: { radius: 0, y: 0 },
+                top: { radius: 0, y: 0 }
+            };
+        }
+        const sinSlope = delta / this.height;
+        const radialNormal = Math.sqrt(1 - sinSlope * sinSlope);
+        const normalY = -sinSlope;
+        return {
+            bottom: {
+                radius: this.radiusBottom * radialNormal,
+                y: -this.height / 2 + this.radiusBottom * normalY
+            },
+            top: { radius: this.radiusTop * radialNormal, y: this.height / 2 + this.radiusTop * normalY }
+        };
+    }
+}
+function intersectSphere(origin, direction, centerY, radius) {
+    return intersectSphereRoots(origin, direction, centerY, radius)[0];
+}
+function intersectSphereRoots(origin, direction, centerY, radius) {
+    if (radius === 0)
+        return [];
+    const relative = [origin[0], origin[1] - centerY, origin[2]];
+    return solveQuadratic(dot(direction, direction), 2 * dot(relative, direction), dot(relative, relative) - radius * radius)
+        .filter(distance => distance >= 0)
+        .map(distance => ({
+        distance,
+        normal: relative.map((value, i) => (value + distance * direction[i]) / radius)
+    }));
+}
+function dot(a, b) {
+    return a[0] * b[0] + a[1] * b[1] + a[2] * b[2];
+}
+//# sourceMappingURL=capsule-shape.js.map

package/dist/lib/shapes/capsule-shape.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"capsule-shape.js","sourceRoot":"","sources":["../../../src/lib/shapes/capsule-shape.ts"],"names":[],"mappings":"AAAA,UAAU;AACV,+BAA+B;AAC/B,oCAAoC;AAEpC,OAAO,EACL,4BAA4B,EAC5B,0BAA0B,EAC1B,iBAAiB,EACjB,cAAc,EAGf,4BAAyB;AAE1B,MAAM,OAAO,GAAG,IAAI,CAAC;AASrB,2EAA2E;AAC3E,MAAM,OAAO,YAAa,SAAQ,iBAAiB;IAMjD,YAAY,QAA2B,EAAE;QACvC,KAAK,CAAC,KAAK,CAAC,MAAM,CAAC,CAAC;QANb,SAAI,GAAG,SAAkB,CAAC;QAOjC,IAAI,CAAC,MAAM,GAAG,KAAK,CAAC,MAAM,IAAI,CAAC,CAAC;QAChC,IAAI,CAAC,YAAY,GAAG,KAAK,CAAC,YAAY,IAAI,GAAG,CAAC;QAC9C,IAAI,CAAC,SAAS,GAAG,KAAK,CAAC,SAAS,IAAI,GAAG,CAAC;QACxC,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,IAAI,CAAC,MAAM,CAAC,IAAI,IAAI,CAAC,MAAM,IAAI,CAAC;YACnD,MAAM,IAAI,KAAK,CAAC,0CAA0C,CAAC,CAAC;QAC9D,IAAI,CAAC,CAAC,IAAI,CAAC,YAAY,EAAE,IAAI,CAAC,SAAS,CAAC,CAAC,KAAK,CAAC,KAAK,CAAC,EAAE,CAAC,MAAM,CAAC,QAAQ,CAAC,KAAK,CAAC,IAAI,KAAK,IAAI,CAAC,CAAC;YAC3F,MAAM,IAAI,KAAK,CAAC,oCAAoC,CAAC,CAAC;QACxD,IAAI,IAAI,CAAC,YAAY,KAAK,CAAC,IAAI,IAAI,CAAC,SAAS,KAAK,CAAC;YACjD,MAAM,IAAI,KAAK,CAAC,uDAAuD,CAAC,CAAC;IAC7E,CAAC;IAED,KAAK;QACH,OAAO,IAAI,YAAY,CAAC;YACtB,MAAM,EAAE,IAAI,CAAC,MAAM;YACnB,YAAY,EAAE,IAAI,CAAC,YAAY;YAC/B,SAAS,EAAE,IAAI,CAAC,SAAS;YACzB,MAAM,EAAE,IAAI,CAAC,MAAM;SACpB,CAAC,CAAC;IACL,CAAC;IACQ,MAAM,CAAC,KAAoB;QAClC,OAAO,CACL,KAAK,YAAY,YAAY;YAC7B,KAAK,CAAC,MAAM,CAAC,KAAK,CAAC;YACnB,IAAI,CAAC,MAAM,KAAK,KAAK,CAAC,MAAM;YAC5B,IAAI,CAAC,YAAY,KAAK,KAAK,CAAC,YAAY;YACxC,IAAI,CAAC,SAAS,KAAK,KAAK,CAAC,SAAS,CACnC,CAAC;IACJ,CAAC;IACD,kBAAkB,CAAC,KAAwB;QACzC,MAAM,OAAO,GAAG,IAAI,CAAC,UAAU,EAAE,CAAC;QAClC,MAAM,MAAM,GAAG,IAAI,CAAC,KAAK,CAAC,KAAK,CAAC,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC;QAC9C,IAAI,OAAO,CAAC,MAAM;YAChB,OAAO,CACL,IAAI,CAAC,KAAK,CAAC,MAAM,EAAE,KAAK,CAAC,CAAC,CAAC,GAAG,OAAO,CAAC,MAAM,CAAC,MAAM,CAAC,IAAI,OAAO,CAAC,MAAM,CAAC,MAAM,GAAG,OAAO,CACxF,CAAC;QACJ,IAAI,KAAK,CAAC,CAAC,CAAC,GAAG,OAAO,CAAC,MAAM,CAAC,CAAC;YAC7B,OAAO,IAAI,CAAC,KAAK,CAAC,MAAM,EAAE,KAAK,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,MAAM,GAAG,CAAC,CAAC,IAAI,IAAI,CAAC,YAAY,GAAG,OAAO,CAAC;QACvF,IAAI,KAAK,CAAC,CAAC,CAAC,GAAG,OAAO,CAAC,GAAG,CAAC,CAAC;YAC1B,OAAO,IAAI,CAAC,KAAK,CAAC,MAAM,EAAE,KAAK,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,MAAM,GAAG,CAAC,CAAC,IAAI,IAAI,CAAC,SAAS,GAAG,OAAO,CAAC;QACpF,MAAM,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,OAAO,CAAC,MAAM,CAAC,CAAC,CAAC,GAAG,CAAC,OAAO,CAAC,GAAG,CAAC,CAAC,GAAG,OAAO,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC;QAC7E,OAAO,CACL,MAAM,IAAI,OAAO,CAAC,MAAM,CAAC,MAAM,GAAG,CAAC,OAAO,CAAC,GAAG,CAAC,MAAM,GAAG,OAAO,CAAC,MAAM,CAAC,MAAM,CAAC,GAAG,CAAC,GAAG,OAAO,CAC7F,CAAC;IACJ,CAAC;IACD,oBAAoB,CAAC,KAAwB;QAC3C,IAAI,IAAI,CAAC,kBAAkB,CAAC,KAAK,CAAC;YAAE,OAAO,CAAC,CAAC;QAC7C,MAAM,OAAO,GAAG,IAAI,CAAC,UAAU,EAAE,CAAC;QAClC,MAAM,MAAM,GAAG,IAAI,CAAC,KAAK,CAAC,KAAK,CAAC,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC;QAC9C,IAAI,OAAO,CAAC,MAAM;YAChB,OAAO,IAAI,CAAC,GAAG,CACb,CAAC,EACD,IAAI,CAAC,KAAK,CAAC,MAAM,EAAE,KAAK,CAAC,CAAC,CAAC,GAAG,OAAO,CAAC,MAAM,CAAC,MAAM,CAAC,GAAG,OAAO,CAAC,MAAM,CAAC,MAAM,CAC7E,CAAC;QACJ,MAAM,KAAK,GAAG,IAAI,CAAC,SAAS,GAAG,IAAI,CAAC,YAAY,CAAC;QACjD,MAAM,WAAW,GAAG,CAAC,KAAK,GAAG,IAAI,CAAC,MAAM,CAAC;QACzC,MAAM,SAAS,GAAG,4BAA4B,CAC5C,MAAM,EACN,KAAK,CAAC,CAAC,CAAC,EACR,CAAC,IAAI,CAAC,MAAM,GAAG,CAAC,EAChB,IAAI,CAAC,YAAY,EACjB,CAAC,CAAC,EACF,WAAW,CACZ,CAAC;QACF,MAAM,MAAM,GAAG,4BAA4B,CACzC,MAAM,EACN,KAAK,CAAC,CAAC,CAAC,EACR,IAAI,CAAC,MAAM,GAAG,CAAC,EACf,IAAI,CAAC,SAAS,EACd,WAAW,EACX,CAAC,CACF,CAAC;QACF,MAAM,IAAI,GAAG,0BAA0B,CACrC,MAAM,EACN,KAAK,CAAC,CAAC,CAAC,EACR,OAAO,CAAC,MAAM,CAAC,MAAM,EACrB,OAAO,CAAC,MAAM,CAAC,CAAC,EAChB,OAAO,CAAC,GAAG,CAAC,MAAM,EAClB,OAAO,CAAC,GAAG,CAAC,CAAC,CACd,CAAC;QACF,OAAO,IAAI,CAAC,GAAG,CAAC,SAAS,EAAE,MAAM,EAAE,IAAI,CAAC,CAAC;IAC3C,CAAC;IACD,YAAY,CAAC,SAA4B;QACvC,MAAM,MAAM,GAAG,IAAI,CAAC,KAAK,CAAC,GAAG,SAAS,CAAC,CAAC;QACxC,MAAM,IAAI,GAAG,MAAM,CAAC,CAAC,CAAC,SAAS,CAAC,GAAG,CAAC,KAAK,CAAC,EAAE,CAAC,KAAK,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;QACzE,MAAM,MAAM,GAAG;YACb,IAAI,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,YAAY;YAC3B,CAAC,IAAI,CAAC,MAAM,GAAG,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,YAAY;YAC9C,IAAI,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,YAAY;SAC5B,CAAC;QACF,MAAM,GAAG,GAAG;YACV,IAAI,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,SAAS;YACxB,IAAI,CAAC,MAAM,GAAG,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,SAAS;YAC1C,IAAI,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,SAAS;SACzB,CAAC;QACF,OAAO,GAAG,CAAC,GAAG,EAAE,SAAS,CAAC,IAAI,GAAG,CAAC,MAAM,EAAE,SAAS,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,MAAM,CAAC;IACtE,CAAC;IACD,iBAAiB,CACf,MAAyB,EACzB,SAA4B;QAE5B,MAAM,OAAO,GAAG,IAAI,CAAC,UAAU,EAAE,CAAC;QAClC,IAAI,OAAO,CAAC,MAAM;YAChB,OAAO,eAAe,CAAC,MAAM,EAAE,SAAS,EAAE,OAAO,CAAC,MAAM,CAAC,MAAM,EAAE,OAAO,CAAC,MAAM,CAAC,MAAM,CAAC,CAAC;QAC1F,MAAM,UAAU,GAA2B,EAAE,CAAC;QAC9C,MAAM,KAAK,GAAG,CAAC,OAAO,CAAC,GAAG,CAAC,MAAM,GAAG,OAAO,CAAC,MAAM,CAAC,MAAM,CAAC,GAAG,CAAC,OAAO,CAAC,GAAG,CAAC,CAAC,GAAG,OAAO,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC;QAChG,MAAM,SAAS,GAAG,OAAO,CAAC,MAAM,CAAC,MAAM,GAAG,KAAK,GAAG,OAAO,CAAC,MAAM,CAAC,CAAC,CAAC;QACnE,MAAM,cAAc,GAAG,SAAS,GAAG,KAAK,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC;QACrD,MAAM,cAAc,GAAG,KAAK,GAAG,SAAS,CAAC,CAAC,CAAC,CAAC;QAC5C,KAAK,MAAM,QAAQ,IAAI,cAAc,CACnC,SAAS,CAAC,CAAC,CAAC,IAAI,CAAC,GAAG,SAAS,CAAC,CAAC,CAAC,IAAI,CAAC,GAAG,cAAc,IAAI,CAAC,EAC3D,CAAC,GAAG,CAAC,MAAM,CAAC,CAAC,CAAC,GAAG,SAAS,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,GAAG,SAAS,CAAC,CAAC,CAAC,GAAG,cAAc,GAAG,cAAc,CAAC,EAC3F,MAAM,CAAC,CAAC,CAAC,IAAI,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,IAAI,CAAC,GAAG,cAAc,IAAI,CAAC,CACtD,EAAE,CAAC;YACF,MAAM,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,GAAG,QAAQ,GAAG,SAAS,CAAC,CAAC,CAAC,CAAC;YAC9C,IAAI,QAAQ,IAAI,CAAC,IAAI,CAAC,IAAI,OAAO,CAAC,MAAM,CAAC,CAAC,GAAG,OAAO,IAAI,CAAC,IAAI,OAAO,CAAC,GAAG,CAAC,CAAC,GAAG,OAAO,EAAE,CAAC;gBACrF,MAAM,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,GAAG,QAAQ,GAAG,SAAS,CAAC,CAAC,CAAC,CAAC;gBAC9C,MAAM,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,GAAG,QAAQ,GAAG,SAAS,CAAC,CAAC,CAAC,CAAC;gBAC9C,MAAM,MAAM,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC,IAAI,CAAC,CAAC;gBACrC,MAAM,MAAM,GAAG,CAAC,CAAC,GAAG,MAAM,EAAE,CAAC,KAAK,EAAE,CAAC,GAAG,MAAM,CAAC,CAAC;gBAChD,MAAM,MAAM,GAAG,IAAI,CAAC,KAAK,CAAC,GAAG,MAAM,CAAC,CAAC;gBACrC,UAAU,CAAC,IAAI,CAAC,EAAC,QAAQ,EAAE,MAAM,EAAE,MAAM,CAAC,GAAG,CAAC,KAAK,CAAC,EAAE,CAAC,KAAK,GAAG,MAAM,CAAC,EAAC,CAAC,CAAC;YAC3E,CAAC;QACH,CAAC;QACD,KAAK,MAAM,SAAS,IAAI,oBAAoB,CAC1C,MAAM,EACN,SAAS,EACT,CAAC,IAAI,CAAC,MAAM,GAAG,CAAC,EAChB,IAAI,CAAC,YAAY,CAClB,EAAE,CAAC;YACF,MAAM,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,GAAG,SAAS,CAAC,QAAQ,GAAG,SAAS,CAAC,CAAC,CAAC,CAAC;YACxD,IAAI,CAAC,IAAI,OAAO,CAAC,MAAM,CAAC,CAAC,GAAG,OAAO;gBAAE,UAAU,CAAC,IAAI,CAAC,SAAS,CAAC,CAAC;QAClE,CAAC;QACD,KAAK,MAAM,MAAM,IAAI,oBAAoB,CAAC,MAAM,EAAE,SAAS,EAAE,IAAI,CAAC,MAAM,GAAG,CAAC,EAAE,IAAI,CAAC,SAAS,CAAC,EAAE,CAAC;YAC9F,MAAM,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC,QAAQ,GAAG,SAAS,CAAC,CAAC,CAAC,CAAC;YACrD,IAAI,CAAC,IAAI,OAAO,CAAC,GAAG,CAAC,CAAC,GAAG,OAAO;gBAAE,UAAU,CAAC,IAAI,CAAC,MAAM,CAAC,CAAC;QAC5D,CAAC;QACD,OAAO,UAAU,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,QAAQ,GAAG,CAAC,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC;IAC/D,CAAC;IAEO,UAAU;QAKhB,MAAM,KAAK,GAAG,IAAI,CAAC,SAAS,GAAG,IAAI,CAAC,YAAY,CAAC;QACjD,IAAI,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,IAAI,IAAI,CAAC,MAAM,EAAE,CAAC;YACnC,MAAM,OAAO,GAAG,IAAI,CAAC,SAAS,IAAI,IAAI,CAAC,YAAY,CAAC;YACpD,OAAO;gBACL,MAAM,EAAE;oBACN,MAAM,EAAE,OAAO,CAAC,CAAC,CAAC,IAAI,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,MAAM,GAAG,CAAC;oBACpD,MAAM,EAAE,OAAO,CAAC,CAAC,CAAC,IAAI,CAAC,SAAS,CAAC,CAAC,CAAC,IAAI,CAAC,YAAY;iBACrD;gBACD,MAAM,EAAE,EAAC,MAAM,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAC;gBACzB,GAAG,EAAE,EAAC,MAAM,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAC;aACvB,CAAC;QACJ,CAAC;QACD,MAAM,QAAQ,GAAG,KAAK,GAAG,IAAI,CAAC,MAAM,CAAC;QACrC,MAAM,YAAY,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,GAAG,QAAQ,GAAG,QAAQ,CAAC,CAAC;QACxD,MAAM,OAAO,GAAG,CAAC,QAAQ,CAAC;QAC1B,OAAO;YACL,MAAM,EAAE;gBACN,MAAM,EAAE,IAAI,CAAC,YAAY,GAAG,YAAY;gBACxC,CAAC,EAAE,CAAC,IAAI,CAAC,MAAM,GAAG,CAAC,GAAG,IAAI,CAAC,YAAY,GAAG,OAAO;aAClD;YACD,GAAG,EAAE,EAAC,MAAM,EAAE,IAAI,CAAC,SAAS,GAAG,YAAY,EAAE,CAAC,EAAE,IAAI,CAAC,MAAM,GAAG,CAAC,GAAG,IAAI,CAAC,SAAS,GAAG,OAAO,EAAC;SAC5F,CAAC;IACJ,CAAC;CACF;AAED,SAAS,eAAe,CACtB,MAAyB,EACzB,SAA4B,EAC5B,OAAe,EACf,MAAc;IAEd,OAAO,oBAAoB,CAAC,MAAM,EAAE,SAAS,EAAE,OAAO,EAAE,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC;AACrE,CAAC;AAED,SAAS,oBAAoB,CAC3B,MAAyB,EACzB,SAA4B,EAC5B,OAAe,EACf,MAAc;IAEd,IAAI,MAAM,KAAK,CAAC;QAAE,OAAO,EAAE,CAAC;IAC5B,MAAM,QAAQ,GAAG,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC,CAAC,GAAG,OAAO,EAAE,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC;IAC7D,OAAO,cAAc,CACnB,GAAG,CAAC,SAAS,EAAE,SAAS,CAAC,EACzB,CAAC,GAAG,GAAG,CAAC,QAAQ,EAAE,SAAS,CAAC,EAC5B,GAAG,CAAC,QAAQ,EAAE,QAAQ,CAAC,GAAG,MAAM,GAAG,MAAM,CAC1C;SACE,MAAM,CAAC,QAAQ,CAAC,EAAE,CAAC,QAAQ,IAAI,CAAC,CAAC;SACjC,GAAG,CAAC,QAAQ,CAAC,EAAE,CAAC,CAAC;QAChB,QAAQ;QACR,MAAM,EAAE,QAAQ,CAAC,GAAG,CAAC,CAAC,KAAK,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,KAAK,GAAG,QAAQ,GAAG,SAAS,CAAC,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC;KAC/E,CAAC,CAAC,CAAC;AACR,CAAC;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/lib/shapes/cylinder-shape.d.ts

@@ -0,0 +1,22 @@
+import { ImplicitShapeBase, type ImplicitShape, type LocalRayIntersection } from "./implicit-shape.js";
+export type CylinderShapeProps = {
+    height?: number;
+    radiusBottom?: number;
+    radiusTop?: number;
+    matrix?: readonly number[];
+};
+/** Analytic closed, Y-aligned glTF cylinder shape, including tapered cylinders. */
+export declare class CylinderShape extends ImplicitShapeBase {
+    readonly type: "cylinder";
+    readonly height: number;
+    readonly radiusBottom: number;
+    readonly radiusTop: number;
+    constructor(props?: CylinderShapeProps);
+    clone(): CylinderShape;
+    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=cylinder-shape.d.ts.map

package/dist/lib/shapes/cylinder-shape.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"cylinder-shape.d.ts","sourceRoot":"","sources":["../../../src/lib/shapes/cylinder-shape.ts"],"names":[],"mappings":"AAIA,OAAO,EAEL,iBAAiB,EAEjB,KAAK,aAAa,EAClB,KAAK,oBAAoB,EAC1B,4BAAyB;AAI1B,MAAM,MAAM,kBAAkB,GAAG;IAC/B,MAAM,CAAC,EAAE,MAAM,CAAC;IAChB,YAAY,CAAC,EAAE,MAAM,CAAC;IACtB,SAAS,CAAC,EAAE,MAAM,CAAC;IACnB,MAAM,CAAC,EAAE,SAAS,MAAM,EAAE,CAAC;CAC5B,CAAC;AAEF,mFAAmF;AACnF,qBAAa,aAAc,SAAQ,iBAAiB;IAClD,QAAQ,CAAC,IAAI,EAAG,UAAU,CAAU;IACpC,QAAQ,CAAC,MAAM,EAAE,MAAM,CAAC;IACxB,QAAQ,CAAC,YAAY,EAAE,MAAM,CAAC;IAC9B,QAAQ,CAAC,SAAS,EAAE,MAAM,CAAC;gBAEf,KAAK,GAAE,kBAAuB;IAa1C,KAAK,IAAI,aAAa;IAQb,MAAM,CAAC,KAAK,EAAE,aAAa,GAAG,OAAO;IAS9C,kBAAkB,CAAC,KAAK,EAAE,SAAS,MAAM,EAAE,GAAG,OAAO;IAOrD,oBAAoB,CAAC,KAAK,EAAE,SAAS,MAAM,EAAE,GAAG,MAAM;IAgBtD,YAAY,CAAC,SAAS,EAAE,SAAS,MAAM,EAAE,GAAG,SAAS,MAAM,EAAE;IAW7D,iBAAiB,CACf,MAAM,EAAE,SAAS,MAAM,EAAE,EACzB,SAAS,EAAE,SAAS,MAAM,EAAE,GAC3B,oBAAoB,GAAG,SAAS;CAqCpC"}

package/dist/lib/shapes/cylinder-shape.js

@@ -0,0 +1,99 @@
+// math.gl
+// SPDX-License-Identifier: MIT
+// Copyright (c) vis.gl contributors
+import { closestDistanceToSegment2D, ImplicitShapeBase, solveQuadratic } from "./implicit-shape.js";
+const EPSILON = 1e-6;
+/** Analytic closed, Y-aligned glTF cylinder shape, including tapered cylinders. */
+export class CylinderShape extends ImplicitShapeBase {
+    constructor(props = {}) {
+        super(props.matrix);
+        this.type = 'cylinder';
+        this.height = props.height ?? 2;
+        this.radiusBottom = props.radiusBottom ?? 0.5;
+        this.radiusTop = props.radiusTop ?? 0.5;
+        if (!Number.isFinite(this.height) || this.height <= 0)
+            throw new Error('Cylinder height must be greater than zero');
+        if (![this.radiusBottom, this.radiusTop].every(value => Number.isFinite(value) && value >= 0))
+            throw new Error('Cylinder radii must be non-negative');
+        if (this.radiusBottom === 0 && this.radiusTop === 0)
+            throw new Error('At least one cylinder radius must be greater than zero');
+    }
+    clone() {
+        return new CylinderShape({
+            height: this.height,
+            radiusBottom: this.radiusBottom,
+            radiusTop: this.radiusTop,
+            matrix: this.matrix
+        });
+    }
+    equals(shape) {
+        return (shape instanceof CylinderShape &&
+            super.equals(shape) &&
+            this.height === shape.height &&
+            this.radiusBottom === shape.radiusBottom &&
+            this.radiusTop === shape.radiusTop);
+    }
+    containsLocalPoint(point) {
+        const half = this.height / 2;
+        if (point[1] < -half - EPSILON || point[1] > half + EPSILON)
+            return false;
+        const t = (point[1] + half) / this.height;
+        const radius = this.radiusBottom + (this.radiusTop - this.radiusBottom) * t;
+        return Math.hypot(point[0], point[2]) <= radius + EPSILON;
+    }
+    distanceToLocalPoint(point) {
+        if (this.containsLocalPoint(point))
+            return 0;
+        const radial = Math.hypot(point[0], point[2]);
+        const half = this.height / 2;
+        const side = closestDistanceToSegment2D(radial, point[1], this.radiusBottom, -half, this.radiusTop, half);
+        const bottom = closestDistanceToSegment2D(radial, point[1], 0, -half, this.radiusBottom, -half);
+        const top = closestDistanceToSegment2D(radial, point[1], 0, half, this.radiusTop, half);
+        return Math.min(side, bottom, top);
+    }
+    supportLocal(direction) {
+        const radialLength = Math.hypot(direction[0], direction[2]);
+        const x = radialLength ? direction[0] / radialLength : 1;
+        const z = radialLength ? direction[2] / radialLength : 0;
+        const half = this.height / 2;
+        const bottomValue = radialLength * this.radiusBottom - direction[1] * half;
+        const topValue = radialLength * this.radiusTop + direction[1] * half;
+        return topValue >= bottomValue
+            ? [x * this.radiusTop, half, z * this.radiusTop]
+            : [x * this.radiusBottom, -half, z * this.radiusBottom];
+    }
+    intersectLocalRay(origin, direction) {
+        const half = this.height / 2;
+        const slope = (this.radiusTop - this.radiusBottom) / this.height;
+        const intercept = (this.radiusTop + this.radiusBottom) / 2;
+        const radialAtOrigin = intercept + slope * origin[1];
+        const radialAlongRay = slope * direction[1];
+        const candidates = [];
+        const roots = solveQuadratic(direction[0] ** 2 + direction[2] ** 2 - radialAlongRay ** 2, 2 * (origin[0] * direction[0] + origin[2] * direction[2] - radialAtOrigin * radialAlongRay), origin[0] ** 2 + origin[2] ** 2 - radialAtOrigin ** 2);
+        for (const distance of roots) {
+            const y = origin[1] + distance * direction[1];
+            if (distance >= 0 && y >= -half - EPSILON && y <= half + EPSILON) {
+                const x = origin[0] + distance * direction[0];
+                const z = origin[2] + distance * direction[2];
+                const radial = Math.hypot(x, z) || 1;
+                const normal = [x / radial, -slope, z / radial];
+                const length = Math.hypot(...normal);
+                candidates.push({ distance, normal: normal.map(value => value / length) });
+            }
+        }
+        if (Math.abs(direction[1]) > 1e-15) {
+            for (const [y, radius, normalY] of [
+                [-half, this.radiusBottom, -1],
+                [half, this.radiusTop, 1]
+            ]) {
+                const distance = (y - origin[1]) / direction[1];
+                const x = origin[0] + distance * direction[0];
+                const z = origin[2] + distance * direction[2];
+                if (distance >= 0 && Math.hypot(x, z) <= radius + EPSILON)
+                    candidates.push({ distance, normal: [0, normalY, 0] });
+            }
+        }
+        return candidates.sort((a, b) => a.distance - b.distance)[0];
+    }
+}
+//# sourceMappingURL=cylinder-shape.js.map

package/dist/lib/shapes/cylinder-shape.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"cylinder-shape.js","sourceRoot":"","sources":["../../../src/lib/shapes/cylinder-shape.ts"],"names":[],"mappings":"AAAA,UAAU;AACV,+BAA+B;AAC/B,oCAAoC;AAEpC,OAAO,EACL,0BAA0B,EAC1B,iBAAiB,EACjB,cAAc,EAGf,4BAAyB;AAE1B,MAAM,OAAO,GAAG,IAAI,CAAC;AASrB,mFAAmF;AACnF,MAAM,OAAO,aAAc,SAAQ,iBAAiB;IAMlD,YAAY,QAA4B,EAAE;QACxC,KAAK,CAAC,KAAK,CAAC,MAAM,CAAC,CAAC;QANb,SAAI,GAAG,UAAmB,CAAC;QAOlC,IAAI,CAAC,MAAM,GAAG,KAAK,CAAC,MAAM,IAAI,CAAC,CAAC;QAChC,IAAI,CAAC,YAAY,GAAG,KAAK,CAAC,YAAY,IAAI,GAAG,CAAC;QAC9C,IAAI,CAAC,SAAS,GAAG,KAAK,CAAC,SAAS,IAAI,GAAG,CAAC;QACxC,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,IAAI,CAAC,MAAM,CAAC,IAAI,IAAI,CAAC,MAAM,IAAI,CAAC;YACnD,MAAM,IAAI,KAAK,CAAC,2CAA2C,CAAC,CAAC;QAC/D,IAAI,CAAC,CAAC,IAAI,CAAC,YAAY,EAAE,IAAI,CAAC,SAAS,CAAC,CAAC,KAAK,CAAC,KAAK,CAAC,EAAE,CAAC,MAAM,CAAC,QAAQ,CAAC,KAAK,CAAC,IAAI,KAAK,IAAI,CAAC,CAAC;YAC3F,MAAM,IAAI,KAAK,CAAC,qCAAqC,CAAC,CAAC;QACzD,IAAI,IAAI,CAAC,YAAY,KAAK,CAAC,IAAI,IAAI,CAAC,SAAS,KAAK,CAAC;YACjD,MAAM,IAAI,KAAK,CAAC,wDAAwD,CAAC,CAAC;IAC9E,CAAC;IAED,KAAK;QACH,OAAO,IAAI,aAAa,CAAC;YACvB,MAAM,EAAE,IAAI,CAAC,MAAM;YACnB,YAAY,EAAE,IAAI,CAAC,YAAY;YAC/B,SAAS,EAAE,IAAI,CAAC,SAAS;YACzB,MAAM,EAAE,IAAI,CAAC,MAAM;SACpB,CAAC,CAAC;IACL,CAAC;IACQ,MAAM,CAAC,KAAoB;QAClC,OAAO,CACL,KAAK,YAAY,aAAa;YAC9B,KAAK,CAAC,MAAM,CAAC,KAAK,CAAC;YACnB,IAAI,CAAC,MAAM,KAAK,KAAK,CAAC,MAAM;YAC5B,IAAI,CAAC,YAAY,KAAK,KAAK,CAAC,YAAY;YACxC,IAAI,CAAC,SAAS,KAAK,KAAK,CAAC,SAAS,CACnC,CAAC;IACJ,CAAC;IACD,kBAAkB,CAAC,KAAwB;QACzC,MAAM,IAAI,GAAG,IAAI,CAAC,MAAM,GAAG,CAAC,CAAC;QAC7B,IAAI,KAAK,CAAC,CAAC,CAAC,GAAG,CAAC,IAAI,GAAG,OAAO,IAAI,KAAK,CAAC,CAAC,CAAC,GAAG,IAAI,GAAG,OAAO;YAAE,OAAO,KAAK,CAAC;QAC1E,MAAM,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,IAAI,CAAC,MAAM,CAAC;QAC1C,MAAM,MAAM,GAAG,IAAI,CAAC,YAAY,GAAG,CAAC,IAAI,CAAC,SAAS,GAAG,IAAI,CAAC,YAAY,CAAC,GAAG,CAAC,CAAC;QAC5E,OAAO,IAAI,CAAC,KAAK,CAAC,KAAK,CAAC,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC,CAAC,IAAI,MAAM,GAAG,OAAO,CAAC;IAC5D,CAAC;IACD,oBAAoB,CAAC,KAAwB;QAC3C,IAAI,IAAI,CAAC,kBAAkB,CAAC,KAAK,CAAC;YAAE,OAAO,CAAC,CAAC;QAC7C,MAAM,MAAM,GAAG,IAAI,CAAC,KAAK,CAAC,KAAK,CAAC,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC;QAC9C,MAAM,IAAI,GAAG,IAAI,CAAC,MAAM,GAAG,CAAC,CAAC;QAC7B,MAAM,IAAI,GAAG,0BAA0B,CACrC,MAAM,EACN,KAAK,CAAC,CAAC,CAAC,EACR,IAAI,CAAC,YAAY,EACjB,CAAC,IAAI,EACL,IAAI,CAAC,SAAS,EACd,IAAI,CACL,CAAC;QACF,MAAM,MAAM,GAAG,0BAA0B,CAAC,MAAM,EAAE,KAAK,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,IAAI,EAAE,IAAI,CAAC,YAAY,EAAE,CAAC,IAAI,CAAC,CAAC;QAChG,MAAM,GAAG,GAAG,0BAA0B,CAAC,MAAM,EAAE,KAAK,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,IAAI,EAAE,IAAI,CAAC,SAAS,EAAE,IAAI,CAAC,CAAC;QACxF,OAAO,IAAI,CAAC,GAAG,CAAC,IAAI,EAAE,MAAM,EAAE,GAAG,CAAC,CAAC;IACrC,CAAC;IACD,YAAY,CAAC,SAA4B;QACvC,MAAM,YAAY,GAAG,IAAI,CAAC,KAAK,CAAC,SAAS,CAAC,CAAC,CAAC,EAAE,SAAS,CAAC,CAAC,CAAC,CAAC,CAAC;QAC5D,MAAM,CAAC,GAAG,YAAY,CAAC,CAAC,CAAC,SAAS,CAAC,CAAC,CAAC,GAAG,YAAY,CAAC,CAAC,CAAC,CAAC,CAAC;QACzD,MAAM,CAAC,GAAG,YAAY,CAAC,CAAC,CAAC,SAAS,CAAC,CAAC,CAAC,GAAG,YAAY,CAAC,CAAC,CAAC,CAAC,CAAC;QACzD,MAAM,IAAI,GAAG,IAAI,CAAC,MAAM,GAAG,CAAC,CAAC;QAC7B,MAAM,WAAW,GAAG,YAAY,GAAG,IAAI,CAAC,YAAY,GAAG,SAAS,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC;QAC3E,MAAM,QAAQ,GAAG,YAAY,GAAG,IAAI,CAAC,SAAS,GAAG,SAAS,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC;QACrE,OAAO,QAAQ,IAAI,WAAW;YAC5B,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,SAAS,EAAE,IAAI,EAAE,CAAC,GAAG,IAAI,CAAC,SAAS,CAAC;YAChD,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,YAAY,EAAE,CAAC,IAAI,EAAE,CAAC,GAAG,IAAI,CAAC,YAAY,CAAC,CAAC;IAC5D,CAAC;IACD,iBAAiB,CACf,MAAyB,EACzB,SAA4B;QAE5B,MAAM,IAAI,GAAG,IAAI,CAAC,MAAM,GAAG,CAAC,CAAC;QAC7B,MAAM,KAAK,GAAG,CAAC,IAAI,CAAC,SAAS,GAAG,IAAI,CAAC,YAAY,CAAC,GAAG,IAAI,CAAC,MAAM,CAAC;QACjE,MAAM,SAAS,GAAG,CAAC,IAAI,CAAC,SAAS,GAAG,IAAI,CAAC,YAAY,CAAC,GAAG,CAAC,CAAC;QAC3D,MAAM,cAAc,GAAG,SAAS,GAAG,KAAK,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC;QACrD,MAAM,cAAc,GAAG,KAAK,GAAG,SAAS,CAAC,CAAC,CAAC,CAAC;QAC5C,MAAM,UAAU,GAA2B,EAAE,CAAC;QAC9C,MAAM,KAAK,GAAG,cAAc,CAC1B,SAAS,CAAC,CAAC,CAAC,IAAI,CAAC,GAAG,SAAS,CAAC,CAAC,CAAC,IAAI,CAAC,GAAG,cAAc,IAAI,CAAC,EAC3D,CAAC,GAAG,CAAC,MAAM,CAAC,CAAC,CAAC,GAAG,SAAS,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,GAAG,SAAS,CAAC,CAAC,CAAC,GAAG,cAAc,GAAG,cAAc,CAAC,EAC3F,MAAM,CAAC,CAAC,CAAC,IAAI,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,IAAI,CAAC,GAAG,cAAc,IAAI,CAAC,CACtD,CAAC;QACF,KAAK,MAAM,QAAQ,IAAI,KAAK,EAAE,CAAC;YAC7B,MAAM,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,GAAG,QAAQ,GAAG,SAAS,CAAC,CAAC,CAAC,CAAC;YAC9C,IAAI,QAAQ,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,IAAI,GAAG,OAAO,IAAI,CAAC,IAAI,IAAI,GAAG,OAAO,EAAE,CAAC;gBACjE,MAAM,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,GAAG,QAAQ,GAAG,SAAS,CAAC,CAAC,CAAC,CAAC;gBAC9C,MAAM,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,GAAG,QAAQ,GAAG,SAAS,CAAC,CAAC,CAAC,CAAC;gBAC9C,MAAM,MAAM,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC,IAAI,CAAC,CAAC;gBACrC,MAAM,MAAM,GAAG,CAAC,CAAC,GAAG,MAAM,EAAE,CAAC,KAAK,EAAE,CAAC,GAAG,MAAM,CAAC,CAAC;gBAChD,MAAM,MAAM,GAAG,IAAI,CAAC,KAAK,CAAC,GAAG,MAAM,CAAC,CAAC;gBACrC,UAAU,CAAC,IAAI,CAAC,EAAC,QAAQ,EAAE,MAAM,EAAE,MAAM,CAAC,GAAG,CAAC,KAAK,CAAC,EAAE,CAAC,KAAK,GAAG,MAAM,CAAC,EAAC,CAAC,CAAC;YAC3E,CAAC;QACH,CAAC;QACD,IAAI,IAAI,CAAC,GAAG,CAAC,SAAS,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,EAAE,CAAC;YACnC,KAAK,MAAM,CAAC,CAAC,EAAE,MAAM,EAAE,OAAO,CAAC,IAAI;gBACjC,CAAC,CAAC,IAAI,EAAE,IAAI,CAAC,YAAY,EAAE,CAAC,CAAC,CAAC;gBAC9B,CAAC,IAAI,EAAE,IAAI,CAAC,SAAS,EAAE,CAAC,CAAC;aAC1B,EAAE,CAAC;gBACF,MAAM,QAAQ,GAAG,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC,GAAG,SAAS,CAAC,CAAC,CAAC,CAAC;gBAChD,MAAM,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,GAAG,QAAQ,GAAG,SAAS,CAAC,CAAC,CAAC,CAAC;gBAC9C,MAAM,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,GAAG,QAAQ,GAAG,SAAS,CAAC,CAAC,CAAC,CAAC;gBAC9C,IAAI,QAAQ,IAAI,CAAC,IAAI,IAAI,CAAC,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC,IAAI,MAAM,GAAG,OAAO;oBACvD,UAAU,CAAC,IAAI,CAAC,EAAC,QAAQ,EAAE,MAAM,EAAE,CAAC,CAAC,EAAE,OAAO,EAAE,CAAC,CAAC,EAAC,CAAC,CAAC;YACzD,CAAC;QACH,CAAC;QACD,OAAO,UAAU,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,QAAQ,GAAG,CAAC,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC;IAC/D,CAAC;CACF"}

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

@@ -0,0 +1,55 @@
+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";
+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): number;
+    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/lib/shapes/implicit-shape.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"implicit-shape.d.ts","sourceRoot":"","sources":["../../../src/lib/shapes/implicit-shape.ts"],"names":[],"mappings":"AAIA,OAAO,EAAC,OAAO,EAAE,OAAO,EAAC,MAAM,eAAe,CAAC;AAE/C,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,MAAM;IAUpC,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/lib/shapes/implicit-shape.js

@@ -0,0 +1,143 @@
+// math.gl
+// SPDX-License-Identifier: MIT
+// Copyright (c) vis.gl contributors
+import { Matrix4, Vector3 } from '@math.gl/core';
+import { INTERSECTION } from "../../constants.js";
+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 INTERSECTION.INSIDE;
+        if (maxDistance < 0)
+            return INTERSECTION.OUTSIDE;
+        return INTERSECTION.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/lib/shapes/implicit-shape.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"implicit-shape.js","sourceRoot":"","sources":["../../../src/lib/shapes/implicit-shape.ts"],"names":[],"mappings":"AAAA,UAAU;AACV,+BAA+B;AAC/B,oCAAoC;AAEpC,OAAO,EAAC,OAAO,EAAE,OAAO,EAAC,MAAM,eAAe,CAAC;AAC/C,OAAO,EAAC,YAAY,EAAC,2BAAwB;AAC7C,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,YAAY,CAAC,MAAM,CAAC;QAChD,IAAI,WAAW,GAAG,CAAC;YAAE,OAAO,YAAY,CAAC,OAAO,CAAC;QACjD,OAAO,YAAY,CAAC,YAAY,CAAC;IACnC,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/lib/shapes/plane-shape.d.ts

@@ -0,0 +1,26 @@
+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): number;
+    getAxisAlignedBoundingBox(): AxisAlignedBoundingBox | undefined;
+    getBoundingSphere(): BoundingSphere | undefined;
+}
+//# sourceMappingURL=plane-shape.d.ts.map

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

@@ -0,0 +1 @@
+{"version":3,"file":"plane-shape.d.ts","sourceRoot":"","sources":["../../../src/lib/shapes/plane-shape.ts"],"names":[],"mappings":"AAMA,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,MAAM;IAWpC,yBAAyB,IAAI,sBAAsB,GAAG,SAAS;IAG/D,iBAAiB,IAAI,cAAc,GAAG,SAAS;CAGzD"}

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

@@ -0,0 +1,72 @@
+// math.gl
+// SPDX-License-Identifier: MIT
+// Copyright (c) vis.gl contributors
+import { Vector3 } from '@math.gl/core';
+import { INTERSECTION } from "../../constants.js";
+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 INTERSECTION.INTERSECTING;
+        const boundaryDistance = plane.getPointDistance(boundaryPoint);
+        if (alignment < 0 && boundaryDistance > 0)
+            return INTERSECTION.INSIDE;
+        if (alignment > 0 && boundaryDistance < 0)
+            return INTERSECTION.OUTSIDE;
+        return INTERSECTION.INTERSECTING;
+    }
+    getAxisAlignedBoundingBox() {
+        return undefined;
+    }
+    getBoundingSphere() {
+        return undefined;
+    }
+}
+//# sourceMappingURL=plane-shape.js.map