hyper-scatter 0.1.0

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

@@ -0,0 +1,321 @@
+/**
+ * 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
+ */
+/**
+ * Create a default hyperbolic view state (centered at origin).
+ */
+export function createHyperbolicView() {
+    return {
+        type: 'poincare',
+        ax: 0,
+        ay: 0,
+        displayZoom: 1,
+    };
+}
+/**
+ * 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 function mobiusTransform(zx, zy, ax, ay) {
+    // Numerator: z - a
+    const numX = zx - ax;
+    const numY = zy - ay;
+    // Denominator: 1 - conj(a) * z = 1 - (ax - i*ay)(zx + i*zy)
+    // = 1 - (ax*zx + ay*zy + i*(ax*zy - ay*zx))
+    const denomX = 1 - (ax * zx + ay * zy);
+    const denomY = -(ax * zy - ay * zx);
+    // Complex division: num / denom
+    const denomNormSq = denomX * denomX + denomY * denomY;
+    if (denomNormSq < 1e-12) {
+        // Degenerate case, return clamped to boundary
+        const norm = Math.sqrt(numX * numX + numY * numY);
+        if (norm < 1e-12)
+            return { x: 0, y: 0 };
+        return { x: (numX / norm) * 0.999, y: (numY / norm) * 0.999 };
+    }
+    const resultX = (numX * denomX + numY * denomY) / denomNormSq;
+    const resultY = (numY * denomX - numX * denomY) / denomNormSq;
+    // Clamp to disk (numerical safety)
+    const rSq = resultX * resultX + resultY * resultY;
+    if (rSq >= 1) {
+        const r = Math.sqrt(rSq);
+        return { x: (resultX / r) * 0.999, y: (resultY / r) * 0.999 };
+    }
+    return { x: resultX, y: resultY };
+}
+/**
+ * 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 function inverseMobiusTransform(wx, wy, ax, ay) {
+    // Numerator: w + a
+    const numX = wx + ax;
+    const numY = wy + ay;
+    // Denominator: 1 + conj(a) * w = 1 + (ax - i*ay)(wx + i*wy)
+    // = 1 + (ax*wx + ay*wy + i*(ax*wy - ay*wx))
+    const denomX = 1 + (ax * wx + ay * wy);
+    const denomY = ax * wy - ay * wx;
+    // Complex division
+    const denomNormSq = denomX * denomX + denomY * denomY;
+    if (denomNormSq < 1e-12) {
+        const norm = Math.sqrt(numX * numX + numY * numY);
+        if (norm < 1e-12)
+            return { x: 0, y: 0 };
+        return { x: (numX / norm) * 0.999, y: (numY / norm) * 0.999 };
+    }
+    const resultX = (numX * denomX + numY * denomY) / denomNormSq;
+    const resultY = (numY * denomX - numX * denomY) / denomNormSq;
+    const rSq = resultX * resultX + resultY * resultY;
+    if (rSq >= 1) {
+        const r = Math.sqrt(rSq);
+        return { x: (resultX / r) * 0.999, y: (resultY / r) * 0.999 };
+    }
+    return { x: resultX, y: resultY };
+}
+/**
+ * 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 function mobiusAdd(ax, ay, bx, by) {
+    const dotAB = ax * bx + ay * by; // ⟨a, b⟩
+    const normASq = ax * ax + ay * ay; // ‖a‖²
+    const normBSq = bx * bx + by * by; // ‖b‖²
+    // Numerator coefficients
+    const coeffA = 1 + 2 * dotAB + normBSq;
+    const coeffB = 1 - normASq;
+    // Denominator
+    const denom = 1 + 2 * dotAB + normASq * normBSq;
+    if (Math.abs(denom) < 1e-12) {
+        return { x: 0, y: 0 };
+    }
+    const resultX = (coeffA * ax + coeffB * bx) / denom;
+    const resultY = (coeffA * ay + coeffB * by) / denom;
+    // Clamp to disk
+    const rSq = resultX * resultX + resultY * resultY;
+    if (rSq >= 1) {
+        const r = Math.sqrt(rSq);
+        return { x: (resultX / r) * 0.99, y: (resultY / r) * 0.99 };
+    }
+    return { x: resultX, y: resultY };
+}
+/**
+ * Mobius negation: -x in the gyrogroup is just negation.
+ */
+export function mobiusNeg(x, y) {
+    return { x: -x, y: -y };
+}
+/**
+ * Mobius subtraction: x ⊖ y = x ⊕ (-y)
+ */
+export function mobiusSub(ax, ay, bx, by) {
+    return mobiusAdd(ax, ay, -bx, -by);
+}
+/**
+ * Project a Poincare disk point to screen coordinates.
+ */
+export function projectPoincare(dataX, dataY, view, width, height) {
+    // First apply the Mobius transformation (camera)
+    const transformed = mobiusTransform(dataX, dataY, view.ax, view.ay);
+    // Then apply display zoom and map to screen
+    // The disk has radius 1, we map it to fit in the canvas
+    const diskRadius = Math.min(width, height) * 0.45 * view.displayZoom;
+    const x = width / 2 + transformed.x * diskRadius;
+    const y = height / 2 - transformed.y * diskRadius; // flip Y
+    return { x, y };
+}
+/**
+ * Unproject screen coordinates to Poincare disk (data space).
+ */
+export function unprojectPoincare(screenX, screenY, view, width, height) {
+    // Map screen to disk coordinates (after camera)
+    const diskRadius = Math.min(width, height) * 0.45 * view.displayZoom;
+    const diskX = (screenX - width / 2) / diskRadius;
+    const diskY = -(screenY - height / 2) / diskRadius; // flip Y
+    // Check if point is outside disk
+    const rSq = diskX * diskX + diskY * diskY;
+    if (rSq >= 1) {
+        // Clamp to disk boundary
+        const r = Math.sqrt(rSq);
+        const clampedX = (diskX / r) * 0.999;
+        const clampedY = (diskY / r) * 0.999;
+        return inverseMobiusTransform(clampedX, clampedY, view.ax, view.ay);
+    }
+    // Apply inverse Mobius to get data coordinates
+    return inverseMobiusTransform(diskX, diskY, view.ax, view.ay);
+}
+/**
+ * Solve for new camera position a' such that T_{a'}(p) = d2
+ *
+ * Given: data point p, target display position d2
+ * Find: camera a' such that (p - a') / (1 - conj(a') * p) = d2
+ *
+ * This is solved as a linear system in the real and imaginary parts of a'.
+ *
+ * The system is:
+ *   (1 - A) * x - B * y = px - d2x
+ *   B * x - (1 + A) * y = d2y - py
+ *
+ * where A = d2x*px - d2y*py, B = d2x*py + d2y*px
+ * and the determinant is A² + B² - 1 = |d2|²|p|² - 1
+ */
+function solveForCamera(px, py, d2x, d2y) {
+    const A = d2x * px - d2y * py;
+    const B = d2x * py + d2y * px;
+    // Determinant = |d2|² * |p|² - 1
+    // For points inside disk, this is always negative
+    const det = A * A + B * B - 1;
+    if (Math.abs(det) < 1e-10) {
+        // Degenerate case: either p or d2 is at origin, or on boundary
+        // Fall back to simple approximation
+        return { x: -d2x, y: -d2y };
+    }
+    const rhsX = px - d2x;
+    const rhsY = d2y - py;
+    // Cramer's rule
+    const x = (-(1 + A) * rhsX + B * rhsY) / det;
+    const y = ((1 - A) * rhsY - B * rhsX) / det;
+    // Clamp to ensure |a'| < 1
+    const rSq = x * x + y * y;
+    if (rSq >= 1) {
+        const r = Math.sqrt(rSq);
+        return { x: (x / r) * 0.99, y: (y / r) * 0.99 };
+    }
+    return { x, y };
+}
+/**
+ * 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 function panPoincare(view, startScreenX, startScreenY, endScreenX, endScreenY, width, height) {
+    const diskRadius = Math.min(width, height) * 0.45 * view.displayZoom;
+    // Convert screen positions to disk coordinates (display space, after camera)
+    const d1x = (startScreenX - width / 2) / diskRadius;
+    const d1y = -(startScreenY - height / 2) / diskRadius;
+    const d2x = (endScreenX - width / 2) / diskRadius;
+    const d2y = -(endScreenY - height / 2) / diskRadius;
+    // Clamp to disk with margin
+    const clampToDisk = (x, y, maxR = 0.95) => {
+        const rSq = x * x + y * y;
+        if (rSq > maxR * maxR) {
+            const r = Math.sqrt(rSq);
+            return { x: (x / r) * maxR, y: (y / r) * maxR };
+        }
+        return { x, y };
+    };
+    const d1 = clampToDisk(d1x, d1y);
+    const d2 = clampToDisk(d2x, d2y);
+    // Get the data point currently displayed at d1
+    // p = T_a^{-1}(d1)
+    const p = inverseMobiusTransform(d1.x, d1.y, view.ax, view.ay);
+    // Solve for new camera a' such that T_{a'}(p) = d2
+    const newA = solveForCamera(p.x, p.y, d2.x, d2.y);
+    return {
+        ...view,
+        ax: newA.x,
+        ay: newA.y,
+    };
+}
+/**
+ * Apply zoom to hyperbolic view (anchor-invariant).
+ *
+ * Display zoom scales the visual representation while keeping the
+ * data point under the cursor stationary.
+ */
+export function zoomPoincare(view, anchorScreenX, anchorScreenY, delta, width, height) {
+    // Calculate new display zoom
+    const zoomFactor = Math.pow(1.1, delta);
+    const newDisplayZoom = Math.max(0.5, Math.min(10, view.displayZoom * zoomFactor));
+    // Get the data point under the anchor
+    const dataPoint = unprojectPoincare(anchorScreenX, anchorScreenY, view, width, height);
+    // Calculate where this data point would appear with new zoom (but same camera)
+    const diskRadiusNew = Math.min(width, height) * 0.45 * newDisplayZoom;
+    // The display position of dataPoint after camera transform
+    const transformed = mobiusTransform(dataPoint.x, dataPoint.y, view.ax, view.ay);
+    // Screen position with new zoom (if we don't adjust camera)
+    const newScreenX = width / 2 + transformed.x * diskRadiusNew;
+    const newScreenY = height / 2 - transformed.y * diskRadiusNew;
+    // If the point moved significantly, adjust camera to compensate
+    const dx = newScreenX - anchorScreenX;
+    const dy = newScreenY - anchorScreenY;
+    if (Math.abs(dx) > 0.5 || Math.abs(dy) > 0.5) {
+        // Pan to bring the anchor point back under the cursor
+        return panPoincare({ ...view, displayZoom: newDisplayZoom }, newScreenX, newScreenY, anchorScreenX, anchorScreenY, width, height);
+    }
+    return {
+        ...view,
+        displayZoom: newDisplayZoom,
+    };
+}
+/**
+ * Hyperbolic distance between two points in the Poincare disk.
+ * d(z1, z2) = 2 * arctanh(|z1 - z2| / |1 - conj(z1) * z2|)
+ */
+export function hyperbolicDistance(x1, y1, x2, y2) {
+    const diffX = x1 - x2;
+    const diffY = y1 - y2;
+    const diffNorm = Math.sqrt(diffX * diffX + diffY * diffY);
+    // 1 - conj(z1) * z2
+    const denomX = 1 - (x1 * x2 + y1 * y2);
+    const denomY = -(x1 * y2 - y1 * x2);
+    const denomNorm = Math.sqrt(denomX * denomX + denomY * denomY);
+    if (denomNorm < 1e-12)
+        return Infinity;
+    // Clamp ratio to avoid atanh(1) = Infinity due to floating-point errors
+    const ratio = Math.min(diffNorm / denomNorm, 1 - 1e-10);
+    return 2 * Math.atanh(ratio);
+}
+/**
+ * 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 function geodesicInterpolate(x1, y1, x2, y2, t) {
+    // Transform so x1 is at origin
+    const transformed = mobiusTransform(x2, y2, x1, y1);
+    // Scale by t (this is exact for geodesics through origin)
+    const scaledNorm = Math.sqrt(transformed.x * transformed.x + transformed.y * transformed.y);
+    if (scaledNorm < 1e-12) {
+        return { x: x1, y: y1 };
+    }
+    // Clamp to avoid atanh(1) = Infinity near boundary
+    const clampedNorm = Math.min(scaledNorm, 1 - 1e-10);
+    // Use arctanh/tanh for proper geodesic scaling
+    const targetNorm = Math.tanh(t * Math.atanh(clampedNorm));
+    const scaled = {
+        x: (transformed.x / scaledNorm) * targetNorm,
+        y: (transformed.y / scaledNorm) * targetNorm,
+    };
+    // Transform back
+    return inverseMobiusTransform(scaled.x, scaled.y, x1, y1);
+}
+//# sourceMappingURL=poincare.js.map

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

@@ -0,0 +1 @@
+{"version":3,"file":"poincare.js","sourceRoot":"","sources":["../../../src/core/math/poincare.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;GAcG;AAIH;;GAEG;AACH,MAAM,UAAU,oBAAoB;IAClC,OAAO;QACL,IAAI,EAAE,UAAU;QAChB,EAAE,EAAE,CAAC;QACL,EAAE,EAAE,CAAC;QACL,WAAW,EAAE,CAAC;KACf,CAAC;AACJ,CAAC;AAED;;;;;;;GAOG;AACH,MAAM,UAAU,eAAe,CAC7B,EAAU,EACV,EAAU,EACV,EAAU,EACV,EAAU;IAEV,mBAAmB;IACnB,MAAM,IAAI,GAAG,EAAE,GAAG,EAAE,CAAC;IACrB,MAAM,IAAI,GAAG,EAAE,GAAG,EAAE,CAAC;IAErB,4DAA4D;IAC5D,4CAA4C;IAC5C,MAAM,MAAM,GAAG,CAAC,GAAG,CAAC,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC,CAAC;IACvC,MAAM,MAAM,GAAG,CAAC,CAAC,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC,CAAC;IAEpC,gCAAgC;IAChC,MAAM,WAAW,GAAG,MAAM,GAAG,MAAM,GAAG,MAAM,GAAG,MAAM,CAAC;IACtD,IAAI,WAAW,GAAG,KAAK,EAAE,CAAC;QACxB,8CAA8C;QAC9C,MAAM,IAAI,GAAG,IAAI,CAAC,IAAI,CAAC,IAAI,GAAG,IAAI,GAAG,IAAI,GAAG,IAAI,CAAC,CAAC;QAClD,IAAI,IAAI,GAAG,KAAK;YAAE,OAAO,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC;QACxC,OAAO,EAAE,CAAC,EAAE,CAAC,IAAI,GAAG,IAAI,CAAC,GAAG,KAAK,EAAE,CAAC,EAAE,CAAC,IAAI,GAAG,IAAI,CAAC,GAAG,KAAK,EAAE,CAAC;IAChE,CAAC;IAED,MAAM,OAAO,GAAG,CAAC,IAAI,GAAG,MAAM,GAAG,IAAI,GAAG,MAAM,CAAC,GAAG,WAAW,CAAC;IAC9D,MAAM,OAAO,GAAG,CAAC,IAAI,GAAG,MAAM,GAAG,IAAI,GAAG,MAAM,CAAC,GAAG,WAAW,CAAC;IAE9D,mCAAmC;IACnC,MAAM,GAAG,GAAG,OAAO,GAAG,OAAO,GAAG,OAAO,GAAG,OAAO,CAAC;IAClD,IAAI,GAAG,IAAI,CAAC,EAAE,CAAC;QACb,MAAM,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC;QACzB,OAAO,EAAE,CAAC,EAAE,CAAC,OAAO,GAAG,CAAC,CAAC,GAAG,KAAK,EAAE,CAAC,EAAE,CAAC,OAAO,GAAG,CAAC,CAAC,GAAG,KAAK,EAAE,CAAC;IAChE,CAAC;IAED,OAAO,EAAE,CAAC,EAAE,OAAO,EAAE,CAAC,EAAE,OAAO,EAAE,CAAC;AACpC,CAAC;AAED;;;;GAIG;AACH,MAAM,UAAU,sBAAsB,CACpC,EAAU,EACV,EAAU,EACV,EAAU,EACV,EAAU;IAEV,mBAAmB;IACnB,MAAM,IAAI,GAAG,EAAE,GAAG,EAAE,CAAC;IACrB,MAAM,IAAI,GAAG,EAAE,GAAG,EAAE,CAAC;IAErB,4DAA4D;IAC5D,4CAA4C;IAC5C,MAAM,MAAM,GAAG,CAAC,GAAG,CAAC,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC,CAAC;IACvC,MAAM,MAAM,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC;IAEjC,mBAAmB;IACnB,MAAM,WAAW,GAAG,MAAM,GAAG,MAAM,GAAG,MAAM,GAAG,MAAM,CAAC;IACtD,IAAI,WAAW,GAAG,KAAK,EAAE,CAAC;QACxB,MAAM,IAAI,GAAG,IAAI,CAAC,IAAI,CAAC,IAAI,GAAG,IAAI,GAAG,IAAI,GAAG,IAAI,CAAC,CAAC;QAClD,IAAI,IAAI,GAAG,KAAK;YAAE,OAAO,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC;QACxC,OAAO,EAAE,CAAC,EAAE,CAAC,IAAI,GAAG,IAAI,CAAC,GAAG,KAAK,EAAE,CAAC,EAAE,CAAC,IAAI,GAAG,IAAI,CAAC,GAAG,KAAK,EAAE,CAAC;IAChE,CAAC;IAED,MAAM,OAAO,GAAG,CAAC,IAAI,GAAG,MAAM,GAAG,IAAI,GAAG,MAAM,CAAC,GAAG,WAAW,CAAC;IAC9D,MAAM,OAAO,GAAG,CAAC,IAAI,GAAG,MAAM,GAAG,IAAI,GAAG,MAAM,CAAC,GAAG,WAAW,CAAC;IAE9D,MAAM,GAAG,GAAG,OAAO,GAAG,OAAO,GAAG,OAAO,GAAG,OAAO,CAAC;IAClD,IAAI,GAAG,IAAI,CAAC,EAAE,CAAC;QACb,MAAM,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC;QACzB,OAAO,EAAE,CAAC,EAAE,CAAC,OAAO,GAAG,CAAC,CAAC,GAAG,KAAK,EAAE,CAAC,EAAE,CAAC,OAAO,GAAG,CAAC,CAAC,GAAG,KAAK,EAAE,CAAC;IAChE,CAAC;IAED,OAAO,EAAE,CAAC,EAAE,OAAO,EAAE,CAAC,EAAE,OAAO,EAAE,CAAC;AACpC,CAAC;AAED;;;;;;;GAOG;AACH,MAAM,UAAU,SAAS,CACvB,EAAU,EACV,EAAU,EACV,EAAU,EACV,EAAU;IAEV,MAAM,KAAK,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC,CAAC,SAAS;IAC1C,MAAM,OAAO,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC,CAAC,OAAO;IAC1C,MAAM,OAAO,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC,CAAC,OAAO;IAE1C,yBAAyB;IACzB,MAAM,MAAM,GAAG,CAAC,GAAG,CAAC,GAAG,KAAK,GAAG,OAAO,CAAC;IACvC,MAAM,MAAM,GAAG,CAAC,GAAG,OAAO,CAAC;IAE3B,cAAc;IACd,MAAM,KAAK,GAAG,CAAC,GAAG,CAAC,GAAG,KAAK,GAAG,OAAO,GAAG,OAAO,CAAC;IAEhD,IAAI,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,GAAG,KAAK,EAAE,CAAC;QAC5B,OAAO,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC;IACxB,CAAC;IAED,MAAM,OAAO,GAAG,CAAC,MAAM,GAAG,EAAE,GAAG,MAAM,GAAG,EAAE,CAAC,GAAG,KAAK,CAAC;IACpD,MAAM,OAAO,GAAG,CAAC,MAAM,GAAG,EAAE,GAAG,MAAM,GAAG,EAAE,CAAC,GAAG,KAAK,CAAC;IAEpD,gBAAgB;IAChB,MAAM,GAAG,GAAG,OAAO,GAAG,OAAO,GAAG,OAAO,GAAG,OAAO,CAAC;IAClD,IAAI,GAAG,IAAI,CAAC,EAAE,CAAC;QACb,MAAM,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC;QACzB,OAAO,EAAE,CAAC,EAAE,CAAC,OAAO,GAAG,CAAC,CAAC,GAAG,IAAI,EAAE,CAAC,EAAE,CAAC,OAAO,GAAG,CAAC,CAAC,GAAG,IAAI,EAAE,CAAC;IAC9D,CAAC;IAED,OAAO,EAAE,CAAC,EAAE,OAAO,EAAE,CAAC,EAAE,OAAO,EAAE,CAAC;AACpC,CAAC;AAED;;GAEG;AACH,MAAM,UAAU,SAAS,CAAC,CAAS,EAAE,CAAS;IAC5C,OAAO,EAAE,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC;AAC1B,CAAC;AAED;;GAEG;AACH,MAAM,UAAU,SAAS,CACvB,EAAU,EACV,EAAU,EACV,EAAU,EACV,EAAU;IAEV,OAAO,SAAS,CAAC,EAAE,EAAE,EAAE,EAAE,CAAC,EAAE,EAAE,CAAC,EAAE,CAAC,CAAC;AACrC,CAAC;AAED;;GAEG;AACH,MAAM,UAAU,eAAe,CAC7B,KAAa,EACb,KAAa,EACb,IAAyB,EACzB,KAAa,EACb,MAAc;IAEd,iDAAiD;IACjD,MAAM,WAAW,GAAG,eAAe,CAAC,KAAK,EAAE,KAAK,EAAE,IAAI,CAAC,EAAE,EAAE,IAAI,CAAC,EAAE,CAAC,CAAC;IAEpE,4CAA4C;IAC5C,wDAAwD;IACxD,MAAM,UAAU,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,EAAE,MAAM,CAAC,GAAG,IAAI,GAAG,IAAI,CAAC,WAAW,CAAC;IACrE,MAAM,CAAC,GAAG,KAAK,GAAG,CAAC,GAAG,WAAW,CAAC,CAAC,GAAG,UAAU,CAAC;IACjD,MAAM,CAAC,GAAG,MAAM,GAAG,CAAC,GAAG,WAAW,CAAC,CAAC,GAAG,UAAU,CAAC,CAAC,SAAS;IAE5D,OAAO,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC;AAClB,CAAC;AAED;;GAEG;AACH,MAAM,UAAU,iBAAiB,CAC/B,OAAe,EACf,OAAe,EACf,IAAyB,EACzB,KAAa,EACb,MAAc;IAEd,gDAAgD;IAChD,MAAM,UAAU,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,EAAE,MAAM,CAAC,GAAG,IAAI,GAAG,IAAI,CAAC,WAAW,CAAC;IACrE,MAAM,KAAK,GAAG,CAAC,OAAO,GAAG,KAAK,GAAG,CAAC,CAAC,GAAG,UAAU,CAAC;IACjD,MAAM,KAAK,GAAG,CAAC,CAAC,OAAO,GAAG,MAAM,GAAG,CAAC,CAAC,GAAG,UAAU,CAAC,CAAC,SAAS;IAE7D,iCAAiC;IACjC,MAAM,GAAG,GAAG,KAAK,GAAG,KAAK,GAAG,KAAK,GAAG,KAAK,CAAC;IAC1C,IAAI,GAAG,IAAI,CAAC,EAAE,CAAC;QACb,yBAAyB;QACzB,MAAM,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC;QACzB,MAAM,QAAQ,GAAG,CAAC,KAAK,GAAG,CAAC,CAAC,GAAG,KAAK,CAAC;QACrC,MAAM,QAAQ,GAAG,CAAC,KAAK,GAAG,CAAC,CAAC,GAAG,KAAK,CAAC;QACrC,OAAO,sBAAsB,CAAC,QAAQ,EAAE,QAAQ,EAAE,IAAI,CAAC,EAAE,EAAE,IAAI,CAAC,EAAE,CAAC,CAAC;IACtE,CAAC;IAED,+CAA+C;IAC/C,OAAO,sBAAsB,CAAC,KAAK,EAAE,KAAK,EAAE,IAAI,CAAC,EAAE,EAAE,IAAI,CAAC,EAAE,CAAC,CAAC;AAChE,CAAC;AAED;;;;;;;;;;;;;;GAcG;AACH,SAAS,cAAc,CACrB,EAAU,EACV,EAAU,EACV,GAAW,EACX,GAAW;IAEX,MAAM,CAAC,GAAG,GAAG,GAAG,EAAE,GAAG,GAAG,GAAG,EAAE,CAAC;IAC9B,MAAM,CAAC,GAAG,GAAG,GAAG,EAAE,GAAG,GAAG,GAAG,EAAE,CAAC;IAE9B,iCAAiC;IACjC,kDAAkD;IAClD,MAAM,GAAG,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;IAE9B,IAAI,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,KAAK,EAAE,CAAC;QAC1B,+DAA+D;QAC/D,oCAAoC;QACpC,OAAO,EAAE,CAAC,EAAE,CAAC,GAAG,EAAE,CAAC,EAAE,CAAC,GAAG,EAAE,CAAC;IAC9B,CAAC;IAED,MAAM,IAAI,GAAG,EAAE,GAAG,GAAG,CAAC;IACtB,MAAM,IAAI,GAAG,GAAG,GAAG,EAAE,CAAC;IAEtB,gBAAgB;IAChB,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,GAAG,CAAC;IAC7C,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,GAAG,CAAC;IAE5C,2BAA2B;IAC3B,MAAM,GAAG,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;IAC1B,IAAI,GAAG,IAAI,CAAC,EAAE,CAAC;QACb,MAAM,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC;QACzB,OAAO,EAAE,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,EAAE,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,EAAE,CAAC;IAClD,CAAC;IAED,OAAO,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC;AAClB,CAAC;AAED;;;;;;;;;;GAUG;AACH,MAAM,UAAU,WAAW,CACzB,IAAyB,EACzB,YAAoB,EACpB,YAAoB,EACpB,UAAkB,EAClB,UAAkB,EAClB,KAAa,EACb,MAAc;IAEd,MAAM,UAAU,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,EAAE,MAAM,CAAC,GAAG,IAAI,GAAG,IAAI,CAAC,WAAW,CAAC;IAErE,6EAA6E;IAC7E,MAAM,GAAG,GAAG,CAAC,YAAY,GAAG,KAAK,GAAG,CAAC,CAAC,GAAG,UAAU,CAAC;IACpD,MAAM,GAAG,GAAG,CAAC,CAAC,YAAY,GAAG,MAAM,GAAG,CAAC,CAAC,GAAG,UAAU,CAAC;IACtD,MAAM,GAAG,GAAG,CAAC,UAAU,GAAG,KAAK,GAAG,CAAC,CAAC,GAAG,UAAU,CAAC;IAClD,MAAM,GAAG,GAAG,CAAC,CAAC,UAAU,GAAG,MAAM,GAAG,CAAC,CAAC,GAAG,UAAU,CAAC;IAEpD,4BAA4B;IAC5B,MAAM,WAAW,GAAG,CAAC,CAAS,EAAE,CAAS,EAAE,OAAe,IAAI,EAAE,EAAE;QAChE,MAAM,GAAG,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;QAC1B,IAAI,GAAG,GAAG,IAAI,GAAG,IAAI,EAAE,CAAC;YACtB,MAAM,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC;YACzB,OAAO,EAAE,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,EAAE,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,EAAE,CAAC;QAClD,CAAC;QACD,OAAO,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC;IAClB,CAAC,CAAC;IAEF,MAAM,EAAE,GAAG,WAAW,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC;IACjC,MAAM,EAAE,GAAG,WAAW,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC;IAEjC,+CAA+C;IAC/C,mBAAmB;IACnB,MAAM,CAAC,GAAG,sBAAsB,CAAC,EAAE,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,EAAE,IAAI,CAAC,EAAE,EAAE,IAAI,CAAC,EAAE,CAAC,CAAC;IAE/D,mDAAmD;IACnD,MAAM,IAAI,GAAG,cAAc,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC;IAElD,OAAO;QACL,GAAG,IAAI;QACP,EAAE,EAAE,IAAI,CAAC,CAAC;QACV,EAAE,EAAE,IAAI,CAAC,CAAC;KACX,CAAC;AACJ,CAAC;AAED;;;;;GAKG;AACH,MAAM,UAAU,YAAY,CAC1B,IAAyB,EACzB,aAAqB,EACrB,aAAqB,EACrB,KAAa,EACb,KAAa,EACb,MAAc;IAEd,6BAA6B;IAC7B,MAAM,UAAU,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,EAAE,KAAK,CAAC,CAAC;IACxC,MAAM,cAAc,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,EAAE,IAAI,CAAC,GAAG,CAAC,EAAE,EAAE,IAAI,CAAC,WAAW,GAAG,UAAU,CAAC,CAAC,CAAC;IAElF,sCAAsC;IACtC,MAAM,SAAS,GAAG,iBAAiB,CAAC,aAAa,EAAE,aAAa,EAAE,IAAI,EAAE,KAAK,EAAE,MAAM,CAAC,CAAC;IAEvF,+EAA+E;IAC/E,MAAM,aAAa,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,EAAE,MAAM,CAAC,GAAG,IAAI,GAAG,cAAc,CAAC;IAEtE,2DAA2D;IAC3D,MAAM,WAAW,GAAG,eAAe,CAAC,SAAS,CAAC,CAAC,EAAE,SAAS,CAAC,CAAC,EAAE,IAAI,CAAC,EAAE,EAAE,IAAI,CAAC,EAAE,CAAC,CAAC;IAEhF,4DAA4D;IAC5D,MAAM,UAAU,GAAG,KAAK,GAAG,CAAC,GAAG,WAAW,CAAC,CAAC,GAAG,aAAa,CAAC;IAC7D,MAAM,UAAU,GAAG,MAAM,GAAG,CAAC,GAAG,WAAW,CAAC,CAAC,GAAG,aAAa,CAAC;IAE9D,gEAAgE;IAChE,MAAM,EAAE,GAAG,UAAU,GAAG,aAAa,CAAC;IACtC,MAAM,EAAE,GAAG,UAAU,GAAG,aAAa,CAAC;IAEtC,IAAI,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,GAAG,IAAI,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,GAAG,EAAE,CAAC;QAC7C,sDAAsD;QACtD,OAAO,WAAW,CAChB,EAAE,GAAG,IAAI,EAAE,WAAW,EAAE,cAAc,EAAE,EACxC,UAAU,EACV,UAAU,EACV,aAAa,EACb,aAAa,EACb,KAAK,EACL,MAAM,CACP,CAAC;IACJ,CAAC;IAED,OAAO;QACL,GAAG,IAAI;QACP,WAAW,EAAE,cAAc;KAC5B,CAAC;AACJ,CAAC;AAED;;;GAGG;AACH,MAAM,UAAU,kBAAkB,CAChC,EAAU,EACV,EAAU,EACV,EAAU,EACV,EAAU;IAEV,MAAM,KAAK,GAAG,EAAE,GAAG,EAAE,CAAC;IACtB,MAAM,KAAK,GAAG,EAAE,GAAG,EAAE,CAAC;IACtB,MAAM,QAAQ,GAAG,IAAI,CAAC,IAAI,CAAC,KAAK,GAAG,KAAK,GAAG,KAAK,GAAG,KAAK,CAAC,CAAC;IAE1D,oBAAoB;IACpB,MAAM,MAAM,GAAG,CAAC,GAAG,CAAC,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC,CAAC;IACvC,MAAM,MAAM,GAAG,CAAC,CAAC,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC,CAAC;IACpC,MAAM,SAAS,GAAG,IAAI,CAAC,IAAI,CAAC,MAAM,GAAG,MAAM,GAAG,MAAM,GAAG,MAAM,CAAC,CAAC;IAE/D,IAAI,SAAS,GAAG,KAAK;QAAE,OAAO,QAAQ,CAAC;IAEvC,wEAAwE;IACxE,MAAM,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,QAAQ,GAAG,SAAS,EAAE,CAAC,GAAG,KAAK,CAAC,CAAC;IAExD,OAAO,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,KAAK,CAAC,CAAC;AAC/B,CAAC;AAED;;;;;;GAMG;AACH,MAAM,UAAU,mBAAmB,CACjC,EAAU,EACV,EAAU,EACV,EAAU,EACV,EAAU,EACV,CAAS;IAET,+BAA+B;IAC/B,MAAM,WAAW,GAAG,eAAe,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC;IAEpD,0DAA0D;IAC1D,MAAM,UAAU,GAAG,IAAI,CAAC,IAAI,CAAC,WAAW,CAAC,CAAC,GAAG,WAAW,CAAC,CAAC,GAAG,WAAW,CAAC,CAAC,GAAG,WAAW,CAAC,CAAC,CAAC,CAAC;IAC5F,IAAI,UAAU,GAAG,KAAK,EAAE,CAAC;QACvB,OAAO,EAAE,CAAC,EAAE,EAAE,EAAE,CAAC,EAAE,EAAE,EAAE,CAAC;IAC1B,CAAC;IAED,mDAAmD;IACnD,MAAM,WAAW,GAAG,IAAI,CAAC,GAAG,CAAC,UAAU,EAAE,CAAC,GAAG,KAAK,CAAC,CAAC;IAEpD,+CAA+C;IAC/C,MAAM,UAAU,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,WAAW,CAAC,CAAC,CAAC;IAC1D,MAAM,MAAM,GAAG;QACb,CAAC,EAAE,CAAC,WAAW,CAAC,CAAC,GAAG,UAAU,CAAC,GAAG,UAAU;QAC5C,CAAC,EAAE,CAAC,WAAW,CAAC,CAAC,GAAG,UAAU,CAAC,GAAG,UAAU;KAC7C,CAAC;IAEF,iBAAiB;IACjB,OAAO,sBAAsB,CAAC,MAAM,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC;AAC5D,CAAC"}

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

@@ -0,0 +1,18 @@
+/**
+ * Seeded pseudo-random number generator for deterministic datasets.
+ * Uses a simple xoshiro128** algorithm.
+ */
+export declare class SeededRNG {
+    private state;
+    constructor(seed: number);
+    /** Returns a float in [0, 1) */
+    random(): number;
+    /** Returns a float in [min, max) */
+    range(min: number, max: number): number;
+    /** Returns an integer in [0, max) */
+    int(max: number): number;
+    /** Returns a normally distributed value (Box-Muller) */
+    gaussian(mean?: number, std?: number): number;
+    private rotl;
+}
+//# sourceMappingURL=rng.d.ts.map

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

@@ -0,0 +1 @@
+{"version":3,"file":"rng.d.ts","sourceRoot":"","sources":["../../src/core/rng.ts"],"names":[],"mappings":"AAAA;;;GAGG;AAEH,qBAAa,SAAS;IACpB,OAAO,CAAC,KAAK,CAAc;gBAEf,IAAI,EAAE,MAAM;IAcxB,gCAAgC;IAChC,MAAM,IAAI,MAAM;IAchB,oCAAoC;IACpC,KAAK,CAAC,GAAG,EAAE,MAAM,EAAE,GAAG,EAAE,MAAM,GAAG,MAAM;IAIvC,qCAAqC;IACrC,GAAG,CAAC,GAAG,EAAE,MAAM,GAAG,MAAM;IAIxB,wDAAwD;IACxD,QAAQ,CAAC,IAAI,SAAI,EAAE,GAAG,SAAI,GAAG,MAAM;IAQnC,OAAO,CAAC,IAAI;CAGb"}

package/dist-lib/core/rng.js

@@ -0,0 +1,52 @@
+/**
+ * Seeded pseudo-random number generator for deterministic datasets.
+ * Uses a simple xoshiro128** algorithm.
+ */
+export class SeededRNG {
+    state;
+    constructor(seed) {
+        // Initialize state from seed using splitmix64
+        this.state = new Uint32Array(4);
+        let s = seed >>> 0;
+        for (let i = 0; i < 4; i++) {
+            s = (s + 0x9e3779b9) >>> 0;
+            let z = s;
+            z = ((z ^ (z >>> 16)) * 0x85ebca6b) >>> 0;
+            z = ((z ^ (z >>> 13)) * 0xc2b2ae35) >>> 0;
+            z = (z ^ (z >>> 16)) >>> 0;
+            this.state[i] = z;
+        }
+    }
+    /** Returns a float in [0, 1) */
+    random() {
+        const result = this.rotl(this.state[1] * 5, 7) * 9;
+        const t = this.state[1] << 9;
+        this.state[2] ^= this.state[0];
+        this.state[3] ^= this.state[1];
+        this.state[1] ^= this.state[2];
+        this.state[0] ^= this.state[3];
+        this.state[2] ^= t;
+        this.state[3] = this.rotl(this.state[3], 11);
+        return (result >>> 0) / 0x100000000;
+    }
+    /** Returns a float in [min, max) */
+    range(min, max) {
+        return min + this.random() * (max - min);
+    }
+    /** Returns an integer in [0, max) */
+    int(max) {
+        return Math.floor(this.random() * max);
+    }
+    /** Returns a normally distributed value (Box-Muller) */
+    gaussian(mean = 0, std = 1) {
+        // Clamp u1 away from 0 to avoid log(0) = -Infinity
+        const u1 = Math.max(this.random(), 1e-10);
+        const u2 = this.random();
+        const z = Math.sqrt(-2 * Math.log(u1)) * Math.cos(2 * Math.PI * u2);
+        return mean + z * std;
+    }
+    rotl(x, k) {
+        return ((x << k) | (x >>> (32 - k))) >>> 0;
+    }
+}
+//# sourceMappingURL=rng.js.map

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

@@ -0,0 +1 @@
+{"version":3,"file":"rng.js","sourceRoot":"","sources":["../../src/core/rng.ts"],"names":[],"mappings":"AAAA;;;GAGG;AAEH,MAAM,OAAO,SAAS;IACZ,KAAK,CAAc;IAE3B,YAAY,IAAY;QACtB,8CAA8C;QAC9C,IAAI,CAAC,KAAK,GAAG,IAAI,WAAW,CAAC,CAAC,CAAC,CAAC;QAChC,IAAI,CAAC,GAAG,IAAI,KAAK,CAAC,CAAC;QACnB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;YAC3B,CAAC,GAAG,CAAC,CAAC,GAAG,UAAU,CAAC,KAAK,CAAC,CAAC;YAC3B,IAAI,CAAC,GAAG,CAAC,CAAC;YACV,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,KAAK,EAAE,CAAC,CAAC,GAAG,UAAU,CAAC,KAAK,CAAC,CAAC;YAC1C,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,KAAK,EAAE,CAAC,CAAC,GAAG,UAAU,CAAC,KAAK,CAAC,CAAC;YAC1C,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,KAAK,EAAE,CAAC,CAAC,KAAK,CAAC,CAAC;YAC3B,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;QACpB,CAAC;IACH,CAAC;IAED,gCAAgC;IAChC,MAAM;QACJ,MAAM,MAAM,GAAG,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC;QACnD,MAAM,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC;QAE7B,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,IAAI,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC;QAC/B,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,IAAI,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC;QAC/B,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,IAAI,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC;QAC/B,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,IAAI,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC;QAC/B,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC;QACnB,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC;QAE7C,OAAO,CAAC,MAAM,KAAK,CAAC,CAAC,GAAG,WAAW,CAAC;IACtC,CAAC;IAED,oCAAoC;IACpC,KAAK,CAAC,GAAW,EAAE,GAAW;QAC5B,OAAO,GAAG,GAAG,IAAI,CAAC,MAAM,EAAE,GAAG,CAAC,GAAG,GAAG,GAAG,CAAC,CAAC;IAC3C,CAAC;IAED,qCAAqC;IACrC,GAAG,CAAC,GAAW;QACb,OAAO,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,MAAM,EAAE,GAAG,GAAG,CAAC,CAAC;IACzC,CAAC;IAED,wDAAwD;IACxD,QAAQ,CAAC,IAAI,GAAG,CAAC,EAAE,GAAG,GAAG,CAAC;QACxB,mDAAmD;QACnD,MAAM,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,MAAM,EAAE,EAAE,KAAK,CAAC,CAAC;QAC1C,MAAM,EAAE,GAAG,IAAI,CAAC,MAAM,EAAE,CAAC;QACzB,MAAM,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,GAAG,EAAE,CAAC,CAAC;QACpE,OAAO,IAAI,GAAG,CAAC,GAAG,GAAG,CAAC;IACxB,CAAC;IAEO,IAAI,CAAC,CAAS,EAAE,CAAS;QAC/B,OAAO,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,KAAK,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC;IAC7C,CAAC;CACF"}

package/dist-lib/core/selection/point_in_polygon.d.ts

@@ -0,0 +1,30 @@
+/**
+ * Point-in-polygon test using ray casting algorithm.
+ * This is the naive, accurate implementation for reference.
+ */
+import type { Bounds2D } from '../types.js';
+/**
+ * Test if a point is inside a polygon using ray casting.
+ * Polygon is given as interleaved x,y coordinates.
+ *
+ * Polygon semantics:
+ * - The polygon is implicitly CLOSED (edge from last vertex to first).
+ * - Vertex order may be CW or CCW (even-odd rule; winding does not matter).
+ * - Boundary rule: points exactly on an edge are considered INSIDE (within tolerance).
+ * - Ray crossing uses a half-open y-interval ((yi > py) !== (yj > py)) to avoid double counts.
+ *
+ * Boundary rule: points exactly on the edge are considered INSIDE.
+ */
+export declare function pointInPolygon(px: number, py: number, polygon: Float32Array): boolean;
+/**
+ * Compute an axis-aligned bounding box (AABB) for an interleaved polygon.
+ *
+ * Returns a finite bounds object. For empty input, returns all zeros.
+ */
+export declare function boundsOfPolygon(polygon: Float32Array): Bounds2D;
+/**
+ * Brute-force lasso selection: test all points against polygon.
+ * Returns indices of points inside the polygon.
+ */
+export declare function lassoSelectBruteForce(positions: Float32Array, polygon: Float32Array): Set<number>;
+//# sourceMappingURL=point_in_polygon.d.ts.map

package/dist-lib/core/selection/point_in_polygon.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"point_in_polygon.d.ts","sourceRoot":"","sources":["../../../src/core/selection/point_in_polygon.ts"],"names":[],"mappings":"AAAA;;;GAGG;AAEH,OAAO,KAAK,EAAE,QAAQ,EAAE,MAAM,aAAa,CAAC;AAE5C;;;;;;;;;;;GAWG;AACH,wBAAgB,cAAc,CAC5B,EAAE,EAAE,MAAM,EACV,EAAE,EAAE,MAAM,EACV,OAAO,EAAE,YAAY,GACpB,OAAO,CA2BT;AAED;;;;GAIG;AACH,wBAAgB,eAAe,CAAC,OAAO,EAAE,YAAY,GAAG,QAAQ,CAoB/D;AA4CD;;;GAGG;AACH,wBAAgB,qBAAqB,CACnC,SAAS,EAAE,YAAY,EACvB,OAAO,EAAE,YAAY,GACpB,GAAG,CAAC,MAAM,CAAC,CAab"}

package/dist-lib/core/selection/point_in_polygon.js

@@ -0,0 +1,112 @@
+/**
+ * Point-in-polygon test using ray casting algorithm.
+ * This is the naive, accurate implementation for reference.
+ */
+/**
+ * Test if a point is inside a polygon using ray casting.
+ * Polygon is given as interleaved x,y coordinates.
+ *
+ * Polygon semantics:
+ * - The polygon is implicitly CLOSED (edge from last vertex to first).
+ * - Vertex order may be CW or CCW (even-odd rule; winding does not matter).
+ * - Boundary rule: points exactly on an edge are considered INSIDE (within tolerance).
+ * - Ray crossing uses a half-open y-interval ((yi > py) !== (yj > py)) to avoid double counts.
+ *
+ * Boundary rule: points exactly on the edge are considered INSIDE.
+ */
+export function pointInPolygon(px, py, polygon) {
+    const n = polygon.length / 2;
+    if (n < 3)
+        return false;
+    let inside = false;
+    for (let i = 0, j = n - 1; i < n; j = i++) {
+        const xi = polygon[i * 2];
+        const yi = polygon[i * 2 + 1];
+        const xj = polygon[j * 2];
+        const yj = polygon[j * 2 + 1];
+        // Check if point is on edge (within tolerance)
+        if (isPointOnSegment(px, py, xi, yi, xj, yj, 1e-9)) {
+            return true; // Boundary counts as inside
+        }
+        // Ray casting
+        if ((yi > py) !== (yj > py)) {
+            const xIntersect = ((xj - xi) * (py - yi)) / (yj - yi) + xi;
+            if (px < xIntersect) {
+                inside = !inside;
+            }
+        }
+    }
+    return inside;
+}
+/**
+ * Compute an axis-aligned bounding box (AABB) for an interleaved polygon.
+ *
+ * Returns a finite bounds object. For empty input, returns all zeros.
+ */
+export function boundsOfPolygon(polygon) {
+    let xMin = Infinity;
+    let yMin = Infinity;
+    let xMax = -Infinity;
+    let yMax = -Infinity;
+    for (let i = 0; i < polygon.length; i += 2) {
+        const x = polygon[i];
+        const y = polygon[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) || !Number.isFinite(xMax) || !Number.isFinite(yMax)) {
+        return { xMin: 0, yMin: 0, xMax: 0, yMax: 0 };
+    }
+    return { xMin, yMin, xMax, yMax };
+}
+/**
+ * Check if a point lies on a line segment (within tolerance).
+ */
+function isPointOnSegment(px, py, x1, y1, x2, y2, tolerance) {
+    // Check if point is within bounding box
+    const minX = Math.min(x1, x2) - tolerance;
+    const maxX = Math.max(x1, x2) + tolerance;
+    const minY = Math.min(y1, y2) - tolerance;
+    const maxY = Math.max(y1, y2) + tolerance;
+    if (px < minX || px > maxX || py < minY || py > maxY) {
+        return false;
+    }
+    // Calculate distance from point to line
+    const dx = x2 - x1;
+    const dy = y2 - y1;
+    const lengthSq = dx * dx + dy * dy;
+    if (lengthSq < tolerance * tolerance) {
+        // Degenerate segment (point)
+        const dist = Math.sqrt((px - x1) * (px - x1) + (py - y1) * (py - y1));
+        return dist <= tolerance;
+    }
+    // Project point onto line
+    const t = Math.max(0, Math.min(1, ((px - x1) * dx + (py - y1) * dy) / lengthSq));
+    const projX = x1 + t * dx;
+    const projY = y1 + t * dy;
+    const dist = Math.sqrt((px - projX) * (px - projX) + (py - projY) * (py - projY));
+    return dist <= tolerance;
+}
+/**
+ * Brute-force lasso selection: test all points against polygon.
+ * Returns indices of points inside the polygon.
+ */
+export function lassoSelectBruteForce(positions, polygon) {
+    const result = new Set();
+    const n = positions.length / 2;
+    for (let i = 0; i < n; i++) {
+        const x = positions[i * 2];
+        const y = positions[i * 2 + 1];
+        if (pointInPolygon(x, y, polygon)) {
+            result.add(i);
+        }
+    }
+    return result;
+}
+//# sourceMappingURL=point_in_polygon.js.map

package/dist-lib/core/selection/point_in_polygon.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"point_in_polygon.js","sourceRoot":"","sources":["../../../src/core/selection/point_in_polygon.ts"],"names":[],"mappings":"AAAA;;;GAGG;AAIH;;;;;;;;;;;GAWG;AACH,MAAM,UAAU,cAAc,CAC5B,EAAU,EACV,EAAU,EACV,OAAqB;IAErB,MAAM,CAAC,GAAG,OAAO,CAAC,MAAM,GAAG,CAAC,CAAC;IAC7B,IAAI,CAAC,GAAG,CAAC;QAAE,OAAO,KAAK,CAAC;IAExB,IAAI,MAAM,GAAG,KAAK,CAAC;IAEnB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,EAAE,CAAC;QAC1C,MAAM,EAAE,GAAG,OAAO,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;QAC1B,MAAM,EAAE,GAAG,OAAO,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC;QAC9B,MAAM,EAAE,GAAG,OAAO,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;QAC1B,MAAM,EAAE,GAAG,OAAO,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC;QAE9B,+CAA+C;QAC/C,IAAI,gBAAgB,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,IAAI,CAAC,EAAE,CAAC;YACnD,OAAO,IAAI,CAAC,CAAC,4BAA4B;QAC3C,CAAC;QAED,cAAc;QACd,IAAI,CAAC,EAAE,GAAG,EAAE,CAAC,KAAK,CAAC,EAAE,GAAG,EAAE,CAAC,EAAE,CAAC;YAC5B,MAAM,UAAU,GAAG,CAAC,CAAC,EAAE,GAAG,EAAE,CAAC,GAAG,CAAC,EAAE,GAAG,EAAE,CAAC,CAAC,GAAG,CAAC,EAAE,GAAG,EAAE,CAAC,GAAG,EAAE,CAAC;YAC5D,IAAI,EAAE,GAAG,UAAU,EAAE,CAAC;gBACpB,MAAM,GAAG,CAAC,MAAM,CAAC;YACnB,CAAC;QACH,CAAC;IACH,CAAC;IAED,OAAO,MAAM,CAAC;AAChB,CAAC;AAED;;;;GAIG;AACH,MAAM,UAAU,eAAe,CAAC,OAAqB;IACnD,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;IAErB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,OAAO,CAAC,MAAM,EAAE,CAAC,IAAI,CAAC,EAAE,CAAC;QAC3C,MAAM,CAAC,GAAG,OAAO,CAAC,CAAC,CAAC,CAAC;QACrB,MAAM,CAAC,GAAG,OAAO,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;QACzB,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;IAED,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,IAAI,CAAC,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,IAAI,CAAC,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,IAAI,CAAC,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,IAAI,CAAC,EAAE,CAAC;QACzG,OAAO,EAAE,IAAI,EAAE,CAAC,EAAE,IAAI,EAAE,CAAC,EAAE,IAAI,EAAE,CAAC,EAAE,IAAI,EAAE,CAAC,EAAE,CAAC;IAChD,CAAC;IAED,OAAO,EAAE,IAAI,EAAE,IAAI,EAAE,IAAI,EAAE,IAAI,EAAE,CAAC;AACpC,CAAC;AAED;;GAEG;AACH,SAAS,gBAAgB,CACvB,EAAU,EACV,EAAU,EACV,EAAU,EACV,EAAU,EACV,EAAU,EACV,EAAU,EACV,SAAiB;IAEjB,wCAAwC;IACxC,MAAM,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,EAAE,EAAE,EAAE,CAAC,GAAG,SAAS,CAAC;IAC1C,MAAM,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,EAAE,EAAE,EAAE,CAAC,GAAG,SAAS,CAAC;IAC1C,MAAM,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,EAAE,EAAE,EAAE,CAAC,GAAG,SAAS,CAAC;IAC1C,MAAM,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,EAAE,EAAE,EAAE,CAAC,GAAG,SAAS,CAAC;IAE1C,IAAI,EAAE,GAAG,IAAI,IAAI,EAAE,GAAG,IAAI,IAAI,EAAE,GAAG,IAAI,IAAI,EAAE,GAAG,IAAI,EAAE,CAAC;QACrD,OAAO,KAAK,CAAC;IACf,CAAC;IAED,wCAAwC;IACxC,MAAM,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC;IACnB,MAAM,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC;IACnB,MAAM,QAAQ,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC;IAEnC,IAAI,QAAQ,GAAG,SAAS,GAAG,SAAS,EAAE,CAAC;QACrC,6BAA6B;QAC7B,MAAM,IAAI,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,EAAE,GAAG,EAAE,CAAC,GAAG,CAAC,EAAE,GAAG,EAAE,CAAC,GAAG,CAAC,EAAE,GAAG,EAAE,CAAC,GAAG,CAAC,EAAE,GAAG,EAAE,CAAC,CAAC,CAAC;QACtE,OAAO,IAAI,IAAI,SAAS,CAAC;IAC3B,CAAC;IAED,0BAA0B;IAC1B,MAAM,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,GAAG,EAAE,CAAC,GAAG,EAAE,GAAG,CAAC,EAAE,GAAG,EAAE,CAAC,GAAG,EAAE,CAAC,GAAG,QAAQ,CAAC,CAAC,CAAC;IACjF,MAAM,KAAK,GAAG,EAAE,GAAG,CAAC,GAAG,EAAE,CAAC;IAC1B,MAAM,KAAK,GAAG,EAAE,GAAG,CAAC,GAAG,EAAE,CAAC;IAE1B,MAAM,IAAI,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,EAAE,GAAG,KAAK,CAAC,GAAG,CAAC,EAAE,GAAG,KAAK,CAAC,GAAG,CAAC,EAAE,GAAG,KAAK,CAAC,GAAG,CAAC,EAAE,GAAG,KAAK,CAAC,CAAC,CAAC;IAClF,OAAO,IAAI,IAAI,SAAS,CAAC;AAC3B,CAAC;AAED;;;GAGG;AACH,MAAM,UAAU,qBAAqB,CACnC,SAAuB,EACvB,OAAqB;IAErB,MAAM,MAAM,GAAG,IAAI,GAAG,EAAU,CAAC;IACjC,MAAM,CAAC,GAAG,SAAS,CAAC,MAAM,GAAG,CAAC,CAAC;IAE/B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;QAC3B,MAAM,CAAC,GAAG,SAAS,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;QAC3B,MAAM,CAAC,GAAG,SAAS,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC;QAC/B,IAAI,cAAc,CAAC,CAAC,EAAE,CAAC,EAAE,OAAO,CAAC,EAAE,CAAC;YAClC,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;QAChB,CAAC;IACH,CAAC;IAED,OAAO,MAAM,CAAC;AAChB,CAAC"}