@johnhenry/math 0.2.0 → 0.3.0

package/dist/EigenGeneral.js

@@ -0,0 +1,327 @@
+import { ComplexNumber } from "./ComplexNumber.js";
+const zeros = (n) => Array.from({ length: n }, () => new Array(n).fill(0));
+export function eigenGeneral(input, maxIterations = 60) {
+    const n = input.length;
+    if (n === 0)
+        return { values: [], vectors: [] };
+    const H = { re: zeros(n), im: zeros(n) };
+    for (let i = 0; i < n; i++) {
+        if (input[i].length !== n)
+            throw new Error("eigen requires a square matrix.");
+        for (let j = 0; j < n; j++) {
+            const v = input[i][j];
+            if (typeof v === "number")
+                H.re[i][j] = v;
+            else {
+                H.re[i][j] = v.re;
+                H.im[i][j] = v.im;
+            }
+            if (!Number.isFinite(H.re[i][j]) || !Number.isFinite(H.im[i][j]))
+                throw new Error("eigen requires finite matrix entries.");
+        }
+    }
+    const Z = { re: zeros(n), im: zeros(n) };
+    for (let i = 0; i < n; i++)
+        Z.re[i][i] = 1;
+    let norm = 0;
+    for (let i = 0; i < n; i++)
+        for (let j = 0; j < n; j++)
+            norm = Math.max(norm, Math.hypot(H.re[i][j], H.im[i][j]));
+    if (norm === 0)
+        norm = 1;
+    hessenberg(H, Z, n);
+    schur(H, Z, n, norm, maxIterations);
+    // eigenvalues & vectors
+    const tNorm = norm * n;
+    const lam = [];
+    for (let k = 0; k < n; k++)
+        lam.push([H.re[k][k], H.im[k][k]]);
+    const realTol = 1e-12 * norm;
+    const realInput = input.every((row) => row.every((v) => typeof v === "number" || v.im === 0));
+    const cols = [];
+    for (let k = 0; k < n; k++) {
+        // back-substitute (T - λ_k I) y = 0 with y_k = 1
+        const yr = new Array(n).fill(0);
+        const yi = new Array(n).fill(0);
+        yr[k] = 1;
+        for (let i = k - 1; i >= 0; i--) {
+            let sr = 0;
+            let si = 0;
+            for (let j = i + 1; j <= k; j++) {
+                sr += H.re[i][j] * yr[j] - H.im[i][j] * yi[j];
+                si += H.re[i][j] * yi[j] + H.im[i][j] * yr[j];
+            }
+            let dr = H.re[i][i] - lam[k][0];
+            let di = H.im[i][i] - lam[k][1];
+            if (Math.hypot(dr, di) < Number.EPSILON * tNorm) {
+                dr = Number.EPSILON * tNorm;
+                di = 0;
+            }
+            const d2 = dr * dr + di * di;
+            yr[i] = -(sr * dr + si * di) / d2;
+            yi[i] = -(si * dr - sr * di) / d2;
+        }
+        // v = Z y
+        const vre = new Array(n).fill(0);
+        const vim = new Array(n).fill(0);
+        for (let i = 0; i < n; i++) {
+            for (let j = 0; j <= k; j++) {
+                vre[i] += Z.re[i][j] * yr[j] - Z.im[i][j] * yi[j];
+                vim[i] += Z.re[i][j] * yi[j] + Z.im[i][j] * yr[j];
+            }
+        }
+        // unit norm, largest component real positive
+        let nrm = 0;
+        let big = 0;
+        let bigIdx = 0;
+        for (let i = 0; i < n; i++) {
+            const m = Math.hypot(vre[i], vim[i]);
+            nrm += m * m;
+            if (m > big * (1 + 1e-12)) {
+                big = m;
+                bigIdx = i;
+            }
+        }
+        nrm = Math.sqrt(nrm);
+        const pr = vre[bigIdx] / big;
+        const pi = -vim[bigIdx] / big; // conj(phase)
+        for (let i = 0; i < n; i++) {
+            const r = (vre[i] * pr - vim[i] * pi) / nrm;
+            const m = (vre[i] * pi + vim[i] * pr) / nrm;
+            vre[i] = r;
+            vim[i] = m;
+        }
+        let lre = lam[k][0];
+        let lim = lam[k][1];
+        if (Math.abs(lim) <= realTol && realInput) {
+            lim = 0;
+            for (let i = 0; i < n; i++)
+                if (Math.abs(vim[i]) <= 1e-10)
+                    vim[i] = 0;
+        }
+        if (lre === 0)
+            lre = 0; // normalise -0
+        cols.push({ lre, lim, vre, vim });
+    }
+    cols.sort((a, b) => (Math.abs(a.lre - b.lre) <= realTol ? b.lim - a.lim : b.lre - a.lre));
+    return {
+        values: cols.map((c) => new ComplexNumber(c.lre, c.lim)),
+        vectors: Array.from({ length: n }, (_, i) => cols.map((c) => new ComplexNumber(c.vre[i], c.vim[i]))),
+    };
+}
+/** Householder reduction to upper Hessenberg form, accumulating the unitary Z (`A = Z H Zᴴ`). */
+function hessenberg(H, Z, n) {
+    for (let k = 0; k < n - 2; k++) {
+        let alpha = 0;
+        for (let i = k + 1; i < n; i++)
+            alpha += H.re[i][k] ** 2 + H.im[i][k] ** 2;
+        alpha = Math.sqrt(alpha);
+        if (alpha === 0)
+            continue;
+        const x0 = Math.hypot(H.re[k + 1][k], H.im[k + 1][k]);
+        // v = x + e^{i arg(x0)} alpha e1
+        const pr = x0 === 0 ? 1 : H.re[k + 1][k] / x0;
+        const pi = x0 === 0 ? 0 : H.im[k + 1][k] / x0;
+        const vr = new Array(n).fill(0);
+        const vi = new Array(n).fill(0);
+        for (let i = k + 1; i < n; i++) {
+            vr[i] = H.re[i][k];
+            vi[i] = H.im[i][k];
+        }
+        vr[k + 1] += pr * alpha;
+        vi[k + 1] += pi * alpha;
+        let vv = 0;
+        for (let i = k + 1; i < n; i++)
+            vv += vr[i] ** 2 + vi[i] ** 2;
+        if (vv === 0)
+            continue;
+        // H <- P H, P = I - 2 v vᴴ / vᴴv   (rows k+1.., all columns)
+        for (let j = 0; j < n; j++) {
+            let sr = 0;
+            let si = 0;
+            for (let i = k + 1; i < n; i++) {
+                sr += vr[i] * H.re[i][j] + vi[i] * H.im[i][j];
+                si += vr[i] * H.im[i][j] - vi[i] * H.re[i][j];
+            }
+            sr = (2 * sr) / vv;
+            si = (2 * si) / vv;
+            for (let i = k + 1; i < n; i++) {
+                H.re[i][j] -= vr[i] * sr - vi[i] * si;
+                H.im[i][j] -= vr[i] * si + vi[i] * sr;
+            }
+        }
+        // H <- H P and Z <- Z P  (columns k+1..)
+        for (const M of [H, Z]) {
+            for (let i = 0; i < n; i++) {
+                let sr = 0;
+                let si = 0;
+                for (let j = k + 1; j < n; j++) {
+                    sr += M.re[i][j] * vr[j] - M.im[i][j] * vi[j];
+                    si += M.re[i][j] * vi[j] + M.im[i][j] * vr[j];
+                }
+                sr = (2 * sr) / vv;
+                si = (2 * si) / vv;
+                for (let j = k + 1; j < n; j++) {
+                    M.re[i][j] -= sr * vr[j] + si * vi[j];
+                    M.im[i][j] -= si * vr[j] - sr * vi[j];
+                }
+            }
+        }
+        for (let i = k + 2; i < n; i++) {
+            H.re[i][k] = 0;
+            H.im[i][k] = 0;
+        }
+    }
+}
+/** Shifted QR iteration on a Hessenberg matrix until upper triangular (complex Schur form). */
+function schur(H, Z, n, norm, maxIterations) {
+    let hi = n - 1;
+    let iter = 0;
+    let total = 0;
+    while (hi > 0) {
+        // find the start of the active unreduced block
+        let l = hi;
+        while (l > 0) {
+            const s = Math.hypot(H.re[l - 1][l - 1], H.im[l - 1][l - 1]) + Math.hypot(H.re[l][l], H.im[l][l]);
+            const sub = Math.hypot(H.re[l][l - 1], H.im[l][l - 1]);
+            if (sub <= Number.EPSILON * (s === 0 ? norm : s)) {
+                H.re[l][l - 1] = 0;
+                H.im[l][l - 1] = 0;
+                break;
+            }
+            l--;
+        }
+        if (l === hi) {
+            hi--;
+            iter = 0;
+            continue;
+        }
+        if (++iter > maxIterations)
+            throw new Error("eigen: QR iteration failed to converge.");
+        if (++total > 1000 * n)
+            throw new Error("eigen: QR iteration failed to converge.");
+        // Wilkinson shift: eigenvalue of the trailing 2x2 closer to H[hi][hi]
+        let mur;
+        let mui;
+        if (iter % 10 === 0) {
+            // exceptional shift
+            mur = H.re[hi][hi] + Math.hypot(H.re[hi][hi - 1], H.im[hi][hi - 1]);
+            mui = H.im[hi][hi];
+        }
+        else {
+            const ar = H.re[hi - 1][hi - 1];
+            const ai = H.im[hi - 1][hi - 1];
+            const br = H.re[hi - 1][hi];
+            const bi = H.im[hi - 1][hi];
+            const cr = H.re[hi][hi - 1];
+            const ci = H.im[hi][hi - 1];
+            const dr = H.re[hi][hi];
+            const di = H.im[hi][hi];
+            // tr/2, disc = sqrt(((a-d)/2)^2 + bc)
+            const hr = (ar - dr) / 2;
+            const hiI = (ai - di) / 2;
+            const qr = hr * hr - hiI * hiI + (br * cr - bi * ci);
+            const qi = 2 * hr * hiI + (br * ci + bi * cr);
+            const qm = Math.hypot(qr, qi);
+            const sr = Math.sqrt((qm + qr) / 2);
+            const si = (qi < 0 ? -1 : 1) * Math.sqrt(Math.max(0, (qm - qr) / 2));
+            const mr = (ar + dr) / 2;
+            const mi = (ai + di) / 2;
+            const e1r = mr + sr;
+            const e1i = mi + si;
+            const e2r = mr - sr;
+            const e2i = mi - si;
+            if (Math.hypot(e1r - dr, e1i - di) <= Math.hypot(e2r - dr, e2i - di)) {
+                mur = e1r;
+                mui = e1i;
+            }
+            else {
+                mur = e2r;
+                mui = e2i;
+            }
+        }
+        for (let i = l; i <= hi; i++) {
+            H.re[i][i] -= mur;
+            H.im[i][i] -= mui;
+        }
+        const rots = []; // [c, sr, si]
+        for (let k = l; k < hi; k++) {
+            const ar = H.re[k][k];
+            const ai = H.im[k][k];
+            const br = H.re[k + 1][k];
+            const bi = H.im[k + 1][k];
+            const am = Math.hypot(ar, ai);
+            const r = Math.hypot(am, Math.hypot(br, bi));
+            let c;
+            let sr;
+            let si;
+            if (r === 0) {
+                c = 1;
+                sr = 0;
+                si = 0;
+            }
+            else if (am === 0) {
+                c = 0;
+                // s = conj(b)/|b|
+                const bm = Math.hypot(br, bi);
+                sr = br / bm;
+                si = -bi / bm;
+            }
+            else {
+                c = am / r;
+                // s = (a/|a|) conj(b) / r
+                const pr = ar / am;
+                const pi = ai / am;
+                sr = (pr * br + pi * bi) / r;
+                si = (pi * br - pr * bi) / r;
+            }
+            rots.push([c, sr, si]);
+            // rows k, k+1 over columns k..n-1: [c s; -conj(s) c]
+            for (let j = k; j < n; j++) {
+                const xr = H.re[k][j];
+                const xi = H.im[k][j];
+                const yr = H.re[k + 1][j];
+                const yi = H.im[k + 1][j];
+                H.re[k][j] = c * xr + (sr * yr - si * yi);
+                H.im[k][j] = c * xi + (sr * yi + si * yr);
+                H.re[k + 1][j] = -(sr * xr + si * xi) + c * yr;
+                H.im[k + 1][j] = -(sr * xi - si * xr) + c * yi;
+            }
+            H.re[k + 1][k] = 0;
+            H.im[k + 1][k] = 0;
+        }
+        for (let idx = 0; idx < rots.length; idx++) {
+            const k = l + idx;
+            const [c, sr, si] = rots[idx];
+            const top = Math.min(hi, k + 1);
+            // columns k, k+1 over rows 0..top (H) and all rows (Z): multiply by Gᴴ = [c -s; conj(s) c]
+            for (const [M, rowEnd] of [
+                [H, top],
+                [Z, n - 1],
+            ]) {
+                for (let i = 0; i <= rowEnd; i++) {
+                    const xr = M.re[i][k];
+                    const xi = M.im[i][k];
+                    const yr = M.re[i][k + 1];
+                    const yi = M.im[i][k + 1];
+                    M.re[i][k] = c * xr + (yr * sr + yi * si);
+                    M.im[i][k] = c * xi + (yi * sr - yr * si);
+                    M.re[i][k + 1] = -(xr * sr - xi * si) + c * yr;
+                    M.im[i][k + 1] = -(xr * si + xi * sr) + c * yi;
+                }
+            }
+        }
+        for (let i = l; i <= hi; i++) {
+            H.re[i][i] += mur;
+            H.im[i][i] += mui;
+        }
+    }
+    // clear numerical dust below the diagonal
+    for (let i = 1; i < n; i++) {
+        for (let j = 0; j < i; j++) {
+            H.re[i][j] = 0;
+            H.im[i][j] = 0;
+        }
+    }
+}
+//# sourceMappingURL=EigenGeneral.js.map

package/dist/EigenGeneral.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"EigenGeneral.js","sourceRoot":"","sources":["../src/EigenGeneral.ts"],"names":[],"mappings":"AAAA,OAAO,EAAE,aAAa,EAAE,MAAM,oBAAoB,CAAC;AAyBnD,MAAM,KAAK,GAAG,CAAC,CAAS,EAAc,EAAE,CAAC,KAAK,CAAC,IAAI,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,EAAE,GAAG,EAAE,CAAC,IAAI,KAAK,CAAS,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC;AAEvG,MAAM,UAAU,YAAY,CAC1B,KAA2D,EAC3D,aAAa,GAAG,EAAE;IAElB,MAAM,CAAC,GAAG,KAAK,CAAC,MAAM,CAAC;IACvB,IAAI,CAAC,KAAK,CAAC;QAAE,OAAO,EAAE,MAAM,EAAE,EAAE,EAAE,OAAO,EAAE,EAAE,EAAE,CAAC;IAChD,MAAM,CAAC,GAAS,EAAE,EAAE,EAAE,KAAK,CAAC,CAAC,CAAC,EAAE,EAAE,EAAE,KAAK,CAAC,CAAC,CAAC,EAAE,CAAC;IAC/C,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;QAC3B,IAAI,KAAK,CAAC,CAAC,CAAC,CAAC,MAAM,KAAK,CAAC;YAAE,MAAM,IAAI,KAAK,CAAC,iCAAiC,CAAC,CAAC;QAC9E,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;YAC3B,MAAM,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;YACtB,IAAI,OAAO,CAAC,KAAK,QAAQ;gBAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;iBACrC,CAAC;gBACJ,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC;gBAClB,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC;YACpB,CAAC;YACD,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;gBAC9D,MAAM,IAAI,KAAK,CAAC,uCAAuC,CAAC,CAAC;QAC7D,CAAC;IACH,CAAC;IACD,MAAM,CAAC,GAAS,EAAE,EAAE,EAAE,KAAK,CAAC,CAAC,CAAC,EAAE,EAAE,EAAE,KAAK,CAAC,CAAC,CAAC,EAAE,CAAC;IAC/C,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE;QAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;IAE3C,IAAI,IAAI,GAAG,CAAC,CAAC;IACb,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,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,EAAE,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IAClH,IAAI,IAAI,KAAK,CAAC;QAAE,IAAI,GAAG,CAAC,CAAC;IAEzB,UAAU,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;IACpB,KAAK,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,IAAI,EAAE,aAAa,CAAC,CAAC;IAEpC,wBAAwB;IACxB,MAAM,KAAK,GAAG,IAAI,GAAG,CAAC,CAAC;IACvB,MAAM,GAAG,GAA4B,EAAE,CAAC;IACxC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE;QAAE,GAAG,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IAC/D,MAAM,OAAO,GAAG,KAAK,GAAG,IAAI,CAAC;IAC7B,MAAM,SAAS,GAAG,KAAK,CAAC,KAAK,CAAC,CAAC,GAAG,EAAE,EAAE,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,OAAO,CAAC,KAAK,QAAQ,IAAI,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC,CAAC;IAE9F,MAAM,IAAI,GAAsE,EAAE,CAAC;IACnF,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;QAC3B,iDAAiD;QACjD,MAAM,EAAE,GAAG,IAAI,KAAK,CAAS,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;QACxC,MAAM,EAAE,GAAG,IAAI,KAAK,CAAS,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;QACxC,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;QACV,KAAK,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;YAChC,IAAI,EAAE,GAAG,CAAC,CAAC;YACX,IAAI,EAAE,GAAG,CAAC,CAAC;YACX,KAAK,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;gBAChC,EAAE,IAAI,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC;gBAC9C,EAAE,IAAI,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC;YAChD,CAAC;YACD,IAAI,EAAE,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;YAChC,IAAI,EAAE,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;YAChC,IAAI,IAAI,CAAC,KAAK,CAAC,EAAE,EAAE,EAAE,CAAC,GAAG,MAAM,CAAC,OAAO,GAAG,KAAK,EAAE,CAAC;gBAChD,EAAE,GAAG,MAAM,CAAC,OAAO,GAAG,KAAK,CAAC;gBAC5B,EAAE,GAAG,CAAC,CAAC;YACT,CAAC;YACD,MAAM,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC;YAC7B,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC,GAAG,EAAE,CAAC;YAClC,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC,GAAG,EAAE,CAAC;QACpC,CAAC;QACD,UAAU;QACV,MAAM,GAAG,GAAG,IAAI,KAAK,CAAS,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;QACzC,MAAM,GAAG,GAAG,IAAI,KAAK,CAAS,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;QACzC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;YAC3B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;gBAC5B,GAAG,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC;gBAClD,GAAG,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC;YACpD,CAAC;QACH,CAAC;QACD,6CAA6C;QAC7C,IAAI,GAAG,GAAG,CAAC,CAAC;QACZ,IAAI,GAAG,GAAG,CAAC,CAAC;QACZ,IAAI,MAAM,GAAG,CAAC,CAAC;QACf,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;YAC3B,MAAM,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;YACrC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC;YACb,IAAI,CAAC,GAAG,GAAG,GAAG,CAAC,CAAC,GAAG,KAAK,CAAC,EAAE,CAAC;gBAC1B,GAAG,GAAG,CAAC,CAAC;gBACR,MAAM,GAAG,CAAC,CAAC;YACb,CAAC;QACH,CAAC;QACD,GAAG,GAAG,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC;QACrB,MAAM,EAAE,GAAG,GAAG,CAAC,MAAM,CAAC,GAAG,GAAG,CAAC;QAC7B,MAAM,EAAE,GAAG,CAAC,GAAG,CAAC,MAAM,CAAC,GAAG,GAAG,CAAC,CAAC,cAAc;QAC7C,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;YAC3B,MAAM,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,EAAE,GAAG,GAAG,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,GAAG,GAAG,CAAC;YAC5C,MAAM,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,EAAE,GAAG,GAAG,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,GAAG,GAAG,CAAC;YAC5C,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;YACX,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;QACb,CAAC;QACD,IAAI,GAAG,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QACpB,IAAI,GAAG,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QACpB,IAAI,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,IAAI,OAAO,IAAI,SAAS,EAAE,CAAC;YAC1C,GAAG,GAAG,CAAC,CAAC;YACR,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE;gBAAE,IAAI,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,IAAI,KAAK;oBAAE,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;QACxE,CAAC;QACD,IAAI,GAAG,KAAK,CAAC;YAAE,GAAG,GAAG,CAAC,CAAC,CAAC,eAAe;QACvC,IAAI,CAAC,IAAI,CAAC,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC,CAAC;IACpC,CAAC;IAED,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,GAAG,CAAC,IAAI,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;IAC1F,OAAO;QACL,MAAM,EAAE,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,IAAI,aAAa,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC;QACxD,OAAO,EAAE,KAAK,CAAC,IAAI,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,IAAI,aAAa,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;KACrG,CAAC;AACJ,CAAC;AAED,iGAAiG;AACjG,SAAS,UAAU,CAAC,CAAO,EAAE,CAAO,EAAE,CAAS;IAC7C,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;QAC/B,IAAI,KAAK,GAAG,CAAC,CAAC;QACd,KAAK,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE;YAAE,KAAK,IAAI,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC;QAC3E,KAAK,GAAG,IAAI,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC;QACzB,IAAI,KAAK,KAAK,CAAC;YAAE,SAAS;QAC1B,MAAM,EAAE,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QACtD,iCAAiC;QACjC,MAAM,EAAE,GAAG,EAAE,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC;QAC9C,MAAM,EAAE,GAAG,EAAE,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC;QAC9C,MAAM,EAAE,GAAG,IAAI,KAAK,CAAS,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;QACxC,MAAM,EAAE,GAAG,IAAI,KAAK,CAAS,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;QACxC,KAAK,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;YAC/B,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;YACnB,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QACrB,CAAC;QACD,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,IAAI,EAAE,GAAG,KAAK,CAAC;QACxB,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,IAAI,EAAE,GAAG,KAAK,CAAC;QACxB,IAAI,EAAE,GAAG,CAAC,CAAC;QACX,KAAK,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE;YAAE,EAAE,IAAI,EAAE,CAAC,CAAC,CAAC,IAAI,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC;QAC9D,IAAI,EAAE,KAAK,CAAC;YAAE,SAAS;QACvB,6DAA6D;QAC7D,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;YAC3B,IAAI,EAAE,GAAG,CAAC,CAAC;YACX,IAAI,EAAE,GAAG,CAAC,CAAC;YACX,KAAK,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;gBAC/B,EAAE,IAAI,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;gBAC9C,EAAE,IAAI,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;YAChD,CAAC;YACD,EAAE,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC,GAAG,EAAE,CAAC;YACnB,EAAE,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC,GAAG,EAAE,CAAC;YACnB,KAAK,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;gBAC/B,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI,EAAE,CAAC,CAAC,CAAC,GAAG,EAAE,GAAG,EAAE,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC;gBACtC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI,EAAE,CAAC,CAAC,CAAC,GAAG,EAAE,GAAG,EAAE,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC;YACxC,CAAC;QACH,CAAC;QACD,yCAAyC;QACzC,KAAK,MAAM,CAAC,IAAI,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC;YACvB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;gBAC3B,IAAI,EAAE,GAAG,CAAC,CAAC;gBACX,IAAI,EAAE,GAAG,CAAC,CAAC;gBACX,KAAK,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;oBAC/B,EAAE,IAAI,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC;oBAC9C,EAAE,IAAI,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC;gBAChD,CAAC;gBACD,EAAE,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC,GAAG,EAAE,CAAC;gBACnB,EAAE,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC,GAAG,EAAE,CAAC;gBACnB,KAAK,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;oBAC/B,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI,EAAE,GAAG,EAAE,CAAC,CAAC,CAAC,GAAG,EAAE,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC;oBACtC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI,EAAE,GAAG,EAAE,CAAC,CAAC,CAAC,GAAG,EAAE,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC;gBACxC,CAAC;YACH,CAAC;QACH,CAAC;QACD,KAAK,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;YAC/B,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;YACf,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;QACjB,CAAC;IACH,CAAC;AACH,CAAC;AAED,+FAA+F;AAC/F,SAAS,KAAK,CAAC,CAAO,EAAE,CAAO,EAAE,CAAS,EAAE,IAAY,EAAE,aAAqB;IAC7E,IAAI,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC;IACf,IAAI,IAAI,GAAG,CAAC,CAAC;IACb,IAAI,KAAK,GAAG,CAAC,CAAC;IACd,OAAO,EAAE,GAAG,CAAC,EAAE,CAAC;QACd,+CAA+C;QAC/C,IAAI,CAAC,GAAG,EAAE,CAAC;QACX,OAAO,CAAC,GAAG,CAAC,EAAE,CAAC;YACb,MAAM,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;YAClG,MAAM,GAAG,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;YACvD,IAAI,GAAG,IAAI,MAAM,CAAC,OAAO,GAAG,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC;gBACjD,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC;gBACnB,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC;gBACnB,MAAM;YACR,CAAC;YACD,CAAC,EAAE,CAAC;QACN,CAAC;QACD,IAAI,CAAC,KAAK,EAAE,EAAE,CAAC;YACb,EAAE,EAAE,CAAC;YACL,IAAI,GAAG,CAAC,CAAC;YACT,SAAS;QACX,CAAC;QACD,IAAI,EAAE,IAAI,GAAG,aAAa;YAAE,MAAM,IAAI,KAAK,CAAC,yCAAyC,CAAC,CAAC;QACvF,IAAI,EAAE,KAAK,GAAG,IAAI,GAAG,CAAC;YAAE,MAAM,IAAI,KAAK,CAAC,yCAAyC,CAAC,CAAC;QAEnF,sEAAsE;QACtE,IAAI,GAAW,CAAC;QAChB,IAAI,GAAW,CAAC;QAChB,IAAI,IAAI,GAAG,EAAE,KAAK,CAAC,EAAE,CAAC;YACpB,oBAAoB;YACpB,GAAG,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC;YACpE,GAAG,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC;QACrB,CAAC;aAAM,CAAC;YACN,MAAM,EAAE,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC;YAChC,MAAM,EAAE,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC;YAChC,MAAM,EAAE,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC;YAC5B,MAAM,EAAE,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC;YAC5B,MAAM,EAAE,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC;YAC5B,MAAM,EAAE,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC;YAC5B,MAAM,EAAE,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC;YACxB,MAAM,EAAE,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC;YACxB,sCAAsC;YACtC,MAAM,EAAE,GAAG,CAAC,EAAE,GAAG,EAAE,CAAC,GAAG,CAAC,CAAC;YACzB,MAAM,GAAG,GAAG,CAAC,EAAE,GAAG,EAAE,CAAC,GAAG,CAAC,CAAC;YAC1B,MAAM,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,GAAG,GAAG,GAAG,GAAG,CAAC,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC,CAAC;YACrD,MAAM,EAAE,GAAG,CAAC,GAAG,EAAE,GAAG,GAAG,GAAG,CAAC,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC,CAAC;YAC9C,MAAM,EAAE,GAAG,IAAI,CAAC,KAAK,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC;YAC9B,MAAM,EAAE,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,EAAE,GAAG,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC;YACpC,MAAM,EAAE,GAAG,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,GAAG,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;YACrE,MAAM,EAAE,GAAG,CAAC,EAAE,GAAG,EAAE,CAAC,GAAG,CAAC,CAAC;YACzB,MAAM,EAAE,GAAG,CAAC,EAAE,GAAG,EAAE,CAAC,GAAG,CAAC,CAAC;YACzB,MAAM,GAAG,GAAG,EAAE,GAAG,EAAE,CAAC;YACpB,MAAM,GAAG,GAAG,EAAE,GAAG,EAAE,CAAC;YACpB,MAAM,GAAG,GAAG,EAAE,GAAG,EAAE,CAAC;YACpB,MAAM,GAAG,GAAG,EAAE,GAAG,EAAE,CAAC;YACpB,IAAI,IAAI,CAAC,KAAK,CAAC,GAAG,GAAG,EAAE,EAAE,GAAG,GAAG,EAAE,CAAC,IAAI,IAAI,CAAC,KAAK,CAAC,GAAG,GAAG,EAAE,EAAE,GAAG,GAAG,EAAE,CAAC,EAAE,CAAC;gBACrE,GAAG,GAAG,GAAG,CAAC;gBACV,GAAG,GAAG,GAAG,CAAC;YACZ,CAAC;iBAAM,CAAC;gBACN,GAAG,GAAG,GAAG,CAAC;gBACV,GAAG,GAAG,GAAG,CAAC;YACZ,CAAC;QACH,CAAC;QAED,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,EAAE,EAAE,CAAC,EAAE,EAAE,CAAC;YAC7B,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI,GAAG,CAAC;YAClB,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI,GAAG,CAAC;QACpB,CAAC;QACD,MAAM,IAAI,GAAoC,EAAE,CAAC,CAAC,cAAc;QAChE,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,EAAE,EAAE,CAAC,EAAE,EAAE,CAAC;YAC5B,MAAM,EAAE,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;YACtB,MAAM,EAAE,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;YACtB,MAAM,EAAE,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;YAC1B,MAAM,EAAE,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;YAC1B,MAAM,EAAE,GAAG,IAAI,CAAC,KAAK,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC;YAC9B,MAAM,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,EAAE,EAAE,IAAI,CAAC,KAAK,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC,CAAC;YAC7C,IAAI,CAAS,CAAC;YACd,IAAI,EAAU,CAAC;YACf,IAAI,EAAU,CAAC;YACf,IAAI,CAAC,KAAK,CAAC,EAAE,CAAC;gBACZ,CAAC,GAAG,CAAC,CAAC;gBACN,EAAE,GAAG,CAAC,CAAC;gBACP,EAAE,GAAG,CAAC,CAAC;YACT,CAAC;iBAAM,IAAI,EAAE,KAAK,CAAC,EAAE,CAAC;gBACpB,CAAC,GAAG,CAAC,CAAC;gBACN,kBAAkB;gBAClB,MAAM,EAAE,GAAG,IAAI,CAAC,KAAK,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC;gBAC9B,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC;gBACb,EAAE,GAAG,CAAC,EAAE,GAAG,EAAE,CAAC;YAChB,CAAC;iBAAM,CAAC;gBACN,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC;gBACX,0BAA0B;gBAC1B,MAAM,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC;gBACnB,MAAM,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC;gBACnB,EAAE,GAAG,CAAC,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC,GAAG,CAAC,CAAC;gBAC7B,EAAE,GAAG,CAAC,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC,GAAG,CAAC,CAAC;YAC/B,CAAC;YACD,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC,CAAC;YACvB,qDAAqD;YACrD,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;gBAC3B,MAAM,EAAE,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;gBACtB,MAAM,EAAE,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;gBACtB,MAAM,EAAE,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;gBAC1B,MAAM,EAAE,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;gBAC1B,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,EAAE,GAAG,CAAC,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC,CAAC;gBAC1C,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,EAAE,GAAG,CAAC,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC,CAAC;gBAC1C,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC,GAAG,CAAC,GAAG,EAAE,CAAC;gBAC/C,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC,GAAG,CAAC,GAAG,EAAE,CAAC;YACjD,CAAC;YACD,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;YACnB,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;QACrB,CAAC;QACD,KAAK,IAAI,GAAG,GAAG,CAAC,EAAE,GAAG,GAAG,IAAI,CAAC,MAAM,EAAE,GAAG,EAAE,EAAE,CAAC;YAC3C,MAAM,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC;YAClB,MAAM,CAAC,CAAC,EAAE,EAAE,EAAE,EAAE,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC;YAC9B,MAAM,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,EAAE,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC;YAChC,2FAA2F;YAC3F,KAAK,MAAM,CAAC,CAAC,EAAE,MAAM,CAAC,IAAI;gBACxB,CAAC,CAAC,EAAE,GAAG,CAAC;gBACR,CAAC,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC;aACc,EAAE,CAAC;gBAC3B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;oBACjC,MAAM,EAAE,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;oBACtB,MAAM,EAAE,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;oBACtB,MAAM,EAAE,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;oBAC1B,MAAM,EAAE,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;oBAC1B,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,EAAE,GAAG,CAAC,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC,CAAC;oBAC1C,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,EAAE,GAAG,CAAC,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC,CAAC;oBAC1C,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC,GAAG,CAAC,GAAG,EAAE,CAAC;oBAC/C,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC,GAAG,CAAC,GAAG,EAAE,CAAC;gBACjD,CAAC;YACH,CAAC;QACH,CAAC;QACD,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,EAAE,EAAE,CAAC,EAAE,EAAE,CAAC;YAC7B,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI,GAAG,CAAC;YAClB,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI,GAAG,CAAC;QACpB,CAAC;IACH,CAAC;IACD,0CAA0C;IAC1C,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;QAC3B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;YAC3B,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;YACf,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;QACjB,CAAC;IACH,CAAC;AACH,CAAC"}

package/dist/Frobenius.d.ts

@@ -0,0 +1,60 @@
+import { Rational } from "./Rational.ts";
+/**
+ * Two "Frobenius" normal forms.
+ *
+ * - {@link frobeniusNormalForm}: the *permutation* normal form of a square
+ *   matrix (Frobenius / Perron-Frobenius block-triangular form). Treat the
+ *   matrix as a weighted digraph (row `i` = node `i`'s out-edges, a non-zero
+ *   entry = an edge); its strongly connected components, put in
+ *   source-to-sink topological order, give a relabelling `P` for which
+ *   `Pᵀ·A·P` is block upper triangular with an irreducible diagonal block per
+ *   component. This is what mallory's matrix panel shows (`MAL-LIN-10`).
+ * - {@link rationalCanonicalForm}: the *similarity* normal form (Frobenius
+ *   normal form in the algebraic sense): `A ~ diag(C(f₁), …, C(f_k))`,
+ *   companion blocks of the invariant factors `f₁ | f₂ | … | f_k`, computed
+ *   exactly over ℚ from the Smith normal form of `xI − A`.
+ */
+export type NumericMatrix = ReadonlyArray<ReadonlyArray<number>>;
+export interface FrobeniusBlock {
+    /** The component's vertices as ORIGINAL matrix indices. */
+    component: number[];
+    /** Row/column range in `permuted` (inclusive start, exclusive end). */
+    start: number;
+    end: number;
+}
+export interface FrobeniusResult {
+    /** `order[i]` = the original index now at permuted position `i`. */
+    order: number[];
+    /** Permutation matrix with `Pᵀ·matrix·P === permuted`. */
+    P: number[][];
+    /** `matrix` reindexed by `order` on both axes: block upper triangular. */
+    permuted: number[][];
+    /** One entry per strongly connected component, source to sink. */
+    blocks: FrobeniusBlock[];
+    /** True when the whole matrix is a single strongly connected component. */
+    irreducible: boolean;
+}
+/**
+ * The block-triangular (permutation) Frobenius normal form of a square matrix.
+ * Diagonal entries are self-loops and never affect the component structure.
+ * Within a component, original index order is the (deterministic) tiebreak.
+ */
+export declare function frobeniusNormalForm(matrix: NumericMatrix): FrobeniusResult;
+export interface RationalCanonicalResult {
+    /** Monic invariant factors `f₁ | f₂ | … | f_k` (non-constant only), each as coefficients lowest degree first. */
+    invariantFactors: Rational[][];
+    /** Block diagonal of the companion matrices of `invariantFactors` (n×n, exact). */
+    form: Rational[][];
+    /** The blocks' sizes (the degrees of the invariant factors). */
+    blockSizes: number[];
+}
+/**
+ * Rational canonical form (Frobenius normal form in the similarity sense) of a
+ * square matrix with integer or {@link Rational} entries, exact over ℚ.
+ * `invariantFactors` multiply to the characteristic polynomial and the last is
+ * the minimal polynomial. Each companion block `C(f)` for
+ * `f = xᵏ + c_{k-1}x^{k-1} + … + c₀` has ones on the subdiagonal and last
+ * column `(-c₀, …, -c_{k-1})ᵀ`. The similarity transform itself is not returned.
+ */
+export declare function rationalCanonicalForm(matrix: ReadonlyArray<ReadonlyArray<number | bigint | Rational>>): RationalCanonicalResult;
+//# sourceMappingURL=Frobenius.d.ts.map

package/dist/Frobenius.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"Frobenius.d.ts","sourceRoot":"","sources":["../src/Frobenius.ts"],"names":[],"mappings":"AACA,OAAO,EAAE,QAAQ,EAAE,MAAM,eAAe,CAAC;AAEzC;;;;;;;;;;;;;;GAcG;AAEH,MAAM,MAAM,aAAa,GAAG,aAAa,CAAC,aAAa,CAAC,MAAM,CAAC,CAAC,CAAC;AAEjE,MAAM,WAAW,cAAc;IAC7B,2DAA2D;IAC3D,SAAS,EAAE,MAAM,EAAE,CAAC;IACpB,uEAAuE;IACvE,KAAK,EAAE,MAAM,CAAC;IACd,GAAG,EAAE,MAAM,CAAC;CACb;AAED,MAAM,WAAW,eAAe;IAC9B,oEAAoE;IACpE,KAAK,EAAE,MAAM,EAAE,CAAC;IAChB,0DAA0D;IAC1D,CAAC,EAAE,MAAM,EAAE,EAAE,CAAC;IACd,0EAA0E;IAC1E,QAAQ,EAAE,MAAM,EAAE,EAAE,CAAC;IACrB,kEAAkE;IAClE,MAAM,EAAE,cAAc,EAAE,CAAC;IACzB,2EAA2E;IAC3E,WAAW,EAAE,OAAO,CAAC;CACtB;AAED;;;;GAIG;AACH,wBAAgB,mBAAmB,CAAC,MAAM,EAAE,aAAa,GAAG,eAAe,CA4B1E;AA2FD,MAAM,WAAW,uBAAuB;IACtC,iHAAiH;IACjH,gBAAgB,EAAE,QAAQ,EAAE,EAAE,CAAC;IAC/B,mFAAmF;IACnF,IAAI,EAAE,QAAQ,EAAE,EAAE,CAAC;IACnB,gEAAgE;IAChE,UAAU,EAAE,MAAM,EAAE,CAAC;CACtB;AAED;;;;;;;GAOG;AACH,wBAAgB,qBAAqB,CACnC,MAAM,EAAE,aAAa,CAAC,aAAa,CAAC,MAAM,GAAG,MAAM,GAAG,QAAQ,CAAC,CAAC,GAC/D,uBAAuB,CAiBzB"}

package/dist/Frobenius.js

@@ -0,0 +1,162 @@
+import { Graph } from "./Graph.js";
+import { Rational } from "./Rational.js";
+/**
+ * The block-triangular (permutation) Frobenius normal form of a square matrix.
+ * Diagonal entries are self-loops and never affect the component structure.
+ * Within a component, original index order is the (deterministic) tiebreak.
+ */
+export function frobeniusNormalForm(matrix) {
+    const n = matrix.length;
+    for (const row of matrix) {
+        if (row.length !== n)
+            throw new Error(`frobeniusNormalForm requires a square matrix (got a ${n}-row matrix with a ${row.length}-entry row).`);
+    }
+    const graph = new Graph(true);
+    for (let i = 0; i < n; i++)
+        graph.addVertex(i);
+    for (let i = 0; i < n; i++)
+        for (let j = 0; j < n; j++)
+            if (matrix[i][j] !== 0)
+                graph.addEdge(i, j, matrix[i][j]);
+    const { graph: dag, components } = graph.condensation();
+    const topo = dag.topologicalSort();
+    if (!topo)
+        throw new Error("frobeniusNormalForm: condensation unexpectedly has a cycle.");
+    const order = [];
+    const blocks = [];
+    for (const c of topo) {
+        const component = [...components[c]].sort((a, b) => a - b);
+        const start = order.length;
+        order.push(...component);
+        blocks.push({ component, start, end: order.length });
+    }
+    const permuted = order.map((oi) => order.map((oj) => matrix[oi][oj]));
+    const P = Array.from({ length: n }, () => new Array(n).fill(0));
+    for (let i = 0; i < n; i++)
+        P[order[i]][i] = 1;
+    return { order, P, permuted, blocks, irreducible: components.length === 1 };
+}
+const trim = (p) => {
+    let n = p.length;
+    while (n > 0 && p[n - 1].isZero())
+        n--;
+    return p.slice(0, n);
+};
+const deg = (p) => p.length - 1;
+const padd = (a, b) => trim(Array.from({ length: Math.max(a.length, b.length) }, (_, i) => (a[i] ?? Rational.Zero).add(b[i] ?? Rational.Zero)));
+const pneg = (a) => a.map((c) => c.negate());
+const pmul = (a, b) => {
+    if (!a.length || !b.length)
+        return [];
+    const out = Array.from({ length: a.length + b.length - 1 }, () => Rational.Zero);
+    for (let i = 0; i < a.length; i++)
+        for (let j = 0; j < b.length; j++)
+            out[i + j] = out[i + j].add(a[i].multiply(b[j]));
+    return trim(out);
+};
+/** Polynomial division with remainder (`b` non-zero). */
+const pdivmod = (a, b) => {
+    let r = [...a];
+    const q = Array.from({ length: Math.max(0, a.length - b.length + 1) }, () => Rational.Zero);
+    const lead = b[b.length - 1];
+    while (r.length >= b.length) {
+        const c = r[r.length - 1].divide(lead);
+        const d = r.length - b.length;
+        q[d] = c;
+        r = trim(r.map((x, i) => (i >= d ? x.subtract(c.multiply(b[i - d])) : x)));
+    }
+    return { q: trim(q), r };
+};
+const monic = (p) => (p.length ? p.map((c) => c.divide(p[p.length - 1])) : p);
+/** Smith normal form diagonal of a polynomial matrix, as monic polynomials (zero entries stay zero). */
+function smithDiagonal(M) {
+    const n = M.length;
+    for (let t = 0; t < n; t++) {
+        for (;;) {
+            // pivot: non-zero entry of least degree in the trailing submatrix
+            let pi = -1;
+            let pj = -1;
+            for (let i = t; i < n; i++) {
+                for (let j = t; j < n; j++) {
+                    if (M[i][j].length && (pi < 0 || deg(M[i][j]) < deg(M[pi][pj]))) {
+                        pi = i;
+                        pj = j;
+                    }
+                }
+            }
+            if (pi < 0)
+                return M.map((row, i) => (i < t ? monic(row[i]) : [])).concat([]);
+            [M[t], M[pi]] = [M[pi], M[t]];
+            for (const row of M)
+                [row[t], row[pj]] = [row[pj], row[t]];
+            const p = M[t][t];
+            let dirty = false;
+            for (let i = t + 1; i < n; i++) {
+                if (!M[i][t].length)
+                    continue;
+                const { q } = pdivmod(M[i][t], p);
+                for (let j = t; j < n; j++)
+                    M[i][j] = padd(M[i][j], pneg(pmul(q, M[t][j])));
+                if (M[i][t].length)
+                    dirty = true;
+            }
+            for (let j = t + 1; j < n; j++) {
+                if (!M[t][j].length)
+                    continue;
+                const { q } = pdivmod(M[t][j], p);
+                for (let i = t; i < n; i++)
+                    M[i][j] = padd(M[i][j], pneg(pmul(q, M[i][t])));
+                if (M[t][j].length)
+                    dirty = true;
+            }
+            if (dirty)
+                continue; // a smaller-degree remainder appeared: re-pivot
+            // pivot must divide everything left; if not, fold the offending row in
+            let bad = -1;
+            for (let i = t + 1; i < n && bad < 0; i++) {
+                for (let j = t + 1; j < n; j++) {
+                    if (M[i][j].length && pdivmod(M[i][j], p).r.length) {
+                        bad = i;
+                        break;
+                    }
+                }
+            }
+            if (bad < 0)
+                break;
+            for (let j = t; j < n; j++)
+                M[t][j] = padd(M[t][j], M[bad][j]);
+        }
+    }
+    return M.map((row, i) => monic(row[i]));
+}
+/**
+ * Rational canonical form (Frobenius normal form in the similarity sense) of a
+ * square matrix with integer or {@link Rational} entries, exact over ℚ.
+ * `invariantFactors` multiply to the characteristic polynomial and the last is
+ * the minimal polynomial. Each companion block `C(f)` for
+ * `f = xᵏ + c_{k-1}x^{k-1} + … + c₀` has ones on the subdiagonal and last
+ * column `(-c₀, …, -c_{k-1})ᵀ`. The similarity transform itself is not returned.
+ */
+export function rationalCanonicalForm(matrix) {
+    const n = matrix.length;
+    for (const row of matrix)
+        if (row.length !== n)
+            throw new Error("rationalCanonicalForm requires a square matrix.");
+    const A = matrix.map((row) => row.map((v) => Rational.from(v)));
+    // xI - A as a polynomial matrix
+    const M = A.map((row, i) => row.map((a, j) => trim(i === j ? [a.negate(), Rational.One] : [a.negate()])));
+    const diag = smithDiagonal(M);
+    const factors = diag.filter((p) => p.length > 1);
+    const form = Array.from({ length: n }, () => Array.from({ length: n }, () => Rational.Zero));
+    let off = 0;
+    for (const f of factors) {
+        const k = deg(f);
+        for (let i = 1; i < k; i++)
+            form[off + i][off + i - 1] = Rational.One;
+        for (let i = 0; i < k; i++)
+            form[off + i][off + k - 1] = f[i].negate();
+        off += k;
+    }
+    return { invariantFactors: factors, form, blockSizes: factors.map(deg) };
+}
+//# sourceMappingURL=Frobenius.js.map

package/dist/Frobenius.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"Frobenius.js","sourceRoot":"","sources":["../src/Frobenius.ts"],"names":[],"mappings":"AAAA,OAAO,EAAE,KAAK,EAAE,MAAM,YAAY,CAAC;AACnC,OAAO,EAAE,QAAQ,EAAE,MAAM,eAAe,CAAC;AAyCzC;;;;GAIG;AACH,MAAM,UAAU,mBAAmB,CAAC,MAAqB;IACvD,MAAM,CAAC,GAAG,MAAM,CAAC,MAAM,CAAC;IACxB,KAAK,MAAM,GAAG,IAAI,MAAM,EAAE,CAAC;QACzB,IAAI,GAAG,CAAC,MAAM,KAAK,CAAC;YAClB,MAAM,IAAI,KAAK,CACb,uDAAuD,CAAC,sBAAsB,GAAG,CAAC,MAAM,cAAc,CACvG,CAAC;IACN,CAAC;IACD,MAAM,KAAK,GAAG,IAAI,KAAK,CAAS,IAAI,CAAC,CAAC;IACtC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE;QAAE,KAAK,CAAC,SAAS,CAAC,CAAC,CAAC,CAAC;IAC/C,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,IAAI,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,KAAK,CAAC;gBAAE,KAAK,CAAC,OAAO,CAAC,CAAC,EAAE,CAAC,EAAE,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IAElH,MAAM,EAAE,KAAK,EAAE,GAAG,EAAE,UAAU,EAAE,GAAG,KAAK,CAAC,YAAY,EAAE,CAAC;IACxD,MAAM,IAAI,GAAG,GAAG,CAAC,eAAe,EAAE,CAAC;IACnC,IAAI,CAAC,IAAI;QAAE,MAAM,IAAI,KAAK,CAAC,6DAA6D,CAAC,CAAC;IAE1F,MAAM,KAAK,GAAa,EAAE,CAAC;IAC3B,MAAM,MAAM,GAAqB,EAAE,CAAC;IACpC,KAAK,MAAM,CAAC,IAAI,IAAI,EAAE,CAAC;QACrB,MAAM,SAAS,GAAG,CAAC,GAAG,UAAU,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;QAC3D,MAAM,KAAK,GAAG,KAAK,CAAC,MAAM,CAAC;QAC3B,KAAK,CAAC,IAAI,CAAC,GAAG,SAAS,CAAC,CAAC;QACzB,MAAM,CAAC,IAAI,CAAC,EAAE,SAAS,EAAE,KAAK,EAAE,GAAG,EAAE,KAAK,CAAC,MAAM,EAAE,CAAC,CAAC;IACvD,CAAC;IACD,MAAM,QAAQ,GAAG,KAAK,CAAC,GAAG,CAAC,CAAC,EAAE,EAAE,EAAE,CAAC,KAAK,CAAC,GAAG,CAAC,CAAC,EAAE,EAAE,EAAE,CAAC,MAAM,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;IACtE,MAAM,CAAC,GAAG,KAAK,CAAC,IAAI,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,EAAE,GAAG,EAAE,CAAC,IAAI,KAAK,CAAS,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC;IACxE,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE;QAAE,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;IAC/C,OAAO,EAAE,KAAK,EAAE,CAAC,EAAE,QAAQ,EAAE,MAAM,EAAE,WAAW,EAAE,UAAU,CAAC,MAAM,KAAK,CAAC,EAAE,CAAC;AAC9E,CAAC;AAOD,MAAM,IAAI,GAAG,CAAC,CAAO,EAAQ,EAAE;IAC7B,IAAI,CAAC,GAAG,CAAC,CAAC,MAAM,CAAC;IACjB,OAAO,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,MAAM,EAAE;QAAE,CAAC,EAAE,CAAC;IACvC,OAAO,CAAC,CAAC,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;AACvB,CAAC,CAAC;AACF,MAAM,GAAG,GAAG,CAAC,CAAO,EAAE,EAAE,CAAC,CAAC,CAAC,MAAM,GAAG,CAAC,CAAC;AACtC,MAAM,IAAI,GAAG,CAAC,CAAO,EAAE,CAAO,EAAQ,EAAE,CACtC,IAAI,CACF,KAAK,CAAC,IAAI,CAAC,EAAE,MAAM,EAAE,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,MAAM,EAAE,CAAC,CAAC,MAAM,CAAC,EAAE,EAAE,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI,QAAQ,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI,QAAQ,CAAC,IAAI,CAAC,CAAC,CACnH,CAAC;AACJ,MAAM,IAAI,GAAG,CAAC,CAAO,EAAQ,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,MAAM,EAAE,CAAC,CAAC;AACzD,MAAM,IAAI,GAAG,CAAC,CAAO,EAAE,CAAO,EAAQ,EAAE;IACtC,IAAI,CAAC,CAAC,CAAC,MAAM,IAAI,CAAC,CAAC,CAAC,MAAM;QAAE,OAAO,EAAE,CAAC;IACtC,MAAM,GAAG,GAAe,KAAK,CAAC,IAAI,CAAC,EAAE,MAAM,EAAE,CAAC,CAAC,MAAM,GAAG,CAAC,CAAC,MAAM,GAAG,CAAC,EAAE,EAAE,GAAG,EAAE,CAAC,QAAQ,CAAC,IAAI,CAAC,CAAC;IAC7F,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,MAAM,EAAE,CAAC,EAAE;QAC/B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,MAAM,EAAE,CAAC,EAAE;YAAE,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IACtF,OAAO,IAAI,CAAC,GAAG,CAAC,CAAC;AACnB,CAAC,CAAC;AACF,yDAAyD;AACzD,MAAM,OAAO,GAAG,CAAC,CAAO,EAAE,CAAO,EAAwB,EAAE;IACzD,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC;IACf,MAAM,CAAC,GAAe,KAAK,CAAC,IAAI,CAAC,EAAE,MAAM,EAAE,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,MAAM,GAAG,CAAC,CAAC,MAAM,GAAG,CAAC,CAAC,EAAE,EAAE,GAAG,EAAE,CAAC,QAAQ,CAAC,IAAI,CAAC,CAAC;IACxG,MAAM,IAAI,GAAG,CAAC,CAAC,CAAC,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC;IAC7B,OAAO,CAAC,CAAC,MAAM,IAAI,CAAC,CAAC,MAAM,EAAE,CAAC;QAC5B,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC,MAAM,CAAC,IAAI,CAAC,CAAC;QACvC,MAAM,CAAC,GAAG,CAAC,CAAC,MAAM,GAAG,CAAC,CAAC,MAAM,CAAC;QAC9B,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;QACT,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,QAAQ,CAAC,CAAC,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IAC7E,CAAC;IACD,OAAO,EAAE,CAAC,EAAE,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC;AAC3B,CAAC,CAAC;AACF,MAAM,KAAK,GAAG,CAAC,CAAO,EAAQ,EAAE,CAAC,CAAC,CAAC,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;AAE1F,wGAAwG;AACxG,SAAS,aAAa,CAAC,CAAW;IAChC,MAAM,CAAC,GAAG,CAAC,CAAC,MAAM,CAAC;IACnB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;QAC3B,SAAS,CAAC;YACR,kEAAkE;YAClE,IAAI,EAAE,GAAG,CAAC,CAAC,CAAC;YACZ,IAAI,EAAE,GAAG,CAAC,CAAC,CAAC;YACZ,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;gBAC3B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;oBAC3B,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,MAAM,IAAI,CAAC,EAAE,GAAG,CAAC,IAAI,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC;wBAChE,EAAE,GAAG,CAAC,CAAC;wBACP,EAAE,GAAG,CAAC,CAAC;oBACT,CAAC;gBACH,CAAC;YACH,CAAC;YACD,IAAI,EAAE,GAAG,CAAC;gBAAE,OAAO,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,KAAK,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,MAAM,CAAC,EAAE,CAAC,CAAC;YAC9E,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;YAC9B,KAAK,MAAM,GAAG,IAAI,CAAC;gBAAE,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;YAC3D,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;YAClB,IAAI,KAAK,GAAG,KAAK,CAAC;YAClB,KAAK,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;gBAC/B,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,MAAM;oBAAE,SAAS;gBAC9B,MAAM,EAAE,CAAC,EAAE,GAAG,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;gBAClC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE;oBAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,IAAI,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;gBAC5E,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,MAAM;oBAAE,KAAK,GAAG,IAAI,CAAC;YACnC,CAAC;YACD,KAAK,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;gBAC/B,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,MAAM;oBAAE,SAAS;gBAC9B,MAAM,EAAE,CAAC,EAAE,GAAG,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;gBAClC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE;oBAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,IAAI,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;gBAC5E,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,MAAM;oBAAE,KAAK,GAAG,IAAI,CAAC;YACnC,CAAC;YACD,IAAI,KAAK;gBAAE,SAAS,CAAC,gDAAgD;YACrE,uEAAuE;YACvE,IAAI,GAAG,GAAG,CAAC,CAAC,CAAC;YACb,KAAK,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,IAAI,GAAG,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;gBAC1C,KAAK,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;oBAC/B,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,MAAM,IAAI,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,MAAM,EAAE,CAAC;wBACnD,GAAG,GAAG,CAAC,CAAC;wBACR,MAAM;oBACR,CAAC;gBACH,CAAC;YACH,CAAC;YACD,IAAI,GAAG,GAAG,CAAC;gBAAE,MAAM;YACnB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE;gBAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QACjE,CAAC;IACH,CAAC;IACD,OAAO,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC,EAAE,EAAE,CAAC,KAAK,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;AAC1C,CAAC;AAWD;;;;;;;GAOG;AACH,MAAM,UAAU,qBAAqB,CACnC,MAAgE;IAEhE,MAAM,CAAC,GAAG,MAAM,CAAC,MAAM,CAAC;IACxB,KAAK,MAAM,GAAG,IAAI,MAAM;QAAE,IAAI,GAAG,CAAC,MAAM,KAAK,CAAC;YAAE,MAAM,IAAI,KAAK,CAAC,iDAAiD,CAAC,CAAC;IACnH,MAAM,CAAC,GAAG,MAAM,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,EAAE,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,QAAQ,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IAChE,gCAAgC;IAChC,MAAM,CAAC,GAAa,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC,EAAE,EAAE,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,IAAI,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,MAAM,EAAE,EAAE,QAAQ,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,MAAM,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;IACpH,MAAM,IAAI,GAAG,aAAa,CAAC,CAAC,CAAC,CAAC;IAC9B,MAAM,OAAO,GAAG,IAAI,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC;IACjD,MAAM,IAAI,GAAiB,KAAK,CAAC,IAAI,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,EAAE,GAAG,EAAE,CAAC,KAAK,CAAC,IAAI,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,EAAE,GAAG,EAAE,CAAC,QAAQ,CAAC,IAAI,CAAC,CAAC,CAAC;IAC3G,IAAI,GAAG,GAAG,CAAC,CAAC;IACZ,KAAK,MAAM,CAAC,IAAI,OAAO,EAAE,CAAC;QACxB,MAAM,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC;QACjB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE;YAAE,IAAI,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,QAAQ,CAAC,GAAG,CAAC;QACtE,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE;YAAE,IAAI,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,MAAM,EAAE,CAAC;QACvE,GAAG,IAAI,CAAC,CAAC;IACX,CAAC;IACD,OAAO,EAAE,gBAAgB,EAAE,OAAO,EAAE,IAAI,EAAE,UAAU,EAAE,OAAO,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,CAAC;AAC3E,CAAC"}

package/dist/Graph.d.ts

@@ -47,6 +47,32 @@ export declare class Graph<T = string | number> {
     hasCycle(): boolean;
     /** Minimum spanning tree (Kruskal) of an undirected graph, as a list of edges. */
     minimumSpanningTree(): Array<Edge<T>>;
+    /**
+     * Minimum spanning tree (or forest, when the graph is disconnected) with
+     * its total weight and the number of connected components. Edge direction
+     * is ignored. `algorithm` picks Kruskal (default, sort + union-find) or
+     * Prim (heap-based, grown from each unvisited vertex in turn); both give
+     * the same total weight, though tie-broken edge choices may differ.
+     */
+    mst(algorithm?: "kruskal" | "prim"): {
+        edges: Array<Edge<T>>;
+        totalWeight: number;
+        components: number;
+    };
+    /** Alias of {@link stronglyConnectedComponents} (Tarjan; components in reverse topological order). */
+    scc(): T[][];
+    /**
+     * The condensation: one vertex per strongly connected component (numbered
+     * by position in {@link scc}'s output), with one weight-1 edge `i -> j`
+     * wherever the original graph has any edge from a member of `i` to a
+     * member of `j`. Always a DAG. Also returns the `components` and the
+     * `componentOf` vertex-to-index map.
+     */
+    condensation(): {
+        graph: Graph<number>;
+        components: T[][];
+        componentOf: Map<T, number>;
+    };
     /** Adjacency matrix (rows/cols follow the returned `order`; `Infinity` = no edge). */
     toAdjacencyMatrix(): {
         matrix: number[][];

package/dist/Graph.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"Graph.d.ts","sourceRoot":"","sources":["../src/Graph.ts"],"names":[],"mappings":"AAAA;;;;;;GAMG;AAEH,MAAM,WAAW,IAAI,CAAC,CAAC;IACrB,IAAI,EAAE,CAAC,CAAC;IACR,EAAE,EAAE,CAAC,CAAC;IACN,MAAM,EAAE,MAAM,CAAC;CAChB;AAoDD,qBAAa,KAAK,CAAC,CAAC,GAAG,MAAM,GAAG,MAAM;IACpC,QAAQ,CAAC,QAAQ,EAAE,OAAO,CAAC;IAC3B,OAAO,CAAC,QAAQ,CAAC,SAAS,CAAgC;gBAE9C,QAAQ,UAAQ;IAI5B,SAAS,CAAC,CAAC,EAAE,CAAC,GAAG,IAAI;IAKrB,OAAO,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,MAAM,SAAI,GAAG,IAAI;IAOrC,QAAQ,IAAI,CAAC,EAAE;IAIf,SAAS,CAAC,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE;IAIpB,KAAK,IAAI,KAAK,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;IAgBvB,kDAAkD;IAClD,GAAG,CAAC,KAAK,EAAE,CAAC,GAAG,CAAC,EAAE;IAgBlB,gDAAgD;IAChD,GAAG,CAAC,KAAK,EAAE,CAAC,GAAG,CAAC,EAAE;IAYlB;;;;;;;OAOG;IACH,OAAO,CAAC,WAAW;IA2BnB,4EAA4E;IAC5E,QAAQ,CAAC,KAAK,EAAE,CAAC,GAAG,GAAG,CAAC,CAAC,EAAE,MAAM,CAAC;IAIlC,uEAAuE;IACvE,YAAY,CAAC,KAAK,EAAE,CAAC,EAAE,GAAG,EAAE,CAAC,GAAG;QAAE,QAAQ,EAAE,MAAM,CAAC;QAAC,IAAI,EAAE,CAAC,EAAE,CAAA;KAAE;IAY/D,gDAAgD;IAChD,mBAAmB,IAAI,CAAC,EAAE,EAAE;IAY5B,uEAAuE;IACvE,eAAe,IAAI,CAAC,EAAE,GAAG,IAAI;IAiB7B,4CAA4C;IAC5C,QAAQ,IAAI,OAAO;IAgBnB,kFAAkF;IAClF,mBAAmB,IAAI,KAAK,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;IAsBrC,sFAAsF;IACtF,iBAAiB,IAAI;QAAE,MAAM,EAAE,MAAM,EAAE,EAAE,CAAC;QAAC,KAAK,EAAE,CAAC,EAAE,CAAA;KAAE;IAavD;;;;;;;;;OASG;IACH,2BAA2B,IAAI,CAAC,EAAE,EAAE;IAwCpC,sDAAsD;IACtD,aAAa,IAAI;QAAE,SAAS,EAAE,MAAM,EAAE,EAAE,CAAC;QAAC,KAAK,EAAE,CAAC,EAAE,CAAA;KAAE;CASvD"}
+{"version":3,"file":"Graph.d.ts","sourceRoot":"","sources":["../src/Graph.ts"],"names":[],"mappings":"AAAA;;;;;;GAMG;AAEH,MAAM,WAAW,IAAI,CAAC,CAAC;IACrB,IAAI,EAAE,CAAC,CAAC;IACR,EAAE,EAAE,CAAC,CAAC;IACN,MAAM,EAAE,MAAM,CAAC;CAChB;AAoDD,qBAAa,KAAK,CAAC,CAAC,GAAG,MAAM,GAAG,MAAM;IACpC,QAAQ,CAAC,QAAQ,EAAE,OAAO,CAAC;IAC3B,OAAO,CAAC,QAAQ,CAAC,SAAS,CAAgC;gBAE9C,QAAQ,UAAQ;IAI5B,SAAS,CAAC,CAAC,EAAE,CAAC,GAAG,IAAI;IAKrB,OAAO,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,MAAM,SAAI,GAAG,IAAI;IAOrC,QAAQ,IAAI,CAAC,EAAE;IAIf,SAAS,CAAC,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE;IAIpB,KAAK,IAAI,KAAK,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;IAgBvB,kDAAkD;IAClD,GAAG,CAAC,KAAK,EAAE,CAAC,GAAG,CAAC,EAAE;IAgBlB,gDAAgD;IAChD,GAAG,CAAC,KAAK,EAAE,CAAC,GAAG,CAAC,EAAE;IAYlB;;;;;;;OAOG;IACH,OAAO,CAAC,WAAW;IA2BnB,4EAA4E;IAC5E,QAAQ,CAAC,KAAK,EAAE,CAAC,GAAG,GAAG,CAAC,CAAC,EAAE,MAAM,CAAC;IAIlC,uEAAuE;IACvE,YAAY,CAAC,KAAK,EAAE,CAAC,EAAE,GAAG,EAAE,CAAC,GAAG;QAAE,QAAQ,EAAE,MAAM,CAAC;QAAC,IAAI,EAAE,CAAC,EAAE,CAAA;KAAE;IAY/D,gDAAgD;IAChD,mBAAmB,IAAI,CAAC,EAAE,EAAE;IAY5B,uEAAuE;IACvE,eAAe,IAAI,CAAC,EAAE,GAAG,IAAI;IAiB7B,4CAA4C;IAC5C,QAAQ,IAAI,OAAO;IAgBnB,kFAAkF;IAClF,mBAAmB,IAAI,KAAK,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;IAsBrC;;;;;;OAMG;IACH,GAAG,CAAC,SAAS,GAAE,SAAS,GAAG,MAAkB,GAAG;QAAE,KAAK,EAAE,KAAK,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC;QAAC,WAAW,EAAE,MAAM,CAAC;QAAC,UAAU,EAAE,MAAM,CAAA;KAAE;IAgClH,sGAAsG;IACtG,GAAG,IAAI,CAAC,EAAE,EAAE;IAIZ;;;;;;OAMG;IACH,YAAY,IAAI;QAAE,KAAK,EAAE,KAAK,CAAC,MAAM,CAAC,CAAC;QAAC,UAAU,EAAE,CAAC,EAAE,EAAE,CAAC;QAAC,WAAW,EAAE,GAAG,CAAC,CAAC,EAAE,MAAM,CAAC,CAAA;KAAE;IAiBxF,sFAAsF;IACtF,iBAAiB,IAAI;QAAE,MAAM,EAAE,MAAM,EAAE,EAAE,CAAC;QAAC,KAAK,EAAE,CAAC,EAAE,CAAA;KAAE;IAavD;;;;;;;;;OASG;IACH,2BAA2B,IAAI,CAAC,EAAE,EAAE;IAwCpC,sDAAsD;IACtD,aAAa,IAAI;QAAE,SAAS,EAAE,MAAM,EAAE,EAAE,CAAC;QAAC,KAAK,EAAE,CAAC,EAAE,CAAA;KAAE;CASvD"}