@johnhenry/math 0.0.0

package/dist/Interval.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"Interval.js","sourceRoot":"","sources":["../src/Interval.ts"],"names":[],"mappings":"AAAA,MAAM,GAAG,GAAG,IAAI,YAAY,CAAC,CAAC,CAAC,CAAC;AAChC,MAAM,GAAG,GAAG,IAAI,cAAc,CAAC,GAAG,CAAC,MAAM,CAAC,CAAC;AAE3C,6NAA6N;AAC7N,SAAS,MAAM,CAAC,KAAa;IAC3B,IAAI,MAAM,CAAC,KAAK,CAAC,KAAK,CAAC,IAAI,KAAK,KAAK,QAAQ;QAAE,OAAO,KAAK,CAAC;IAC5D,IAAI,KAAK,KAAK,CAAC;QAAE,OAAO,MAAM,CAAC,SAAS,CAAC;IACzC,GAAG,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC;IACf,GAAG,CAAC,CAAC,CAAC,GAAI,GAAG,CAAC,CAAC,CAAY,GAAG,CAAC,KAAK,GAAG,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC;IACrD,OAAO,GAAG,CAAC,CAAC,CAAW,CAAC;AAC1B,CAAC;AAED,oLAAoL;AACpL,SAAS,QAAQ,CAAC,KAAa;IAC7B,OAAO,CAAC,MAAM,CAAC,CAAC,KAAK,CAAC,CAAC;AACzB,CAAC;AAED;;;;;;;;;;;;GAYG;AACH,SAAS,OAAO,CAAC,CAAW;IAC1B,OAAO,IAAI,QAAQ,CAAC,QAAQ,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,MAAM,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;AACpD,CAAC;AAED;;;;;;;;GAQG;AACH,MAAM,OAAO,QAAQ;IACV,EAAE,CAAS;IACX,EAAE,CAAS;IAEpB,YAAY,EAAU,EAAE,KAAa,EAAE;QACrC,IAAI,EAAE,GAAG,EAAE;YAAE,MAAM,IAAI,KAAK,CAAC,sBAAsB,EAAE,KAAK,EAAE,GAAG,CAAC,CAAC;QACjE,IAAI,CAAC,EAAE,GAAG,EAAE,CAAC;QACb,IAAI,CAAC,EAAE,GAAG,EAAE,CAAC;IACf,CAAC;IAED,MAAM,CAAC,EAAE,CAAC,EAAU,EAAE,EAAU;QAC9B,OAAO,IAAI,QAAQ,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC;IAC9B,CAAC;IAED,sCAAsC;IACtC,MAAM,CAAC,KAAK,CAAC,CAAS;QACpB,OAAO,IAAI,QAAQ,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;IAC5B,CAAC;IAED,IAAI,KAAK;QACP,OAAO,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,EAAE,CAAC;IAC3B,CAAC;IACD,IAAI,QAAQ;QACV,OAAO,CAAC,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC;IACjC,CAAC;IACD,IAAI,MAAM;QACR,OAAO,IAAI,CAAC,KAAK,GAAG,CAAC,CAAC;IACxB,CAAC;IAED,QAAQ,CAAC,CAAS;QAChB,OAAO,IAAI,CAAC,EAAE,IAAI,CAAC,IAAI,CAAC,IAAI,IAAI,CAAC,EAAE,CAAC;IACtC,CAAC;IACD,QAAQ,CAAC,CAAW;QAClB,OAAO,IAAI,CAAC,EAAE,IAAI,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,EAAE,IAAI,IAAI,CAAC,EAAE,CAAC;IAC5C,CAAC;IAED,GAAG,CAAC,CAAW;QACb,OAAO,OAAO,CAAC,IAAI,QAAQ,CAAC,IAAI,CAAC,EAAE,GAAG,CAAC,CAAC,EAAE,EAAE,IAAI,CAAC,EAAE,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;IAC/D,CAAC;IACD,QAAQ,CAAC,CAAW;QAClB,OAAO,OAAO,CAAC,IAAI,QAAQ,CAAC,IAAI,CAAC,EAAE,GAAG,CAAC,CAAC,EAAE,EAAE,IAAI,CAAC,EAAE,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;IAC/D,CAAC;IACD,qEAAqE;IACrE,MAAM;QACJ,OAAO,IAAI,QAAQ,CAAC,CAAC,IAAI,CAAC,EAAE,EAAE,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC;IAC1C,CAAC;IAED,QAAQ,CAAC,CAAW;QAClB,MAAM,CAAC,GAAG,CAAC,IAAI,CAAC,EAAE,GAAG,CAAC,CAAC,EAAE,EAAE,IAAI,CAAC,EAAE,GAAG,CAAC,CAAC,EAAE,EAAE,IAAI,CAAC,EAAE,GAAG,CAAC,CAAC,EAAE,EAAE,IAAI,CAAC,EAAE,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC;QAC3E,OAAO,OAAO,CAAC,IAAI,QAAQ,CAAC,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;IAC/D,CAAC;IAED,MAAM,CAAC,CAAW;QAChB,IAAI,CAAC,CAAC,QAAQ,CAAC,CAAC,CAAC;YAAE,MAAM,IAAI,KAAK,CAAC,mDAAmD,CAAC,CAAC;QACxF,sEAAsE;QACtE,kEAAkE;QAClE,uEAAuE;QACvE,SAAS;QACT,MAAM,UAAU,GAAG,OAAO,CAAC,IAAI,QAAQ,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,EAAE,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;QAC7D,OAAO,IAAI,CAAC,QAAQ,CAAC,UAAU,CAAC,CAAC;IACnC,CAAC;IAED,uEAAuE;IACvE,IAAI,CAAC,CAAW;QACd,OAAO,IAAI,QAAQ,CAAC,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,EAAE,EAAE,CAAC,CAAC,EAAE,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,EAAE,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;IACxE,CAAC;IAED,oEAAoE;IACpE,SAAS,CAAC,CAAW;QACnB,MAAM,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,EAAE,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC;QACnC,MAAM,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,EAAE,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC;QACnC,OAAO,EAAE,IAAI,EAAE,CAAC,CAAC,CAAC,IAAI,QAAQ,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC;IAChD,CAAC;IAED,GAAG;QACD,IAAI,IAAI,CAAC,EAAE,IAAI,CAAC;YAAE,OAAO,IAAI,CAAC;QAC9B,IAAI,IAAI,CAAC,EAAE,IAAI,CAAC;YAAE,OAAO,IAAI,CAAC,MAAM,EAAE,CAAC;QACvC,OAAO,IAAI,QAAQ,CAAC,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,CAAC,IAAI,CAAC,EAAE,EAAE,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC;IACtD,CAAC;IAED,IAAI;QACF,IAAI,IAAI,CAAC,EAAE,GAAG,CAAC;YAAE,MAAM,IAAI,KAAK,CAAC,2CAA2C,CAAC,CAAC;QAC9E,OAAO,OAAO,CAAC,IAAI,QAAQ,CAAC,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,EAAE,CAAC,EAAE,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;IACvE,CAAC;IACD,GAAG;QACD,OAAO,OAAO,CAAC,IAAI,QAAQ,CAAC,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,EAAE,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;IACrE,CAAC;IACD,GAAG;QACD,IAAI,IAAI,CAAC,EAAE,IAAI,CAAC;YAAE,MAAM,IAAI,KAAK,CAAC,8CAA8C,CAAC,CAAC;QAClF,OAAO,OAAO,CAAC,IAAI,QAAQ,CAAC,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,EAAE,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;IACrE,CAAC;IAED,qJAAqJ;IACrJ,GAAG,CAAC,CAAS;QACX,IAAI,CAAC,MAAM,CAAC,SAAS,CAAC,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC;YAAE,MAAM,IAAI,KAAK,CAAC,+CAA+C,CAAC,CAAC;QACpG,IAAI,CAAC,KAAK,CAAC;YAAE,OAAO,QAAQ,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC;QACtC,MAAM,CAAC,GAAG,IAAI,CAAC,EAAE,IAAI,CAAC,CAAC;QACvB,MAAM,CAAC,GAAG,IAAI,CAAC,EAAE,IAAI,CAAC,CAAC;QACvB,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,IAAI,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC;YAAE,OAAO,OAAO,CAAC,IAAI,QAAQ,CAAC,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;QACrF,OAAO,OAAO,CAAC,IAAI,QAAQ,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;IAC/D,CAAC;IAEO,OAAO,CAAC,CAAwB,EAAE,mBAA2B,EAAE,MAAc;QACnF,MAAM,UAAU,GAAG,CAAC,CAAC,CAAC,IAAI,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC;QAC5C,gFAAgF;QAChF,MAAM,MAAM,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,IAAI,CAAC,EAAE,GAAG,mBAAmB,CAAC,GAAG,MAAM,CAAC,CAAC;QACnE,MAAM,IAAI,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,IAAI,CAAC,EAAE,GAAG,mBAAmB,CAAC,GAAG,MAAM,CAAC,CAAC;QAClE,KAAK,IAAI,CAAC,GAAG,MAAM,EAAE,CAAC,IAAI,IAAI,EAAE,CAAC,EAAE;YAAE,UAAU,CAAC,IAAI,CAAC,CAAC,CAAC,mBAAmB,GAAG,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC;QAC1F,OAAO,OAAO,CAAC,IAAI,QAAQ,CAAC,IAAI,CAAC,GAAG,CAAC,GAAG,UAAU,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,GAAG,UAAU,CAAC,CAAC,CAAC,CAAC;IACjF,CAAC;IAED,GAAG;QACD,IAAI,IAAI,CAAC,KAAK,IAAI,CAAC,GAAG,IAAI,CAAC,EAAE;YAAE,OAAO,IAAI,QAAQ,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;QAC1D,OAAO,IAAI,CAAC,OAAO,CAAC,IAAI,CAAC,GAAG,EAAE,IAAI,CAAC,EAAE,GAAG,CAAC,EAAE,IAAI,CAAC,EAAE,CAAC,CAAC;IACtD,CAAC;IACD,GAAG;QACD,IAAI,IAAI,CAAC,KAAK,IAAI,CAAC,GAAG,IAAI,CAAC,EAAE;YAAE,OAAO,IAAI,QAAQ,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;QAC1D,OAAO,IAAI,CAAC,OAAO,CAAC,IAAI,CAAC,GAAG,EAAE,CAAC,EAAE,IAAI,CAAC,EAAE,CAAC,CAAC;IAC5C,CAAC;IAED,MAAM,CAAC,CAAW;QAChB,OAAO,IAAI,CAAC,EAAE,KAAK,CAAC,CAAC,EAAE,IAAI,IAAI,CAAC,EAAE,KAAK,CAAC,CAAC,EAAE,CAAC;IAC9C,CAAC;IAED,QAAQ;QACN,OAAO,IAAI,IAAI,CAAC,EAAE,KAAK,IAAI,CAAC,EAAE,GAAG,CAAC;IACpC,CAAC;CACF"}

package/dist/Lattice.d.ts

@@ -0,0 +1,48 @@
+/**
+ * Lattice — axial-coordinate neighbor helpers for hexagonal and triangular
+ * grids (part of johnhenry/math#30's "lattice-geometry helpers" item,
+ * upstream for the generalized Wang tile laboratory,
+ * johnhenry/mallory#92). Both the solver (adjacency for constraint
+ * propagation) and the app (drawing) need a shared, tested "direction d's
+ * neighbor" definition per lattice — this is that definition, kept purely
+ * combinatorial (integer coordinate offsets), independent of any pixel
+ * geometry (see {@link Polygon.regular} for that side).
+ */
+/** One of the 6 axial directions on a hex grid, in the standard pointy-top ordering (E, NE, NW, W, SW, SE). */
+export type HexDirection = 0 | 1 | 2 | 3 | 4 | 5;
+/**
+ * The 6 axial-coordinate offsets, indexed by {@link HexDirection}. Opposite
+ * directions are `d` and `(d+3) % 6` — verified: `hexNeighbor` in direction
+ * `d` followed by `hexNeighbor` in direction `(d+3) % 6` returns to the
+ * start cell, for every `d`.
+ */
+export declare const HEX_AXIAL_DIRECTIONS: ReadonlyArray<readonly [number, number]>;
+/** The axial-coordinate neighbor of `(q, r)` in direction `d`. */
+export declare function hexNeighbor(q: number, r: number, d: HexDirection): readonly [number, number];
+/** Every axial neighbor of `(q, r)`, in direction order. */
+export declare function hexNeighbors(q: number, r: number): Array<readonly [number, number]>;
+/**
+ * A triangular lattice cell's orientation: `(x, y)` with `x + y` even is
+ * "up" (a triangle pointing up, vertex at the top edge), odd is "down" —
+ * a lattice of alternating up/down triangles tiling the plane, indexed by
+ * integer `(x, y)` (NOT pixel coordinates; see {@link Polygon.regular} for
+ * that side).
+ */
+export type TriOrientation = "up" | "down";
+export declare function triOrientation(x: number, y: number): TriOrientation;
+/** An "up" triangle's 3 edge-sharing neighbors (all "down"): its left/right neighbors in the same row, and the one above sharing its top edge. */
+export type UpTriDirection = "left" | "right" | "top";
+/** A "down" triangle's 3 edge-sharing neighbors (all "up"): its left/right neighbors in the same row, and the one below sharing its bottom edge. */
+export type DownTriDirection = "left" | "right" | "bottom";
+export type TriDirection = UpTriDirection | DownTriDirection;
+/**
+ * The neighbor of triangular cell `(x, y)` in `direction` — reciprocal
+ * pairs: an "up" cell's `"right"` neighbor's `"left"` is itself, and an
+ * "up" cell's `"top"` neighbor's `"bottom"` is itself (verified before
+ * writing tests). `direction` must match `(x, y)`'s own orientation
+ * (`"top"` is only valid on an "up" cell, `"bottom"` only on "down").
+ */
+export declare function triNeighbor(x: number, y: number, direction: TriDirection): readonly [number, number];
+/** Every edge-sharing neighbor of triangular cell `(x, y)`, in the direction order valid for its own orientation. */
+export declare function triNeighbors(x: number, y: number): Array<readonly [number, number]>;
+//# sourceMappingURL=Lattice.d.ts.map

package/dist/Lattice.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"Lattice.d.ts","sourceRoot":"","sources":["../src/Lattice.ts"],"names":[],"mappings":"AAAA;;;;;;;;;GASG;AAEH,+GAA+G;AAC/G,MAAM,MAAM,YAAY,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;AAEjD;;;;;GAKG;AACH,eAAO,MAAM,oBAAoB,EAAE,aAAa,CAAC,SAAS,CAAC,MAAM,EAAE,MAAM,CAAC,CAOzE,CAAC;AAEF,kEAAkE;AAClE,wBAAgB,WAAW,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,YAAY,GAAG,SAAS,CAAC,MAAM,EAAE,MAAM,CAAC,CAG5F;AAED,4DAA4D;AAC5D,wBAAgB,YAAY,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,GAAG,KAAK,CAAC,SAAS,CAAC,MAAM,EAAE,MAAM,CAAC,CAAC,CAEnF;AAED;;;;;;GAMG;AACH,MAAM,MAAM,cAAc,GAAG,IAAI,GAAG,MAAM,CAAC;AAE3C,wBAAgB,cAAc,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,GAAG,cAAc,CAEnE;AAED,kJAAkJ;AAClJ,MAAM,MAAM,cAAc,GAAG,MAAM,GAAG,OAAO,GAAG,KAAK,CAAC;AACtD,oJAAoJ;AACpJ,MAAM,MAAM,gBAAgB,GAAG,MAAM,GAAG,OAAO,GAAG,QAAQ,CAAC;AAC3D,MAAM,MAAM,YAAY,GAAG,cAAc,GAAG,gBAAgB,CAAC;AAE7D;;;;;;GAMG;AACH,wBAAgB,WAAW,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,SAAS,EAAE,YAAY,GAAG,SAAS,CAAC,MAAM,EAAE,MAAM,CAAC,CAYpG;AAED,qHAAqH;AACrH,wBAAgB,YAAY,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,GAAG,KAAK,CAAC,SAAS,CAAC,MAAM,EAAE,MAAM,CAAC,CAAC,CAKnF"}

package/dist/Lattice.js

@@ -0,0 +1,65 @@
+/**
+ * Lattice — axial-coordinate neighbor helpers for hexagonal and triangular
+ * grids (part of johnhenry/math#30's "lattice-geometry helpers" item,
+ * upstream for the generalized Wang tile laboratory,
+ * johnhenry/mallory#92). Both the solver (adjacency for constraint
+ * propagation) and the app (drawing) need a shared, tested "direction d's
+ * neighbor" definition per lattice — this is that definition, kept purely
+ * combinatorial (integer coordinate offsets), independent of any pixel
+ * geometry (see {@link Polygon.regular} for that side).
+ */
+/**
+ * The 6 axial-coordinate offsets, indexed by {@link HexDirection}. Opposite
+ * directions are `d` and `(d+3) % 6` — verified: `hexNeighbor` in direction
+ * `d` followed by `hexNeighbor` in direction `(d+3) % 6` returns to the
+ * start cell, for every `d`.
+ */
+export const HEX_AXIAL_DIRECTIONS = [
+    [1, 0],
+    [1, -1],
+    [0, -1],
+    [-1, 0],
+    [-1, 1],
+    [0, 1],
+];
+/** The axial-coordinate neighbor of `(q, r)` in direction `d`. */
+export function hexNeighbor(q, r, d) {
+    const [dq, dr] = HEX_AXIAL_DIRECTIONS[d];
+    return [q + dq, r + dr];
+}
+/** Every axial neighbor of `(q, r)`, in direction order. */
+export function hexNeighbors(q, r) {
+    return [0, 1, 2, 3, 4, 5].map((d) => hexNeighbor(q, r, d));
+}
+export function triOrientation(x, y) {
+    return (((x + y) % 2) + 2) % 2 === 0 ? "up" : "down";
+}
+/**
+ * The neighbor of triangular cell `(x, y)` in `direction` — reciprocal
+ * pairs: an "up" cell's `"right"` neighbor's `"left"` is itself, and an
+ * "up" cell's `"top"` neighbor's `"bottom"` is itself (verified before
+ * writing tests). `direction` must match `(x, y)`'s own orientation
+ * (`"top"` is only valid on an "up" cell, `"bottom"` only on "down").
+ */
+export function triNeighbor(x, y, direction) {
+    const orientation = triOrientation(x, y);
+    if (direction === "left")
+        return [x - 1, y];
+    if (direction === "right")
+        return [x + 1, y];
+    if (direction === "top") {
+        if (orientation !== "up")
+            throw new Error(`triNeighbor: "top" is only valid on an "up" cell, (${x},${y}) is "down".`);
+        return [x, y + 1];
+    }
+    if (orientation !== "down")
+        throw new Error(`triNeighbor: "bottom" is only valid on a "down" cell, (${x},${y}) is "up".`);
+    return [x, y - 1];
+}
+/** Every edge-sharing neighbor of triangular cell `(x, y)`, in the direction order valid for its own orientation. */
+export function triNeighbors(x, y) {
+    const orientation = triOrientation(x, y);
+    const directions = orientation === "up" ? ["left", "right", "top"] : ["left", "right", "bottom"];
+    return directions.map((d) => triNeighbor(x, y, d));
+}
+//# sourceMappingURL=Lattice.js.map

package/dist/Lattice.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"Lattice.js","sourceRoot":"","sources":["../src/Lattice.ts"],"names":[],"mappings":"AAAA;;;;;;;;;GASG;AAKH;;;;;GAKG;AACH,MAAM,CAAC,MAAM,oBAAoB,GAA6C;IAC5E,CAAC,CAAC,EAAE,CAAC,CAAC;IACN,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;IACP,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;IACP,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC;IACP,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC;IACP,CAAC,CAAC,EAAE,CAAC,CAAC;CACP,CAAC;AAEF,kEAAkE;AAClE,MAAM,UAAU,WAAW,CAAC,CAAS,EAAE,CAAS,EAAE,CAAe;IAC/D,MAAM,CAAC,EAAE,EAAE,EAAE,CAAC,GAAG,oBAAoB,CAAC,CAAC,CAA8B,CAAC;IACtE,OAAO,CAAC,CAAC,GAAG,EAAE,EAAE,CAAC,GAAG,EAAE,CAAC,CAAC;AAC1B,CAAC;AAED,4DAA4D;AAC5D,MAAM,UAAU,YAAY,CAAC,CAAS,EAAE,CAAS;IAC/C,OAAO,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,WAAW,CAAC,CAAC,EAAE,CAAC,EAAE,CAAiB,CAAC,CAAC,CAAC;AAC7E,CAAC;AAWD,MAAM,UAAU,cAAc,CAAC,CAAS,EAAE,CAAS;IACjD,OAAO,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,MAAM,CAAC;AACvD,CAAC;AAQD;;;;;;GAMG;AACH,MAAM,UAAU,WAAW,CAAC,CAAS,EAAE,CAAS,EAAE,SAAuB;IACvE,MAAM,WAAW,GAAG,cAAc,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;IACzC,IAAI,SAAS,KAAK,MAAM;QAAE,OAAO,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,CAAC;IAC5C,IAAI,SAAS,KAAK,OAAO;QAAE,OAAO,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,CAAC;IAC7C,IAAI,SAAS,KAAK,KAAK,EAAE,CAAC;QACxB,IAAI,WAAW,KAAK,IAAI;YACtB,MAAM,IAAI,KAAK,CAAC,sDAAsD,CAAC,IAAI,CAAC,cAAc,CAAC,CAAC;QAC9F,OAAO,CAAC,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC;IACpB,CAAC;IACD,IAAI,WAAW,KAAK,MAAM;QACxB,MAAM,IAAI,KAAK,CAAC,0DAA0D,CAAC,IAAI,CAAC,YAAY,CAAC,CAAC;IAChG,OAAO,CAAC,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC;AACpB,CAAC;AAED,qHAAqH;AACrH,MAAM,UAAU,YAAY,CAAC,CAAS,EAAE,CAAS;IAC/C,MAAM,WAAW,GAAG,cAAc,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;IACzC,MAAM,UAAU,GACd,WAAW,KAAK,IAAI,CAAC,CAAC,CAAC,CAAC,MAAM,EAAE,OAAO,EAAE,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,MAAM,EAAE,OAAO,EAAE,QAAQ,CAAC,CAAC;IAChF,OAAO,UAAU,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,WAAW,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;AACrD,CAAC"}

package/dist/Logic.d.ts

@@ -0,0 +1,22 @@
+/**
+ * Logic — boolean helper predicates.
+ *
+ * Ported from Mallory's ActionScript `Logic`. The AS3 source did not actually
+ * compile (`And(a*;b:*)` and friends are syntactically invalid); the intent is
+ * clear from the comments and reproduced here as ordinary boolean functions.
+ */
+export declare class Logic {
+    /** Always true, whatever the arguments. */
+    static truth(..._anything: unknown[]): boolean;
+    /** Always false, whatever the arguments. */
+    static falsehood(..._anything: unknown[]): boolean;
+    /** Logical AND. */
+    static and(a: unknown, b: unknown): boolean;
+    /** Logical (inclusive) OR. */
+    static or(a: unknown, b: unknown): boolean;
+    /** Logical exclusive OR. */
+    static xor(a: unknown, b: unknown): boolean;
+    /** Logical NOT. */
+    static not(a: unknown): boolean;
+}
+//# sourceMappingURL=Logic.d.ts.map

package/dist/Logic.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"Logic.d.ts","sourceRoot":"","sources":["../src/Logic.ts"],"names":[],"mappings":"AAAA;;;;;;GAMG;AACH,qBAAa,KAAK;IAChB,2CAA2C;IAC3C,MAAM,CAAC,KAAK,CAAC,GAAG,SAAS,EAAE,OAAO,EAAE,GAAG,OAAO;IAI9C,4CAA4C;IAC5C,MAAM,CAAC,SAAS,CAAC,GAAG,SAAS,EAAE,OAAO,EAAE,GAAG,OAAO;IAIlD,mBAAmB;IACnB,MAAM,CAAC,GAAG,CAAC,CAAC,EAAE,OAAO,EAAE,CAAC,EAAE,OAAO,GAAG,OAAO;IAI3C,8BAA8B;IAC9B,MAAM,CAAC,EAAE,CAAC,CAAC,EAAE,OAAO,EAAE,CAAC,EAAE,OAAO,GAAG,OAAO;IAI1C,4BAA4B;IAC5B,MAAM,CAAC,GAAG,CAAC,CAAC,EAAE,OAAO,EAAE,CAAC,EAAE,OAAO,GAAG,OAAO;IAI3C,mBAAmB;IACnB,MAAM,CAAC,GAAG,CAAC,CAAC,EAAE,OAAO,GAAG,OAAO;CAGhC"}

package/dist/Logic.js

@@ -0,0 +1,34 @@
+/**
+ * Logic — boolean helper predicates.
+ *
+ * Ported from Mallory's ActionScript `Logic`. The AS3 source did not actually
+ * compile (`And(a*;b:*)` and friends are syntactically invalid); the intent is
+ * clear from the comments and reproduced here as ordinary boolean functions.
+ */
+export class Logic {
+    /** Always true, whatever the arguments. */
+    static truth(..._anything) {
+        return true;
+    }
+    /** Always false, whatever the arguments. */
+    static falsehood(..._anything) {
+        return false;
+    }
+    /** Logical AND. */
+    static and(a, b) {
+        return Boolean(a) && Boolean(b);
+    }
+    /** Logical (inclusive) OR. */
+    static or(a, b) {
+        return Boolean(a) || Boolean(b);
+    }
+    /** Logical exclusive OR. */
+    static xor(a, b) {
+        return Boolean(a) !== Boolean(b);
+    }
+    /** Logical NOT. */
+    static not(a) {
+        return !a;
+    }
+}
+//# sourceMappingURL=Logic.js.map

package/dist/Logic.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"Logic.js","sourceRoot":"","sources":["../src/Logic.ts"],"names":[],"mappings":"AAAA;;;;;;GAMG;AACH,MAAM,OAAO,KAAK;IAChB,2CAA2C;IAC3C,MAAM,CAAC,KAAK,CAAC,GAAG,SAAoB;QAClC,OAAO,IAAI,CAAC;IACd,CAAC;IAED,4CAA4C;IAC5C,MAAM,CAAC,SAAS,CAAC,GAAG,SAAoB;QACtC,OAAO,KAAK,CAAC;IACf,CAAC;IAED,mBAAmB;IACnB,MAAM,CAAC,GAAG,CAAC,CAAU,EAAE,CAAU;QAC/B,OAAO,OAAO,CAAC,CAAC,CAAC,IAAI,OAAO,CAAC,CAAC,CAAC,CAAC;IAClC,CAAC;IAED,8BAA8B;IAC9B,MAAM,CAAC,EAAE,CAAC,CAAU,EAAE,CAAU;QAC9B,OAAO,OAAO,CAAC,CAAC,CAAC,IAAI,OAAO,CAAC,CAAC,CAAC,CAAC;IAClC,CAAC;IAED,4BAA4B;IAC5B,MAAM,CAAC,GAAG,CAAC,CAAU,EAAE,CAAU;QAC/B,OAAO,OAAO,CAAC,CAAC,CAAC,KAAK,OAAO,CAAC,CAAC,CAAC,CAAC;IACnC,CAAC;IAED,mBAAmB;IACnB,MAAM,CAAC,GAAG,CAAC,CAAU;QACnB,OAAO,CAAC,CAAC,CAAC;IACZ,CAAC;CACF"}

package/dist/MatrixMath.d.ts

@@ -0,0 +1,71 @@
+import { Vector } from "./Vector.ts";
+import type { Matrix } from "./VectorUtils.ts";
+/**
+ * MatrixMath — numerical linear algebra: decompositions (LU, QR, Cholesky,
+ * eigen, SVD), direct solvers, rank/RREF/null space, least squares, the
+ * pseudo-inverse, and matrix norms. Complements the elementary matrix ops in
+ * {@link RealMath} (determinant, inverse, multiply) with the higher-level
+ * factorisations they lacked.
+ *
+ * Matrices are the library's `Matrix<number>` (a {@link Vector} of row Vectors);
+ * plain `number[][]` is also accepted as input.
+ */
+type Mat = Matrix<number>;
+type MatrixInput = Mat | number[][];
+export interface LUResult {
+    L: Mat;
+    U: Mat;
+    P: Mat;
+    sign: number;
+}
+export interface QRResult {
+    Q: Mat;
+    R: Mat;
+}
+export interface EigenResult {
+    values: Vector<number>;
+    vectors: Mat;
+}
+export interface SVDResult {
+    U: Mat;
+    S: Vector<number>;
+    V: Mat;
+}
+export declare class MatrixMath {
+    /** LU decomposition with partial pivoting: `P·A = L·U`. */
+    static lu(A: MatrixInput): LUResult;
+    /** Solve the linear system `A·x = b` via LU decomposition. */
+    static solve(A: MatrixInput, b: number[] | Vector<number>): Vector<number>;
+    /** Reduced row echelon form (Gauss-Jordan elimination). */
+    static rref(A: MatrixInput): Mat;
+    /** Rank of a matrix (number of nonzero rows in its RREF). */
+    static rank(A: MatrixInput): number;
+    /** A basis for the null space (kernel) of `A`, as a matrix whose rows are basis vectors. */
+    static nullSpace(A: MatrixInput): Mat;
+    /** QR decomposition via modified Gram-Schmidt: `A = Q·R` (Q orthonormal). */
+    static qr(A: MatrixInput): QRResult;
+    /** Cholesky decomposition of a symmetric positive-definite matrix: `A = L·Lᵀ`. */
+    static cholesky(A: MatrixInput): Mat;
+    /**
+     * Eigenvalues and orthonormal eigenvectors of a symmetric matrix via the
+     * cyclic Jacobi rotation algorithm. Values are returned in descending order.
+     */
+    static eigenSymmetric(A: MatrixInput, maxSweeps?: number): EigenResult;
+    /**
+     * Singular value decomposition `A = U·Σ·Vᵀ` (via the eigendecomposition of
+     * `AᵀA`). Returns `U`, the singular values `S` (descending), and `V`.
+     */
+    static svd(A: MatrixInput): SVDResult;
+    /** Least-squares solution of an overdetermined system `A·x ≈ b` (normal equations). */
+    static leastSquares(A: MatrixInput, b: number[] | Vector<number>): Vector<number>;
+    /** The Moore-Penrose pseudo-inverse `A⁺` (via SVD). */
+    static pseudoInverse(A: MatrixInput): Mat;
+    /** The Frobenius norm (root of the sum of squared entries). */
+    static frobeniusNorm(A: MatrixInput): number;
+    /** The spectral (2-)norm: the largest singular value. */
+    static spectralNorm(A: MatrixInput): number;
+    /** The 2-norm condition number `σ_max / σ_min`. */
+    static conditionNumber(A: MatrixInput): number;
+}
+export {};
+//# sourceMappingURL=MatrixMath.d.ts.map

package/dist/MatrixMath.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"MatrixMath.d.ts","sourceRoot":"","sources":["../src/MatrixMath.ts"],"names":[],"mappings":"AAAA,OAAO,EAAE,MAAM,EAAE,MAAM,aAAa,CAAC;AACrC,OAAO,KAAK,EAAE,MAAM,EAAE,MAAM,kBAAkB,CAAC;AAE/C;;;;;;;;;GASG;AAEH,KAAK,GAAG,GAAG,MAAM,CAAC,MAAM,CAAC,CAAC;AAC1B,KAAK,WAAW,GAAG,GAAG,GAAG,MAAM,EAAE,EAAE,CAAC;AAoBpC,MAAM,WAAW,QAAQ;IACvB,CAAC,EAAE,GAAG,CAAC;IACP,CAAC,EAAE,GAAG,CAAC;IACP,CAAC,EAAE,GAAG,CAAC;IACP,IAAI,EAAE,MAAM,CAAC;CACd;AAED,MAAM,WAAW,QAAQ;IACvB,CAAC,EAAE,GAAG,CAAC;IACP,CAAC,EAAE,GAAG,CAAC;CACR;AAED,MAAM,WAAW,WAAW;IAC1B,MAAM,EAAE,MAAM,CAAC,MAAM,CAAC,CAAC;IACvB,OAAO,EAAE,GAAG,CAAC;CACd;AAED,MAAM,WAAW,SAAS;IACxB,CAAC,EAAE,GAAG,CAAC;IACP,CAAC,EAAE,MAAM,CAAC,MAAM,CAAC,CAAC;IAClB,CAAC,EAAE,GAAG,CAAC;CACR;AAED,qBAAa,UAAU;IACrB,2DAA2D;IAC3D,MAAM,CAAC,EAAE,CAAC,CAAC,EAAE,WAAW,GAAG,QAAQ;IAyBnC,8DAA8D;IAC9D,MAAM,CAAC,KAAK,CAAC,CAAC,EAAE,WAAW,EAAE,CAAC,EAAE,MAAM,EAAE,GAAG,MAAM,CAAC,MAAM,CAAC,GAAG,MAAM,CAAC,MAAM,CAAC;IA2B1E,2DAA2D;IAC3D,MAAM,CAAC,IAAI,CAAC,CAAC,EAAE,WAAW,GAAG,GAAG;IA2BhC,6DAA6D;IAC7D,MAAM,CAAC,IAAI,CAAC,CAAC,EAAE,WAAW,GAAG,MAAM;IAKnC,4FAA4F;IAC5F,MAAM,CAAC,SAAS,CAAC,CAAC,EAAE,WAAW,GAAG,GAAG;IAsBrC,6EAA6E;IAC7E,MAAM,CAAC,EAAE,CAAC,CAAC,EAAE,WAAW,GAAG,QAAQ;IAuBnC,kFAAkF;IAClF,MAAM,CAAC,QAAQ,CAAC,CAAC,EAAE,WAAW,GAAG,GAAG;IAkBpC;;;OAGG;IACH,MAAM,CAAC,cAAc,CAAC,CAAC,EAAE,WAAW,EAAE,SAAS,SAAM,GAAG,WAAW;IA4CnE;;;OAGG;IACH,MAAM,CAAC,GAAG,CAAC,CAAC,EAAE,WAAW,GAAG,SAAS;IAoBrC,uFAAuF;IACvF,MAAM,CAAC,YAAY,CAAC,CAAC,EAAE,WAAW,EAAE,CAAC,EAAE,MAAM,EAAE,GAAG,MAAM,CAAC,MAAM,CAAC,GAAG,MAAM,CAAC,MAAM,CAAC;IASjF,uDAAuD;IACvD,MAAM,CAAC,aAAa,CAAC,CAAC,EAAE,WAAW,GAAG,GAAG;IAazC,+DAA+D;IAC/D,MAAM,CAAC,aAAa,CAAC,CAAC,EAAE,WAAW,GAAG,MAAM;IAI5C,yDAAyD;IACzD,MAAM,CAAC,YAAY,CAAC,CAAC,EAAE,WAAW,GAAG,MAAM;IAI3C,mDAAmD;IACnD,MAAM,CAAC,eAAe,CAAC,CAAC,EAAE,WAAW,GAAG,MAAM;CAK/C"}

package/dist/MatrixMath.js

@@ -0,0 +1,302 @@
+import { Vector } from "./Vector.js";
+const EPS = 1e-12;
+const raw = (m) => [...m].map((r) => [...r]);
+const wrap = (a) => Vector.fromArray(a.map((r) => Vector.fromArray(r)));
+const vec = (a) => Vector.fromArray(a);
+const identity = (n) => Array.from({ length: n }, (_, i) => Array.from({ length: n }, (_, j) => (i === j ? 1 : 0)));
+const zeros = (m, n) => Array.from({ length: m }, () => new Array(n).fill(0));
+const transposeRaw = (a) => a[0].map((_, j) => a.map((row) => row[j]));
+const matmul = (a, b) => {
+    const m = a.length;
+    const n = b[0].length;
+    const k = b.length;
+    const out = zeros(m, n);
+    for (let i = 0; i < m; i++)
+        for (let p = 0; p < k; p++)
+            for (let j = 0; j < n; j++)
+                out[i][j] += a[i][p] * b[p][j];
+    return out;
+};
+export class MatrixMath {
+    /** LU decomposition with partial pivoting: `P·A = L·U`. */
+    static lu(A) {
+        const U = raw(A);
+        const n = U.length;
+        const L = identity(n);
+        const P = identity(n);
+        let sign = 1;
+        for (let k = 0; k < n; k++) {
+            let pivot = k;
+            for (let i = k + 1; i < n; i++)
+                if (Math.abs(U[i][k]) > Math.abs(U[pivot][k]))
+                    pivot = i;
+            if (pivot !== k) {
+                [U[k], U[pivot]] = [U[pivot], U[k]];
+                [P[k], P[pivot]] = [P[pivot], P[k]];
+                for (let j = 0; j < k; j++)
+                    [L[k][j], L[pivot][j]] = [L[pivot][j], L[k][j]];
+                sign = -sign;
+            }
+            if (Math.abs(U[k][k]) < EPS)
+                continue; // singular column
+            for (let i = k + 1; i < n; i++) {
+                const f = U[i][k] / U[k][k];
+                L[i][k] = f;
+                for (let j = k; j < n; j++)
+                    U[i][j] -= f * U[k][j];
+            }
+        }
+        return { L: wrap(L), U: wrap(U), P: wrap(P), sign };
+    }
+    /** Solve the linear system `A·x = b` via LU decomposition. */
+    static solve(A, b) {
+        const a = raw(A);
+        const n = a.length;
+        const rhs = [...b];
+        const { L, U, P } = MatrixMath.lu(a);
+        const Lr = raw(L);
+        const Ur = raw(U);
+        const Pr = raw(P);
+        // Pb
+        const pb = Pr.map((row) => row.reduce((s, v, j) => s + v * rhs[j], 0));
+        // forward solve L y = Pb
+        const y = new Array(n).fill(0);
+        for (let i = 0; i < n; i++) {
+            let s = pb[i];
+            for (let j = 0; j < i; j++)
+                s -= Lr[i][j] * y[j];
+            y[i] = s; // L has unit diagonal
+        }
+        // back solve U x = y
+        const x = new Array(n).fill(0);
+        for (let i = n - 1; i >= 0; i--) {
+            let s = y[i];
+            for (let j = i + 1; j < n; j++)
+                s -= Ur[i][j] * x[j];
+            x[i] = Math.abs(Ur[i][i]) < EPS ? 0 : s / Ur[i][i];
+        }
+        return vec(x);
+    }
+    /** Reduced row echelon form (Gauss-Jordan elimination). */
+    static rref(A) {
+        const a = raw(A);
+        const rows = a.length;
+        const cols = a[0].length;
+        let lead = 0;
+        for (let r = 0; r < rows && lead < cols; r++, lead++) {
+            let i = r;
+            while (lead < cols && Math.abs(a[i][lead]) < EPS) {
+                i++;
+                if (i === rows) {
+                    i = r;
+                    lead++;
+                    if (lead === cols)
+                        return wrap(a);
+                }
+            }
+            [a[i], a[r]] = [a[r], a[i]];
+            const pivot = a[r][lead];
+            for (let j = 0; j < cols; j++)
+                a[r][j] /= pivot;
+            for (let k = 0; k < rows; k++) {
+                if (k === r)
+                    continue;
+                const factor = a[k][lead];
+                for (let j = 0; j < cols; j++)
+                    a[k][j] -= factor * a[r][j];
+            }
+        }
+        return wrap(a);
+    }
+    /** Rank of a matrix (number of nonzero rows in its RREF). */
+    static rank(A) {
+        const r = raw(MatrixMath.rref(A));
+        return r.filter((row) => row.some((v) => Math.abs(v) > 1e-9)).length;
+    }
+    /** A basis for the null space (kernel) of `A`, as a matrix whose rows are basis vectors. */
+    static nullSpace(A) {
+        const R = raw(MatrixMath.rref(A));
+        const rows = R.length;
+        const cols = R[0].length;
+        const pivotCols = [];
+        for (let r = 0; r < rows; r++) {
+            const lead = R[r].findIndex((v) => Math.abs(v) > 1e-9);
+            if (lead !== -1)
+                pivotCols.push(lead);
+        }
+        const freeCols = [];
+        for (let c = 0; c < cols; c++)
+            if (!pivotCols.includes(c))
+                freeCols.push(c);
+        const basis = [];
+        for (const free of freeCols) {
+            const v = new Array(cols).fill(0);
+            v[free] = 1;
+            for (let r = 0; r < pivotCols.length; r++)
+                v[pivotCols[r]] = -R[r][free];
+            basis.push(v);
+        }
+        return wrap(basis.length ? basis : [new Array(cols).fill(0)]);
+    }
+    /** QR decomposition via modified Gram-Schmidt: `A = Q·R` (Q orthonormal). */
+    static qr(A) {
+        const a = raw(A);
+        const m = a.length;
+        const n = a[0].length;
+        const Q = zeros(m, n);
+        const R = zeros(n, n);
+        const cols = transposeRaw(a); // columns of A
+        for (let j = 0; j < n; j++) {
+            const v = [...cols[j]];
+            for (let i = 0; i < j; i++) {
+                let dot = 0;
+                for (let r = 0; r < m; r++)
+                    dot += Q[r][i] * v[r];
+                R[i][j] = dot;
+                for (let r = 0; r < m; r++)
+                    v[r] -= dot * Q[r][i];
+            }
+            let norm = Math.hypot(...v);
+            R[j][j] = norm;
+            if (norm < EPS)
+                norm = 1;
+            for (let r = 0; r < m; r++)
+                Q[r][j] = v[r] / norm;
+        }
+        return { Q: wrap(Q), R: wrap(R) };
+    }
+    /** Cholesky decomposition of a symmetric positive-definite matrix: `A = L·Lᵀ`. */
+    static cholesky(A) {
+        const a = raw(A);
+        const n = a.length;
+        const L = zeros(n, n);
+        for (let j = 0; j < n; j++) {
+            let d = a[j][j];
+            for (let k = 0; k < j; k++)
+                d -= L[j][k] * L[j][k];
+            if (d <= 0)
+                throw new Error("Matrix is not positive-definite.");
+            L[j][j] = Math.sqrt(d);
+            for (let i = j + 1; i < n; i++) {
+                let s = a[i][j];
+                for (let k = 0; k < j; k++)
+                    s -= L[i][k] * L[j][k];
+                L[i][j] = s / L[j][j];
+            }
+        }
+        return wrap(L);
+    }
+    /**
+     * Eigenvalues and orthonormal eigenvectors of a symmetric matrix via the
+     * cyclic Jacobi rotation algorithm. Values are returned in descending order.
+     */
+    static eigenSymmetric(A, maxSweeps = 100) {
+        const a = raw(A);
+        const n = a.length;
+        const V = identity(n);
+        for (let sweep = 0; sweep < maxSweeps; sweep++) {
+            let off = 0;
+            for (let p = 0; p < n; p++)
+                for (let q = p + 1; q < n; q++)
+                    off += a[p][q] * a[p][q];
+            if (off < EPS)
+                break;
+            for (let p = 0; p < n; p++) {
+                for (let q = p + 1; q < n; q++) {
+                    if (Math.abs(a[p][q]) < EPS)
+                        continue;
+                    const theta = (a[q][q] - a[p][p]) / (2 * a[p][q]);
+                    const t = Math.sign(theta || 1) / (Math.abs(theta) + Math.sqrt(theta * theta + 1));
+                    const c = 1 / Math.sqrt(t * t + 1);
+                    const s = t * c;
+                    for (let i = 0; i < n; i++) {
+                        const aip = a[i][p];
+                        const aiq = a[i][q];
+                        a[i][p] = c * aip - s * aiq;
+                        a[i][q] = s * aip + c * aiq;
+                    }
+                    for (let i = 0; i < n; i++) {
+                        const api = a[p][i];
+                        const aqi = a[q][i];
+                        a[p][i] = c * api - s * aqi;
+                        a[q][i] = s * api + c * aqi;
+                    }
+                    for (let i = 0; i < n; i++) {
+                        const vip = V[i][p];
+                        const viq = V[i][q];
+                        V[i][p] = c * vip - s * viq;
+                        V[i][q] = s * vip + c * viq;
+                    }
+                }
+            }
+        }
+        const pairs = a.map((row, i) => ({ value: row[i], vec: V.map((r) => r[i]) }));
+        pairs.sort((x, y) => y.value - x.value);
+        return {
+            values: vec(pairs.map((p) => p.value)),
+            vectors: wrap(transposeRaw(pairs.map((p) => p.vec))),
+        };
+    }
+    /**
+     * Singular value decomposition `A = U·Σ·Vᵀ` (via the eigendecomposition of
+     * `AᵀA`). Returns `U`, the singular values `S` (descending), and `V`.
+     */
+    static svd(A) {
+        const a = raw(A);
+        const m = a.length;
+        const n = a[0].length;
+        const at = transposeRaw(a);
+        const ata = matmul(at, a); // n×n
+        const { values, vectors } = MatrixMath.eigenSymmetric(ata);
+        const V = raw(vectors); // columns are eigenvectors of AᵀA
+        const sv = [...values].map((v) => Math.sqrt(Math.max(v, 0)));
+        // U columns: A·v_i / σ_i
+        const U = zeros(m, n);
+        for (let i = 0; i < n; i++) {
+            const vi = V.map((r) => r[i]);
+            const avi = a.map((row) => row.reduce((s, val, j) => s + val * vi[j], 0));
+            const sigma = sv[i];
+            for (let r = 0; r < m; r++)
+                U[r][i] = sigma > EPS ? avi[r] / sigma : 0;
+        }
+        return { U: wrap(U), S: vec(sv), V: wrap(V) };
+    }
+    /** Least-squares solution of an overdetermined system `A·x ≈ b` (normal equations). */
+    static leastSquares(A, b) {
+        const a = raw(A);
+        const at = transposeRaw(a);
+        const ata = matmul(at, a);
+        const rhs = [...b];
+        const atb = at.map((row) => row.reduce((s, v, j) => s + v * rhs[j], 0));
+        return MatrixMath.solve(wrap(ata), atb);
+    }
+    /** The Moore-Penrose pseudo-inverse `A⁺` (via SVD). */
+    static pseudoInverse(A) {
+        const { U, S, V } = MatrixMath.svd(A);
+        const Ur = raw(U);
+        const Vr = raw(V);
+        const s = [...S];
+        const n = s.length;
+        // A⁺ = V · Σ⁺ · Uᵀ
+        const sInv = s.map((v) => (v > 1e-9 ? 1 / v : 0));
+        const sigmaPlus = zeros(n, n);
+        for (let i = 0; i < n; i++)
+            sigmaPlus[i][i] = sInv[i];
+        return wrap(matmul(matmul(Vr, sigmaPlus), transposeRaw(Ur)));
+    }
+    /** The Frobenius norm (root of the sum of squared entries). */
+    static frobeniusNorm(A) {
+        return Math.hypot(...raw(A).flat());
+    }
+    /** The spectral (2-)norm: the largest singular value. */
+    static spectralNorm(A) {
+        return Math.max(...MatrixMath.svd(A).S);
+    }
+    /** The 2-norm condition number `σ_max / σ_min`. */
+    static conditionNumber(A) {
+        const s = [...MatrixMath.svd(A).S];
+        const min = Math.min(...s);
+        return min < EPS ? Number.POSITIVE_INFINITY : Math.max(...s) / min;
+    }
+}
+//# sourceMappingURL=MatrixMath.js.map

package/dist/MatrixMath.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"MatrixMath.js","sourceRoot":"","sources":["../src/MatrixMath.ts"],"names":[],"mappings":"AAAA,OAAO,EAAE,MAAM,EAAE,MAAM,aAAa,CAAC;AAiBrC,MAAM,GAAG,GAAG,KAAK,CAAC;AAElB,MAAM,GAAG,GAAG,CAAC,CAAc,EAAc,EAAE,CAAC,CAAC,GAAI,CAAgC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;AACtG,MAAM,IAAI,GAAG,CAAC,CAAa,EAAO,EAAE,CAAC,MAAM,CAAC,SAAS,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,MAAM,CAAC,SAAS,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;AACzF,MAAM,GAAG,GAAG,CAAC,CAAW,EAAkB,EAAE,CAAC,MAAM,CAAC,SAAS,CAAC,CAAC,CAAC,CAAC;AACjE,MAAM,QAAQ,GAAG,CAAC,CAAS,EAAc,EAAE,CACzC,KAAK,CAAC,IAAI,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,KAAK,CAAC,IAAI,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;AAC9F,MAAM,KAAK,GAAG,CAAC,CAAS,EAAE,CAAS,EAAc,EAAE,CAAC,KAAK,CAAC,IAAI,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,EAAE,GAAG,EAAE,CAAC,IAAI,KAAK,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC;AAC1G,MAAM,YAAY,GAAG,CAAC,CAAa,EAAc,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;AAC/F,MAAM,MAAM,GAAG,CAAC,CAAa,EAAE,CAAa,EAAc,EAAE;IAC1D,MAAM,CAAC,GAAG,CAAC,CAAC,MAAM,CAAC;IACnB,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,MAAM,CAAC;IACtB,MAAM,CAAC,GAAG,CAAC,CAAC,MAAM,CAAC;IACnB,MAAM,GAAG,GAAG,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;IACxB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE;QAAE,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE;YAAE,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE;gBAAE,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IACnH,OAAO,GAAG,CAAC;AACb,CAAC,CAAC;AAyBF,MAAM,OAAO,UAAU;IACrB,2DAA2D;IAC3D,MAAM,CAAC,EAAE,CAAC,CAAc;QACtB,MAAM,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC;QACjB,MAAM,CAAC,GAAG,CAAC,CAAC,MAAM,CAAC;QACnB,MAAM,CAAC,GAAG,QAAQ,CAAC,CAAC,CAAC,CAAC;QACtB,MAAM,CAAC,GAAG,QAAQ,CAAC,CAAC,CAAC,CAAC;QACtB,IAAI,IAAI,GAAG,CAAC,CAAC;QACb,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;YAC3B,IAAI,KAAK,GAAG,CAAC,CAAC;YACd,KAAK,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE;gBAAE,IAAI,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC;oBAAE,KAAK,GAAG,CAAC,CAAC;YACzF,IAAI,KAAK,KAAK,CAAC,EAAE,CAAC;gBAChB,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,KAAK,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,KAAK,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;gBACpC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,KAAK,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,KAAK,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;gBACpC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE;oBAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;gBAC5E,IAAI,GAAG,CAAC,IAAI,CAAC;YACf,CAAC;YACD,IAAI,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,GAAG;gBAAE,SAAS,CAAC,kBAAkB;YACzD,KAAK,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;gBAC/B,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;gBAC5B,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;gBACZ,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE;oBAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;YACrD,CAAC;QACH,CAAC;QACD,OAAO,EAAE,CAAC,EAAE,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,IAAI,CAAC,CAAC,CAAC,EAAE,IAAI,EAAE,CAAC;IACtD,CAAC;IAED,8DAA8D;IAC9D,MAAM,CAAC,KAAK,CAAC,CAAc,EAAE,CAA4B;QACvD,MAAM,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC;QACjB,MAAM,CAAC,GAAG,CAAC,CAAC,MAAM,CAAC;QACnB,MAAM,GAAG,GAAG,CAAC,GAAI,CAAsB,CAAC,CAAC;QACzC,MAAM,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,GAAG,UAAU,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;QACrC,MAAM,EAAE,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC;QAClB,MAAM,EAAE,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC;QAClB,MAAM,EAAE,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC;QAClB,KAAK;QACL,MAAM,EAAE,GAAG,EAAE,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,EAAE,CAAC,GAAG,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;QACvE,yBAAyB;QACzB,MAAM,CAAC,GAAG,IAAI,KAAK,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;QAC/B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;YAC3B,IAAI,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC;YACd,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE;gBAAE,CAAC,IAAI,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;YACjD,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,sBAAsB;QAClC,CAAC;QACD,qBAAqB;QACrB,MAAM,CAAC,GAAG,IAAI,KAAK,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;QAC/B,KAAK,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;YAChC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;YACb,KAAK,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE;gBAAE,CAAC,IAAI,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;YACrD,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QACrD,CAAC;QACD,OAAO,GAAG,CAAC,CAAC,CAAC,CAAC;IAChB,CAAC;IAED,2DAA2D;IAC3D,MAAM,CAAC,IAAI,CAAC,CAAc;QACxB,MAAM,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC;QACjB,MAAM,IAAI,GAAG,CAAC,CAAC,MAAM,CAAC;QACtB,MAAM,IAAI,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,MAAM,CAAC;QACzB,IAAI,IAAI,GAAG,CAAC,CAAC;QACb,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,IAAI,IAAI,GAAG,IAAI,EAAE,CAAC,EAAE,EAAE,IAAI,EAAE,EAAE,CAAC;YACrD,IAAI,CAAC,GAAG,CAAC,CAAC;YACV,OAAO,IAAI,GAAG,IAAI,IAAI,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,GAAG,GAAG,EAAE,CAAC;gBACjD,CAAC,EAAE,CAAC;gBACJ,IAAI,CAAC,KAAK,IAAI,EAAE,CAAC;oBACf,CAAC,GAAG,CAAC,CAAC;oBACN,IAAI,EAAE,CAAC;oBACP,IAAI,IAAI,KAAK,IAAI;wBAAE,OAAO,IAAI,CAAC,CAAC,CAAC,CAAC;gBACpC,CAAC;YACH,CAAC;YACD,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;YAC5B,MAAM,KAAK,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC;YACzB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,EAAE,CAAC,EAAE;gBAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI,KAAK,CAAC;YAChD,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,EAAE,CAAC,EAAE,EAAE,CAAC;gBAC9B,IAAI,CAAC,KAAK,CAAC;oBAAE,SAAS;gBACtB,MAAM,MAAM,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC;gBAC1B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,EAAE,CAAC,EAAE;oBAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI,MAAM,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;YAC7D,CAAC;QACH,CAAC;QACD,OAAO,IAAI,CAAC,CAAC,CAAC,CAAC;IACjB,CAAC;IAED,6DAA6D;IAC7D,MAAM,CAAC,IAAI,CAAC,CAAc;QACxB,MAAM,CAAC,GAAG,GAAG,CAAC,UAAU,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC;QAClC,OAAO,CAAC,CAAC,MAAM,CAAC,CAAC,GAAG,EAAE,EAAE,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,MAAM,CAAC;IACvE,CAAC;IAED,4FAA4F;IAC5F,MAAM,CAAC,SAAS,CAAC,CAAc;QAC7B,MAAM,CAAC,GAAG,GAAG,CAAC,UAAU,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC;QAClC,MAAM,IAAI,GAAG,CAAC,CAAC,MAAM,CAAC;QACtB,MAAM,IAAI,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,MAAM,CAAC;QACzB,MAAM,SAAS,GAAa,EAAE,CAAC;QAC/B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,EAAE,CAAC,EAAE,EAAE,CAAC;YAC9B,MAAM,IAAI,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,SAAS,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC;YACvD,IAAI,IAAI,KAAK,CAAC,CAAC;gBAAE,SAAS,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC;QACxC,CAAC;QACD,MAAM,QAAQ,GAAa,EAAE,CAAC;QAC9B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,EAAE,CAAC,EAAE;YAAE,IAAI,CAAC,SAAS,CAAC,QAAQ,CAAC,CAAC,CAAC;gBAAE,QAAQ,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;QAE5E,MAAM,KAAK,GAAe,EAAE,CAAC;QAC7B,KAAK,MAAM,IAAI,IAAI,QAAQ,EAAE,CAAC;YAC5B,MAAM,CAAC,GAAG,IAAI,KAAK,CAAC,IAAI,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;YAClC,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC;YACZ,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,SAAS,CAAC,MAAM,EAAE,CAAC,EAAE;gBAAE,CAAC,CAAC,SAAS,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC;YACzE,KAAK,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;QAChB,CAAC;QACD,OAAO,IAAI,CAAC,KAAK,CAAC,MAAM,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,IAAI,KAAK,CAAC,IAAI,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IAChE,CAAC;IAED,6EAA6E;IAC7E,MAAM,CAAC,EAAE,CAAC,CAAc;QACtB,MAAM,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC;QACjB,MAAM,CAAC,GAAG,CAAC,CAAC,MAAM,CAAC;QACnB,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,MAAM,CAAC;QACtB,MAAM,CAAC,GAAG,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;QACtB,MAAM,CAAC,GAAG,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;QACtB,MAAM,IAAI,GAAG,YAAY,CAAC,CAAC,CAAC,CAAC,CAAC,eAAe;QAC7C,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;YAC3B,MAAM,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC;YACvB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;gBAC3B,IAAI,GAAG,GAAG,CAAC,CAAC;gBACZ,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE;oBAAE,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;gBAClD,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC;gBACd,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE;oBAAE,CAAC,CAAC,CAAC,CAAC,IAAI,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;YACpD,CAAC;YACD,IAAI,IAAI,GAAG,IAAI,CAAC,KAAK,CAAC,GAAG,CAAC,CAAC,CAAC;YAC5B,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC;YACf,IAAI,IAAI,GAAG,GAAG;gBAAE,IAAI,GAAG,CAAC,CAAC;YACzB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE;gBAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC;QACpD,CAAC;QACD,OAAO,EAAE,CAAC,EAAE,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC;IACpC,CAAC;IAED,kFAAkF;IAClF,MAAM,CAAC,QAAQ,CAAC,CAAc;QAC5B,MAAM,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC;QACjB,MAAM,CAAC,GAAG,CAAC,CAAC,MAAM,CAAC;QACnB,MAAM,CAAC,GAAG,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;QACtB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;YAC3B,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;YAChB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE;gBAAE,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;YACnD,IAAI,CAAC,IAAI,CAAC;gBAAE,MAAM,IAAI,KAAK,CAAC,kCAAkC,CAAC,CAAC;YAChE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;YACvB,KAAK,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;gBAC/B,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;gBAChB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE;oBAAE,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;gBACnD,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;YACxB,CAAC;QACH,CAAC;QACD,OAAO,IAAI,CAAC,CAAC,CAAC,CAAC;IACjB,CAAC;IAED;;;OAGG;IACH,MAAM,CAAC,cAAc,CAAC,CAAc,EAAE,SAAS,GAAG,GAAG;QACnD,MAAM,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC;QACjB,MAAM,CAAC,GAAG,CAAC,CAAC,MAAM,CAAC;QACnB,MAAM,CAAC,GAAG,QAAQ,CAAC,CAAC,CAAC,CAAC;QACtB,KAAK,IAAI,KAAK,GAAG,CAAC,EAAE,KAAK,GAAG,SAAS,EAAE,KAAK,EAAE,EAAE,CAAC;YAC/C,IAAI,GAAG,GAAG,CAAC,CAAC;YACZ,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE;gBAAE,KAAK,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE;oBAAE,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;YACrF,IAAI,GAAG,GAAG,GAAG;gBAAE,MAAM;YACrB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;gBAC3B,KAAK,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;oBAC/B,IAAI,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,GAAG;wBAAE,SAAS;oBACtC,MAAM,KAAK,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;oBAClD,MAAM,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,KAAK,IAAI,CAAC,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,KAAK,GAAG,KAAK,GAAG,CAAC,CAAC,CAAC,CAAC;oBACnF,MAAM,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC;oBACnC,MAAM,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;oBAChB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;wBAC3B,MAAM,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;wBACpB,MAAM,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;wBACpB,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,GAAG,GAAG,CAAC,GAAG,GAAG,CAAC;wBAC5B,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,GAAG,GAAG,CAAC,GAAG,GAAG,CAAC;oBAC9B,CAAC;oBACD,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;wBAC3B,MAAM,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;wBACpB,MAAM,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;wBACpB,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,GAAG,GAAG,CAAC,GAAG,GAAG,CAAC;wBAC5B,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,GAAG,GAAG,CAAC,GAAG,GAAG,CAAC;oBAC9B,CAAC;oBACD,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;wBAC3B,MAAM,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;wBACpB,MAAM,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;wBACpB,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,GAAG,GAAG,CAAC,GAAG,GAAG,CAAC;wBAC5B,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,GAAG,GAAG,CAAC,GAAG,GAAG,CAAC;oBAC9B,CAAC;gBACH,CAAC;YACH,CAAC;QACH,CAAC;QACD,MAAM,KAAK,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,EAAE,KAAK,EAAE,GAAG,CAAC,CAAC,CAAC,EAAE,GAAG,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;QAC9E,KAAK,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,KAAK,GAAG,CAAC,CAAC,KAAK,CAAC,CAAC;QACxC,OAAO;YACL,MAAM,EAAE,GAAG,CAAC,KAAK,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC;YACtC,OAAO,EAAE,IAAI,CAAC,YAAY,CAAC,KAAK,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;SACrD,CAAC;IACJ,CAAC;IAED;;;OAGG;IACH,MAAM,CAAC,GAAG,CAAC,CAAc;QACvB,MAAM,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC;QACjB,MAAM,CAAC,GAAG,CAAC,CAAC,MAAM,CAAC;QACnB,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,MAAM,CAAC;QACtB,MAAM,EAAE,GAAG,YAAY,CAAC,CAAC,CAAC,CAAC;QAC3B,MAAM,GAAG,GAAG,MAAM,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,MAAM;QACjC,MAAM,EAAE,MAAM,EAAE,OAAO,EAAE,GAAG,UAAU,CAAC,cAAc,CAAC,GAAG,CAAC,CAAC;QAC3D,MAAM,CAAC,GAAG,GAAG,CAAC,OAAO,CAAC,CAAC,CAAC,kCAAkC;QAC1D,MAAM,EAAE,GAAG,CAAC,GAAG,MAAM,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;QAC7D,yBAAyB;QACzB,MAAM,CAAC,GAAG,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;QACtB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;YAC3B,MAAM,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;YAC9B,MAAM,GAAG,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,EAAE,CAAC,GAAG,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,GAAG,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,GAAG,GAAG,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;YAC1E,MAAM,KAAK,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC;YACpB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE;gBAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,GAAG,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC;QACzE,CAAC;QACD,OAAO,EAAE,CAAC,EAAE,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,GAAG,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC;IAChD,CAAC;IAED,uFAAuF;IACvF,MAAM,CAAC,YAAY,CAAC,CAAc,EAAE,CAA4B;QAC9D,MAAM,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC;QACjB,MAAM,EAAE,GAAG,YAAY,CAAC,CAAC,CAAC,CAAC;QAC3B,MAAM,GAAG,GAAG,MAAM,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC;QAC1B,MAAM,GAAG,GAAG,CAAC,GAAI,CAAsB,CAAC,CAAC;QACzC,MAAM,GAAG,GAAG,EAAE,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,EAAE,CAAC,GAAG,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;QACxE,OAAO,UAAU,CAAC,KAAK,CAAC,IAAI,CAAC,GAAG,CAAC,EAAE,GAAG,CAAC,CAAC;IAC1C,CAAC;IAED,uDAAuD;IACvD,MAAM,CAAC,aAAa,CAAC,CAAc;QACjC,MAAM,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,GAAG,UAAU,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;QACtC,MAAM,EAAE,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC;QAClB,MAAM,EAAE,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC;QAClB,MAAM,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC;QACjB,MAAM,CAAC,GAAG,CAAC,CAAC,MAAM,CAAC;QACnB,mBAAmB;QACnB,MAAM,IAAI,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QAClD,MAAM,SAAS,GAAG,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;QAC9B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE;YAAE,SAAS,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC;QACtD,OAAO,IAAI,CAAC,MAAM,CAAC,MAAM,CAAC,EAAE,EAAE,SAAS,CAAC,EAAE,YAAY,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;IAC/D,CAAC;IAED,+DAA+D;IAC/D,MAAM,CAAC,aAAa,CAAC,CAAc;QACjC,OAAO,IAAI,CAAC,KAAK,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,IAAI,EAAE,CAAC,CAAC;IACtC,CAAC;IAED,yDAAyD;IACzD,MAAM,CAAC,YAAY,CAAC,CAAc;QAChC,OAAO,IAAI,CAAC,GAAG,CAAC,GAAG,UAAU,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IAC1C,CAAC;IAED,mDAAmD;IACnD,MAAM,CAAC,eAAe,CAAC,CAAc;QACnC,MAAM,CAAC,GAAG,CAAC,GAAG,UAAU,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QACnC,MAAM,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC;QAC3B,OAAO,GAAG,GAAG,GAAG,CAAC,CAAC,CAAC,MAAM,CAAC,iBAAiB,CAAC,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,CAAC;IACrE,CAAC;CACF"}