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

package/dist/discontinuous/myriahedral/unfold.js

@@ -0,0 +1,191 @@
+// math.gl
+// SPDX-License-Identifier: MIT
+// SPDX-FileCopyrightText: Copyright (c) vis.gl contributors
+import { Matrix4 } from '@math.gl/core';
+import { HingeTree } from '@math.gl/geometry';
+import { center, cross, dot, edgeKey, faceNormal, normalize, subtract } from "./mesh.js";
+/** Kruskal maximum spanning tree, with stable input-index tie breaking. */
+export function maximumSpanningTree(count, edges) {
+    const roots = Array.from({ length: count }, (_, i) => i), ranks = new Uint32Array(count);
+    const find = (i) => {
+        while (roots[i] !== i) {
+            roots[i] = roots[roots[i]];
+            i = roots[i];
+        }
+        return i;
+    };
+    const order = edges
+        .map((edge, i) => {
+        if (edge.faces.length !== 2 ||
+            edge.faces.some(f => !Number.isInteger(f) || f < 0 || f >= count) ||
+            !Number.isFinite(edge.weight))
+            throw new Error('Invalid dual edge');
+        return i;
+    })
+        .sort((a, b) => edges[b].weight - edges[a].weight || a - b);
+    const selected = [];
+    for (const i of order) {
+        let a = find(edges[i].faces[0]), b = find(edges[i].faces[1]);
+        if (a === b)
+            continue;
+        if (ranks[a] < ranks[b])
+            [a, b] = [b, a];
+        roots[b] = a;
+        if (ranks[a] === ranks[b])
+            ranks[a]++;
+        selected.push(i);
+        if (selected.length === count - 1)
+            break;
+    }
+    if (selected.length !== count - 1)
+        throw new Error('Disconnected face graph');
+    return selected;
+}
+/** Internal layout builder for closed, consistently wound convex meshes.
+ * Each panel uses the original solid coordinate frame, so rest transforms are identity.
+ */
+export function createLayout(mesh, seed = 1, options = {}) {
+    const normals = mesh.faces.map(f => faceNormal(mesh, f));
+    const map = new Map();
+    mesh.faces.forEach((face, i) => face.forEach((a, j) => {
+        const b = face[(j + 1) % face.length], key = edgeKey(a, b);
+        if (!map.has(key))
+            map.set(key, { a: Math.min(a, b), b: Math.max(a, b), faces: [], weight: 0, hinge: false });
+        map.get(key).faces.push(i);
+    }));
+    const edges = [...map.values()];
+    for (const edge of edges) {
+        if (edge.faces.length !== 2)
+            throw new Error('Expected a closed two-manifold mesh');
+        // Favor shorter joins, with seeded variation to explore different unfoldings.
+        // This is a layout heuristic, not an overlap-minimizing objective.
+        let hash = Math.imul(edge.a + 1, 73856093) ^ Math.imul(edge.b + 1, 19349663) ^ Math.imul(seed, 83492791);
+        hash = Math.imul(hash ^ (hash >>> 16), 0x45d9f3b);
+        hash ^= hash >>> 16;
+        edge.weight =
+            -Math.hypot(...subtract(mesh.vertices[edge.a], mesh.vertices[edge.b])) +
+                (0.25 * (hash >>> 0)) / 4294967296;
+    }
+    for (const edge of edges)
+        if (options.weight)
+            edge.weight = options.weight(edge, edge.weight);
+    const allowed = options.joins && new Set(options.joins.map(([a, b]) => edgeKey(a, b)));
+    const candidates = allowed
+        ? edges.filter(e => allowed.has(edgeKey(e.faces[0], e.faces[1])))
+        : edges;
+    if (allowed && candidates.length !== allowed.size)
+        throw new Error('Preset contains a non-adjacent join');
+    const selected = maximumSpanningTree(mesh.faces.length, candidates);
+    const adjacency = mesh.faces.map(() => []);
+    for (const i of selected) {
+        const edge = candidates[i];
+        edge.hinge = true;
+        const [a, b] = edge.faces;
+        adjacency[a].push({ face: b, edge });
+        adjacency[b].push({ face: a, edge });
+    }
+    const panels = mesh.faces.map(() => ({ parent: -1 }));
+    const targetAngles = new Float64Array(panels.length);
+    const order = [0], seen = new Set([0]);
+    for (let i = 0; i < order.length; i++) {
+        const parent = order[i];
+        for (const { face, edge } of adjacency[parent]) {
+            if (seen.has(face))
+                continue;
+            seen.add(face);
+            order.push(face);
+            const axis = normalize(subtract(mesh.vertices[edge.b], mesh.vertices[edge.a]));
+            panels[face] = { parent, pivot: mesh.vertices[edge.a], axis };
+            targetAngles[face] = Math.atan2(dot(axis, cross(normals[face], normals[parent])), dot(normals[face], normals[parent]));
+        }
+    }
+    const origin = center(mesh.faces[0].map(i => mesh.vertices[i]));
+    const x = normalize(subtract(mesh.vertices[mesh.faces[0][1]], mesh.vertices[mesh.faces[0][0]]));
+    const z = normals[0], y = cross(z, x);
+    const rootTransform = new Matrix4([
+        x[0],
+        y[0],
+        z[0],
+        0,
+        x[1],
+        y[1],
+        z[1],
+        0,
+        x[2],
+        y[2],
+        z[2],
+        0,
+        -dot(x, origin),
+        -dot(y, origin),
+        -dot(z, origin),
+        1
+    ]);
+    const layout = restoreLayout({
+        mesh,
+        edges,
+        panels,
+        targetAngles,
+        rootTransform: Array.from(rootTransform),
+        overlaps: 0
+    });
+    layout.overlaps = countOverlaps(layout.evaluate(1));
+    return layout;
+}
+/** Restore pose evaluation without repeating spanning-tree selection or overlap counting. */
+export function restoreLayout(data) {
+    const { mesh, panels, targetAngles } = data;
+    if (panels.length !== mesh.faces.length ||
+        targetAngles.length !== panels.length ||
+        !targetAngles.every(Number.isFinite) ||
+        data.rootTransform.length !== 16 ||
+        !data.rootTransform.every(Number.isFinite))
+        throw new Error('Invalid prepared hinge layout');
+    const tree = new HingeTree(panels);
+    const rootTransform = new Matrix4(data.rootTransform);
+    const poses = new Float64Array(panels.length * 16), angles = new Float64Array(panels.length);
+    const evaluate = (amount) => {
+        if (!Number.isFinite(amount) || amount < 0 || amount > 1)
+            throw new RangeError('Unfold amount must be within [0, 1]');
+        for (let i = 0; i < angles.length; i++)
+            angles[i] = targetAngles[i] * amount;
+        tree.evaluate(angles, poses, rootTransform);
+        const matrix = new Matrix4();
+        return mesh.faces.map((face, i) => {
+            matrix.copy(poses.subarray(i * 16, i * 16 + 16));
+            return face.map(v => Array.from(matrix.transformAsPoint(mesh.vertices[v])));
+        });
+    };
+    return { ...data, tree, rootTransform, evaluate };
+}
+/** 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 function countOverlaps(faces) {
+    const bounds = faces
+        .map((points, i) => ({
+        i,
+        minX: Math.min(...points.map(p => p[0])),
+        maxX: Math.max(...points.map(p => p[0])),
+        minY: Math.min(...points.map(p => p[1])),
+        maxY: Math.max(...points.map(p => p[1]))
+    }))
+        .sort((a, b) => a.minX - b.minX);
+    let count = 0;
+    for (let i = 0; i < bounds.length; i++)
+        for (let j = i + 1; j < bounds.length && bounds[j].minX < bounds[i].maxX - 1e-8; j++) {
+            const a = bounds[i], b = bounds[j];
+            if (b.minY >= a.maxY - 1e-8 || a.minY >= b.maxY - 1e-8)
+                continue;
+            const p = faces[a.i], q = faces[b.i];
+            const separate = [p, q].some(poly => poly.some((v, k) => {
+                const next = poly[(k + 1) % poly.length], dx = next[0] - v[0], dy = next[1] - v[1], length = Math.hypot(dx, dy);
+                const project = (point) => (point[0] * dy - point[1] * dx) / length;
+                const pp = p.map(project), qq = q.map(project);
+                return (Math.max(...pp) <= Math.min(...qq) + 1e-8 || Math.max(...qq) <= Math.min(...pp) + 1e-8);
+            }));
+            if (!separate)
+                count++;
+        }
+    return count;
+}
+//# sourceMappingURL=unfold.js.map

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

@@ -0,0 +1 @@
+{"version":3,"file":"unfold.js","sourceRoot":"","sources":["../../../src/discontinuous/myriahedral/unfold.ts"],"names":[],"mappings":"AAAA,UAAU;AACV,+BAA+B;AAC/B,4DAA4D;AAE5D,OAAO,EAAC,OAAO,EAAC,MAAM,eAAe,CAAC;AACtC,OAAO,EAAC,SAAS,EAAkB,MAAM,mBAAmB,CAAC;AAC7D,OAAO,EACL,MAAM,EACN,KAAK,EACL,GAAG,EACH,OAAO,EACP,UAAU,EACV,SAAS,EACT,QAAQ,EAGT,kBAAe;AAGhB,2EAA2E;AAC3E,MAAM,UAAU,mBAAmB,CACjC,KAAa,EACb,KAA0C;IAE1C,MAAM,KAAK,GAAG,KAAK,CAAC,IAAI,CAAC,EAAC,MAAM,EAAE,KAAK,EAAC,EAAE,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,EACpD,KAAK,GAAG,IAAI,WAAW,CAAC,KAAK,CAAC,CAAC;IACjC,MAAM,IAAI,GAAG,CAAC,CAAS,EAAU,EAAE;QACjC,OAAO,KAAK,CAAC,CAAC,CAAC,KAAK,CAAC,EAAE,CAAC;YACtB,KAAK,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC;YAC3B,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC;QACf,CAAC;QACD,OAAO,CAAC,CAAC;IACX,CAAC,CAAC;IACF,MAAM,KAAK,GAAG,KAAK;SAChB,GAAG,CAAC,CAAC,IAAI,EAAE,CAAC,EAAE,EAAE;QACf,IACE,IAAI,CAAC,KAAK,CAAC,MAAM,KAAK,CAAC;YACvB,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,MAAM,CAAC,SAAS,CAAC,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,IAAI,KAAK,CAAC;YACjE,CAAC,MAAM,CAAC,QAAQ,CAAC,IAAI,CAAC,MAAM,CAAC;YAE7B,MAAM,IAAI,KAAK,CAAC,mBAAmB,CAAC,CAAC;QACvC,OAAO,CAAC,CAAC;IACX,CAAC,CAAC;SACD,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,MAAM,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC,MAAM,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC;IAC9D,MAAM,QAAQ,GAAa,EAAE,CAAC;IAC9B,KAAK,MAAM,CAAC,IAAI,KAAK,EAAE,CAAC;QACtB,IAAI,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,EAC7B,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC;QAC9B,IAAI,CAAC,KAAK,CAAC;YAAE,SAAS;QACtB,IAAI,KAAK,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC;YAAE,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;QACzC,KAAK,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;QACb,IAAI,KAAK,CAAC,CAAC,CAAC,KAAK,KAAK,CAAC,CAAC,CAAC;YAAE,KAAK,CAAC,CAAC,CAAC,EAAE,CAAC;QACtC,QAAQ,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;QACjB,IAAI,QAAQ,CAAC,MAAM,KAAK,KAAK,GAAG,CAAC;YAAE,MAAM;IAC3C,CAAC;IACD,IAAI,QAAQ,CAAC,MAAM,KAAK,KAAK,GAAG,CAAC;QAAE,MAAM,IAAI,KAAK,CAAC,yBAAyB,CAAC,CAAC;IAC9E,OAAO,QAAQ,CAAC;AAClB,CAAC;AAED;;GAEG;AACH,MAAM,UAAU,YAAY,CAC1B,IAAgB,EAChB,IAAI,GAAG,CAAC,EACR,UAGI,EAAE;IAEN,MAAM,OAAO,GAAG,IAAI,CAAC,KAAK,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,UAAU,CAAC,IAAI,EAAE,CAAC,CAAC,CAAC,CAAC;IACzD,MAAM,GAAG,GAAG,IAAI,GAAG,EAAoB,CAAC;IACxC,IAAI,CAAC,KAAK,CAAC,OAAO,CAAC,CAAC,IAAI,EAAE,CAAC,EAAE,EAAE,CAC7B,IAAI,CAAC,OAAO,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE;QACpB,MAAM,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,MAAM,CAAC,EACnC,GAAG,GAAG,OAAO,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;QACtB,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC;YACf,GAAG,CAAC,GAAG,CAAC,GAAG,EAAE,EAAC,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,KAAK,EAAE,EAAE,EAAE,MAAM,EAAE,CAAC,EAAE,KAAK,EAAE,KAAK,EAAC,CAAC,CAAC;QAC3F,GAAG,CAAC,GAAG,CAAC,GAAG,CAAE,CAAC,KAAK,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;IAC9B,CAAC,CAAC,CACH,CAAC;IACF,MAAM,KAAK,GAAG,CAAC,GAAG,GAAG,CAAC,MAAM,EAAE,CAAC,CAAC;IAChC,KAAK,MAAM,IAAI,IAAI,KAAK,EAAE,CAAC;QACzB,IAAI,IAAI,CAAC,KAAK,CAAC,MAAM,KAAK,CAAC;YAAE,MAAM,IAAI,KAAK,CAAC,qCAAqC,CAAC,CAAC;QACpF,8EAA8E;QAC9E,mEAAmE;QACnE,IAAI,IAAI,GACN,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,EAAE,QAAQ,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,EAAE,QAAQ,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,IAAI,EAAE,QAAQ,CAAC,CAAC;QAChG,IAAI,GAAG,IAAI,CAAC,IAAI,CAAC,IAAI,GAAG,CAAC,IAAI,KAAK,EAAE,CAAC,EAAE,SAAS,CAAC,CAAC;QAClD,IAAI,IAAI,IAAI,KAAK,EAAE,CAAC;QACpB,IAAI,CAAC,MAAM;YACT,CAAC,IAAI,CAAC,KAAK,CAAC,GAAG,QAAQ,CAAC,IAAI,CAAC,QAAQ,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,QAAQ,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC;gBACtE,CAAC,IAAI,GAAG,CAAC,IAAI,KAAK,CAAC,CAAC,CAAC,GAAG,UAAU,CAAC;IACvC,CAAC;IACD,KAAK,MAAM,IAAI,IAAI,KAAK;QAAE,IAAI,OAAO,CAAC,MAAM;YAAE,IAAI,CAAC,MAAM,GAAG,OAAO,CAAC,MAAM,CAAC,IAAI,EAAE,IAAI,CAAC,MAAM,CAAC,CAAC;IAC9F,MAAM,OAAO,GAAG,OAAO,CAAC,KAAK,IAAI,IAAI,GAAG,CAAC,OAAO,CAAC,KAAK,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,EAAE,CAAC,OAAO,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;IACvF,MAAM,UAAU,GAAG,OAAO;QACxB,CAAC,CAAC,KAAK,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,CAAC,OAAO,CAAC,GAAG,CAAC,OAAO,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QACjE,CAAC,CAAC,KAAK,CAAC;IACV,IAAI,OAAO,IAAI,UAAU,CAAC,MAAM,KAAK,OAAO,CAAC,IAAI;QAC/C,MAAM,IAAI,KAAK,CAAC,qCAAqC,CAAC,CAAC;IACzD,MAAM,QAAQ,GAAG,mBAAmB,CAAC,IAAI,CAAC,KAAK,CAAC,MAAM,EAAE,UAAU,CAAC,CAAC;IACpE,MAAM,SAAS,GAAuC,IAAI,CAAC,KAAK,CAAC,GAAG,CAAC,GAAG,EAAE,CAAC,EAAE,CAAC,CAAC;IAC/E,KAAK,MAAM,CAAC,IAAI,QAAQ,EAAE,CAAC;QACzB,MAAM,IAAI,GAAG,UAAU,CAAC,CAAC,CAAC,CAAC;QAC3B,IAAI,CAAC,KAAK,GAAG,IAAI,CAAC;QAClB,MAAM,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC;QAC1B,SAAS,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,EAAC,IAAI,EAAE,CAAC,EAAE,IAAI,EAAC,CAAC,CAAC;QACnC,SAAS,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,EAAC,IAAI,EAAE,CAAC,EAAE,IAAI,EAAC,CAAC,CAAC;IACrC,CAAC;IACD,MAAM,MAAM,GAAiB,IAAI,CAAC,KAAK,CAAC,GAAG,CAAC,GAAG,EAAE,CAAC,CAAC,EAAC,MAAM,EAAE,CAAC,CAAC,EAAC,CAAC,CAAC,CAAC;IAClE,MAAM,YAAY,GAAG,IAAI,YAAY,CAAC,MAAM,CAAC,MAAM,CAAC,CAAC;IACrD,MAAM,KAAK,GAAG,CAAC,CAAC,CAAC,EACf,IAAI,GAAG,IAAI,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IACtB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,KAAK,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;QACtC,MAAM,MAAM,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC;QACxB,KAAK,MAAM,EAAC,IAAI,EAAE,IAAI,EAAC,IAAI,SAAS,CAAC,MAAM,CAAC,EAAE,CAAC;YAC7C,IAAI,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC;gBAAE,SAAS;YAC7B,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC;YACf,KAAK,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC;YACjB,MAAM,IAAI,GAAG,SAAS,CAAC,QAAQ,CAAC,IAAI,CAAC,QAAQ,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,QAAQ,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;YAC/E,MAAM,CAAC,IAAI,CAAC,GAAG,EAAC,MAAM,EAAE,KAAK,EAAE,IAAI,CAAC,QAAQ,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,IAAI,EAAC,CAAC;YAC5D,YAAY,CAAC,IAAI,CAAC,GAAG,IAAI,CAAC,KAAK,CAC7B,GAAG,CAAC,IAAI,EAAE,KAAK,CAAC,OAAO,CAAC,IAAI,CAAC,EAAE,OAAO,CAAC,MAAM,CAAC,CAAC,CAAC,EAChD,GAAG,CAAC,OAAO,CAAC,IAAI,CAAC,EAAE,OAAO,CAAC,MAAM,CAAC,CAAC,CACpC,CAAC;QACJ,CAAC;IACH,CAAC;IACD,MAAM,MAAM,GAAG,MAAM,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IAChE,MAAM,CAAC,GAAG,SAAS,CAAC,QAAQ,CAAC,IAAI,CAAC,QAAQ,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,QAAQ,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IAChG,MAAM,CAAC,GAAG,OAAO,CAAC,CAAC,CAAC,EAClB,CAAC,GAAG,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;IAClB,MAAM,aAAa,GAAG,IAAI,OAAO,CAAC;QAChC,CAAC,CAAC,CAAC,CAAC;QACJ,CAAC,CAAC,CAAC,CAAC;QACJ,CAAC,CAAC,CAAC,CAAC;QACJ,CAAC;QACD,CAAC,CAAC,CAAC,CAAC;QACJ,CAAC,CAAC,CAAC,CAAC;QACJ,CAAC,CAAC,CAAC,CAAC;QACJ,CAAC;QACD,CAAC,CAAC,CAAC,CAAC;QACJ,CAAC,CAAC,CAAC,CAAC;QACJ,CAAC,CAAC,CAAC,CAAC;QACJ,CAAC;QACD,CAAC,GAAG,CAAC,CAAC,EAAE,MAAM,CAAC;QACf,CAAC,GAAG,CAAC,CAAC,EAAE,MAAM,CAAC;QACf,CAAC,GAAG,CAAC,CAAC,EAAE,MAAM,CAAC;QACf,CAAC;KACF,CAAC,CAAC;IACH,MAAM,MAAM,GAAG,aAAa,CAAC;QAC3B,IAAI;QACJ,KAAK;QACL,MAAM;QACN,YAAY;QACZ,aAAa,EAAE,KAAK,CAAC,IAAI,CAAC,aAAa,CAAC;QACxC,QAAQ,EAAE,CAAC;KACZ,CAAC,CAAC;IACH,MAAM,CAAC,QAAQ,GAAG,aAAa,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC;IACpD,OAAO,MAAM,CAAC;AAChB,CAAC;AAYD,6FAA6F;AAC7F,MAAM,UAAU,aAAa,CAAC,IAAyB;IACrD,MAAM,EAAC,IAAI,EAAE,MAAM,EAAE,YAAY,EAAC,GAAG,IAAI,CAAC;IAC1C,IACE,MAAM,CAAC,MAAM,KAAK,IAAI,CAAC,KAAK,CAAC,MAAM;QACnC,YAAY,CAAC,MAAM,KAAK,MAAM,CAAC,MAAM;QACrC,CAAC,YAAY,CAAC,KAAK,CAAC,MAAM,CAAC,QAAQ,CAAC;QACpC,IAAI,CAAC,aAAa,CAAC,MAAM,KAAK,EAAE;QAChC,CAAC,IAAI,CAAC,aAAa,CAAC,KAAK,CAAC,MAAM,CAAC,QAAQ,CAAC;QAE1C,MAAM,IAAI,KAAK,CAAC,+BAA+B,CAAC,CAAC;IACnD,MAAM,IAAI,GAAG,IAAI,SAAS,CAAC,MAAM,CAAC,CAAC;IACnC,MAAM,aAAa,GAAG,IAAI,OAAO,CAAC,IAAI,CAAC,aAAa,CAAC,CAAC;IACtD,MAAM,KAAK,GAAG,IAAI,YAAY,CAAC,MAAM,CAAC,MAAM,GAAG,EAAE,CAAC,EAChD,MAAM,GAAG,IAAI,YAAY,CAAC,MAAM,CAAC,MAAM,CAAC,CAAC;IAC3C,MAAM,QAAQ,GAAG,CAAC,MAAc,EAAY,EAAE;QAC5C,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,MAAM,CAAC,IAAI,MAAM,GAAG,CAAC,IAAI,MAAM,GAAG,CAAC;YACtD,MAAM,IAAI,UAAU,CAAC,qCAAqC,CAAC,CAAC;QAC9D,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,MAAM,CAAC,MAAM,EAAE,CAAC,EAAE;YAAE,MAAM,CAAC,CAAC,CAAC,GAAG,YAAY,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC;QAC7E,IAAI,CAAC,QAAQ,CAAC,MAAM,EAAE,KAAK,EAAE,aAAa,CAAC,CAAC;QAC5C,MAAM,MAAM,GAAG,IAAI,OAAO,EAAE,CAAC;QAC7B,OAAO,IAAI,CAAC,KAAK,CAAC,GAAG,CAAC,CAAC,IAAI,EAAE,CAAC,EAAE,EAAE;YAChC,MAAM,CAAC,IAAI,CAAC,KAAK,CAAC,QAAQ,CAAC,CAAC,GAAG,EAAE,EAAE,CAAC,GAAG,EAAE,GAAG,EAAE,CAAC,CAAC,CAAC;YACjD,OAAO,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,KAAK,CAAC,IAAI,CAAC,MAAM,CAAC,gBAAgB,CAAC,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAS,CAAC,CAAC;QACtF,CAAC,CAAC,CAAC;IACL,CAAC,CAAC;IACF,OAAO,EAAC,GAAG,IAAI,EAAE,IAAI,EAAE,aAAa,EAAE,QAAQ,EAAC,CAAC;AAClD,CAAC;AAED;;GAEG;AACH,MAAM,UAAU,aAAa,CAAC,KAAe;IAC3C,MAAM,MAAM,GAAG,KAAK;SACjB,GAAG,CAAC,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC;QACnB,CAAC;QACD,IAAI,EAAE,IAAI,CAAC,GAAG,CAAC,GAAG,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QACxC,IAAI,EAAE,IAAI,CAAC,GAAG,CAAC,GAAG,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QACxC,IAAI,EAAE,IAAI,CAAC,GAAG,CAAC,GAAG,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QACxC,IAAI,EAAE,IAAI,CAAC,GAAG,CAAC,GAAG,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;KACzC,CAAC,CAAC;SACF,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,IAAI,GAAG,CAAC,CAAC,IAAI,CAAC,CAAC;IACnC,IAAI,KAAK,GAAG,CAAC,CAAC;IACd,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,MAAM,CAAC,MAAM,EAAE,CAAC,EAAE;QACpC,KAAK,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,MAAM,CAAC,MAAM,IAAI,MAAM,CAAC,CAAC,CAAC,CAAC,IAAI,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC,IAAI,GAAG,IAAI,EAAE,CAAC,EAAE,EAAE,CAAC;YACrF,MAAM,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,EACjB,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC;YAChB,IAAI,CAAC,CAAC,IAAI,IAAI,CAAC,CAAC,IAAI,GAAG,IAAI,IAAI,CAAC,CAAC,IAAI,IAAI,CAAC,CAAC,IAAI,GAAG,IAAI;gBAAE,SAAS;YACjE,MAAM,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,EAClB,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;YACjB,MAAM,QAAQ,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,IAAI,CAAC,IAAI,CAAC,EAAE,CAClC,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE;gBACjB,MAAM,IAAI,GAAG,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,MAAM,CAAC,EACtC,EAAE,GAAG,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,EACnB,EAAE,GAAG,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,EACnB,MAAM,GAAG,IAAI,CAAC,KAAK,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC;gBAC9B,MAAM,OAAO,GAAG,CAAC,KAAW,EAAE,EAAE,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,EAAE,GAAG,KAAK,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,GAAG,MAAM,CAAC;gBAC1E,MAAM,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,OAAO,CAAC,EACvB,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,OAAO,CAAC,CAAC;gBACtB,OAAO,CACL,IAAI,CAAC,GAAG,CAAC,GAAG,EAAE,CAAC,IAAI,IAAI,CAAC,GAAG,CAAC,GAAG,EAAE,CAAC,GAAG,IAAI,IAAI,IAAI,CAAC,GAAG,CAAC,GAAG,EAAE,CAAC,IAAI,IAAI,CAAC,GAAG,CAAC,GAAG,EAAE,CAAC,GAAG,IAAI,CACvF,CAAC;YACJ,CAAC,CAAC,CACH,CAAC;YACF,IAAI,CAAC,QAAQ;gBAAE,KAAK,EAAE,CAAC;QACzB,CAAC;IACH,OAAO,KAAK,CAAC;AACf,CAAC"}

package/dist/discontinuous/myriahedral.d.ts

@@ -0,0 +1,62 @@
+import { DiscontinuousMap, type GeographicPosition } from "./map-layout.js";
+import { createLayout, type PreparedHingeLayout } from "./myriahedral/unfold.js";
+import { type Polyhedron, type Solid } from "./myriahedral/mesh.js";
+export type MyriahedralMesh = Polyhedron;
+export type MyriahedralSolid = Solid;
+export type MyriahedralPreset = 'polyconic' | 'cylindric' | 'azimuthal' | 'hemispheres' | 'sealine' | 'coastline' | 'polyhedral';
+export type MyriahedralOptions = {
+    /** Independent cut strategy, not an implementation of an analytic namesake. Default cylindric. */
+    preset?: MyriahedralPreset;
+    /** Multiple of four, 12–96. Default 48. Graticule presets only. */
+    columns?: number;
+    /** Polyhedral preset only. Default icosahedron. */
+    solid?: MyriahedralSolid;
+    /** Polyhedral preset only, 0–3. Default 1. Zero keeps the original polygonal faces. */
+    subdivisions?: number;
+    /** Optional closed, outward-wound, convex planar panels on the unit sphere. */
+    mesh?: MyriahedralMesh;
+    /** Longitude in degrees for the backbone/secondary join preference. Default 20. */
+    centralMeridian?: number;
+    /** Positive output scale; default 1. Input and texture coordinates remain spherical. */
+    radius?: number;
+    /** Required for sealine/coastline. Called with degrees; no bundled geography or network access. */
+    isLand?: (position: GeographicPosition) => boolean;
+    /** Maximum angular interval for land sampling in degrees; default .25, range [.05, 10]. */
+    landSampleStep?: number;
+    /** Final retained-join weight override. Higher values are less likely to be cut. */
+    edgeWeight?: (edge: Readonly<MyriahedralEdge>, presetWeight: number) => number;
+};
+export type MyriahedralEdge = {
+    vertices: readonly [number, number];
+    panels: readonly [number, number];
+    midpoint: GeographicPosition;
+    arcLength: number;
+    landFraction?: number;
+};
+export declare const MYRIAHEDRAL_PRESETS: readonly MyriahedralPreset[];
+/** Build an independent spherical panel net and its rigid unfolding animation.
+ * Each call owns its mesh. Geographical costs are sampled once per edge.
+ */
+export declare function createMyriahedralLayout(options?: MyriahedralOptions): MyriahedralLayout;
+export type MyriahedralLayout = ReturnType<typeof createLayout> & {
+    preset: MyriahedralPreset;
+    radius: number;
+    /** Undefined without an isLand callback. */
+    cutLandRadians?: number;
+};
+export type PreparedMyriahedralLayout = PreparedHingeLayout & {
+    preset: MyriahedralPreset;
+    radius: number;
+    cutLandRadians?: number;
+};
+export declare function restoreMyriahedralLayout(data: PreparedMyriahedralLayout): MyriahedralLayout;
+/** A static seam-aware projection with the same panel IDs and layout as its animation.
+ * Radial/gnomonic mapping is locally approximate, not equal-area or conformal.
+ */
+export declare class MyriahedralMap extends DiscontinuousMap {
+    readonly layout: MyriahedralLayout;
+    constructor(options?: MyriahedralOptions, prepared?: MyriahedralLayout);
+}
+/** Create a reusable static projection; use map.layout.evaluate(t) for unfolding. */
+export declare function createMyriahedralMap(options?: MyriahedralOptions): MyriahedralMap;
+//# sourceMappingURL=myriahedral.d.ts.map

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

@@ -0,0 +1 @@
+{"version":3,"file":"myriahedral.d.ts","sourceRoot":"","sources":["../../src/discontinuous/myriahedral.ts"],"names":[],"mappings":"AAIA,OAAO,EAAC,gBAAgB,EAAE,KAAK,kBAAkB,EAAiB,wBAAqB;AAGvF,OAAO,EACL,YAAY,EAEZ,KAAK,mBAAmB,EAEzB,gCAA6B;AAC9B,OAAO,EAQL,KAAK,UAAU,EACf,KAAK,KAAK,EAEX,8BAA2B;AAE5B,MAAM,MAAM,eAAe,GAAG,UAAU,CAAC;AACzC,MAAM,MAAM,gBAAgB,GAAG,KAAK,CAAC;AACrC,MAAM,MAAM,iBAAiB,GACzB,WAAW,GACX,WAAW,GACX,WAAW,GACX,aAAa,GACb,SAAS,GACT,WAAW,GACX,YAAY,CAAC;AACjB,MAAM,MAAM,kBAAkB,GAAG;IAC/B,kGAAkG;IAClG,MAAM,CAAC,EAAE,iBAAiB,CAAC;IAC3B,mEAAmE;IACnE,OAAO,CAAC,EAAE,MAAM,CAAC;IACjB,mDAAmD;IACnD,KAAK,CAAC,EAAE,gBAAgB,CAAC;IACzB,uFAAuF;IACvF,YAAY,CAAC,EAAE,MAAM,CAAC;IACtB,+EAA+E;IAC/E,IAAI,CAAC,EAAE,eAAe,CAAC;IACvB,mFAAmF;IACnF,eAAe,CAAC,EAAE,MAAM,CAAC;IACzB,wFAAwF;IACxF,MAAM,CAAC,EAAE,MAAM,CAAC;IAChB,mGAAmG;IACnG,MAAM,CAAC,EAAE,CAAC,QAAQ,EAAE,kBAAkB,KAAK,OAAO,CAAC;IACnD,2FAA2F;IAC3F,cAAc,CAAC,EAAE,MAAM,CAAC;IACxB,oFAAoF;IACpF,UAAU,CAAC,EAAE,CAAC,IAAI,EAAE,QAAQ,CAAC,eAAe,CAAC,EAAE,YAAY,EAAE,MAAM,KAAK,MAAM,CAAC;CAChF,CAAC;AACF,MAAM,MAAM,eAAe,GAAG;IAC5B,QAAQ,EAAE,SAAS,CAAC,MAAM,EAAE,MAAM,CAAC,CAAC;IACpC,MAAM,EAAE,SAAS,CAAC,MAAM,EAAE,MAAM,CAAC,CAAC;IAClC,QAAQ,EAAE,kBAAkB,CAAC;IAC7B,SAAS,EAAE,MAAM,CAAC;IAClB,YAAY,CAAC,EAAE,MAAM,CAAC;CACvB,CAAC;AAEF,eAAO,MAAM,mBAAmB,EAAE,SAAS,iBAAiB,EAQ1D,CAAC;AAEH;;GAEG;AACH,wBAAgB,uBAAuB,CAAC,OAAO,GAAE,kBAAuB,GAAG,iBAAiB,CAgH3F;AACD,MAAM,MAAM,iBAAiB,GAAG,UAAU,CAAC,OAAO,YAAY,CAAC,GAAG;IAChE,MAAM,EAAE,iBAAiB,CAAC;IAC1B,MAAM,EAAE,MAAM,CAAC;IACf,4CAA4C;IAC5C,cAAc,CAAC,EAAE,MAAM,CAAC;CACzB,CAAC;AAEF,MAAM,MAAM,yBAAyB,GAAG,mBAAmB,GAAG;IAC5D,MAAM,EAAE,iBAAiB,CAAC;IAC1B,MAAM,EAAE,MAAM,CAAC;IACf,cAAc,CAAC,EAAE,MAAM,CAAC;CACzB,CAAC;AAEF,wBAAgB,wBAAwB,CAAC,IAAI,EAAE,yBAAyB,GAAG,iBAAiB,CAa3F;AAED;;GAEG;AACH,qBAAa,cAAe,SAAQ,gBAAgB;IAClD,QAAQ,CAAC,MAAM,EAAE,iBAAiB,CAAC;gBACvB,OAAO,GAAE,kBAAuB,EAAE,QAAQ,CAAC,EAAE,iBAAiB;CAK3E;AAED,qFAAqF;AACrF,wBAAgB,oBAAoB,CAAC,OAAO,GAAE,kBAAuB,GAAG,cAAc,CAErF"}

package/dist/discontinuous/myriahedral.js

@@ -0,0 +1,217 @@
+// math.gl
+// SPDX-License-Identifier: MIT
+// SPDX-FileCopyrightText: Copyright (c) vis.gl contributors
+import { DiscontinuousMap } from "./map-layout.js";
+import { coordinates, geographic, angle, RAD } from "./math.js";
+import { createGraticule } from "./myriahedral/graticule.js";
+import { createLayout, restoreLayout } from "./myriahedral/unfold.js";
+import { createSolid, center, cross, dot, faceNormal, normalize, subtract } from "./myriahedral/mesh.js";
+export const MYRIAHEDRAL_PRESETS = Object.freeze([
+    'polyconic',
+    'cylindric',
+    'azimuthal',
+    'hemispheres',
+    'sealine',
+    'coastline',
+    'polyhedral'
+]);
+/** Build an independent spherical panel net and its rigid unfolding animation.
+ * Each call owns its mesh. Geographical costs are sampled once per edge.
+ */
+export function createMyriahedralLayout(options = {}) {
+    const preset = options.preset ?? 'cylindric';
+    if (!MYRIAHEDRAL_PRESETS.includes(preset))
+        throw new RangeError('Unknown myriahedral preset');
+    const radius = options.radius ?? 1;
+    const longitude = options.centralMeridian ?? 20;
+    const step = options.landSampleStep ?? 0.25;
+    if (!Number.isFinite(radius) || radius <= 0)
+        throw new RangeError('Radius must be finite and positive');
+    if (!Number.isFinite(longitude))
+        throw new RangeError('Central meridian must be finite');
+    if (!Number.isFinite(step) || step < 0.05 || step > 10)
+        throw new RangeError('Land sample step must be within [.05, 10]');
+    if ((preset === 'sealine' || preset === 'coastline') && !options.isLand)
+        throw new Error('Sealine and coastline require an isLand callback');
+    const central = (((((longitude + 180) % 360) + 360) % 360) - 180) * RAD;
+    const mesh = options.mesh
+        ? snapshotMesh(options.mesh)
+        : preset === 'polyhedral'
+            ? createSolid(options.solid ?? 'icosahedron', options.subdivisions ?? 1)
+            : createGraticule(options.columns ?? 48);
+    // Select a stable root near the requested center (north pole for radial nets).
+    const target = geographic([
+        longitude,
+        preset === 'azimuthal' || preset === 'hemispheres' ? 90 : 0
+    ]);
+    let root = 0;
+    for (let i = 1; i < mesh.faces.length; i++)
+        if (dot(normalize(center(mesh.faces[i].map(v => mesh.vertices[v]))), target) >
+            dot(normalize(center(mesh.faces[root].map(v => mesh.vertices[v]))), target) + 1e-12)
+            root = i;
+    [mesh.faces[0], mesh.faces[root]] = [mesh.faces[root], mesh.faces[0]];
+    const samples = new Map();
+    const data = (edge) => {
+        const key = `${edge.a}:${edge.b}`;
+        let value = samples.get(key);
+        if (value)
+            return value;
+        const a = mesh.vertices[edge.a], b = mesh.vertices[edge.b];
+        const arcLength = angle(a, b);
+        let landFraction;
+        if (options.isLand) {
+            const count = Math.max(1, Math.ceil(arcLength / (step * RAD)));
+            let hits = 0;
+            for (let i = 0; i < count; i++) {
+                const t = (i + 0.5) / count;
+                const p = a.map((x, k) => x * Math.sin((1 - t) * arcLength) + b[k] * Math.sin(t * arcLength));
+                if (options.isLand(coordinates(p)))
+                    hits++;
+            }
+            landFraction = hits / count;
+        }
+        value = Object.freeze({
+            vertices: Object.freeze([edge.a, edge.b]),
+            panels: Object.freeze([...edge.faces]),
+            midpoint: Object.freeze(coordinates(normalize(center([a, b])))),
+            arcLength,
+            landFraction
+        });
+        samples.set(key, value);
+        return value;
+    };
+    const weight = (edge) => {
+        const a = mesh.vertices[edge.a], b = mesh.vertices[edge.b];
+        const d = data(edge);
+        const lat = d.midpoint[1] * RAD;
+        const lon = d.midpoint[0] * RAD;
+        const delta = Math.abs(Math.atan2(Math.sin(lon - central), Math.cos(lon - central)));
+        const parallel = Math.abs(a[2] - b[2]) < 1e-10;
+        const seam = !parallel && Math.abs(Math.abs(lon) - Math.PI) < 1e-8;
+        let result;
+        if (preset === 'polyhedral') {
+            // Geographic continuity, then short local joins. Stable input order resolves ties.
+            result = (d.landFraction ?? 0) * 10000 - d.arcLength;
+        }
+        else if (preset === 'polyconic') {
+            // Latitude strips connected along a central meridian.
+            result = (parallel ? 0 : 100) - delta - Math.abs(lat) * 0.01;
+        }
+        else {
+            // Longitude strips connected at the equator, north cap, or both polar caps.
+            const anchor = preset === 'azimuthal'
+                ? Math.PI / 2
+                : preset === 'hemispheres'
+                    ? (Math.sign(lat) * Math.PI) / 2
+                    : 0;
+            result = (parallel ? 100 : 0) - Math.abs(lat - anchor) - delta * 0.01;
+            if (preset === 'hemispheres' && parallel && Math.abs(a[2]) < 1e-10)
+                result = -100 - delta;
+            if (preset === 'sealine')
+                result += d.landFraction * 10000;
+            if (preset === 'coastline')
+                result += (1 - d.landFraction) * 10000;
+            if ((preset === 'cylindric' || preset === 'sealine' || preset === 'coastline') &&
+                !parallel &&
+                Math.abs(Math.min(a[2], b[2])) < 1e-10)
+                result = 1e6;
+        }
+        if (seam && preset !== 'polyhedral')
+            result = -1e9;
+        return options.edgeWeight ? options.edgeWeight(d, result) : result;
+    };
+    const layout = createLayout(mesh, 1, { weight });
+    const unscaled = layout.evaluate;
+    const evaluate = (amount) => unscaled(amount).map(face => face.map(p => p.map(v => v * radius)));
+    const cutLandRadians = options.isLand
+        ? layout.edges.reduce((sum, edge) => sum + (edge.hinge ? 0 : data(edge).arcLength * data(edge).landFraction), 0)
+        : undefined;
+    return { ...layout, evaluate, preset, radius, cutLandRadians };
+}
+export function restoreMyriahedralLayout(data) {
+    if (!Number.isFinite(data.radius) || data.radius <= 0)
+        throw new Error('Invalid prepared map radius');
+    const layout = restoreLayout({ ...data, mesh: snapshotMesh(data.mesh) });
+    const evaluate = (amount) => layout.evaluate(amount).map(face => face.map(p => p.map(v => v * data.radius)));
+    return {
+        ...layout,
+        preset: data.preset,
+        radius: data.radius,
+        cutLandRadians: data.cutLandRadians,
+        evaluate
+    };
+}
+/** A static seam-aware projection with the same panel IDs and layout as its animation.
+ * Radial/gnomonic mapping is locally approximate, not equal-area or conformal.
+ */
+export class MyriahedralMap extends DiscontinuousMap {
+    constructor(options = {}, prepared) {
+        const layout = prepared ?? createMyriahedralLayout(options);
+        super(regions(layout));
+        this.layout = layout;
+    }
+}
+/** Create a reusable static projection; use map.layout.evaluate(t) for unfolding. */
+export function createMyriahedralMap(options = {}) {
+    return new MyriahedralMap(options);
+}
+function regions(layout) {
+    const flat = layout.evaluate(1);
+    return layout.mesh.faces.map((face, i) => {
+        const origin = [...layout.mesh.vertices[face[0]]];
+        const x = normalize(subtract(layout.mesh.vertices[face[1]], origin));
+        const normal = faceNormal(layout.mesh, face), y = cross(normal, x);
+        const height = dot(normal, origin);
+        const fx = normalize(subtract(flat[i][1], flat[i][0]));
+        const radius = layout.radius;
+        return {
+            id: `panel-${i}`,
+            boundary: face.map(v => coordinates(layout.mesh.vertices[v])),
+            interior: coordinates(normalize(center(face.map(v => layout.mesh.vertices[v])))),
+            placement: [
+                radius * fx[0],
+                radius * fx[1],
+                -radius * fx[1],
+                radius * fx[0],
+                flat[i][0][0],
+                flat[i][0][1]
+            ],
+            project: p => {
+                const direction = geographic(p);
+                const point = direction.map(v => (v * height) / dot(normal, direction));
+                return [dot(subtract(point, origin), x), dot(subtract(point, origin), y)];
+            },
+            unproject: ([u, v]) => coordinates(normalize(origin.map((c, k) => c + u * x[k] + v * y[k])))
+        };
+    });
+}
+function snapshotMesh(input) {
+    if (!input.vertices?.length || !input.faces?.length || input.faces.length > 20000)
+        throw new Error('Expected a nonempty mesh with at most 20000 panels');
+    const mesh = {
+        vertices: input.vertices.map(v => [...v]),
+        faces: input.faces.map(f => [...f])
+    };
+    if (mesh.vertices.some(v => v.length !== 3 || !v.every(Number.isFinite) || Math.abs(Math.hypot(...v) - 1) > 1e-8))
+        throw new Error('Mesh vertices must lie on the unit sphere');
+    for (const face of mesh.faces) {
+        if (face.length < 3 ||
+            new Set(face).size !== face.length ||
+            face.some(v => !Number.isInteger(v) || v < 0 || v >= mesh.vertices.length))
+            throw new Error('Invalid panel indices');
+        const normal = faceNormal(mesh, face), origin = mesh.vertices[face[0]];
+        if (!normal.every(Number.isFinite) ||
+            dot(normal, origin) <= 1e-8 ||
+            face.some(v => Math.abs(dot(normal, subtract(mesh.vertices[v], origin))) > 1e-8))
+            throw new Error('Panels must be planar and outward-wound');
+        for (let i = 0; i < face.length; i++) {
+            const a = mesh.vertices[face[i]], b = mesh.vertices[face[(i + 1) % face.length]];
+            if (angle(a, b) < 1e-10 ||
+                angle(a, b) >= Math.PI - 1e-8 ||
+                face.some(v => dot(cross(a, b), mesh.vertices[v]) < -1e-8))
+                throw new Error('Panels must be convex with nondegenerate minor-arc edges');
+        }
+    }
+    return mesh;
+}
+//# sourceMappingURL=myriahedral.js.map

package/dist/discontinuous/myriahedral.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"myriahedral.js","sourceRoot":"","sources":["../../src/discontinuous/myriahedral.ts"],"names":[],"mappings":"AAAA,UAAU;AACV,+BAA+B;AAC/B,4DAA4D;AAE5D,OAAO,EAAC,gBAAgB,EAA0C,wBAAqB;AACvF,OAAO,EAAC,WAAW,EAAE,UAAU,EAAE,KAAK,EAAE,GAAG,EAAC,kBAAe;AAC3D,OAAO,EAAC,eAAe,EAAC,mCAAgC;AACxD,OAAO,EACL,YAAY,EACZ,aAAa,EAGd,gCAA6B;AAC9B,OAAO,EACL,WAAW,EACX,MAAM,EACN,KAAK,EACL,GAAG,EACH,UAAU,EACV,SAAS,EACT,QAAQ,EAIT,8BAA2B;AA0C5B,MAAM,CAAC,MAAM,mBAAmB,GAAiC,MAAM,CAAC,MAAM,CAAC;IAC7E,WAAW;IACX,WAAW;IACX,WAAW;IACX,aAAa;IACb,SAAS;IACT,WAAW;IACX,YAAY;CACb,CAAC,CAAC;AAEH;;GAEG;AACH,MAAM,UAAU,uBAAuB,CAAC,UAA8B,EAAE;IACtE,MAAM,MAAM,GAAG,OAAO,CAAC,MAAM,IAAI,WAAW,CAAC;IAC7C,IAAI,CAAC,mBAAmB,CAAC,QAAQ,CAAC,MAAM,CAAC;QAAE,MAAM,IAAI,UAAU,CAAC,4BAA4B,CAAC,CAAC;IAC9F,MAAM,MAAM,GAAG,OAAO,CAAC,MAAM,IAAI,CAAC,CAAC;IACnC,MAAM,SAAS,GAAG,OAAO,CAAC,eAAe,IAAI,EAAE,CAAC;IAChD,MAAM,IAAI,GAAG,OAAO,CAAC,cAAc,IAAI,IAAI,CAAC;IAC5C,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,MAAM,CAAC,IAAI,MAAM,IAAI,CAAC;QACzC,MAAM,IAAI,UAAU,CAAC,oCAAoC,CAAC,CAAC;IAC7D,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,SAAS,CAAC;QAAE,MAAM,IAAI,UAAU,CAAC,iCAAiC,CAAC,CAAC;IACzF,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,IAAI,CAAC,IAAI,IAAI,GAAG,IAAI,IAAI,IAAI,GAAG,EAAE;QACpD,MAAM,IAAI,UAAU,CAAC,2CAA2C,CAAC,CAAC;IACpE,IAAI,CAAC,MAAM,KAAK,SAAS,IAAI,MAAM,KAAK,WAAW,CAAC,IAAI,CAAC,OAAO,CAAC,MAAM;QACrE,MAAM,IAAI,KAAK,CAAC,kDAAkD,CAAC,CAAC;IACtE,MAAM,OAAO,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,SAAS,GAAG,GAAG,CAAC,GAAG,GAAG,CAAC,GAAG,GAAG,CAAC,GAAG,GAAG,CAAC,GAAG,GAAG,CAAC,GAAG,GAAG,CAAC;IACxE,MAAM,IAAI,GAAG,OAAO,CAAC,IAAI;QACvB,CAAC,CAAC,YAAY,CAAC,OAAO,CAAC,IAAI,CAAC;QAC5B,CAAC,CAAC,MAAM,KAAK,YAAY;YACvB,CAAC,CAAC,WAAW,CAAC,OAAO,CAAC,KAAK,IAAI,aAAa,EAAE,OAAO,CAAC,YAAY,IAAI,CAAC,CAAC;YACxE,CAAC,CAAC,eAAe,CAAC,OAAO,CAAC,OAAO,IAAI,EAAE,CAAC,CAAC;IAC7C,+EAA+E;IAC/E,MAAM,MAAM,GAAG,UAAU,CAAC;QACxB,SAAS;QACT,MAAM,KAAK,WAAW,IAAI,MAAM,KAAK,aAAa,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;KAC5D,CAAC,CAAC;IACH,IAAI,IAAI,GAAG,CAAC,CAAC;IACb,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,MAAM,EAAE,CAAC,EAAE;QACxC,IACE,GAAG,CAAC,SAAS,CAAC,MAAM,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,MAAM,CAAC;YACxE,GAAG,CAAC,SAAS,CAAC,MAAM,CAAC,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,MAAM,CAAC,GAAG,KAAK;YAEnF,IAAI,GAAG,CAAC,CAAC;IACb,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,EAAE,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC;IACtE,MAAM,OAAO,GAAG,IAAI,GAAG,EAA2B,CAAC;IACnD,MAAM,IAAI,GAAG,CAAC,IAAc,EAAmB,EAAE;QAC/C,MAAM,GAAG,GAAG,GAAG,IAAI,CAAC,CAAC,IAAI,IAAI,CAAC,CAAC,EAAE,CAAC;QAClC,IAAI,KAAK,GAAG,OAAO,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;QAC7B,IAAI,KAAK;YAAE,OAAO,KAAK,CAAC;QACxB,MAAM,CAAC,GAAG,IAAI,CAAC,QAAQ,CAAC,IAAI,CAAC,CAAC,CAAC,EAC7B,CAAC,GAAG,IAAI,CAAC,QAAQ,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;QAC5B,MAAM,SAAS,GAAG,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;QAC9B,IAAI,YAAgC,CAAC;QACrC,IAAI,OAAO,CAAC,MAAM,EAAE,CAAC;YACnB,MAAM,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,IAAI,CAAC,IAAI,CAAC,SAAS,GAAG,CAAC,IAAI,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC;YAC/D,IAAI,IAAI,GAAG,CAAC,CAAC;YACb,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,KAAK,EAAE,CAAC,EAAE,EAAE,CAAC;gBAC/B,MAAM,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,CAAC,GAAG,KAAK,CAAC;gBAC5B,MAAM,CAAC,GAAG,CAAC,CAAC,GAAG,CACb,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,SAAS,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,SAAS,CAAC,CACrE,CAAC;gBACV,IAAI,OAAO,CAAC,MAAM,CAAC,WAAW,CAAC,CAAC,CAAC,CAAC;oBAAE,IAAI,EAAE,CAAC;YAC7C,CAAC;YACD,YAAY,GAAG,IAAI,GAAG,KAAK,CAAC;QAC9B,CAAC;QACD,KAAK,GAAG,MAAM,CAAC,MAAM,CAAC;YACpB,QAAQ,EAAE,MAAM,CAAC,MAAM,CAAC,CAAC,IAAI,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,CAAC,CAA8B;YACtE,MAAM,EAAE,MAAM,CAAC,MAAM,CAAC,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,CAA8B;YACnE,QAAQ,EAAE,MAAM,CAAC,MAAM,CAAC,WAAW,CAAC,SAAS,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;YAC/D,SAAS;YACT,YAAY;SACb,CAAC,CAAC;QACH,OAAO,CAAC,GAAG,CAAC,GAAG,EAAE,KAAK,CAAC,CAAC;QACxB,OAAO,KAAK,CAAC;IACf,CAAC,CAAC;IACF,MAAM,MAAM,GAAG,CAAC,IAAc,EAAE,EAAE;QAChC,MAAM,CAAC,GAAG,IAAI,CAAC,QAAQ,CAAC,IAAI,CAAC,CAAC,CAAC,EAC7B,CAAC,GAAG,IAAI,CAAC,QAAQ,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;QAC5B,MAAM,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC;QACrB,MAAM,GAAG,GAAG,CAAC,CAAC,QAAQ,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC;QAChC,MAAM,GAAG,GAAG,CAAC,CAAC,QAAQ,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC;QAChC,MAAM,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,GAAG,CAAC,GAAG,GAAG,OAAO,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,GAAG,GAAG,OAAO,CAAC,CAAC,CAAC,CAAC;QACrF,MAAM,QAAQ,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC;QAC/C,MAAM,IAAI,GAAG,CAAC,QAAQ,IAAI,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC;QACnE,IAAI,MAAc,CAAC;QACnB,IAAI,MAAM,KAAK,YAAY,EAAE,CAAC;YAC5B,mFAAmF;YACnF,MAAM,GAAG,CAAC,CAAC,CAAC,YAAY,IAAI,CAAC,CAAC,GAAG,KAAK,GAAG,CAAC,CAAC,SAAS,CAAC;QACvD,CAAC;aAAM,IAAI,MAAM,KAAK,WAAW,EAAE,CAAC;YAClC,sDAAsD;YACtD,MAAM,GAAG,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC;QAC/D,CAAC;aAAM,CAAC;YACN,4EAA4E;YAC5E,MAAM,MAAM,GACV,MAAM,KAAK,WAAW;gBACpB,CAAC,CAAC,IAAI,CAAC,EAAE,GAAG,CAAC;gBACb,CAAC,CAAC,MAAM,KAAK,aAAa;oBACxB,CAAC,CAAC,CAAC,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC,GAAG,CAAC;oBAChC,CAAC,CAAC,CAAC,CAAC;YACV,MAAM,GAAG,CAAC,QAAQ,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,GAAG,MAAM,CAAC,GAAG,KAAK,GAAG,IAAI,CAAC;YACtE,IAAI,MAAM,KAAK,aAAa,IAAI,QAAQ,IAAI,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK;gBAAE,MAAM,GAAG,CAAC,GAAG,GAAG,KAAK,CAAC;YAC1F,IAAI,MAAM,KAAK,SAAS;gBAAE,MAAM,IAAI,CAAC,CAAC,YAAa,GAAG,KAAK,CAAC;YAC5D,IAAI,MAAM,KAAK,WAAW;gBAAE,MAAM,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,YAAa,CAAC,GAAG,KAAK,CAAC;YACpE,IACE,CAAC,MAAM,KAAK,WAAW,IAAI,MAAM,KAAK,SAAS,IAAI,MAAM,KAAK,WAAW,CAAC;gBAC1E,CAAC,QAAQ;gBACT,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK;gBAEtC,MAAM,GAAG,GAAG,CAAC;QACjB,CAAC;QACD,IAAI,IAAI,IAAI,MAAM,KAAK,YAAY;YAAE,MAAM,GAAG,CAAC,GAAG,CAAC;QACnD,OAAO,OAAO,CAAC,UAAU,CAAC,CAAC,CAAC,OAAO,CAAC,UAAU,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC,CAAC,CAAC,MAAM,CAAC;IACrE,CAAC,CAAC;IACF,MAAM,MAAM,GAAG,YAAY,CAAC,IAAI,EAAE,CAAC,EAAE,EAAC,MAAM,EAAC,CAAC,CAAC;IAC/C,MAAM,QAAQ,GAAG,MAAM,CAAC,QAAQ,CAAC;IACjC,MAAM,QAAQ,GAAG,CAAC,MAAc,EAAE,EAAE,CAClC,QAAQ,CAAC,MAAM,CAAC,CAAC,GAAG,CAAC,IAAI,CAAC,EAAE,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,MAAM,CAAS,CAAC,CAAC,CAAC;IAC9E,MAAM,cAAc,GAAG,OAAO,CAAC,MAAM;QACnC,CAAC,CAAC,MAAM,CAAC,KAAK,CAAC,MAAM,CACjB,CAAC,GAAG,EAAE,IAAI,EAAE,EAAE,CAAC,GAAG,GAAG,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,SAAS,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,YAAa,CAAC,EACvF,CAAC,CACF;QACH,CAAC,CAAC,SAAS,CAAC;IACd,OAAO,EAAC,GAAG,MAAM,EAAE,QAAQ,EAAE,MAAM,EAAE,MAAM,EAAE,cAAc,EAAC,CAAC;AAC/D,CAAC;AAcD,MAAM,UAAU,wBAAwB,CAAC,IAA+B;IACtE,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,IAAI,CAAC,MAAM,CAAC,IAAI,IAAI,CAAC,MAAM,IAAI,CAAC;QACnD,MAAM,IAAI,KAAK,CAAC,6BAA6B,CAAC,CAAC;IACjD,MAAM,MAAM,GAAG,aAAa,CAAC,EAAC,GAAG,IAAI,EAAE,IAAI,EAAE,YAAY,CAAC,IAAI,CAAC,IAAI,CAAC,EAAC,CAAC,CAAC;IACvE,MAAM,QAAQ,GAAG,CAAC,MAAc,EAAE,EAAE,CAClC,MAAM,CAAC,QAAQ,CAAC,MAAM,CAAC,CAAC,GAAG,CAAC,IAAI,CAAC,EAAE,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,IAAI,CAAC,MAAM,CAAS,CAAC,CAAC,CAAC;IAC1F,OAAO;QACL,GAAG,MAAM;QACT,MAAM,EAAE,IAAI,CAAC,MAAM;QACnB,MAAM,EAAE,IAAI,CAAC,MAAM;QACnB,cAAc,EAAE,IAAI,CAAC,cAAc;QACnC,QAAQ;KACT,CAAC;AACJ,CAAC;AAED;;GAEG;AACH,MAAM,OAAO,cAAe,SAAQ,gBAAgB;IAElD,YAAY,UAA8B,EAAE,EAAE,QAA4B;QACxE,MAAM,MAAM,GAAG,QAAQ,IAAI,uBAAuB,CAAC,OAAO,CAAC,CAAC;QAC5D,KAAK,CAAC,OAAO,CAAC,MAAM,CAAC,CAAC,CAAC;QACvB,IAAI,CAAC,MAAM,GAAG,MAAM,CAAC;IACvB,CAAC;CACF;AAED,qFAAqF;AACrF,MAAM,UAAU,oBAAoB,CAAC,UAA8B,EAAE;IACnE,OAAO,IAAI,cAAc,CAAC,OAAO,CAAC,CAAC;AACrC,CAAC;AAED,SAAS,OAAO,CAAC,MAAyB;IACxC,MAAM,IAAI,GAAG,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC;IAChC,OAAO,MAAM,CAAC,IAAI,CAAC,KAAK,CAAC,GAAG,CAAC,CAAC,IAAI,EAAE,CAAC,EAAE,EAAE;QACvC,MAAM,MAAM,GAAG,CAAC,GAAG,MAAM,CAAC,IAAI,CAAC,QAAQ,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAS,CAAC;QAC1D,MAAM,CAAC,GAAG,SAAS,CAAC,QAAQ,CAAC,MAAM,CAAC,IAAI,CAAC,QAAQ,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC,CAAC;QACrE,MAAM,MAAM,GAAG,UAAU,CAAC,MAAM,CAAC,IAAI,EAAE,IAAI,CAAC,EAC1C,CAAC,GAAG,KAAK,CAAC,MAAM,EAAE,CAAC,CAAC,CAAC;QACvB,MAAM,MAAM,GAAG,GAAG,CAAC,MAAM,EAAE,MAAM,CAAC,CAAC;QACnC,MAAM,EAAE,GAAG,SAAS,CAAC,QAAQ,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QACvD,MAAM,MAAM,GAAG,MAAM,CAAC,MAAM,CAAC;QAC7B,OAAO;YACL,EAAE,EAAE,SAAS,CAAC,EAAE;YAChB,QAAQ,EAAE,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,WAAW,CAAC,MAAM,CAAC,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC;YAC7D,QAAQ,EAAE,WAAW,CAAC,SAAS,CAAC,MAAM,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,MAAM,CAAC,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;YAChF,SAAS,EAAE;gBACT,MAAM,GAAG,EAAE,CAAC,CAAC,CAAC;gBACd,MAAM,GAAG,EAAE,CAAC,CAAC,CAAC;gBACd,CAAC,MAAM,GAAG,EAAE,CAAC,CAAC,CAAC;gBACf,MAAM,GAAG,EAAE,CAAC,CAAC,CAAC;gBACd,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;gBACb,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;aACd;YACD,OAAO,EAAE,CAAC,CAAC,EAAE;gBACX,MAAM,SAAS,GAAG,UAAU,CAAC,CAAC,CAAC,CAAC;gBAChC,MAAM,KAAK,GAAG,SAAS,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC,GAAG,GAAG,CAAC,MAAM,EAAE,SAAS,CAAC,CAAS,CAAC;gBAChF,OAAO,CAAC,GAAG,CAAC,QAAQ,CAAC,KAAK,EAAE,MAAM,CAAC,EAAE,CAAC,CAAC,EAAE,GAAG,CAAC,QAAQ,CAAC,KAAK,EAAE,MAAM,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;YAC5E,CAAC;YACD,SAAS,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,EAAE,CACpB,WAAW,CAAC,SAAS,CAAC,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAS,CAAC,CAAC;SAChF,CAAC;IACJ,CAAC,CAAC,CAAC;AACL,CAAC;AAED,SAAS,YAAY,CAAC,KAAsB;IAC1C,IAAI,CAAC,KAAK,CAAC,QAAQ,EAAE,MAAM,IAAI,CAAC,KAAK,CAAC,KAAK,EAAE,MAAM,IAAI,KAAK,CAAC,KAAK,CAAC,MAAM,GAAG,KAAK;QAC/E,MAAM,IAAI,KAAK,CAAC,oDAAoD,CAAC,CAAC;IACxE,MAAM,IAAI,GAAG;QACX,QAAQ,EAAE,KAAK,CAAC,QAAQ,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,CAAS,CAAC;QACjD,KAAK,EAAE,KAAK,CAAC,KAAK,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;KACpC,CAAC;IACF,IACE,IAAI,CAAC,QAAQ,CAAC,IAAI,CAChB,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,MAAM,KAAK,CAAC,IAAI,CAAC,CAAC,CAAC,KAAK,CAAC,MAAM,CAAC,QAAQ,CAAC,IAAI,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,KAAK,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAC1F;QAED,MAAM,IAAI,KAAK,CAAC,2CAA2C,CAAC,CAAC;IAC/D,KAAK,MAAM,IAAI,IAAI,IAAI,CAAC,KAAK,EAAE,CAAC;QAC9B,IACE,IAAI,CAAC,MAAM,GAAG,CAAC;YACf,IAAI,GAAG,CAAC,IAAI,CAAC,CAAC,IAAI,KAAK,IAAI,CAAC,MAAM;YAClC,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,MAAM,CAAC,SAAS,CAAC,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,IAAI,IAAI,CAAC,QAAQ,CAAC,MAAM,CAAC;YAE1E,MAAM,IAAI,KAAK,CAAC,uBAAuB,CAAC,CAAC;QAC3C,MAAM,MAAM,GAAG,UAAU,CAAC,IAAI,EAAE,IAAI,CAAC,EACnC,MAAM,GAAG,IAAI,CAAC,QAAQ,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC;QAClC,IACE,CAAC,MAAM,CAAC,KAAK,CAAC,MAAM,CAAC,QAAQ,CAAC;YAC9B,GAAG,CAAC,MAAM,EAAE,MAAM,CAAC,IAAI,IAAI;YAC3B,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,MAAM,EAAE,QAAQ,CAAC,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC;YAEhF,MAAM,IAAI,KAAK,CAAC,yCAAyC,CAAC,CAAC;QAC7D,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;YACrC,MAAM,CAAC,GAAG,IAAI,CAAC,QAAQ,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,EAC9B,CAAC,GAAG,IAAI,CAAC,QAAQ,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,MAAM,CAAC,CAAC,CAAC;YACjD,IACE,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,KAAK;gBACnB,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC,IAAI,IAAI,CAAC,EAAE,GAAG,IAAI;gBAC7B,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,IAAI,CAAC;gBAE1D,MAAM,IAAI,KAAK,CAAC,0DAA0D,CAAC,CAAC;QAChF,CAAC;IACH,CAAC;IACD,OAAO,IAAI,CAAC;AACd,CAAC"}

package/dist/discontinuous/petroff.d.ts

@@ -0,0 +1,8 @@
+import { DiscontinuousMap } from "./map-layout.js";
+import { type Position } from "./math.js";
+import type { QuincuncialMapOptions } from "./quincuncial-star.js";
+/** Published polar quincuncial aspect; local longitude is measured from a quadrant's center. */
+export declare function petroffForward(lambda: number, phi: number, quadrant: number, south: boolean): Position;
+/** Petroff's 2021 equal-area square, north pole at center and south pole at corners. */
+export declare function createPetroffMap(options?: QuincuncialMapOptions): DiscontinuousMap;
+//# sourceMappingURL=petroff.d.ts.map

package/dist/discontinuous/petroff.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"petroff.d.ts","sourceRoot":"","sources":["../../src/discontinuous/petroff.ts"],"names":[],"mappings":"AAMA,OAAO,EAAC,gBAAgB,EAAiB,wBAAqB;AAC9D,OAAO,EAAiD,KAAK,QAAQ,EAAC,kBAAe;AACrF,OAAO,KAAK,EAAC,qBAAqB,EAAC,8BAA2B;AA0C9D,gGAAgG;AAChG,wBAAgB,cAAc,CAC5B,MAAM,EAAE,MAAM,EACd,GAAG,EAAE,MAAM,EACX,QAAQ,EAAE,MAAM,EAChB,KAAK,EAAE,OAAO,GACb,QAAQ,CA+CV;AAuCD,wFAAwF;AACxF,wBAAgB,gBAAgB,CAAC,OAAO,GAAE,qBAA0B,oBAkCnE"}