hyper-scatter 0.1.0

package/dist-lib/core/dataset.js

@@ -0,0 +1,152 @@
+/**
+ * Deterministic dataset generation for testing.
+ */
+import { SeededRNG } from './rng.js';
+// For very large-N benchmarks we want determinism but also *reasonable* setup time.
+// Dataset generation is not the performance target of this project; rendering and
+// interaction are. For N in the multi-millions, avoid per-point trig/log-heavy
+// sampling where possible.
+const FAST_PATH_THRESHOLD_N = 2_000_000;
+function gaussianFastApprox(rng, mean = 0, std = 1) {
+    // Irwin–Hall approximation: sum of uniforms approximates a normal.
+    // Using 6 uniforms is a good speed/quality tradeoff.
+    // Output has mean 0, variance 0.5 -> std ~0.707; we rescale to unit std.
+    const u = rng.random() + rng.random() + rng.random() +
+        rng.random() + rng.random() + rng.random();
+    const z = (u - 3.0) * 1.4142135623730951; // sqrt(2)
+    return mean + z * std;
+}
+function fillPoincareClusterCentersRandom(rng, labelCount, outCx, outCy, outJitterStd) {
+    // Random cluster centers inside disk with area-uniform radius.
+    // We still scale jitter by (1-|center|^2)/2 to keep clusters visually
+    // comparable in hyperbolic units.
+    const hyperbolicStd = 0.30;
+    for (let label = 0; label < labelCount; label++) {
+        const r = Math.sqrt(rng.random()) * 0.92;
+        const theta = rng.range(0, 2 * Math.PI);
+        const cx = r * Math.cos(theta);
+        const cy = r * Math.sin(theta);
+        outCx[label] = cx;
+        outCy[label] = cy;
+        const r2 = cx * cx + cy * cy;
+        outJitterStd[label] = hyperbolicStd * 0.5 * (1 - r2);
+    }
+}
+/**
+ * Generate a deterministic dataset.
+ */
+export function generateDataset(config) {
+    const rng = new SeededRNG(config.seed);
+    const positions = new Float32Array(config.n * 2);
+    const labels = new Uint16Array(config.n);
+    const fast = config.n >= FAST_PATH_THRESHOLD_N;
+    const distribution = config.distribution ?? 'default';
+    if (config.geometry === 'euclidean') {
+        // Generate clustered Euclidean data
+        const clusterCenters = [];
+        for (let i = 0; i < config.labelCount; i++) {
+            clusterCenters.push({
+                x: rng.range(-1, 1),
+                y: rng.range(-1, 1),
+            });
+        }
+        for (let i = 0; i < config.n; i++) {
+            const label = rng.int(config.labelCount);
+            const center = clusterCenters[label];
+            // For huge N, avoid Box–Muller (log/cos/sqrt) per point.
+            // The approximation keeps overall visual structure and determinism.
+            const noiseX = fast ? gaussianFastApprox(rng, 0, 0.15) : rng.gaussian(0, 0.15);
+            const noiseY = fast ? gaussianFastApprox(rng, 0, 0.15) : rng.gaussian(0, 0.15);
+            positions[i * 2] = center.x + noiseX;
+            positions[i * 2 + 1] = center.y + noiseY;
+            labels[i] = label;
+        }
+    }
+    else {
+        // Generate Poincare disk data
+        // Points must be inside the unit disk (|z| < 1)
+        if (distribution === 'clustered') {
+            const cx = new Float32Array(config.labelCount);
+            const cy = new Float32Array(config.labelCount);
+            const jitterStd = new Float32Array(config.labelCount);
+            fillPoincareClusterCentersRandom(rng, config.labelCount, cx, cy, jitterStd);
+            for (let i = 0; i < config.n; i++) {
+                const label = rng.int(config.labelCount);
+                const j = jitterStd[label];
+                const noiseX = fast ? gaussianFastApprox(rng, 0, j) : rng.gaussian(0, j);
+                const noiseY = fast ? gaussianFastApprox(rng, 0, j) : rng.gaussian(0, j);
+                let x = cx[label] + noiseX;
+                let y = cy[label] + noiseY;
+                const rSq = x * x + y * y;
+                if (rSq >= 1) {
+                    const invR = 0.999 / Math.sqrt(rSq);
+                    x *= invR;
+                    y *= invR;
+                }
+                positions[i * 2] = x;
+                positions[i * 2 + 1] = y;
+                labels[i] = label;
+            }
+            return {
+                n: config.n,
+                positions,
+                labels,
+                geometry: config.geometry,
+            };
+        }
+        // For huge N, avoid sin/cos per point by sampling angle from a lookup table.
+        // Deterministic (seeded) and close enough for stress/perf testing.
+        const tableSize = fast ? 4096 : 0;
+        const cosTable = tableSize ? new Float32Array(tableSize) : null;
+        const sinTable = tableSize ? new Float32Array(tableSize) : null;
+        if (tableSize && cosTable && sinTable) {
+            for (let t = 0; t < tableSize; t++) {
+                const theta = (t / tableSize) * 2 * Math.PI;
+                cosTable[t] = Math.cos(theta);
+                sinTable[t] = Math.sin(theta);
+            }
+        }
+        for (let i = 0; i < config.n; i++) {
+            const label = rng.int(config.labelCount);
+            let r;
+            if (config.boundaryStress) {
+                // Stress test: many points near boundary
+                r = rng.range(0.9, 0.999);
+            }
+            else {
+                // Normal distribution: use sqrt for uniform area distribution
+                // Then compress slightly to avoid exact boundary
+                r = Math.sqrt(rng.random()) * 0.95;
+            }
+            if (tableSize && cosTable && sinTable) {
+                const t = rng.int(tableSize);
+                positions[i * 2] = r * cosTable[t];
+                positions[i * 2 + 1] = r * sinTable[t];
+            }
+            else {
+                const theta = rng.range(0, 2 * Math.PI);
+                positions[i * 2] = r * Math.cos(theta);
+                positions[i * 2 + 1] = r * Math.sin(theta);
+            }
+            labels[i] = label;
+        }
+    }
+    return {
+        n: config.n,
+        positions,
+        labels,
+        geometry: config.geometry,
+    };
+}
+/**
+ * Standard test datasets.
+ */
+export const STANDARD_DATASETS = {
+    euclidean_10k: { seed: 42, n: 10_000, labelCount: 10, geometry: 'euclidean' },
+    euclidean_100k: { seed: 42, n: 100_000, labelCount: 10, geometry: 'euclidean' },
+    euclidean_1m: { seed: 42, n: 1_000_000, labelCount: 10, geometry: 'euclidean' },
+    poincare_10k: { seed: 42, n: 10_000, labelCount: 10, geometry: 'poincare' },
+    poincare_100k: { seed: 42, n: 100_000, labelCount: 10, geometry: 'poincare' },
+    poincare_boundary_stress: { seed: 42, n: 100_000, labelCount: 10, geometry: 'poincare', boundaryStress: true },
+};
+//# sourceMappingURL=dataset.js.map

package/dist-lib/core/dataset.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"dataset.js","sourceRoot":"","sources":["../../src/core/dataset.ts"],"names":[],"mappings":"AAAA;;GAEG;AAGH,OAAO,EAAE,SAAS,EAAE,MAAM,UAAU,CAAC;AAiBrC,oFAAoF;AACpF,kFAAkF;AAClF,+EAA+E;AAC/E,2BAA2B;AAC3B,MAAM,qBAAqB,GAAG,SAAS,CAAC;AAExC,SAAS,kBAAkB,CAAC,GAAc,EAAE,IAAI,GAAG,CAAC,EAAE,GAAG,GAAG,CAAC;IAC3D,mEAAmE;IACnE,qDAAqD;IACrD,yEAAyE;IACzE,MAAM,CAAC,GACL,GAAG,CAAC,MAAM,EAAE,GAAG,GAAG,CAAC,MAAM,EAAE,GAAG,GAAG,CAAC,MAAM,EAAE;QAC1C,GAAG,CAAC,MAAM,EAAE,GAAG,GAAG,CAAC,MAAM,EAAE,GAAG,GAAG,CAAC,MAAM,EAAE,CAAC;IAC7C,MAAM,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,CAAC,GAAG,kBAAkB,CAAC,CAAC,UAAU;IACpD,OAAO,IAAI,GAAG,CAAC,GAAG,GAAG,CAAC;AACxB,CAAC;AAED,SAAS,gCAAgC,CACvC,GAAc,EACd,UAAkB,EAClB,KAAmB,EACnB,KAAmB,EACnB,YAA0B;IAE1B,+DAA+D;IAC/D,sEAAsE;IACtE,kCAAkC;IAClC,MAAM,aAAa,GAAG,IAAI,CAAC;IAE3B,KAAK,IAAI,KAAK,GAAG,CAAC,EAAE,KAAK,GAAG,UAAU,EAAE,KAAK,EAAE,EAAE,CAAC;QAChD,MAAM,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,MAAM,EAAE,CAAC,GAAG,IAAI,CAAC;QACzC,MAAM,KAAK,GAAG,GAAG,CAAC,KAAK,CAAC,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC,CAAC;QACxC,MAAM,EAAE,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC;QAC/B,MAAM,EAAE,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC;QAC/B,KAAK,CAAC,KAAK,CAAC,GAAG,EAAE,CAAC;QAClB,KAAK,CAAC,KAAK,CAAC,GAAG,EAAE,CAAC;QAClB,MAAM,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC;QAC7B,YAAY,CAAC,KAAK,CAAC,GAAG,aAAa,GAAG,GAAG,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC;IACvD,CAAC;AACH,CAAC;AAED;;GAEG;AACH,MAAM,UAAU,eAAe,CAAC,MAAqB;IACnD,MAAM,GAAG,GAAG,IAAI,SAAS,CAAC,MAAM,CAAC,IAAI,CAAC,CAAC;IACvC,MAAM,SAAS,GAAG,IAAI,YAAY,CAAC,MAAM,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;IACjD,MAAM,MAAM,GAAG,IAAI,WAAW,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC;IAEzC,MAAM,IAAI,GAAG,MAAM,CAAC,CAAC,IAAI,qBAAqB,CAAC;IAE/C,MAAM,YAAY,GAAwB,MAAM,CAAC,YAAY,IAAI,SAAS,CAAC;IAE3E,IAAI,MAAM,CAAC,QAAQ,KAAK,WAAW,EAAE,CAAC;QACpC,oCAAoC;QACpC,MAAM,cAAc,GAAoC,EAAE,CAAC;QAC3D,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,MAAM,CAAC,UAAU,EAAE,CAAC,EAAE,EAAE,CAAC;YAC3C,cAAc,CAAC,IAAI,CAAC;gBAClB,CAAC,EAAE,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC;gBACnB,CAAC,EAAE,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC;aACpB,CAAC,CAAC;QACL,CAAC;QAED,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,MAAM,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;YAClC,MAAM,KAAK,GAAG,GAAG,CAAC,GAAG,CAAC,MAAM,CAAC,UAAU,CAAC,CAAC;YACzC,MAAM,MAAM,GAAG,cAAc,CAAC,KAAK,CAAC,CAAC;YAErC,yDAAyD;YACzD,oEAAoE;YACpE,MAAM,MAAM,GAAG,IAAI,CAAC,CAAC,CAAC,kBAAkB,CAAC,GAAG,EAAE,CAAC,EAAE,IAAI,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,QAAQ,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC;YAC/E,MAAM,MAAM,GAAG,IAAI,CAAC,CAAC,CAAC,kBAAkB,CAAC,GAAG,EAAE,CAAC,EAAE,IAAI,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,QAAQ,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC;YAC/E,SAAS,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,GAAG,MAAM,CAAC;YACrC,SAAS,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,GAAG,MAAM,CAAC;YACzC,MAAM,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC;QACpB,CAAC;IACH,CAAC;SAAM,CAAC;QACN,8BAA8B;QAC9B,gDAAgD;QAEhD,IAAI,YAAY,KAAK,WAAW,EAAE,CAAC;YACjC,MAAM,EAAE,GAAG,IAAI,YAAY,CAAC,MAAM,CAAC,UAAU,CAAC,CAAC;YAC/C,MAAM,EAAE,GAAG,IAAI,YAAY,CAAC,MAAM,CAAC,UAAU,CAAC,CAAC;YAC/C,MAAM,SAAS,GAAG,IAAI,YAAY,CAAC,MAAM,CAAC,UAAU,CAAC,CAAC;YACtD,gCAAgC,CAAC,GAAG,EAAE,MAAM,CAAC,UAAU,EAAE,EAAE,EAAE,EAAE,EAAE,SAAS,CAAC,CAAC;YAE5E,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,MAAM,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;gBAClC,MAAM,KAAK,GAAG,GAAG,CAAC,GAAG,CAAC,MAAM,CAAC,UAAU,CAAC,CAAC;gBACzC,MAAM,CAAC,GAAG,SAAS,CAAC,KAAK,CAAC,CAAC;gBAC3B,MAAM,MAAM,GAAG,IAAI,CAAC,CAAC,CAAC,kBAAkB,CAAC,GAAG,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,QAAQ,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;gBACzE,MAAM,MAAM,GAAG,IAAI,CAAC,CAAC,CAAC,kBAAkB,CAAC,GAAG,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,QAAQ,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;gBAEzE,IAAI,CAAC,GAAG,EAAE,CAAC,KAAK,CAAC,GAAG,MAAM,CAAC;gBAC3B,IAAI,CAAC,GAAG,EAAE,CAAC,KAAK,CAAC,GAAG,MAAM,CAAC;gBAC3B,MAAM,GAAG,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;gBAC1B,IAAI,GAAG,IAAI,CAAC,EAAE,CAAC;oBACb,MAAM,IAAI,GAAG,KAAK,GAAG,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC;oBACpC,CAAC,IAAI,IAAI,CAAC;oBACV,CAAC,IAAI,IAAI,CAAC;gBACZ,CAAC;gBAED,SAAS,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC;gBACrB,SAAS,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC;gBACzB,MAAM,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC;YACpB,CAAC;YAED,OAAO;gBACL,CAAC,EAAE,MAAM,CAAC,CAAC;gBACX,SAAS;gBACT,MAAM;gBACN,QAAQ,EAAE,MAAM,CAAC,QAAQ;aAC1B,CAAC;QACJ,CAAC;QAED,6EAA6E;QAC7E,mEAAmE;QACnE,MAAM,SAAS,GAAG,IAAI,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC;QAClC,MAAM,QAAQ,GAAG,SAAS,CAAC,CAAC,CAAC,IAAI,YAAY,CAAC,SAAS,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC;QAChE,MAAM,QAAQ,GAAG,SAAS,CAAC,CAAC,CAAC,IAAI,YAAY,CAAC,SAAS,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC;QAChE,IAAI,SAAS,IAAI,QAAQ,IAAI,QAAQ,EAAE,CAAC;YACtC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,SAAS,EAAE,CAAC,EAAE,EAAE,CAAC;gBACnC,MAAM,KAAK,GAAG,CAAC,CAAC,GAAG,SAAS,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC;gBAC5C,QAAQ,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC;gBAC9B,QAAQ,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC;YAChC,CAAC;QACH,CAAC;QAED,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,MAAM,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;YAClC,MAAM,KAAK,GAAG,GAAG,CAAC,GAAG,CAAC,MAAM,CAAC,UAAU,CAAC,CAAC;YACzC,IAAI,CAAS,CAAC;YAEd,IAAI,MAAM,CAAC,cAAc,EAAE,CAAC;gBAC1B,yCAAyC;gBACzC,CAAC,GAAG,GAAG,CAAC,KAAK,CAAC,GAAG,EAAE,KAAK,CAAC,CAAC;YAC5B,CAAC;iBAAM,CAAC;gBACN,8DAA8D;gBAC9D,iDAAiD;gBACjD,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,MAAM,EAAE,CAAC,GAAG,IAAI,CAAC;YACrC,CAAC;YAED,IAAI,SAAS,IAAI,QAAQ,IAAI,QAAQ,EAAE,CAAC;gBACtC,MAAM,CAAC,GAAG,GAAG,CAAC,GAAG,CAAC,SAAS,CAAC,CAAC;gBAC7B,SAAS,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,QAAQ,CAAC,CAAC,CAAC,CAAC;gBACnC,SAAS,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,QAAQ,CAAC,CAAC,CAAC,CAAC;YACzC,CAAC;iBAAM,CAAC;gBACN,MAAM,KAAK,GAAG,GAAG,CAAC,KAAK,CAAC,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC,CAAC;gBACxC,SAAS,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC;gBACvC,SAAS,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC;YAC7C,CAAC;YACD,MAAM,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC;QACpB,CAAC;IACH,CAAC;IAED,OAAO;QACL,CAAC,EAAE,MAAM,CAAC,CAAC;QACX,SAAS;QACT,MAAM;QACN,QAAQ,EAAE,MAAM,CAAC,QAAQ;KAC1B,CAAC;AACJ,CAAC;AAED;;GAEG;AACH,MAAM,CAAC,MAAM,iBAAiB,GAAG;IAC/B,aAAa,EAAE,EAAE,IAAI,EAAE,EAAE,EAAE,CAAC,EAAE,MAAM,EAAE,UAAU,EAAE,EAAE,EAAE,QAAQ,EAAE,WAAoB,EAAE;IACtF,cAAc,EAAE,EAAE,IAAI,EAAE,EAAE,EAAE,CAAC,EAAE,OAAO,EAAE,UAAU,EAAE,EAAE,EAAE,QAAQ,EAAE,WAAoB,EAAE;IACxF,YAAY,EAAE,EAAE,IAAI,EAAE,EAAE,EAAE,CAAC,EAAE,SAAS,EAAE,UAAU,EAAE,EAAE,EAAE,QAAQ,EAAE,WAAoB,EAAE;IACxF,YAAY,EAAE,EAAE,IAAI,EAAE,EAAE,EAAE,CAAC,EAAE,MAAM,EAAE,UAAU,EAAE,EAAE,EAAE,QAAQ,EAAE,UAAmB,EAAE;IACpF,aAAa,EAAE,EAAE,IAAI,EAAE,EAAE,EAAE,CAAC,EAAE,OAAO,EAAE,UAAU,EAAE,EAAE,EAAE,QAAQ,EAAE,UAAmB,EAAE;IACtF,wBAAwB,EAAE,EAAE,IAAI,EAAE,EAAE,EAAE,CAAC,EAAE,OAAO,EAAE,UAAU,EAAE,EAAE,EAAE,QAAQ,EAAE,UAAmB,EAAE,cAAc,EAAE,IAAI,EAAE;CACxH,CAAC"}

package/dist-lib/core/lasso_simplify.d.ts

@@ -0,0 +1,16 @@
+/**
+ * Lasso polygon simplification utilities.
+ *
+ * These are UI-level helpers (not benchmark-critical math) and are intended to
+ * match the "feel" of Embedding Atlas' lasso:
+ * - sample points while dragging
+ * - smooth (Chaikin)
+ * - simplify until under vertex budget
+ */
+/**
+ * Simplify a closed polygon represented as an interleaved XY list.
+ *
+ * Returns a new `Float32Array` of interleaved XY coordinates.
+ */
+export declare function simplifyPolygonData(points: ArrayLike<number>, maxVerts: number): Float32Array;
+//# sourceMappingURL=lasso_simplify.d.ts.map

package/dist-lib/core/lasso_simplify.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"lasso_simplify.d.ts","sourceRoot":"","sources":["../../src/core/lasso_simplify.ts"],"names":[],"mappings":"AAAA;;;;;;;;GAQG;AAqDH;;;;GAIG;AACH,wBAAgB,mBAAmB,CAAC,MAAM,EAAE,SAAS,CAAC,MAAM,CAAC,EAAE,QAAQ,EAAE,MAAM,GAAG,YAAY,CAiH7F"}

package/dist-lib/core/lasso_simplify.js

@@ -0,0 +1,173 @@
+/**
+ * Lasso polygon simplification utilities.
+ *
+ * These are UI-level helpers (not benchmark-critical math) and are intended to
+ * match the "feel" of Embedding Atlas' lasso:
+ * - sample points while dragging
+ * - smooth (Chaikin)
+ * - simplify until under vertex budget
+ */
+function boundingRectXY(points) {
+    let xMin = Infinity;
+    let yMin = Infinity;
+    let xMax = -Infinity;
+    let yMax = -Infinity;
+    for (let i = 0; i < points.length; i += 2) {
+        const x = points[i];
+        const y = points[i + 1];
+        if (x < xMin)
+            xMin = x;
+        if (x > xMax)
+            xMax = x;
+        if (y < yMin)
+            yMin = y;
+        if (y > yMax)
+            yMax = y;
+    }
+    if (!Number.isFinite(xMin) || !Number.isFinite(yMin)) {
+        return { xMin: 0, yMin: 0, xMax: 0, yMax: 0 };
+    }
+    return { xMin, yMin, xMax, yMax };
+}
+function chaikinClosedXY(points, iterations) {
+    // Chaikin smoothing on a closed polygon.
+    let pts = Array.from(points);
+    let n = Math.floor(pts.length / 2);
+    if (n < 3)
+        return pts;
+    for (let it = 0; it < iterations; it++) {
+        n = Math.floor(pts.length / 2);
+        if (n < 3)
+            break;
+        const out = new Array(n * 4);
+        let w = 0;
+        for (let i = 0; i < n; i++) {
+            const i0 = i;
+            const i1 = (i + 1) % n;
+            const x0 = pts[i0 * 2];
+            const y0 = pts[i0 * 2 + 1];
+            const x1 = pts[i1 * 2];
+            const y1 = pts[i1 * 2 + 1];
+            // Q = 0.75*p0 + 0.25*p1
+            out[w++] = 0.75 * x0 + 0.25 * x1;
+            out[w++] = 0.75 * y0 + 0.25 * y1;
+            // R = 0.25*p0 + 0.75*p1
+            out[w++] = 0.25 * x0 + 0.75 * x1;
+            out[w++] = 0.25 * y0 + 0.75 * y1;
+        }
+        pts = out;
+    }
+    return pts;
+}
+/**
+ * Simplify a closed polygon represented as an interleaved XY list.
+ *
+ * Returns a new `Float32Array` of interleaved XY coordinates.
+ */
+export function simplifyPolygonData(points, maxVerts) {
+    // Embedding Atlas strategy (ported conceptually):
+    // 1) Chaikin smoothing on the closed polygon.
+    // 2) Simplify with a tolerance derived from bounding box; increase until within budget.
+    const n0 = Math.floor(points.length / 2);
+    if (n0 <= maxVerts)
+        return new Float32Array(points);
+    // Safety: avoid pathological blow-ups if the user drags for a long time.
+    // If we have too many points, pre-subsample before smoothing.
+    const MAX_RAW = 2048;
+    let base = points;
+    if (n0 > MAX_RAW) {
+        const out = new Float32Array(MAX_RAW * 2);
+        for (let i = 0; i < MAX_RAW; i++) {
+            const src = Math.floor((i * n0) / MAX_RAW);
+            out[i * 2] = points[src * 2];
+            out[i * 2 + 1] = points[src * 2 + 1];
+        }
+        base = out;
+    }
+    const smoothed = chaikinClosedXY(base, 5);
+    const rect = boundingRectXY(smoothed);
+    const span = Math.max(rect.xMax - rect.xMin, rect.yMax - rect.yMin);
+    const distSqToSegment = (px, py, ax, ay, bx, by) => {
+        const abx = bx - ax;
+        const aby = by - ay;
+        const apx = px - ax;
+        const apy = py - ay;
+        const abLenSq = abx * abx + aby * aby;
+        let t = 0;
+        if (abLenSq > 1e-12) {
+            t = (apx * abx + apy * aby) / abLenSq;
+            if (t < 0)
+                t = 0;
+            else if (t > 1)
+                t = 1;
+        }
+        const qx = ax + t * abx;
+        const qy = ay + t * aby;
+        const dx = px - qx;
+        const dy = py - qy;
+        return dx * dx + dy * dy;
+    };
+    const rdp = (eps) => {
+        const m = Math.floor(smoothed.length / 2);
+        if (m <= 2)
+            return smoothed.slice();
+        const epsSq = eps * eps;
+        const keep = new Uint8Array(m);
+        keep[0] = 1;
+        keep[m - 1] = 1;
+        const stack = [0, m - 1];
+        while (stack.length > 0) {
+            const end = stack.pop();
+            const start = stack.pop();
+            const ax = smoothed[start * 2];
+            const ay = smoothed[start * 2 + 1];
+            const bx = smoothed[end * 2];
+            const by = smoothed[end * 2 + 1];
+            let maxD = -1;
+            let maxI = -1;
+            for (let i = start + 1; i < end; i++) {
+                const px = smoothed[i * 2];
+                const py = smoothed[i * 2 + 1];
+                const d = distSqToSegment(px, py, ax, ay, bx, by);
+                if (d > maxD) {
+                    maxD = d;
+                    maxI = i;
+                }
+            }
+            if (maxD > epsSq && maxI >= 0) {
+                keep[maxI] = 1;
+                stack.push(start, maxI);
+                stack.push(maxI, end);
+            }
+        }
+        const out = [];
+        for (let i = 0; i < m; i++) {
+            if (!keep[i])
+                continue;
+            out.push(smoothed[i * 2], smoothed[i * 2 + 1]);
+        }
+        return out;
+    };
+    // Match Embedding Atlas: start with tolerance ~= bbox/100, then increase slowly.
+    let tol = Math.max(1e-12, span / 100);
+    let simplified = rdp(tol);
+    let it = 0;
+    while (Math.floor(simplified.length / 2) > maxVerts && it < 20) {
+        tol *= 1.1;
+        simplified = rdp(tol);
+        it++;
+    }
+    // Final safety: if still too large, fall back to uniform subsampling.
+    const n = Math.floor(simplified.length / 2);
+    if (n > maxVerts) {
+        const out = new Float32Array(maxVerts * 2);
+        for (let i = 0; i < maxVerts; i++) {
+            const src = Math.floor((i * n) / maxVerts);
+            out[i * 2] = simplified[src * 2];
+            out[i * 2 + 1] = simplified[src * 2 + 1];
+        }
+        return out;
+    }
+    return new Float32Array(simplified);
+}
+//# sourceMappingURL=lasso_simplify.js.map

package/dist-lib/core/lasso_simplify.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"lasso_simplify.js","sourceRoot":"","sources":["../../src/core/lasso_simplify.ts"],"names":[],"mappings":"AAAA;;;;;;;;GAQG;AAEH,SAAS,cAAc,CAAC,MAAyB;IAC/C,IAAI,IAAI,GAAG,QAAQ,CAAC;IACpB,IAAI,IAAI,GAAG,QAAQ,CAAC;IACpB,IAAI,IAAI,GAAG,CAAC,QAAQ,CAAC;IACrB,IAAI,IAAI,GAAG,CAAC,QAAQ,CAAC;IACrB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,MAAM,CAAC,MAAM,EAAE,CAAC,IAAI,CAAC,EAAE,CAAC;QAC1C,MAAM,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC;QACpB,MAAM,CAAC,GAAG,MAAM,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;QACxB,IAAI,CAAC,GAAG,IAAI;YAAE,IAAI,GAAG,CAAC,CAAC;QACvB,IAAI,CAAC,GAAG,IAAI;YAAE,IAAI,GAAG,CAAC,CAAC;QACvB,IAAI,CAAC,GAAG,IAAI;YAAE,IAAI,GAAG,CAAC,CAAC;QACvB,IAAI,CAAC,GAAG,IAAI;YAAE,IAAI,GAAG,CAAC,CAAC;IACzB,CAAC;IACD,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,IAAI,CAAC,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,IAAI,CAAC,EAAE,CAAC;QACrD,OAAO,EAAE,IAAI,EAAE,CAAC,EAAE,IAAI,EAAE,CAAC,EAAE,IAAI,EAAE,CAAC,EAAE,IAAI,EAAE,CAAC,EAAE,CAAC;IAChD,CAAC;IACD,OAAO,EAAE,IAAI,EAAE,IAAI,EAAE,IAAI,EAAE,IAAI,EAAE,CAAC;AACpC,CAAC;AAED,SAAS,eAAe,CAAC,MAAyB,EAAE,UAAkB;IACpE,yCAAyC;IACzC,IAAI,GAAG,GAAG,KAAK,CAAC,IAAI,CAAC,MAAM,CAAC,CAAC;IAC7B,IAAI,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,GAAG,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC;IACnC,IAAI,CAAC,GAAG,CAAC;QAAE,OAAO,GAAG,CAAC;IAEtB,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,UAAU,EAAE,EAAE,EAAE,EAAE,CAAC;QACvC,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,GAAG,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC;QAC/B,IAAI,CAAC,GAAG,CAAC;YAAE,MAAM;QACjB,MAAM,GAAG,GAAa,IAAI,KAAK,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;QACvC,IAAI,CAAC,GAAG,CAAC,CAAC;QACV,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;YAC3B,MAAM,EAAE,GAAG,CAAC,CAAC;YACb,MAAM,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC;YACvB,MAAM,EAAE,GAAG,GAAG,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC;YACvB,MAAM,EAAE,GAAG,GAAG,CAAC,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC;YAC3B,MAAM,EAAE,GAAG,GAAG,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC;YACvB,MAAM,EAAE,GAAG,GAAG,CAAC,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC;YAE3B,wBAAwB;YACxB,GAAG,CAAC,CAAC,EAAE,CAAC,GAAG,IAAI,GAAG,EAAE,GAAG,IAAI,GAAG,EAAE,CAAC;YACjC,GAAG,CAAC,CAAC,EAAE,CAAC,GAAG,IAAI,GAAG,EAAE,GAAG,IAAI,GAAG,EAAE,CAAC;YACjC,wBAAwB;YACxB,GAAG,CAAC,CAAC,EAAE,CAAC,GAAG,IAAI,GAAG,EAAE,GAAG,IAAI,GAAG,EAAE,CAAC;YACjC,GAAG,CAAC,CAAC,EAAE,CAAC,GAAG,IAAI,GAAG,EAAE,GAAG,IAAI,GAAG,EAAE,CAAC;QACnC,CAAC;QACD,GAAG,GAAG,GAAG,CAAC;IACZ,CAAC;IAED,OAAO,GAAG,CAAC;AACb,CAAC;AAED;;;;GAIG;AACH,MAAM,UAAU,mBAAmB,CAAC,MAAyB,EAAE,QAAgB;IAC7E,kDAAkD;IAClD,8CAA8C;IAC9C,wFAAwF;IACxF,MAAM,EAAE,GAAG,IAAI,CAAC,KAAK,CAAC,MAAM,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC;IACzC,IAAI,EAAE,IAAI,QAAQ;QAAE,OAAO,IAAI,YAAY,CAAC,MAAa,CAAC,CAAC;IAE3D,yEAAyE;IACzE,8DAA8D;IAC9D,MAAM,OAAO,GAAG,IAAI,CAAC;IACrB,IAAI,IAAI,GAAsB,MAAM,CAAC;IACrC,IAAI,EAAE,GAAG,OAAO,EAAE,CAAC;QACjB,MAAM,GAAG,GAAG,IAAI,YAAY,CAAC,OAAO,GAAG,CAAC,CAAC,CAAC;QAC1C,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,OAAO,EAAE,CAAC,EAAE,EAAE,CAAC;YACjC,MAAM,GAAG,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,GAAG,OAAO,CAAC,CAAC;YAC3C,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,MAAM,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC;YAC7B,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,MAAM,CAAC,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC;QACvC,CAAC;QACD,IAAI,GAAG,GAAG,CAAC;IACb,CAAC;IAED,MAAM,QAAQ,GAAG,eAAe,CAAC,IAAI,EAAE,CAAC,CAAC,CAAC;IAC1C,MAAM,IAAI,GAAG,cAAc,CAAC,QAAQ,CAAC,CAAC;IACtC,MAAM,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,IAAI,GAAG,IAAI,CAAC,IAAI,EAAE,IAAI,CAAC,IAAI,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC;IAEpE,MAAM,eAAe,GAAG,CAAC,EAAU,EAAE,EAAU,EAAE,EAAU,EAAE,EAAU,EAAE,EAAU,EAAE,EAAU,EAAU,EAAE;QACzG,MAAM,GAAG,GAAG,EAAE,GAAG,EAAE,CAAC;QACpB,MAAM,GAAG,GAAG,EAAE,GAAG,EAAE,CAAC;QACpB,MAAM,GAAG,GAAG,EAAE,GAAG,EAAE,CAAC;QACpB,MAAM,GAAG,GAAG,EAAE,GAAG,EAAE,CAAC;QACpB,MAAM,OAAO,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,CAAC;QACtC,IAAI,CAAC,GAAG,CAAC,CAAC;QACV,IAAI,OAAO,GAAG,KAAK,EAAE,CAAC;YACpB,CAAC,GAAG,CAAC,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,CAAC,GAAG,OAAO,CAAC;YACtC,IAAI,CAAC,GAAG,CAAC;gBAAE,CAAC,GAAG,CAAC,CAAC;iBACZ,IAAI,CAAC,GAAG,CAAC;gBAAE,CAAC,GAAG,CAAC,CAAC;QACxB,CAAC;QACD,MAAM,EAAE,GAAG,EAAE,GAAG,CAAC,GAAG,GAAG,CAAC;QACxB,MAAM,EAAE,GAAG,EAAE,GAAG,CAAC,GAAG,GAAG,CAAC;QACxB,MAAM,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC;QACnB,MAAM,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC;QACnB,OAAO,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC;IAC3B,CAAC,CAAC;IAEF,MAAM,GAAG,GAAG,CAAC,GAAW,EAAY,EAAE;QACpC,MAAM,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,QAAQ,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC;QAC1C,IAAI,CAAC,IAAI,CAAC;YAAE,OAAO,QAAQ,CAAC,KAAK,EAAE,CAAC;QAEpC,MAAM,KAAK,GAAG,GAAG,GAAG,GAAG,CAAC;QACxB,MAAM,IAAI,GAAG,IAAI,UAAU,CAAC,CAAC,CAAC,CAAC;QAC/B,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;QACZ,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC;QAEhB,MAAM,KAAK,GAAa,CAAC,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC;QACnC,OAAO,KAAK,CAAC,MAAM,GAAG,CAAC,EAAE,CAAC;YACxB,MAAM,GAAG,GAAG,KAAK,CAAC,GAAG,EAAG,CAAC;YACzB,MAAM,KAAK,GAAG,KAAK,CAAC,GAAG,EAAG,CAAC;YAE3B,MAAM,EAAE,GAAG,QAAQ,CAAC,KAAK,GAAG,CAAC,CAAC,CAAC;YAC/B,MAAM,EAAE,GAAG,QAAQ,CAAC,KAAK,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC;YACnC,MAAM,EAAE,GAAG,QAAQ,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC;YAC7B,MAAM,EAAE,GAAG,QAAQ,CAAC,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC;YAEjC,IAAI,IAAI,GAAG,CAAC,CAAC,CAAC;YACd,IAAI,IAAI,GAAG,CAAC,CAAC,CAAC;YACd,KAAK,IAAI,CAAC,GAAG,KAAK,GAAG,CAAC,EAAE,CAAC,GAAG,GAAG,EAAE,CAAC,EAAE,EAAE,CAAC;gBACrC,MAAM,EAAE,GAAG,QAAQ,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;gBAC3B,MAAM,EAAE,GAAG,QAAQ,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC;gBAC/B,MAAM,CAAC,GAAG,eAAe,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC;gBAClD,IAAI,CAAC,GAAG,IAAI,EAAE,CAAC;oBACb,IAAI,GAAG,CAAC,CAAC;oBACT,IAAI,GAAG,CAAC,CAAC;gBACX,CAAC;YACH,CAAC;YAED,IAAI,IAAI,GAAG,KAAK,IAAI,IAAI,IAAI,CAAC,EAAE,CAAC;gBAC9B,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC;gBACf,KAAK,CAAC,IAAI,CAAC,KAAK,EAAE,IAAI,CAAC,CAAC;gBACxB,KAAK,CAAC,IAAI,CAAC,IAAI,EAAE,GAAG,CAAC,CAAC;YACxB,CAAC;QACH,CAAC;QAED,MAAM,GAAG,GAAa,EAAE,CAAC;QACzB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;YAC3B,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC;gBAAE,SAAS;YACvB,GAAG,CAAC,IAAI,CAAC,QAAQ,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,QAAQ,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;QACjD,CAAC;QACD,OAAO,GAAG,CAAC;IACb,CAAC,CAAC;IAEF,iFAAiF;IACjF,IAAI,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,EAAE,IAAI,GAAG,GAAG,CAAC,CAAC;IACtC,IAAI,UAAU,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC;IAC1B,IAAI,EAAE,GAAG,CAAC,CAAC;IACX,OAAO,IAAI,CAAC,KAAK,CAAC,UAAU,CAAC,MAAM,GAAG,CAAC,CAAC,GAAG,QAAQ,IAAI,EAAE,GAAG,EAAE,EAAE,CAAC;QAC/D,GAAG,IAAI,GAAG,CAAC;QACX,UAAU,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC;QACtB,EAAE,EAAE,CAAC;IACP,CAAC;IAED,sEAAsE;IACtE,MAAM,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,UAAU,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC;IAC5C,IAAI,CAAC,GAAG,QAAQ,EAAE,CAAC;QACjB,MAAM,GAAG,GAAG,IAAI,YAAY,CAAC,QAAQ,GAAG,CAAC,CAAC,CAAC;QAC3C,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,QAAQ,EAAE,CAAC,EAAE,EAAE,CAAC;YAClC,MAAM,GAAG,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,QAAQ,CAAC,CAAC;YAC3C,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,UAAU,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC;YACjC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,UAAU,CAAC,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC;QAC3C,CAAC;QACD,OAAO,GAAG,CAAC;IACb,CAAC;IAED,OAAO,IAAI,YAAY,CAAC,UAAU,CAAC,CAAC;AACtC,CAAC"}

package/dist-lib/core/math/euclidean.d.ts

@@ -0,0 +1,31 @@
+/**
+ * Euclidean geometry utilities.
+ */
+import { EuclideanViewState } from '../types.js';
+/**
+ * Create a default Euclidean view state.
+ */
+export declare function createEuclideanView(): EuclideanViewState;
+/**
+ * Project a data point to screen coordinates.
+ */
+export declare function projectEuclidean(dataX: number, dataY: number, view: EuclideanViewState, width: number, height: number): {
+    x: number;
+    y: number;
+};
+/**
+ * Unproject screen coordinates to data space.
+ */
+export declare function unprojectEuclidean(screenX: number, screenY: number, view: EuclideanViewState, width: number, height: number): {
+    x: number;
+    y: number;
+};
+/**
+ * Apply pan to view (anchor-invariant: point under cursor stays under cursor).
+ */
+export declare function panEuclidean(view: EuclideanViewState, deltaScreenX: number, deltaScreenY: number, width: number, height: number): EuclideanViewState;
+/**
+ * Apply zoom to view (anchor-invariant: point under cursor stays under cursor).
+ */
+export declare function zoomEuclidean(view: EuclideanViewState, anchorScreenX: number, anchorScreenY: number, delta: number, width: number, height: number): EuclideanViewState;
+//# sourceMappingURL=euclidean.d.ts.map

package/dist-lib/core/math/euclidean.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"euclidean.d.ts","sourceRoot":"","sources":["../../../src/core/math/euclidean.ts"],"names":[],"mappings":"AAAA;;GAEG;AAEH,OAAO,EAAE,kBAAkB,EAAE,MAAM,aAAa,CAAC;AAEjD;;GAEG;AACH,wBAAgB,mBAAmB,IAAI,kBAAkB,CAOxD;AAED;;GAEG;AACH,wBAAgB,gBAAgB,CAC9B,KAAK,EAAE,MAAM,EACb,KAAK,EAAE,MAAM,EACb,IAAI,EAAE,kBAAkB,EACxB,KAAK,EAAE,MAAM,EACb,MAAM,EAAE,MAAM,GACb;IAAE,CAAC,EAAE,MAAM,CAAC;IAAC,CAAC,EAAE,MAAM,CAAA;CAAE,CAK1B;AAED;;GAEG;AACH,wBAAgB,kBAAkB,CAChC,OAAO,EAAE,MAAM,EACf,OAAO,EAAE,MAAM,EACf,IAAI,EAAE,kBAAkB,EACxB,KAAK,EAAE,MAAM,EACb,MAAM,EAAE,MAAM,GACb;IAAE,CAAC,EAAE,MAAM,CAAC;IAAC,CAAC,EAAE,MAAM,CAAA;CAAE,CAK1B;AAED;;GAEG;AACH,wBAAgB,YAAY,CAC1B,IAAI,EAAE,kBAAkB,EACxB,YAAY,EAAE,MAAM,EACpB,YAAY,EAAE,MAAM,EACpB,KAAK,EAAE,MAAM,EACb,MAAM,EAAE,MAAM,GACb,kBAAkB,CAOpB;AAED;;GAEG;AACH,wBAAgB,aAAa,CAC3B,IAAI,EAAE,kBAAkB,EACxB,aAAa,EAAE,MAAM,EACrB,aAAa,EAAE,MAAM,EACrB,KAAK,EAAE,MAAM,EACb,KAAK,EAAE,MAAM,EACb,MAAM,EAAE,MAAM,GACb,kBAAkB,CAmBpB"}

package/dist-lib/core/math/euclidean.js

@@ -0,0 +1,64 @@
+/**
+ * Euclidean geometry utilities.
+ */
+/**
+ * Create a default Euclidean view state.
+ */
+export function createEuclideanView() {
+    return {
+        type: 'euclidean',
+        centerX: 0,
+        centerY: 0,
+        zoom: 1,
+    };
+}
+/**
+ * Project a data point to screen coordinates.
+ */
+export function projectEuclidean(dataX, dataY, view, width, height) {
+    const scale = Math.min(width, height) * 0.4 * view.zoom;
+    const x = width / 2 + (dataX - view.centerX) * scale;
+    const y = height / 2 - (dataY - view.centerY) * scale; // flip Y
+    return { x, y };
+}
+/**
+ * Unproject screen coordinates to data space.
+ */
+export function unprojectEuclidean(screenX, screenY, view, width, height) {
+    const scale = Math.min(width, height) * 0.4 * view.zoom;
+    const x = view.centerX + (screenX - width / 2) / scale;
+    const y = view.centerY - (screenY - height / 2) / scale; // flip Y
+    return { x, y };
+}
+/**
+ * Apply pan to view (anchor-invariant: point under cursor stays under cursor).
+ */
+export function panEuclidean(view, deltaScreenX, deltaScreenY, width, height) {
+    const scale = Math.min(width, height) * 0.4 * view.zoom;
+    return {
+        ...view,
+        centerX: view.centerX - deltaScreenX / scale,
+        centerY: view.centerY + deltaScreenY / scale, // flip Y
+    };
+}
+/**
+ * Apply zoom to view (anchor-invariant: point under cursor stays under cursor).
+ */
+export function zoomEuclidean(view, anchorScreenX, anchorScreenY, delta, width, height) {
+    // Get data point under anchor before zoom
+    const anchorData = unprojectEuclidean(anchorScreenX, anchorScreenY, view, width, height);
+    // Apply zoom (delta > 0 = zoom in)
+    const zoomFactor = Math.pow(1.1, delta);
+    const newZoom = Math.max(0.1, Math.min(100, view.zoom * zoomFactor));
+    // Calculate new center so anchor point stays in place
+    const newScale = Math.min(width, height) * 0.4 * newZoom;
+    const newCenterX = anchorData.x - (anchorScreenX - width / 2) / newScale;
+    const newCenterY = anchorData.y + (anchorScreenY - height / 2) / newScale;
+    return {
+        ...view,
+        centerX: newCenterX,
+        centerY: newCenterY,
+        zoom: newZoom,
+    };
+}
+//# sourceMappingURL=euclidean.js.map

package/dist-lib/core/math/euclidean.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"euclidean.js","sourceRoot":"","sources":["../../../src/core/math/euclidean.ts"],"names":[],"mappings":"AAAA;;GAEG;AAIH;;GAEG;AACH,MAAM,UAAU,mBAAmB;IACjC,OAAO;QACL,IAAI,EAAE,WAAW;QACjB,OAAO,EAAE,CAAC;QACV,OAAO,EAAE,CAAC;QACV,IAAI,EAAE,CAAC;KACR,CAAC;AACJ,CAAC;AAED;;GAEG;AACH,MAAM,UAAU,gBAAgB,CAC9B,KAAa,EACb,KAAa,EACb,IAAwB,EACxB,KAAa,EACb,MAAc;IAEd,MAAM,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,EAAE,MAAM,CAAC,GAAG,GAAG,GAAG,IAAI,CAAC,IAAI,CAAC;IACxD,MAAM,CAAC,GAAG,KAAK,GAAG,CAAC,GAAG,CAAC,KAAK,GAAG,IAAI,CAAC,OAAO,CAAC,GAAG,KAAK,CAAC;IACrD,MAAM,CAAC,GAAG,MAAM,GAAG,CAAC,GAAG,CAAC,KAAK,GAAG,IAAI,CAAC,OAAO,CAAC,GAAG,KAAK,CAAC,CAAC,SAAS;IAChE,OAAO,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC;AAClB,CAAC;AAED;;GAEG;AACH,MAAM,UAAU,kBAAkB,CAChC,OAAe,EACf,OAAe,EACf,IAAwB,EACxB,KAAa,EACb,MAAc;IAEd,MAAM,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,EAAE,MAAM,CAAC,GAAG,GAAG,GAAG,IAAI,CAAC,IAAI,CAAC;IACxD,MAAM,CAAC,GAAG,IAAI,CAAC,OAAO,GAAG,CAAC,OAAO,GAAG,KAAK,GAAG,CAAC,CAAC,GAAG,KAAK,CAAC;IACvD,MAAM,CAAC,GAAG,IAAI,CAAC,OAAO,GAAG,CAAC,OAAO,GAAG,MAAM,GAAG,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,SAAS;IAClE,OAAO,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC;AAClB,CAAC;AAED;;GAEG;AACH,MAAM,UAAU,YAAY,CAC1B,IAAwB,EACxB,YAAoB,EACpB,YAAoB,EACpB,KAAa,EACb,MAAc;IAEd,MAAM,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,EAAE,MAAM,CAAC,GAAG,GAAG,GAAG,IAAI,CAAC,IAAI,CAAC;IACxD,OAAO;QACL,GAAG,IAAI;QACP,OAAO,EAAE,IAAI,CAAC,OAAO,GAAG,YAAY,GAAG,KAAK;QAC5C,OAAO,EAAE,IAAI,CAAC,OAAO,GAAG,YAAY,GAAG,KAAK,EAAE,SAAS;KACxD,CAAC;AACJ,CAAC;AAED;;GAEG;AACH,MAAM,UAAU,aAAa,CAC3B,IAAwB,EACxB,aAAqB,EACrB,aAAqB,EACrB,KAAa,EACb,KAAa,EACb,MAAc;IAEd,0CAA0C;IAC1C,MAAM,UAAU,GAAG,kBAAkB,CAAC,aAAa,EAAE,aAAa,EAAE,IAAI,EAAE,KAAK,EAAE,MAAM,CAAC,CAAC;IAEzF,mCAAmC;IACnC,MAAM,UAAU,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,EAAE,KAAK,CAAC,CAAC;IACxC,MAAM,OAAO,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,EAAE,IAAI,CAAC,GAAG,CAAC,GAAG,EAAE,IAAI,CAAC,IAAI,GAAG,UAAU,CAAC,CAAC,CAAC;IAErE,sDAAsD;IACtD,MAAM,QAAQ,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,EAAE,MAAM,CAAC,GAAG,GAAG,GAAG,OAAO,CAAC;IACzD,MAAM,UAAU,GAAG,UAAU,CAAC,CAAC,GAAG,CAAC,aAAa,GAAG,KAAK,GAAG,CAAC,CAAC,GAAG,QAAQ,CAAC;IACzE,MAAM,UAAU,GAAG,UAAU,CAAC,CAAC,GAAG,CAAC,aAAa,GAAG,MAAM,GAAG,CAAC,CAAC,GAAG,QAAQ,CAAC;IAE1E,OAAO;QACL,GAAG,IAAI;QACP,OAAO,EAAE,UAAU;QACnB,OAAO,EAAE,UAAU;QACnB,IAAI,EAAE,OAAO;KACd,CAAC;AACJ,CAAC"}

package/dist-lib/core/math/poincare.d.ts

@@ -0,0 +1,117 @@
+/**
+ * Poincare disk (hyperbolic) geometry utilities.
+ *
+ * The Poincare disk model represents hyperbolic space inside the unit disk.
+ * Points are represented as complex numbers z with |z| < 1.
+ *
+ * Isometries (distance-preserving transformations) are Mobius transformations
+ * of the form: f(z) = (z - a) / (1 - conj(a) * z) where |a| < 1
+ *
+ * This maps point 'a' to the origin, which is our "camera translation".
+ *
+ * References:
+ * - https://math.libretexts.org/Bookshelves/Geometry/Geometry_with_an_Introduction_to_Cosmic_Topology_(Hitchman)/05:_Hyperbolic_Geometry/5.01:_The_Poincare_Disk_Model
+ * - https://geoopt.readthedocs.io/en/latest/extended/stereographic.html
+ */
+import { HyperbolicViewState } from '../types.js';
+/**
+ * Create a default hyperbolic view state (centered at origin).
+ */
+export declare function createHyperbolicView(): HyperbolicViewState;
+/**
+ * Apply Mobius transformation: T_a(z) = (z - a) / (1 - conj(a) * z)
+ * This is the hyperbolic isometry that translates 'a' to the origin.
+ *
+ * Mathematical derivation:
+ * - Numerator: z - a (complex subtraction)
+ * - Denominator: 1 - conj(a) * z where conj(a) = ax - i*ay
+ */
+export declare function mobiusTransform(zx: number, zy: number, ax: number, ay: number): {
+    x: number;
+    y: number;
+};
+/**
+ * Inverse Mobius transformation.
+ * If forward is T_a(z) = (z - a) / (1 - conj(a) * z)
+ * Then inverse is T_a^{-1}(w) = (w + a) / (1 + conj(a) * w)
+ */
+export declare function inverseMobiusTransform(wx: number, wy: number, ax: number, ay: number): {
+    x: number;
+    y: number;
+};
+/**
+ * Mobius addition (gyroaddition) in the Poincare ball.
+ * Formula from geoopt for curvature k = -1:
+ *
+ * x ⊕ y = ((1 + 2⟨x,y⟩ + ‖y‖²)x + (1 - ‖x‖²)y) / (1 + 2⟨x,y⟩ + ‖x‖²‖y‖²)
+ *
+ * where ⟨x,y⟩ is the real inner product (dot product).
+ */
+export declare function mobiusAdd(ax: number, ay: number, bx: number, by: number): {
+    x: number;
+    y: number;
+};
+/**
+ * Mobius negation: -x in the gyrogroup is just negation.
+ */
+export declare function mobiusNeg(x: number, y: number): {
+    x: number;
+    y: number;
+};
+/**
+ * Mobius subtraction: x ⊖ y = x ⊕ (-y)
+ */
+export declare function mobiusSub(ax: number, ay: number, bx: number, by: number): {
+    x: number;
+    y: number;
+};
+/**
+ * Project a Poincare disk point to screen coordinates.
+ */
+export declare function projectPoincare(dataX: number, dataY: number, view: HyperbolicViewState, width: number, height: number): {
+    x: number;
+    y: number;
+};
+/**
+ * Unproject screen coordinates to Poincare disk (data space).
+ */
+export declare function unprojectPoincare(screenX: number, screenY: number, view: HyperbolicViewState, width: number, height: number): {
+    x: number;
+    y: number;
+};
+/**
+ * Apply pan to hyperbolic view (anchor-invariant).
+ *
+ * The point under the cursor at drag start should stay under the cursor
+ * throughout the drag. This is the correct "natural" drag behavior.
+ *
+ * Algorithm:
+ * 1. Find the data point p under the start screen position
+ * 2. Get the target disk position d2 from end screen position
+ * 3. Solve for new camera a' such that T_{a'}(p) = d2
+ */
+export declare function panPoincare(view: HyperbolicViewState, startScreenX: number, startScreenY: number, endScreenX: number, endScreenY: number, width: number, height: number): HyperbolicViewState;
+/**
+ * Apply zoom to hyperbolic view (anchor-invariant).
+ *
+ * Display zoom scales the visual representation while keeping the
+ * data point under the cursor stationary.
+ */
+export declare function zoomPoincare(view: HyperbolicViewState, anchorScreenX: number, anchorScreenY: number, delta: number, width: number, height: number): HyperbolicViewState;
+/**
+ * Hyperbolic distance between two points in the Poincare disk.
+ * d(z1, z2) = 2 * arctanh(|z1 - z2| / |1 - conj(z1) * z2|)
+ */
+export declare function hyperbolicDistance(x1: number, y1: number, x2: number, y2: number): number;
+/**
+ * Geodesic interpolation between two points.
+ * gamma(t) = x ⊕ (t ⊗ ((-x) ⊕ y))
+ *
+ * For simplicity, we use a linear approximation in disk space
+ * which is accurate for small distances.
+ */
+export declare function geodesicInterpolate(x1: number, y1: number, x2: number, y2: number, t: number): {
+    x: number;
+    y: number;
+};
+//# sourceMappingURL=poincare.d.ts.map

package/dist-lib/core/math/poincare.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"poincare.d.ts","sourceRoot":"","sources":["../../../src/core/math/poincare.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;GAcG;AAEH,OAAO,EAAE,mBAAmB,EAAE,MAAM,aAAa,CAAC;AAElD;;GAEG;AACH,wBAAgB,oBAAoB,IAAI,mBAAmB,CAO1D;AAED;;;;;;;GAOG;AACH,wBAAgB,eAAe,CAC7B,EAAE,EAAE,MAAM,EACV,EAAE,EAAE,MAAM,EACV,EAAE,EAAE,MAAM,EACV,EAAE,EAAE,MAAM,GACT;IAAE,CAAC,EAAE,MAAM,CAAC;IAAC,CAAC,EAAE,MAAM,CAAA;CAAE,CA8B1B;AAED;;;;GAIG;AACH,wBAAgB,sBAAsB,CACpC,EAAE,EAAE,MAAM,EACV,EAAE,EAAE,MAAM,EACV,EAAE,EAAE,MAAM,EACV,EAAE,EAAE,MAAM,GACT;IAAE,CAAC,EAAE,MAAM,CAAC;IAAC,CAAC,EAAE,MAAM,CAAA;CAAE,CA4B1B;AAED;;;;;;;GAOG;AACH,wBAAgB,SAAS,CACvB,EAAE,EAAE,MAAM,EACV,EAAE,EAAE,MAAM,EACV,EAAE,EAAE,MAAM,EACV,EAAE,EAAE,MAAM,GACT;IAAE,CAAC,EAAE,MAAM,CAAC;IAAC,CAAC,EAAE,MAAM,CAAA;CAAE,CA2B1B;AAED;;GAEG;AACH,wBAAgB,SAAS,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,GAAG;IAAE,CAAC,EAAE,MAAM,CAAC;IAAC,CAAC,EAAE,MAAM,CAAA;CAAE,CAExE;AAED;;GAEG;AACH,wBAAgB,SAAS,CACvB,EAAE,EAAE,MAAM,EACV,EAAE,EAAE,MAAM,EACV,EAAE,EAAE,MAAM,EACV,EAAE,EAAE,MAAM,GACT;IAAE,CAAC,EAAE,MAAM,CAAC;IAAC,CAAC,EAAE,MAAM,CAAA;CAAE,CAE1B;AAED;;GAEG;AACH,wBAAgB,eAAe,CAC7B,KAAK,EAAE,MAAM,EACb,KAAK,EAAE,MAAM,EACb,IAAI,EAAE,mBAAmB,EACzB,KAAK,EAAE,MAAM,EACb,MAAM,EAAE,MAAM,GACb;IAAE,CAAC,EAAE,MAAM,CAAC;IAAC,CAAC,EAAE,MAAM,CAAA;CAAE,CAW1B;AAED;;GAEG;AACH,wBAAgB,iBAAiB,CAC/B,OAAO,EAAE,MAAM,EACf,OAAO,EAAE,MAAM,EACf,IAAI,EAAE,mBAAmB,EACzB,KAAK,EAAE,MAAM,EACb,MAAM,EAAE,MAAM,GACb;IAAE,CAAC,EAAE,MAAM,CAAC;IAAC,CAAC,EAAE,MAAM,CAAA;CAAE,CAkB1B;AAqDD;;;;;;;;;;GAUG;AACH,wBAAgB,WAAW,CACzB,IAAI,EAAE,mBAAmB,EACzB,YAAY,EAAE,MAAM,EACpB,YAAY,EAAE,MAAM,EACpB,UAAU,EAAE,MAAM,EAClB,UAAU,EAAE,MAAM,EAClB,KAAK,EAAE,MAAM,EACb,MAAM,EAAE,MAAM,GACb,mBAAmB,CAkCrB;AAED;;;;;GAKG;AACH,wBAAgB,YAAY,CAC1B,IAAI,EAAE,mBAAmB,EACzB,aAAa,EAAE,MAAM,EACrB,aAAa,EAAE,MAAM,EACrB,KAAK,EAAE,MAAM,EACb,KAAK,EAAE,MAAM,EACb,MAAM,EAAE,MAAM,GACb,mBAAmB,CAuCrB;AAED;;;GAGG;AACH,wBAAgB,kBAAkB,CAChC,EAAE,EAAE,MAAM,EACV,EAAE,EAAE,MAAM,EACV,EAAE,EAAE,MAAM,EACV,EAAE,EAAE,MAAM,GACT,MAAM,CAgBR;AAED;;;;;;GAMG;AACH,wBAAgB,mBAAmB,CACjC,EAAE,EAAE,MAAM,EACV,EAAE,EAAE,MAAM,EACV,EAAE,EAAE,MAAM,EACV,EAAE,EAAE,MAAM,EACV,CAAC,EAAE,MAAM,GACR;IAAE,CAAC,EAAE,MAAM,CAAC;IAAC,CAAC,EAAE,MAAM,CAAA;CAAE,CAsB1B"}