@math.gl/geospatial 5.0.0-alpha.9 → 5.0.0-beta.2

package/dist/discontinuous/mercator-map.js

@@ -0,0 +1,246 @@
+// math.gl
+// SPDX-License-Identifier: MIT
+// SPDX-FileCopyrightText: Copyright (c) vis.gl contributors
+import polygonClipping from 'polygon-clipping';
+import { tracePolygonConstraints, composePolygonEdge } from "./polygon-constraints.js";
+import { sourceArc } from "./source-arc.js";
+import { DiscontinuousMap } from "./map-layout.js";
+import { angle, coordinates, geographic, RAD, unit } from "./math.js";
+/** Spherical Mercator, clipped to a finite latitude band after seam splitting. */
+export class MercatorMap extends DiscontinuousMap {
+    constructor(radius = 1, centralMeridian = 0, maxLatitude = 85.0511287798066) {
+        if (!Number.isFinite(radius) ||
+            radius <= 0 ||
+            !Number.isFinite(centralMeridian) ||
+            !Number.isFinite(maxLatitude) ||
+            maxLatitude <= 0 ||
+            maxLatitude >= 90)
+            throw new Error('Invalid Mercator radius, central meridian or latitude limit');
+        const regions = [];
+        for (const sign of [1, -1])
+            for (let i = 0; i < 4; i++) {
+                const west = -180 + i * 90, east = west + 90, center = west + 45;
+                regions.push({
+                    id: `${sign > 0 ? 'north' : 'south'}-${i}`,
+                    layout: { kind: 'main' },
+                    boundary: [
+                        [west + centralMeridian, sign * 90],
+                        [west + centralMeridian, 0],
+                        [east + centralMeridian, 0]
+                    ],
+                    interior: [center + centralMeridian, sign * 45],
+                    project: ([lon, lat]) => {
+                        const delta = ((((lon - centralMeridian - center + 180) % 360) + 360) % 360) - 180;
+                        // Keep longitude unchanged in the finite polar extension. Tapering
+                        // X toward zero moves latitude-band intersections into other charts.
+                        // Directional pole limits are supplied below; the extension is clipped.
+                        const outside = Math.max(0, Math.abs(lat) - maxLatitude);
+                        const phi = Math.min(Math.abs(lat), maxLatitude) * RAD;
+                        return [
+                            radius * (center + delta) * RAD,
+                            radius * Math.sign(lat) * (Math.asinh(Math.tan(phi)) + outside * RAD)
+                        ];
+                    },
+                    unproject: ([x, y]) => [
+                        x / radius / RAD + centralMeridian,
+                        Math.atan(Math.sinh(y / radius)) / RAD
+                    ]
+                });
+            }
+        super(regions);
+        this.radius = radius;
+        this.centralMeridian = centralMeridian;
+        this.limit = maxLatitude;
+        this.ymax = radius * Math.asinh(Math.tan(maxLatitude * RAD));
+    }
+    resolveSourceLimit(point, towards) {
+        if (Math.hypot(point[0], point[1]) < 1e-14 && Math.hypot(towards[0], towards[1]) > 1e-14)
+            return geographic([Math.atan2(towards[1], towards[0]) / RAD, point[2] > 0 ? 90 : -90]);
+        return point;
+    }
+    containsConvexPolygon() {
+        return false;
+    }
+    findRegions(p) {
+        geographic(p);
+        return Math.abs(p[1]) <= this.limit + 1e-10 ? super.findRegions(p) : [];
+    }
+    projectPoint(p, id) {
+        geographic(p);
+        return Math.abs(p[1]) <= this.limit + 1e-10 ? super.projectPoint(p, id) : null;
+    }
+    unprojectPoint(p, id) {
+        return super
+            .unprojectPoint(p, id)
+            .filter(hit => Math.abs(hit.position[1]) <= this.limit + 1e-10);
+    }
+    source([x, y]) {
+        return [
+            x / this.radius / RAD + this.centralMeridian,
+            Math.atan(Math.sinh(y / this.radius)) / RAD
+        ];
+    }
+    projectLine(line, options) {
+        const result = [];
+        for (const part of super.projectLine(line, options)) {
+            let current;
+            for (let i = 1; i < part.coordinates.length; i++) {
+                const a = part.coordinates[i - 1], b = part.coordinates[i];
+                const dy = b[1] - a[1];
+                const t0 = dy
+                    ? Math.max(0, Math.min((-this.ymax - a[1]) / dy, (this.ymax - a[1]) / dy))
+                    : 0;
+                const t1 = dy
+                    ? Math.min(1, Math.max((-this.ymax - a[1]) / dy, (this.ymax - a[1]) / dy))
+                    : 1;
+                if (t0 > t1 || (!dy && Math.abs(a[1]) > this.ymax)) {
+                    current = undefined;
+                    continue;
+                }
+                const at = (t) => [a[0] + (b[0] - a[0]) * t, a[1] + dy * t];
+                // Locate the latitude cut on the sampled source great-circle arc, not
+                // on its projected chord. Retain the parent segment throughout this stage.
+                const endpoint = (t) => {
+                    if (t === 0 || t === 1)
+                        return { xy: at(t), source: part.sourceCoordinates[i - 1 + t], t };
+                    const va = geographic(part.sourceCoordinates[i - 1]), vb = geographic(part.sourceCoordinates[i]), arc = angle(va, vb);
+                    const sourceAt = (t) => {
+                        if (!arc)
+                            return part.sourceCoordinates[i - 1];
+                        const a = Math.sin((1 - t) * arc) / Math.sin(arc), b = Math.sin(t * arc) / Math.sin(arc);
+                        return coordinates(unit(va.map((v, k) => a * v + b * vb[k])));
+                    };
+                    const latitude = Math.sign(at(t)[1]) * this.limit;
+                    let lo = 0, hi = 1;
+                    for (let j = 0; j < 52; j++) {
+                        const mid = (lo + hi) / 2;
+                        if (sourceAt(mid)[1] < latitude === part.sourceCoordinates[i - 1][1] < latitude)
+                            lo = mid;
+                        else
+                            hi = mid;
+                    }
+                    const f = (lo + hi) / 2, source = sourceAt(f);
+                    const xy = super.projectPoint(source, part.regionId).position;
+                    return { xy, source, t: f };
+                };
+                const start = endpoint(t0), end = endpoint(t1);
+                const append = (point, index) => {
+                    current.coordinates.push(point.xy);
+                    current.sourceCoordinates.push(point.source);
+                    const source = part.provenance, target = current.provenance;
+                    if (!source || !target)
+                        return;
+                    const segment = source.segments[i];
+                    const from = source.segments[i - 1] === segment ? source.parameters[i - 1] : 0;
+                    target.segments.push(point.t === 0 || point.t === 1 ? source.segments[index] : segment);
+                    target.parameters.push(point.t === 0 || point.t === 1
+                        ? source.parameters[index]
+                        : from + point.t * (source.parameters[i] - from));
+                    target.vertices.push(point.t === 0 || point.t === 1 ? source.vertices[index] : 0xffffffff);
+                    if (options?.pathDistances)
+                        target.angularDistances.push(source.angularDistances[i - 1] +
+                            point.t * (source.angularDistances[i] - source.angularDistances[i - 1]));
+                };
+                if (!current || t0 > 0) {
+                    current = {
+                        regionId: part.regionId,
+                        coordinates: [],
+                        sourceCoordinates: [],
+                        ...(part.provenance
+                            ? {
+                                provenance: {
+                                    segments: [],
+                                    parameters: [],
+                                    vertices: [],
+                                    angularDistances: [],
+                                    ...(options?.pathBoundaries
+                                        ? { startReason: t0 > 0 ? 'coverage-boundary' : part.provenance.startReason }
+                                        : {})
+                                }
+                            }
+                            : {})
+                    };
+                    result.push(current);
+                    append(start, i - 1);
+                }
+                append(end, i);
+                if (current.provenance && options?.pathBoundaries)
+                    current.provenance.endReason = t1 < 1 ? 'coverage-boundary' : part.provenance.endReason;
+                if (t1 < 1)
+                    current = undefined;
+            }
+        }
+        return result;
+    }
+    projectPolygon(rings, options) {
+        const x = Math.PI * this.radius, y = this.ymax;
+        const box = [
+            [
+                [-x, -y],
+                [x, -y],
+                [x, y],
+                [-x, y],
+                [-x, -y]
+            ]
+        ];
+        return super.projectPolygon(rings, options).flatMap(part => polygonClipping.intersection(part.coordinates, box).map(polygon => {
+            const traced = tracePolygonConstraints(part.coordinates, polygon, 1e-10, part.edgeProvenance);
+            const sources = traced.rings.map(r => r.map(p => this.source(p)));
+            traced.owners.forEach((owners, r) => owners.forEach((owner, i) => {
+                if (!owner)
+                    return;
+                const [parentRing, parentEdge] = owner, a = part.sourceCoordinates[parentRing][parentEdge], b = part.sourceCoordinates[parentRing][parentEdge + 1], at = sourceArc(a, b);
+                for (const j of [i, i + 1]) {
+                    const q = traced.rings[r][j], first = part.coordinates[parentRing][parentEdge], last = part.coordinates[parentRing][parentEdge + 1];
+                    if (Math.hypot(q[0] - first[0], q[1] - first[1]) < 1e-10 * this.radius)
+                        sources[r][j] = [...a];
+                    else if (Math.hypot(q[0] - last[0], q[1] - last[1]) < 1e-10 * this.radius)
+                        sources[r][j] = [...b];
+                    else {
+                        const latitude = Math.sign(q[1]) * this.limit;
+                        let lo = 0, hi = 1;
+                        if (a[1] < latitude === b[1] < latitude)
+                            throw new Error('Unresolved Mercator source constraint intersection');
+                        for (let k = 0; k < 52; k++) {
+                            const mid = (lo + hi) / 2;
+                            if (at(mid)[1] < latitude === a[1] < latitude)
+                                lo = mid;
+                            else
+                                hi = mid;
+                        }
+                        sources[r][j] = at((lo + hi) / 2);
+                    }
+                }
+            }));
+            for (const r of sources)
+                r[r.length - 1] = [...r[0]];
+            return {
+                regionId: part.regionId,
+                sourceCoordinates: sources,
+                coordinates: sources.map(r => r.map(p => super.projectPoint(p, part.regionId).position)),
+                ...(options?.polygonEdgeProvenance
+                    ? {
+                        edgeProvenance: traced.owners.map((owners, r) => owners.map((owner, i) => {
+                            if (owner) {
+                                const [pr, pe] = owner;
+                                return composePolygonEdge(part.edgeProvenance[pr][pe], part.sourceCoordinates[pr][pe], part.sourceCoordinates[pr][pe + 1], sources[r][i], sources[r][i + 1]);
+                            }
+                            if (!sources[r]
+                                .slice(i, i + 2)
+                                .every(p => Math.abs(Math.abs(p[1]) - this.limit) < 1e-7))
+                                throw new Error('Unknown Mercator clipping constraint');
+                            return {
+                                kind: 'coverage-boundary',
+                                ring: 0xffffffff,
+                                segment: 0xffffffff,
+                                t0: 0,
+                                t1: 0
+                            };
+                        }))
+                    }
+                    : {})
+            };
+        }));
+    }
+}
+//# sourceMappingURL=mercator-map.js.map

package/dist/discontinuous/mercator-map.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"mercator-map.js","sourceRoot":"","sources":["../../src/discontinuous/mercator-map.ts"],"names":[],"mappings":"AAAA,UAAU;AACV,+BAA+B;AAC/B,4DAA4D;AAE5D,OAAO,eAAe,MAAM,kBAAkB,CAAC;AAC/C,OAAO,EAAC,uBAAuB,EAAE,kBAAkB,EAAC,iCAA8B;AAClF,OAAO,EAAC,SAAS,EAAC,wBAAqB;AACvC,OAAO,EACL,gBAAgB,EAKjB,wBAAqB;AACtB,OAAO,EAAC,KAAK,EAAE,WAAW,EAAE,UAAU,EAAE,GAAG,EAAE,IAAI,EAAc,kBAAe;AAE9E,kFAAkF;AAClF,MAAM,OAAO,WAAY,SAAQ,gBAAgB;IAG/C,YACW,SAAS,CAAC,EACV,kBAAkB,CAAC,EAC5B,WAAW,GAAG,gBAAgB;QAE9B,IACE,CAAC,MAAM,CAAC,QAAQ,CAAC,MAAM,CAAC;YACxB,MAAM,IAAI,CAAC;YACX,CAAC,MAAM,CAAC,QAAQ,CAAC,eAAe,CAAC;YACjC,CAAC,MAAM,CAAC,QAAQ,CAAC,WAAW,CAAC;YAC7B,WAAW,IAAI,CAAC;YAChB,WAAW,IAAI,EAAE;YAEjB,MAAM,IAAI,KAAK,CAAC,6DAA6D,CAAC,CAAC;QACjF,MAAM,OAAO,GAAgB,EAAE,CAAC;QAChC,KAAK,MAAM,IAAI,IAAI,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;YACxB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;gBAC3B,MAAM,IAAI,GAAG,CAAC,GAAG,GAAG,CAAC,GAAG,EAAE,EACxB,IAAI,GAAG,IAAI,GAAG,EAAE,EAChB,MAAM,GAAG,IAAI,GAAG,EAAE,CAAC;gBACrB,OAAO,CAAC,IAAI,CAAC;oBACX,EAAE,EAAE,GAAG,IAAI,GAAG,CAAC,CAAC,CAAC,CAAC,OAAO,CAAC,CAAC,CAAC,OAAO,IAAI,CAAC,EAAE;oBAC1C,MAAM,EAAE,EAAC,IAAI,EAAE,MAAM,EAAC;oBACtB,QAAQ,EAAE;wBACR,CAAC,IAAI,GAAG,eAAe,EAAE,IAAI,GAAG,EAAE,CAAC;wBACnC,CAAC,IAAI,GAAG,eAAe,EAAE,CAAC,CAAC;wBAC3B,CAAC,IAAI,GAAG,eAAe,EAAE,CAAC,CAAC;qBAC5B;oBACD,QAAQ,EAAE,CAAC,MAAM,GAAG,eAAe,EAAE,IAAI,GAAG,EAAE,CAAC;oBAC/C,OAAO,EAAE,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC,EAAE,EAAE;wBACtB,MAAM,KAAK,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,GAAG,eAAe,GAAG,MAAM,GAAG,GAAG,CAAC,GAAG,GAAG,CAAC,GAAG,GAAG,CAAC,GAAG,GAAG,CAAC,GAAG,GAAG,CAAC;wBACnF,mEAAmE;wBACnE,qEAAqE;wBACrE,wEAAwE;wBACxE,MAAM,OAAO,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,WAAW,CAAC,CAAC;wBAEzD,MAAM,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,WAAW,CAAC,GAAG,GAAG,CAAC;wBACvD,OAAO;4BACL,MAAM,GAAG,CAAC,MAAM,GAAG,KAAK,CAAC,GAAG,GAAG;4BAC/B,MAAM,GAAG,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,OAAO,GAAG,GAAG,CAAC;yBACtE,CAAC;oBACJ,CAAC;oBACD,SAAS,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,EAAE,CAAC;wBACrB,CAAC,GAAG,MAAM,GAAG,GAAG,GAAG,eAAe;wBAClC,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,GAAG,GAAG;qBACvC;iBACF,CAAC,CAAC;YACL,CAAC;QACH,KAAK,CAAC,OAAO,CAAC,CAAC;QA/CN,WAAM,GAAN,MAAM,CAAI;QACV,oBAAe,GAAf,eAAe,CAAI;QA+C5B,IAAI,CAAC,KAAK,GAAG,WAAW,CAAC;QACzB,IAAI,CAAC,IAAI,GAAG,MAAM,GAAG,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,GAAG,CAAC,WAAW,GAAG,GAAG,CAAC,CAAC,CAAC;IAC/D,CAAC;IACkB,kBAAkB,CAAC,KAAa,EAAE,OAAe;QAClE,IAAI,IAAI,CAAC,KAAK,CAAC,KAAK,CAAC,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,IAAI,IAAI,CAAC,KAAK,CAAC,OAAO,CAAC,CAAC,CAAC,EAAE,OAAO,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK;YACtF,OAAO,UAAU,CAAC,CAAC,IAAI,CAAC,KAAK,CAAC,OAAO,CAAC,CAAC,CAAC,EAAE,OAAO,CAAC,CAAC,CAAC,CAAC,GAAG,GAAG,EAAE,KAAK,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;QACzF,OAAO,KAAK,CAAC;IACf,CAAC;IACQ,qBAAqB;QAC5B,OAAO,KAAK,CAAC;IACf,CAAC;IACQ,WAAW,CAAC,CAAqB;QACxC,UAAU,CAAC,CAAC,CAAC,CAAC;QACd,OAAO,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI,IAAI,CAAC,KAAK,GAAG,KAAK,CAAC,CAAC,CAAC,KAAK,CAAC,WAAW,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC;IAC1E,CAAC;IACQ,YAAY,CAAC,CAAqB,EAAE,EAAW;QACtD,UAAU,CAAC,CAAC,CAAC,CAAC;QACd,OAAO,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI,IAAI,CAAC,KAAK,GAAG,KAAK,CAAC,CAAC,CAAC,KAAK,CAAC,YAAY,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC;IACjF,CAAC;IACQ,cAAc,CAAC,CAAqB,EAAE,EAAW;QACxD,OAAO,KAAK;aACT,cAAc,CAAC,CAAC,EAAE,EAAE,CAAC;aACrB,MAAM,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,IAAI,IAAI,CAAC,KAAK,GAAG,KAAK,CAAC,CAAC;IACpE,CAAC;IACO,MAAM,CAAC,CAAC,CAAC,EAAE,CAAC,CAAqB;QACvC,OAAO;YACL,CAAC,GAAG,IAAI,CAAC,MAAM,GAAG,GAAG,GAAG,IAAI,CAAC,eAAe;YAC5C,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,GAAG,IAAI,CAAC,MAAM,CAAC,CAAC,GAAG,GAAG;SAC5C,CAAC;IACJ,CAAC;IACQ,WAAW,CAClB,IAAmC,EACnC,OAA4B;QAE5B,MAAM,MAAM,GAAc,EAAE,CAAC;QAC7B,KAAK,MAAM,IAAI,IAAI,KAAK,CAAC,WAAW,CAAC,IAAI,EAAE,OAAO,CAAC,EAAE,CAAC;YACpD,IAAI,OAA4B,CAAC;YACjC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,WAAW,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;gBACjD,MAAM,CAAC,GAAG,IAAI,CAAC,WAAW,CAAC,CAAC,GAAG,CAAC,CAAC,EAC/B,CAAC,GAAG,IAAI,CAAC,WAAW,CAAC,CAAC,CAAC,CAAC;gBAC1B,MAAM,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;gBACvB,MAAM,EAAE,GAAG,EAAE;oBACX,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,IAAI,CAAC,IAAI,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,EAAE,EAAE,CAAC,IAAI,CAAC,IAAI,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC;oBAC1E,CAAC,CAAC,CAAC,CAAC;gBACN,MAAM,EAAE,GAAG,EAAE;oBACX,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,IAAI,CAAC,IAAI,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,EAAE,EAAE,CAAC,IAAI,CAAC,IAAI,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC;oBAC1E,CAAC,CAAC,CAAC,CAAC;gBACN,IAAI,EAAE,GAAG,EAAE,IAAI,CAAC,CAAC,EAAE,IAAI,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,EAAE,CAAC;oBACnD,OAAO,GAAG,SAAS,CAAC;oBACpB,SAAS;gBACX,CAAC;gBACD,MAAM,EAAE,GAAG,CAAC,CAAS,EAAoB,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC,CAAC;gBACtF,sEAAsE;gBACtE,2EAA2E;gBAC3E,MAAM,QAAQ,GAAG,CAAC,CAAS,EAAE,EAAE;oBAC7B,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,KAAK,CAAC;wBAAE,OAAO,EAAC,EAAE,EAAE,EAAE,CAAC,CAAC,CAAC,EAAE,MAAM,EAAE,IAAI,CAAC,iBAAiB,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,EAAC,CAAC;oBACzF,MAAM,EAAE,GAAG,UAAU,CAAC,IAAI,CAAC,iBAAiB,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAClD,EAAE,GAAG,UAAU,CAAC,IAAI,CAAC,iBAAiB,CAAC,CAAC,CAAC,CAAC,EAC1C,GAAG,GAAG,KAAK,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC;oBACtB,MAAM,QAAQ,GAAG,CAAC,CAAS,EAAE,EAAE;wBAC7B,IAAI,CAAC,GAAG;4BAAE,OAAO,IAAI,CAAC,iBAAiB,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;wBAC/C,MAAM,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,EAC/C,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;wBACxC,OAAO,WAAW,CAAC,IAAI,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,CAAW,CAAC,CAAC,CAAC;oBAC1E,CAAC,CAAC;oBACF,MAAM,QAAQ,GAAG,IAAI,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC;oBAClD,IAAI,EAAE,GAAG,CAAC,EACR,EAAE,GAAG,CAAC,CAAC;oBACT,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,EAAE,EAAE,CAAC,EAAE,EAAE,CAAC;wBAC5B,MAAM,GAAG,GAAG,CAAC,EAAE,GAAG,EAAE,CAAC,GAAG,CAAC,CAAC;wBAC1B,IAAI,QAAQ,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,QAAQ,KAAK,IAAI,CAAC,iBAAiB,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,QAAQ;4BAC7E,EAAE,GAAG,GAAG,CAAC;;4BACN,EAAE,GAAG,GAAG,CAAC;oBAChB,CAAC;oBACD,MAAM,CAAC,GAAG,CAAC,EAAE,GAAG,EAAE,CAAC,GAAG,CAAC,EACrB,MAAM,GAAG,QAAQ,CAAC,CAAC,CAAC,CAAC;oBACvB,MAAM,EAAE,GAAG,KAAK,CAAC,YAAY,CAAC,MAAM,EAAE,IAAI,CAAC,QAAQ,CAAE,CAAC,QAAQ,CAAC;oBAC/D,OAAO,EAAC,EAAE,EAAE,MAAM,EAAE,CAAC,EAAE,CAAC,EAAC,CAAC;gBAC5B,CAAC,CAAC;gBACF,MAAM,KAAK,GAAG,QAAQ,CAAC,EAAE,CAAC,EACxB,GAAG,GAAG,QAAQ,CAAC,EAAE,CAAC,CAAC;gBACrB,MAAM,MAAM,GAAG,CAAC,KAAmB,EAAE,KAAa,EAAE,EAAE;oBACpD,OAAQ,CAAC,WAAW,CAAC,IAAI,CAAC,KAAK,CAAC,EAAE,CAAC,CAAC;oBACpC,OAAQ,CAAC,iBAAiB,CAAC,IAAI,CAAC,KAAK,CAAC,MAAM,CAAC,CAAC;oBAC9C,MAAM,MAAM,GAAG,IAAI,CAAC,UAAU,EAC5B,MAAM,GAAG,OAAQ,CAAC,UAAU,CAAC;oBAC/B,IAAI,CAAC,MAAM,IAAI,CAAC,MAAM;wBAAE,OAAO;oBAC/B,MAAM,OAAO,GAAG,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC;oBACnC,MAAM,IAAI,GAAG,MAAM,CAAC,QAAQ,CAAC,CAAC,GAAG,CAAC,CAAC,KAAK,OAAO,CAAC,CAAC,CAAC,MAAM,CAAC,UAAU,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;oBAC/E,MAAM,CAAC,QAAQ,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC,KAAK,CAAC,IAAI,KAAK,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,MAAM,CAAC,QAAQ,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,OAAO,CAAC,CAAC;oBACxF,MAAM,CAAC,UAAU,CAAC,IAAI,CACpB,KAAK,CAAC,CAAC,KAAK,CAAC,IAAI,KAAK,CAAC,CAAC,KAAK,CAAC;wBAC5B,CAAC,CAAC,MAAM,CAAC,UAAU,CAAC,KAAK,CAAC;wBAC1B,CAAC,CAAC,IAAI,GAAG,KAAK,CAAC,CAAC,GAAG,CAAC,MAAM,CAAC,UAAU,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,CACnD,CAAC;oBACF,MAAM,CAAC,QAAQ,CAAC,IAAI,CAClB,KAAK,CAAC,CAAC,KAAK,CAAC,IAAI,KAAK,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,MAAM,CAAC,QAAQ,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,UAAU,CACrE,CAAC;oBACF,IAAI,OAAO,EAAE,aAAa;wBACxB,MAAM,CAAC,gBAAgB,CAAC,IAAI,CAC1B,MAAM,CAAC,gBAAgB,CAAC,CAAC,GAAG,CAAC,CAAC;4BAC5B,KAAK,CAAC,CAAC,GAAG,CAAC,MAAM,CAAC,gBAAgB,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC,gBAAgB,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAC1E,CAAC;gBACN,CAAC,CAAC;gBACF,IAAI,CAAC,OAAO,IAAI,EAAE,GAAG,CAAC,EAAE,CAAC;oBACvB,OAAO,GAAG;wBACR,QAAQ,EAAE,IAAI,CAAC,QAAQ;wBACvB,WAAW,EAAE,EAAE;wBACf,iBAAiB,EAAE,EAAE;wBACrB,GAAG,CAAC,IAAI,CAAC,UAAU;4BACjB,CAAC,CAAC;gCACE,UAAU,EAAE;oCACV,QAAQ,EAAE,EAAE;oCACZ,UAAU,EAAE,EAAE;oCACd,QAAQ,EAAE,EAAE;oCACZ,gBAAgB,EAAE,EAAE;oCACpB,GAAG,CAAC,OAAO,EAAE,cAAc;wCACzB,CAAC,CAAC,EAAC,WAAW,EAAE,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,mBAAmB,CAAC,CAAC,CAAC,IAAI,CAAC,UAAU,CAAC,WAAW,EAAC;wCAC3E,CAAC,CAAC,EAAE,CAAC;iCACR;6BACF;4BACH,CAAC,CAAC,EAAE,CAAC;qBACR,CAAC;oBACF,MAAM,CAAC,IAAI,CAAC,OAAO,CAAC,CAAC;oBACrB,MAAM,CAAC,KAAK,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC;gBACvB,CAAC;gBACD,MAAM,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC;gBACf,IAAI,OAAO,CAAC,UAAU,IAAI,OAAO,EAAE,cAAc;oBAC/C,OAAO,CAAC,UAAU,CAAC,SAAS,GAAG,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,mBAAmB,CAAC,CAAC,CAAC,IAAI,CAAC,UAAW,CAAC,SAAS,CAAC;gBAC3F,IAAI,EAAE,GAAG,CAAC;oBAAE,OAAO,GAAG,SAAS,CAAC;YAClC,CAAC;QACH,CAAC;QACD,OAAO,MAAM,CAAC;IAChB,CAAC;IACQ,cAAc,CACrB,KAAiD,EACjD,OAA4B;QAE5B,MAAM,CAAC,GAAG,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,MAAM,EAC7B,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC;QAChB,MAAM,GAAG,GAAyB;YAChC;gBACE,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;gBACR,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;gBACP,CAAC,CAAC,EAAE,CAAC,CAAC;gBACN,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC;gBACP,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;aACT;SACF,CAAC;QACF,OAAO,KAAK,CAAC,cAAc,CAAC,KAAK,EAAE,OAAO,CAAC,CAAC,OAAO,CAAC,IAAI,CAAC,EAAE,CACzD,eAAe,CAAC,YAAY,CAAC,IAAI,CAAC,WAAW,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,OAAO,CAAC,EAAE;YAChE,MAAM,MAAM,GAAG,uBAAuB,CACpC,IAAI,CAAC,WAAW,EAChB,OAAO,EACP,KAAK,EACL,IAAI,CAAC,cAAc,CACpB,CAAC;YACF,MAAM,OAAO,GAAG,MAAM,CAAC,KAAK,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,IAAI,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;YAClE,MAAM,CAAC,MAAM,CAAC,OAAO,CAAC,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE,CAClC,MAAM,CAAC,OAAO,CAAC,CAAC,KAAK,EAAE,CAAC,EAAE,EAAE;gBAC1B,IAAI,CAAC,KAAK;oBAAE,OAAO;gBACnB,MAAM,CAAC,UAAU,EAAE,UAAU,CAAC,GAAG,KAAK,EACpC,CAAC,GAAG,IAAI,CAAC,iBAAiB,CAAC,UAAU,CAAC,CAAC,UAAU,CAAC,EAClD,CAAC,GAAG,IAAI,CAAC,iBAAiB,CAAC,UAAU,CAAC,CAAC,UAAU,GAAG,CAAC,CAAC,EACtD,EAAE,GAAG,SAAS,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;gBACvB,KAAK,MAAM,CAAC,IAAI,CAAC,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC;oBAC3B,MAAM,CAAC,GAAG,MAAM,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAC1B,KAAK,GAAG,IAAI,CAAC,WAAW,CAAC,UAAU,CAAC,CAAC,UAAU,CAAC,EAChD,IAAI,GAAG,IAAI,CAAC,WAAW,CAAC,UAAU,CAAC,CAAC,UAAU,GAAG,CAAC,CAAC,CAAC;oBACtD,IAAI,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,GAAG,IAAI,CAAC,MAAM;wBACpE,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC;yBACpB,IAAI,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,GAAG,IAAI,CAAC,MAAM;wBACvE,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC;yBACpB,CAAC;wBACJ,MAAM,QAAQ,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC;wBAC9C,IAAI,EAAE,GAAG,CAAC,EACR,EAAE,GAAG,CAAC,CAAC;wBACT,IAAI,CAAC,CAAC,CAAC,CAAC,GAAG,QAAQ,KAAK,CAAC,CAAC,CAAC,CAAC,GAAG,QAAQ;4BACrC,MAAM,IAAI,KAAK,CAAC,oDAAoD,CAAC,CAAC;wBACxE,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,EAAE,EAAE,CAAC,EAAE,EAAE,CAAC;4BAC5B,MAAM,GAAG,GAAG,CAAC,EAAE,GAAG,EAAE,CAAC,GAAG,CAAC,CAAC;4BAC1B,IAAI,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,QAAQ,KAAK,CAAC,CAAC,CAAC,CAAC,GAAG,QAAQ;gCAAE,EAAE,GAAG,GAAG,CAAC;;gCACnD,EAAE,GAAG,GAAG,CAAC;wBAChB,CAAC;wBACD,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,EAAE,GAAG,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC;oBACpC,CAAC;gBACH,CAAC;YACH,CAAC,CAAC,CACH,CAAC;YACF,KAAK,MAAM,CAAC,IAAI,OAAO;gBAAE,CAAC,CAAC,CAAC,CAAC,MAAM,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;YACrD,OAAO;gBACL,QAAQ,EAAE,IAAI,CAAC,QAAQ;gBACvB,iBAAiB,EAAE,OAAO;gBAC1B,WAAW,EAAE,OAAO,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,KAAK,CAAC,YAAY,CAAC,CAAC,EAAE,IAAI,CAAC,QAAQ,CAAE,CAAC,QAAQ,CAAC,CAAC;gBACzF,GAAG,CAAC,OAAO,EAAE,qBAAqB;oBAChC,CAAC,CAAC;wBACE,cAAc,EAAE,MAAM,CAAC,MAAM,CAAC,GAAG,CAAC,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE,CAC9C,MAAM,CAAC,GAAG,CAAC,CAAC,KAAK,EAAE,CAAC,EAAE,EAAE;4BACtB,IAAI,KAAK,EAAE,CAAC;gCACV,MAAM,CAAC,EAAE,EAAE,EAAE,CAAC,GAAG,KAAK,CAAC;gCACvB,OAAO,kBAAkB,CACvB,IAAI,CAAC,cAAe,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,EAC5B,IAAI,CAAC,iBAAiB,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,EAC9B,IAAI,CAAC,iBAAiB,CAAC,EAAE,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,EAClC,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EACb,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAClB,CAAC;4BACJ,CAAC;4BACD,IACE,CAAC,OAAO,CAAC,CAAC,CAAC;iCACR,KAAK,CAAC,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC;iCACf,KAAK,CAAC,CAAC,CAAC,EAAE,CAAC,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,GAAG,IAAI,CAAC;gCAE3D,MAAM,IAAI,KAAK,CAAC,sCAAsC,CAAC,CAAC;4BAC1D,OAAO;gCACL,IAAI,EAAE,mBAA4B;gCAClC,IAAI,EAAE,UAAU;gCAChB,OAAO,EAAE,UAAU;gCACnB,EAAE,EAAE,CAAC;gCACL,EAAE,EAAE,CAAC;6BACN,CAAC;wBACJ,CAAC,CAAC,CACH;qBACF;oBACH,CAAC,CAAC,EAAE,CAAC;aACR,CAAC;QACJ,CAAC,CAAC,CACH,CAAC;IACJ,CAAC;CACF"}

package/dist/discontinuous/myriahedral/graticule.d.ts

@@ -0,0 +1,6 @@
+import type { Polyhedron } from "./mesh.js";
+/** Independent, closed graticule mesh. Equal-latitude trapezoids are planar.
+ * Pole cells are triangles; shared longitude seams use the same vertex IDs.
+ */
+export declare function createGraticule(columns: number): Polyhedron;
+//# sourceMappingURL=graticule.d.ts.map

package/dist/discontinuous/myriahedral/graticule.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"graticule.d.ts","sourceRoot":"","sources":["../../../src/discontinuous/myriahedral/graticule.ts"],"names":[],"mappings":"AAIA,OAAO,KAAK,EAAC,UAAU,EAAO,kBAAe;AAE7C;;GAEG;AACH,wBAAgB,eAAe,CAAC,OAAO,EAAE,MAAM,GAAG,UAAU,CA8B3D"}

package/dist/discontinuous/myriahedral/graticule.js

@@ -0,0 +1,37 @@
+// math.gl
+// SPDX-License-Identifier: MIT
+// SPDX-FileCopyrightText: Copyright (c) vis.gl contributors
+/** Independent, closed graticule mesh. Equal-latitude trapezoids are planar.
+ * Pole cells are triangles; shared longitude seams use the same vertex IDs.
+ */
+export function createGraticule(columns) {
+    if (!Number.isInteger(columns) || columns < 12 || columns > 96 || columns % 4)
+        throw new RangeError('Use a multiple of four between 12 and 96 columns');
+    const rows = columns / 2;
+    const vertices = [[0, 0, -1]];
+    for (let row = 1; row < rows; row++) {
+        const latitude = -Math.PI / 2 + (row * Math.PI) / rows;
+        for (let col = 0; col < columns; col++) {
+            const longitude = -Math.PI + (col * 2 * Math.PI) / columns;
+            vertices.push([
+                Math.cos(latitude) * Math.cos(longitude),
+                Math.cos(latitude) * Math.sin(longitude),
+                Math.sin(latitude)
+            ]);
+        }
+    }
+    const north = vertices.push([0, 0, 1]) - 1;
+    const at = (row, col) => row === 0 ? 0 : row === rows ? north : 1 + (row - 1) * columns + (col % columns);
+    const faces = [];
+    for (let row = 0; row < rows; row++)
+        for (let col = 0; col < columns; col++) {
+            const cell = [at(row, col), at(row, col + 1), at(row + 1, col + 1), at(row + 1, col)];
+            faces.push([...new Set(cell)]);
+        }
+    // Root just north of the equator in Africa. Its southern edge points east,
+    // establishing an east-right, north-up flat frame without a camera spin.
+    const root = (rows / 2) * columns + Math.floor((columns * 200) / 360);
+    [faces[0], faces[root]] = [faces[root], faces[0]];
+    return { vertices, faces };
+}
+//# sourceMappingURL=graticule.js.map

package/dist/discontinuous/myriahedral/graticule.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"graticule.js","sourceRoot":"","sources":["../../../src/discontinuous/myriahedral/graticule.ts"],"names":[],"mappings":"AAAA,UAAU;AACV,+BAA+B;AAC/B,4DAA4D;AAI5D;;GAEG;AACH,MAAM,UAAU,eAAe,CAAC,OAAe;IAC7C,IAAI,CAAC,MAAM,CAAC,SAAS,CAAC,OAAO,CAAC,IAAI,OAAO,GAAG,EAAE,IAAI,OAAO,GAAG,EAAE,IAAI,OAAO,GAAG,CAAC;QAC3E,MAAM,IAAI,UAAU,CAAC,kDAAkD,CAAC,CAAC;IAC3E,MAAM,IAAI,GAAG,OAAO,GAAG,CAAC,CAAC;IACzB,MAAM,QAAQ,GAAW,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;IACtC,KAAK,IAAI,GAAG,GAAG,CAAC,EAAE,GAAG,GAAG,IAAI,EAAE,GAAG,EAAE,EAAE,CAAC;QACpC,MAAM,QAAQ,GAAG,CAAC,IAAI,CAAC,EAAE,GAAG,CAAC,GAAG,CAAC,GAAG,GAAG,IAAI,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC;QACvD,KAAK,IAAI,GAAG,GAAG,CAAC,EAAE,GAAG,GAAG,OAAO,EAAE,GAAG,EAAE,EAAE,CAAC;YACvC,MAAM,SAAS,GAAG,CAAC,IAAI,CAAC,EAAE,GAAG,CAAC,GAAG,GAAG,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC,GAAG,OAAO,CAAC;YAC3D,QAAQ,CAAC,IAAI,CAAC;gBACZ,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,SAAS,CAAC;gBACxC,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,SAAS,CAAC;gBACxC,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC;aACnB,CAAC,CAAC;QACL,CAAC;IACH,CAAC;IACD,MAAM,KAAK,GAAG,QAAQ,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;IAC3C,MAAM,EAAE,GAAG,CAAC,GAAW,EAAE,GAAW,EAAE,EAAE,CACtC,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,IAAI,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC,CAAC,GAAG,OAAO,GAAG,CAAC,GAAG,GAAG,OAAO,CAAC,CAAC;IACnF,MAAM,KAAK,GAAe,EAAE,CAAC;IAC7B,KAAK,IAAI,GAAG,GAAG,CAAC,EAAE,GAAG,GAAG,IAAI,EAAE,GAAG,EAAE;QACjC,KAAK,IAAI,GAAG,GAAG,CAAC,EAAE,GAAG,GAAG,OAAO,EAAE,GAAG,EAAE,EAAE,CAAC;YACvC,MAAM,IAAI,GAAG,CAAC,EAAE,CAAC,GAAG,EAAE,GAAG,CAAC,EAAE,EAAE,CAAC,GAAG,EAAE,GAAG,GAAG,CAAC,CAAC,EAAE,EAAE,CAAC,GAAG,GAAG,CAAC,EAAE,GAAG,GAAG,CAAC,CAAC,EAAE,EAAE,CAAC,GAAG,GAAG,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC;YACtF,KAAK,CAAC,IAAI,CAAC,CAAC,GAAG,IAAI,GAAG,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;QACjC,CAAC;IACH,2EAA2E;IAC3E,yEAAyE;IACzE,MAAM,IAAI,GAAG,CAAC,IAAI,GAAG,CAAC,CAAC,GAAG,OAAO,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,OAAO,GAAG,GAAG,CAAC,GAAG,GAAG,CAAC,CAAC;IACtE,CAAC,KAAK,CAAC,CAAC,CAAC,EAAE,KAAK,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,KAAK,CAAC,IAAI,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC;IAClD,OAAO,EAAC,QAAQ,EAAE,KAAK,EAAC,CAAC;AAC3B,CAAC"}

package/dist/discontinuous/myriahedral/mesh.d.ts

@@ -0,0 +1,18 @@
+export type Vec3 = [number, number, number];
+export type Polyhedron = {
+    vertices: Vec3[];
+    faces: number[][];
+};
+export type Solid = 'tetrahedron' | 'cube' | 'octahedron' | 'dodecahedron' | 'icosahedron';
+export declare const subtract: (a: Vec3, b: Vec3) => Vec3;
+export declare const dot: (a: Vec3, b: Vec3) => number;
+export declare const cross: (a: Vec3, b: Vec3) => Vec3;
+export declare const normalize: (v: Vec3) => Vec3;
+export declare const center: (points: Vec3[]) => Vec3;
+/** Independent construction from symmetric coordinates, with a small convex hull
+ * used only for the 4–20 vertices of the regular solids. No imported mesh tables.
+ */
+export declare function createSolid(solid: Solid, subdivisions?: number): Polyhedron;
+export declare function edgeKey(a: number, b: number): string;
+export declare function faceNormal(mesh: Polyhedron, face: number[]): Vec3;
+//# sourceMappingURL=mesh.d.ts.map

package/dist/discontinuous/myriahedral/mesh.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"mesh.d.ts","sourceRoot":"","sources":["../../../src/discontinuous/myriahedral/mesh.ts"],"names":[],"mappings":"AAIA,MAAM,MAAM,IAAI,GAAG,CAAC,MAAM,EAAE,MAAM,EAAE,MAAM,CAAC,CAAC;AAC5C,MAAM,MAAM,UAAU,GAAG;IAAC,QAAQ,EAAE,IAAI,EAAE,CAAC;IAAC,KAAK,EAAE,MAAM,EAAE,EAAE,CAAA;CAAC,CAAC;AAC/D,MAAM,MAAM,KAAK,GAAG,aAAa,GAAG,MAAM,GAAG,YAAY,GAAG,cAAc,GAAG,aAAa,CAAC;AAC3F,eAAO,MAAM,QAAQ,GAAI,GAAG,IAAI,EAAE,GAAG,IAAI,KAAG,IAA+C,CAAC;AAC5F,eAAO,MAAM,GAAG,GAAI,GAAG,IAAI,EAAE,GAAG,IAAI,KAAG,MAAiD,CAAC;AACzF,eAAO,MAAM,KAAK,GAAI,GAAG,IAAI,EAAE,GAAG,IAAI,KAAG,IAIxC,CAAC;AACF,eAAO,MAAM,SAAS,GAAI,GAAG,IAAI,KAAG,IAGnC,CAAC;AACF,eAAO,MAAM,MAAM,GAAI,QAAQ,IAAI,EAAE,KAAG,IAC0C,CAAC;AAEnF;;GAEG;AACH,wBAAgB,WAAW,CAAC,KAAK,EAAE,KAAK,EAAE,YAAY,SAAI,GAAG,UAAU,CAyEtE;AACD,wBAAgB,OAAO,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,GAAG,MAAM,CAEpD;AACD,wBAAgB,UAAU,CAAC,IAAI,EAAE,UAAU,EAAE,IAAI,EAAE,MAAM,EAAE,GAAG,IAAI,CAOjE"}

package/dist/discontinuous/myriahedral/mesh.js

@@ -0,0 +1,131 @@
+// math.gl
+// SPDX-License-Identifier: MIT
+// SPDX-FileCopyrightText: Copyright (c) vis.gl contributors
+export const subtract = (a, b) => [a[0] - b[0], a[1] - b[1], a[2] - b[2]];
+export const dot = (a, b) => a[0] * b[0] + a[1] * b[1] + a[2] * b[2];
+export const cross = (a, b) => [
+    a[1] * b[2] - a[2] * b[1],
+    a[2] * b[0] - a[0] * b[2],
+    a[0] * b[1] - a[1] * b[0]
+];
+export const normalize = (v) => {
+    const n = Math.hypot(...v);
+    return v.map(x => x / n);
+};
+export const center = (points) => [0, 1, 2].map(i => points.reduce((s, p) => s + p[i], 0) / points.length);
+/** Independent construction from symmetric coordinates, with a small convex hull
+ * used only for the 4–20 vertices of the regular solids. No imported mesh tables.
+ */
+export function createSolid(solid, subdivisions = 0) {
+    if (!Number.isInteger(subdivisions) || subdivisions < 0 || subdivisions > 3)
+        throw new RangeError('Subdivision must be 0–3');
+    let vertices = [];
+    switch (solid) {
+        case 'tetrahedron':
+            vertices = [
+                [1, 1, 1],
+                [1, -1, -1],
+                [-1, 1, -1],
+                [-1, -1, 1]
+            ];
+            break;
+        case 'cube':
+            for (const x of [-1, 1])
+                for (const y of [-1, 1])
+                    for (const z of [-1, 1])
+                        vertices.push([x, y, z]);
+            break;
+        case 'octahedron':
+            vertices = [
+                [1, 0, 0],
+                [-1, 0, 0],
+                [0, 1, 0],
+                [0, -1, 0],
+                [0, 0, 1],
+                [0, 0, -1]
+            ];
+            break;
+        case 'icosahedron': {
+            const phi = (1 + Math.sqrt(5)) / 2;
+            for (const a of [-1, 1])
+                for (const b of [-phi, phi])
+                    vertices.push([0, a, b], [a, b, 0], [b, 0, a]);
+            break;
+        }
+        case 'dodecahedron': {
+            // The dual of the icosahedron: its face centers are dodecahedron vertices.
+            const ico = createSolid('icosahedron');
+            vertices = ico.faces.map(face => center(face.map(i => ico.vertices[i])));
+            break;
+        }
+        default:
+            throw new RangeError('Unknown solid');
+    }
+    vertices = vertices.map(normalize);
+    let faces = convexFaces(vertices);
+    // Triangulate polygonal faces around a shared spherical face center before
+    // recursive subdivision; triangle edges and all subsequent midpoints are shared.
+    if (subdivisions)
+        faces = faces.flatMap(face => {
+            if (face.length === 3)
+                return [face];
+            const middle = vertices.push(normalize(center(face.map(i => vertices[i])))) - 1;
+            return face.map((v, i) => [v, face[(i + 1) % face.length], middle]);
+        });
+    for (let level = 0; level < subdivisions; level++) {
+        const midpoints = new Map();
+        const midpoint = (a, b) => {
+            const key = edgeKey(a, b);
+            if (!midpoints.has(key))
+                midpoints.set(key, vertices.push(normalize(center([vertices[a], vertices[b]]))) - 1);
+            return midpoints.get(key);
+        };
+        faces = faces.flatMap(([a, b, c]) => {
+            const ab = midpoint(a, b), bc = midpoint(b, c), ca = midpoint(c, a);
+            return [
+                [a, ab, ca],
+                [ab, b, bc],
+                [ca, bc, c],
+                [ab, bc, ca]
+            ];
+        });
+    }
+    return { vertices, faces };
+}
+export function edgeKey(a, b) {
+    return a < b ? `${a}:${b}` : `${b}:${a}`;
+}
+export function faceNormal(mesh, face) {
+    return normalize(cross(subtract(mesh.vertices[face[1]], mesh.vertices[face[0]]), subtract(mesh.vertices[face[2]], mesh.vertices[face[0]])));
+}
+function convexFaces(vertices) {
+    const found = new Map();
+    for (let a = 0; a < vertices.length; a++)
+        for (let b = a + 1; b < vertices.length; b++)
+            for (let c = b + 1; c < vertices.length; c++) {
+                let n = cross(subtract(vertices[b], vertices[a]), subtract(vertices[c], vertices[a]));
+                if (Math.hypot(...n) < 1e-8)
+                    continue;
+                n = normalize(n);
+                let distances = vertices.map(v => dot(n, subtract(v, vertices[a])));
+                if (distances.some(d => d > 1e-8) && distances.some(d => d < -1e-8))
+                    continue;
+                if (dot(n, vertices[a]) < 0) {
+                    n = n.map(x => -x);
+                    distances = distances.map(x => -x);
+                }
+                const ids = distances.flatMap((d, i) => (Math.abs(d) < 1e-8 ? [i] : []));
+                const key = ids.join(':');
+                if (found.has(key))
+                    continue;
+                const origin = center(ids.map(i => vertices[i]));
+                const x = normalize(subtract(vertices[ids[0]], origin)), y = cross(n, x);
+                ids.sort((i, j) => {
+                    const p = subtract(vertices[i], origin), q = subtract(vertices[j], origin);
+                    return Math.atan2(dot(p, y), dot(p, x)) - Math.atan2(dot(q, y), dot(q, x));
+                });
+                found.set(key, ids);
+            }
+    return [...found.values()];
+}
+//# sourceMappingURL=mesh.js.map

package/dist/discontinuous/myriahedral/mesh.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"mesh.js","sourceRoot":"","sources":["../../../src/discontinuous/myriahedral/mesh.ts"],"names":[],"mappings":"AAAA,UAAU;AACV,+BAA+B;AAC/B,4DAA4D;AAK5D,MAAM,CAAC,MAAM,QAAQ,GAAG,CAAC,CAAO,EAAE,CAAO,EAAQ,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;AAC5F,MAAM,CAAC,MAAM,GAAG,GAAG,CAAC,CAAO,EAAE,CAAO,EAAU,EAAE,CAAC,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;AACzF,MAAM,CAAC,MAAM,KAAK,GAAG,CAAC,CAAO,EAAE,CAAO,EAAQ,EAAE,CAAC;IAC/C,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;IACzB,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;IACzB,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;CAC1B,CAAC;AACF,MAAM,CAAC,MAAM,SAAS,GAAG,CAAC,CAAO,EAAQ,EAAE;IACzC,MAAM,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,GAAG,CAAC,CAAC,CAAC;IAC3B,OAAO,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,CAAS,CAAC;AACnC,CAAC,CAAC;AACF,MAAM,CAAC,MAAM,MAAM,GAAG,CAAC,MAAc,EAAQ,EAAE,CAC7C,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,MAAM,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,MAAM,CAAC,MAAM,CAAS,CAAC;AAEnF;;GAEG;AACH,MAAM,UAAU,WAAW,CAAC,KAAY,EAAE,YAAY,GAAG,CAAC;IACxD,IAAI,CAAC,MAAM,CAAC,SAAS,CAAC,YAAY,CAAC,IAAI,YAAY,GAAG,CAAC,IAAI,YAAY,GAAG,CAAC;QACzE,MAAM,IAAI,UAAU,CAAC,yBAAyB,CAAC,CAAC;IAClD,IAAI,QAAQ,GAAW,EAAE,CAAC;IAC1B,QAAQ,KAAK,EAAE,CAAC;QACd,KAAK,aAAa;YAChB,QAAQ,GAAG;gBACT,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;gBACT,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;gBACX,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;gBACX,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC;aACZ,CAAC;YACF,MAAM;QACR,KAAK,MAAM;YACT,KAAK,MAAM,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC;gBACrB,KAAK,MAAM,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC;oBAAE,KAAK,MAAM,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC;wBAAE,QAAQ,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;YAC7E,MAAM;QACR,KAAK,YAAY;YACf,QAAQ,GAAG;gBACT,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;gBACT,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;gBACV,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;gBACT,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC;gBACV,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;gBACT,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;aACX,CAAC;YACF,MAAM;QACR,KAAK,aAAa,CAAC,CAAC,CAAC;YACnB,MAAM,GAAG,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;YACnC,KAAK,MAAM,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC;gBACrB,KAAK,MAAM,CAAC,IAAI,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC;oBAAE,QAAQ,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;YAC9E,MAAM;QACR,CAAC;QACD,KAAK,cAAc,CAAC,CAAC,CAAC;YACpB,2EAA2E;YAC3E,MAAM,GAAG,GAAG,WAAW,CAAC,aAAa,CAAC,CAAC;YACvC,QAAQ,GAAG,GAAG,CAAC,KAAK,CAAC,GAAG,CAAC,IAAI,CAAC,EAAE,CAAC,MAAM,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,GAAG,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;YACzE,MAAM;QACR,CAAC;QACD;YACE,MAAM,IAAI,UAAU,CAAC,eAAe,CAAC,CAAC;IAC1C,CAAC;IACD,QAAQ,GAAG,QAAQ,CAAC,GAAG,CAAC,SAAS,CAAC,CAAC;IACnC,IAAI,KAAK,GAAG,WAAW,CAAC,QAAQ,CAAC,CAAC;IAClC,2EAA2E;IAC3E,iFAAiF;IACjF,IAAI,YAAY;QACd,KAAK,GAAG,KAAK,CAAC,OAAO,CAAC,IAAI,CAAC,EAAE;YAC3B,IAAI,IAAI,CAAC,MAAM,KAAK,CAAC;gBAAE,OAAO,CAAC,IAAI,CAAC,CAAC;YACrC,MAAM,MAAM,GAAG,QAAQ,CAAC,IAAI,CAAC,SAAS,CAAC,MAAM,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;YAChF,OAAO,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,MAAM,CAAC,EAAE,MAAM,CAAC,CAAC,CAAC;QACtE,CAAC,CAAC,CAAC;IACL,KAAK,IAAI,KAAK,GAAG,CAAC,EAAE,KAAK,GAAG,YAAY,EAAE,KAAK,EAAE,EAAE,CAAC;QAClD,MAAM,SAAS,GAAG,IAAI,GAAG,EAAkB,CAAC;QAC5C,MAAM,QAAQ,GAAG,CAAC,CAAS,EAAE,CAAS,EAAU,EAAE;YAChD,MAAM,GAAG,GAAG,OAAO,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;YAC1B,IAAI,CAAC,SAAS,CAAC,GAAG,CAAC,GAAG,CAAC;gBACrB,SAAS,CAAC,GAAG,CAAC,GAAG,EAAE,QAAQ,CAAC,IAAI,CAAC,SAAS,CAAC,MAAM,CAAC,CAAC,QAAQ,CAAC,CAAC,CAAC,EAAE,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;YACvF,OAAO,SAAS,CAAC,GAAG,CAAC,GAAG,CAAE,CAAC;QAC7B,CAAC,CAAC;QACF,KAAK,GAAG,KAAK,CAAC,OAAO,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,EAAE,EAAE;YAClC,MAAM,EAAE,GAAG,QAAQ,CAAC,CAAC,EAAE,CAAC,CAAC,EACvB,EAAE,GAAG,QAAQ,CAAC,CAAC,EAAE,CAAC,CAAC,EACnB,EAAE,GAAG,QAAQ,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;YACtB,OAAO;gBACL,CAAC,CAAC,EAAE,EAAE,EAAE,EAAE,CAAC;gBACX,CAAC,EAAE,EAAE,CAAC,EAAE,EAAE,CAAC;gBACX,CAAC,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC;gBACX,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC;aACb,CAAC;QACJ,CAAC,CAAC,CAAC;IACL,CAAC;IACD,OAAO,EAAC,QAAQ,EAAE,KAAK,EAAC,CAAC;AAC3B,CAAC;AACD,MAAM,UAAU,OAAO,CAAC,CAAS,EAAE,CAAS;IAC1C,OAAO,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,IAAI,CAAC,EAAE,CAAC;AAC3C,CAAC;AACD,MAAM,UAAU,UAAU,CAAC,IAAgB,EAAE,IAAc;IACzD,OAAO,SAAS,CACd,KAAK,CACH,QAAQ,CAAC,IAAI,CAAC,QAAQ,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,QAAQ,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,EACxD,QAAQ,CAAC,IAAI,CAAC,QAAQ,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,QAAQ,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CACzD,CACF,CAAC;AACJ,CAAC;AACD,SAAS,WAAW,CAAC,QAAgB;IACnC,MAAM,KAAK,GAAG,IAAI,GAAG,EAAoB,CAAC;IAC1C,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,QAAQ,CAAC,MAAM,EAAE,CAAC,EAAE;QACtC,KAAK,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,QAAQ,CAAC,MAAM,EAAE,CAAC,EAAE;YAC1C,KAAK,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,QAAQ,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;gBAC7C,IAAI,CAAC,GAAG,KAAK,CAAC,QAAQ,CAAC,QAAQ,CAAC,CAAC,CAAC,EAAE,QAAQ,CAAC,CAAC,CAAC,CAAC,EAAE,QAAQ,CAAC,QAAQ,CAAC,CAAC,CAAC,EAAE,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;gBACtF,IAAI,IAAI,CAAC,KAAK,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI;oBAAE,SAAS;gBACtC,CAAC,GAAG,SAAS,CAAC,CAAC,CAAC,CAAC;gBACjB,IAAI,SAAS,GAAG,QAAQ,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,EAAE,QAAQ,CAAC,CAAC,EAAE,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;gBACpE,IAAI,SAAS,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,IAAI,CAAC,IAAI,SAAS,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,IAAI,CAAC;oBAAE,SAAS;gBAC9E,IAAI,GAAG,CAAC,CAAC,EAAE,QAAQ,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC;oBAC5B,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAS,CAAC;oBAC3B,SAAS,GAAG,SAAS,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;gBACrC,CAAC;gBACD,MAAM,GAAG,GAAG,SAAS,CAAC,OAAO,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;gBACzE,MAAM,GAAG,GAAG,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC;gBAC1B,IAAI,KAAK,CAAC,GAAG,CAAC,GAAG,CAAC;oBAAE,SAAS;gBAC7B,MAAM,MAAM,GAAG,MAAM,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;gBACjD,MAAM,CAAC,GAAG,SAAS,CAAC,QAAQ,CAAC,QAAQ,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC,EACrD,CAAC,GAAG,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;gBAClB,GAAG,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE;oBAChB,MAAM,CAAC,GAAG,QAAQ,CAAC,QAAQ,CAAC,CAAC,CAAC,EAAE,MAAM,CAAC,EACrC,CAAC,GAAG,QAAQ,CAAC,QAAQ,CAAC,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC;oBACpC,OAAO,IAAI,CAAC,KAAK,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;gBAC7E,CAAC,CAAC,CAAC;gBACH,KAAK,CAAC,GAAG,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC;YACtB,CAAC;IACL,OAAO,CAAC,GAAG,KAAK,CAAC,MAAM,EAAE,CAAC,CAAC;AAC7B,CAAC"}

package/dist/discontinuous/myriahedral/unfold.d.ts

@@ -0,0 +1,56 @@
+import { Matrix4 } from '@math.gl/core';
+import { HingeTree, type HingePanel } from '@math.gl/geometry';
+import { type Polyhedron, type Vec3 } from "./mesh.js";
+export type MeshEdge = {
+    a: number;
+    b: number;
+    faces: number[];
+    weight: number;
+    hinge: boolean;
+};
+/** Kruskal maximum spanning tree, with stable input-index tie breaking. */
+export declare function maximumSpanningTree(count: number, edges: {
+    faces: number[];
+    weight: number;
+}[]): number[];
+/** Internal layout builder for closed, consistently wound convex meshes.
+ * Each panel uses the original solid coordinate frame, so rest transforms are identity.
+ */
+export declare function createLayout(mesh: Polyhedron, seed?: number, options?: {
+    weight?: (edge: MeshEdge, legacyWeight: number) => number;
+    joins?: number[][];
+}): {
+    tree: HingeTree;
+    rootTransform: Matrix4;
+    evaluate: (amount: number) => Vec3[][];
+    mesh: Polyhedron;
+    edges: MeshEdge[];
+    panels: HingePanel[];
+    targetAngles: Float64Array;
+    overlaps: number;
+};
+/** Data-only layout for worker transport; topology and geographic weights are already resolved. */
+export type PreparedHingeLayout = {
+    mesh: Polyhedron;
+    edges: MeshEdge[];
+    panels: HingePanel[];
+    targetAngles: Float64Array;
+    rootTransform: number[];
+    overlaps: number;
+};
+/** Restore pose evaluation without repeating spanning-tree selection or overlap counting. */
+export declare function restoreLayout(data: PreparedHingeLayout): {
+    tree: HingeTree;
+    rootTransform: Matrix4;
+    evaluate: (amount: number) => Vec3[][];
+    mesh: Polyhedron;
+    edges: MeshEdge[];
+    panels: HingePanel[];
+    targetAngles: Float64Array;
+    overlaps: number;
+};
+/** Count positive-area overlaps of convex faces in the final XY net.
+ * Sweep by X bounds, then separating axes; mere shared edges do not overlap.
+ */
+export declare function countOverlaps(faces: Vec3[][]): number;
+//# sourceMappingURL=unfold.d.ts.map

package/dist/discontinuous/myriahedral/unfold.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"unfold.d.ts","sourceRoot":"","sources":["../../../src/discontinuous/myriahedral/unfold.ts"],"names":[],"mappings":"AAIA,OAAO,EAAC,OAAO,EAAC,MAAM,eAAe,CAAC;AACtC,OAAO,EAAC,SAAS,EAAE,KAAK,UAAU,EAAC,MAAM,mBAAmB,CAAC;AAC7D,OAAO,EAQL,KAAK,UAAU,EACf,KAAK,IAAI,EACV,kBAAe;AAEhB,MAAM,MAAM,QAAQ,GAAG;IAAC,CAAC,EAAE,MAAM,CAAC;IAAC,CAAC,EAAE,MAAM,CAAC;IAAC,KAAK,EAAE,MAAM,EAAE,CAAC;IAAC,MAAM,EAAE,MAAM,CAAC;IAAC,KAAK,EAAE,OAAO,CAAA;CAAC,CAAC;AAC/F,2EAA2E;AAC3E,wBAAgB,mBAAmB,CACjC,KAAK,EAAE,MAAM,EACb,KAAK,EAAE;IAAC,KAAK,EAAE,MAAM,EAAE,CAAC;IAAC,MAAM,EAAE,MAAM,CAAA;CAAC,EAAE,GACzC,MAAM,EAAE,CAkCV;AAED;;GAEG;AACH,wBAAgB,YAAY,CAC1B,IAAI,EAAE,UAAU,EAChB,IAAI,SAAI,EACR,OAAO,GAAE;IACP,MAAM,CAAC,EAAE,CAAC,IAAI,EAAE,QAAQ,EAAE,YAAY,EAAE,MAAM,KAAK,MAAM,CAAC;IAC1D,KAAK,CAAC,EAAE,MAAM,EAAE,EAAE,CAAC;CACf;;;uBAuHoB,MAAM,KAAG,IAAI,EAAE,EAAE;UAvBrC,UAAU;WACT,QAAQ,EAAE;YACT,UAAU,EAAE;kBACN,YAAY;cAEhB,MAAM;EATjB;AAED,mGAAmG;AACnG,MAAM,MAAM,mBAAmB,GAAG;IAChC,IAAI,EAAE,UAAU,CAAC;IACjB,KAAK,EAAE,QAAQ,EAAE,CAAC;IAClB,MAAM,EAAE,UAAU,EAAE,CAAC;IACrB,YAAY,EAAE,YAAY,CAAC;IAC3B,aAAa,EAAE,MAAM,EAAE,CAAC;IACxB,QAAQ,EAAE,MAAM,CAAC;CAClB,CAAC;AAEF,6FAA6F;AAC7F,wBAAgB,aAAa,CAAC,IAAI,EAAE,mBAAmB;;;uBAc3B,MAAM,KAAG,IAAI,EAAE,EAAE;UAvBrC,UAAU;WACT,QAAQ,EAAE;YACT,UAAU,EAAE;kBACN,YAAY;cAEhB,MAAM;EA8BjB;AAED;;GAEG;AACH,wBAAgB,aAAa,CAAC,KAAK,EAAE,IAAI,EAAE,EAAE,GAAG,MAAM,CAmCrD"}