universal-physics-tensor 0.4.6 → 0.5.1

package/dist/numerical/curvature-lowering-helpers.js

@@ -0,0 +1,412 @@
+/**
+ * Numerical helpers for Riemann-curvature lowering (Task 6 [U] / v0.5.0 1c-ii).
+ *
+ * Mirrors `connection-lowering-helpers.ts`: produces Γ and ∂Γ tensors from
+ * coordinate-dependent metric closures, then assembles the Riemann tensor via
+ * the corrected (Carroll-Ch.3 §3.4 / Adam+Eve F4-S3) formula:
+ *
+ *   R^ρ_{σμν} = ∂_μ Γ^ρ_{σν} − ∂_ν Γ^ρ_{σμ} + Σ_λ (Γ^ρ_{λμ} Γ^λ_{σν} − Γ^ρ_{λν} Γ^λ_{σμ})
+ *
+ * Index conventions (pinned via JSDoc + a runtime assert):
+ *   - Γ[ρ][σ][ν] = Γ^ρ_{σν}  (matches computeChristoffelTensor's [α][μ][ν] with
+ *     σ ↔ μ — i.e., the SECOND axis is σ, the THIRD is ν, per Carroll Ch.5).
+ *   - dGamma[λ][ρ][σ][ν] = ∂_λ Γ^ρ_{σν}     (I3)
+ *   - R[ρ][σ][μ][ν] = R^ρ_{σμν}              (Carroll)
+ *
+ * NOTE on the σ/μ swap vs `computeChristoffelTensor`:
+ *   `computeChristoffelTensor` was authored for v0.4.0 covariant-derivative
+ *   lowering, which sums Γ^α_{μλ} T^λ over the middle index `μ` (the wrt index
+ *   of the cov-deriv). v0.5.0 Riemann needs σ in the SECOND-lower slot of every
+ *   Γ. Because Γ is symmetric in its two lower indices (Christoffel symmetry),
+ *   the two storage conventions are numerically identical — we use the
+ *   existing [α][middle][last] storage and access it with σ in the middle slot.
+ *
+ * The fast path uses the v0.4.0 `computeChristoffelTensor` infrastructure
+ * (no new tensor-engine code) and the v0.4.0 `pderivNumericalFn` pattern
+ * (centered finite differences). M11: ∂Γ is computed by repeated centered
+ * finite-differences on `christoffel()`, not a new AST node type or AD pass.
+ *
+ * @module numerical/curvature-lowering-helpers
+ */
+import { NumericalBackendError } from './errors.js';
+import { computeChristoffelTensor, flattenNA } from './connection-lowering-helpers.js';
+import { pderivNumericalFn } from './pderiv.js';
+// ---------------------------------------------------------------------------
+// Finite-difference step
+// ---------------------------------------------------------------------------
+/**
+ * Step size for finite differences on the Christoffel function.
+ *
+ * Inner step (∂g for one Γ evaluation) uses the v0.4.0 default
+ * `1e-6·max(|x|,1)` via `pderivNumericalFn`-style centered FD. This gives ∂g
+ * with relative precision ≈ 1e-11 on Schwarzschild components (verified by
+ * inspection at r=3·r_s).
+ *
+ * Outer step (∂Γ for one R evaluation) is `1e-4·max(|x|,1)` and we use a
+ * **4th-order centered stencil**:
+ *
+ *   f'(x) ≈ [−f(x+2h) + 8 f(x+h) − 8 f(x-h) + f(x-2h)] / (12 h)
+ *
+ * Truncation O(h⁴) instead of O(h²). At h=0.886 m (r=3·r_s, M=M_sun),
+ * truncation ≈ h⁴·Γ⁽⁵⁾/Γ ~ 1e-13 relative, round-off ε·|Γ|/h ~ 1e-15 relative.
+ * Both sit well below the 1e-9 gate.
+ *
+ * The 2nd-order outer FD that v0.4.0 covariant-derivative lowering uses
+ * leaves ~3e-6 relative error on R^t_{rtr} due to noise propagation through
+ * the c²-scaled g_{tt} component — 4th-order outer FD recovers the precision.
+ */
+export function outerStep(x) {
+    return 1e-4 * Math.max(Math.abs(x), 1);
+}
+/**
+ * Build a flat ∂_μ g sampler around point x using a **4th-order centered
+ * stencil**:
+ *
+ *   ∂_μ g ≈ [−g(x+2h) + 8 g(x+h) − 8 g(x−h) + g(x−2h)] / (12 h)
+ *
+ * Returned function matches the `getMetricDeriv` signature expected by
+ * `computeChristoffelTensor`.
+ *
+ * Inner step `h_i = 1e-3·max(|x|,1)` is intentionally LARGER than the v0.4.0
+ * 1e-6 default to compensate for the c²·g_{tt} ~ 6e16 cancellation noise: a
+ * 4th-order stencil at h_i ≈ 8.86 m around r=3·r_s gives truncation
+ * O(h_i⁴·g⁽⁵⁾) ~ 1e-13 relative on each ∂g component while keeping round-off
+ * ε|g|/h_i ~ 1e-15 relative — orders below the 1e-9 component-match gate.
+ */
+function makeInnerGradFn(gFn, x0, N) {
+    // v0.5.1 PD-7: delegates to the shared `pderivNumericalFn(..., {order: 4})`
+    // 4th-order centered stencil. The explicit `h = 1e-3·max(|x|,1)` override
+    // is load-bearing: the c²·g_{tt} cancellation noise on Schwarzschild
+    // demands a larger step than the v0.5.0-default 1e-4 to balance truncation
+    // vs round-off at the ≤1e-9 component-match gate. (See module JSDoc.)
+    return (mu) => {
+        const xc = x0[mu];
+        const h = 1e-3 * Math.max(Math.abs(xc), 1);
+        const d = pderivNumericalFn(gFn, x0, mu, { order: 4, h });
+        const flat = Array.isArray(d) ? d : [d];
+        if (flat.length !== N * N) {
+            throw new NumericalBackendError(`curvature-lowering: metric closure returned wrong shape — expected ${N * N}, ` +
+                `got ${flat.length}`);
+        }
+        return flat;
+    };
+}
+/**
+ * Compute Γ^ρ_{σν} at coordinate x. Internally builds a centered-FD ∂g sampler
+ * around x and delegates to `computeChristoffelTensor` (the v0.4.0 helper).
+ *
+ * Returns a fully-materialised 3-deep nested array Γ[ρ][σ][ν]. (We do NOT keep
+ * the EngineTensor wrapper because we immediately rebuild dGamma and R as
+ * plain nested arrays — round-tripping through the engine on every FD sample
+ * would dominate the cost.)
+ */
+export function christoffelAt(x, gFn, gInverseFn, N, engine) {
+    const gInvFlat = flattenNA(gInverseFn(x));
+    if (gInvFlat.length !== N * N) {
+        throw new NumericalBackendError(`curvature-lowering: gInverseFn returned wrong shape — expected ${N * N}, got ${gInvFlat.length}`);
+    }
+    const gradFn = makeInnerGradFn(gFn, x, N);
+    const tensor = computeChristoffelTensor(gInvFlat, gradFn, N, engine);
+    const nested = engine.toNested(tensor);
+    engine.dispose?.(tensor);
+    // Cast to the expected nested shape; computeChristoffelTensor produces
+    // shape [N,N,N] so this is structurally [ρ][middle][last]. We treat
+    // [middle] as σ and [last] as ν (Christoffel symmetry makes the two
+    // lower-slot conventions numerically interchangeable).
+    return nested;
+}
+/**
+ * dGamma[λ][ρ][σ][ν] = ∂_λ Γ^ρ_{σν} via centered FD on `christoffelAt`. (I3)
+ *
+ * Includes a runtime index-assert verifying `dGamma[1][1][1][1]` is finite
+ * before returning (catches silent shape-permutation bugs introduced by future
+ * edits to `computeChristoffelTensor`'s output ordering).
+ */
+export function dGammaAt(x, gFn, gInverseFn, N, engine) {
+    // Allocate dGamma[λ][ρ][σ][ν]
+    const dGamma = Array.from({ length: N }, () => Array.from({ length: N }, () => Array.from({ length: N }, () => new Array(N).fill(0))));
+    // v0.5.1 PD-7: 4th-order centered stencil via the shared
+    // `pderivNumericalFn(..., {order: 4})` (h = 1e-4·max(|x|,1) default matches
+    // the v0.5.0 `outerStep`). Each axis returns a flat number[] of length N³;
+    // unflatten into the `dGamma[λ][ρ][σ][ν]` axis-0 slab.
+    const christFn = (xs) => christoffelAt(xs, gFn, gInverseFn, N, engine);
+    for (let lam = 0; lam < N; lam++) {
+        const d = pderivNumericalFn(christFn, x, lam, { order: 4 });
+        const flat = Array.isArray(d) ? d : [d];
+        if (flat.length !== N * N * N) {
+            throw new NumericalBackendError(`dGammaAt: christoffelAt derivative shape mismatch — expected ${N * N * N}, got ${flat.length}`);
+        }
+        const N2 = N * N;
+        for (let rho = 0; rho < N; rho++) {
+            for (let sigma = 0; sigma < N; sigma++) {
+                for (let nu = 0; nu < N; nu++) {
+                    dGamma[lam][rho][sigma][nu] = flat[rho * N2 + sigma * N + nu];
+                }
+            }
+        }
+    }
+    // I3 runtime index-assert: verify the container is populated as expected.
+    // dGamma[1][1][1][1] is ∂_r Γ^r_{rr} for Schwarzschild — a non-trivial real
+    // number, so a `Number.isFinite` check catches both shape regressions and
+    // numerical blow-ups.
+    const probe = dGamma[1]?.[1]?.[1]?.[1];
+    if (typeof probe !== 'number' || !Number.isFinite(probe)) {
+        throw new NumericalBackendError(`curvature-lowering: dGamma index-assert failed — dGamma[1][1][1][1] is ` +
+            `${probe} (expected finite number). Shape regression in christoffelAt?`);
+    }
+    return dGamma;
+}
+/**
+ * Build the full 4×4×4×4 Riemann tensor R^ρ_{σμν} from Γ and ∂Γ per the
+ * Carroll formula (Adam+Eve F4-S3):
+ *
+ *   R[ρ][σ][μ][ν] = ∂_μ Γ[ρ][σ][ν] − ∂_ν Γ[ρ][σ][μ]
+ *                 + Σ_λ (Γ[ρ][λ][μ] Γ[λ][σ][ν] − Γ[ρ][λ][ν] Γ[λ][σ][μ])
+ *
+ * σ in the SECOND lower slot of each Γ (F4/S3 correction).
+ */
+export function buildRiemann(gamma, dGamma, N) {
+    const R = Array.from({ length: N }, () => Array.from({ length: N }, () => Array.from({ length: N }, () => new Array(N).fill(0))));
+    for (let rho = 0; rho < N; rho++) {
+        for (let sigma = 0; sigma < N; sigma++) {
+            for (let mu = 0; mu < N; mu++) {
+                for (let nu = 0; nu < N; nu++) {
+                    // ∂_μ Γ^ρ_{σν} − ∂_ν Γ^ρ_{σμ}
+                    let value = dGamma[mu][rho][sigma][nu] - dGamma[nu][rho][sigma][mu];
+                    // Σ_λ (Γ^ρ_{λμ} Γ^λ_{σν} − Γ^ρ_{λν} Γ^λ_{σμ})
+                    for (let lam = 0; lam < N; lam++) {
+                        value += gamma[rho][lam][mu] * gamma[lam][sigma][nu]
+                            - gamma[rho][lam][nu] * gamma[lam][sigma][mu];
+                    }
+                    R[rho][sigma][mu][nu] = value;
+                }
+            }
+        }
+    }
+    return R;
+}
+// ---------------------------------------------------------------------------
+// v0.5.0 Task 9: helpers for bianchiResidual()
+// ---------------------------------------------------------------------------
+/**
+ * Compute the upper-Riemann tensor R^ρ_{σμν} at coordinate x using the same
+ * Γ + ∂Γ pipeline as the riemann-tensor lowering case. Encapsulates the
+ * full christoffelAt + dGammaAt + buildRiemann sequence so callers (Task 9
+ * Bianchi residual) can sample R at perturbed coordinates without
+ * re-implementing the FD machinery.
+ */
+export function riemannUpperAt(x, gFn, gInverseFn, N, engine) {
+    const gamma = christoffelAt(x, gFn, gInverseFn, N, engine);
+    const dGamma = dGammaAt(x, gFn, gInverseFn, N, engine);
+    return buildRiemann(gamma, dGamma, N);
+}
+/**
+ * Lower the upper-ρ index of R^ρ_{σμν} via the covariant metric g_{aρ}:
+ *
+ *   R_{aσμν} = Σ_ρ g_{aρ} R^ρ_{σμν}
+ *
+ * Output index order: `[a][σ][μ][ν]` — all four lower. (a is the freshly
+ * lowered index in the first slot.)
+ */
+export function lowerFirstIndex(R, gLowerFlat, N) {
+    const Rlow = Array.from({ length: N }, () => Array.from({ length: N }, () => Array.from({ length: N }, () => new Array(N).fill(0))));
+    for (let a = 0; a < N; a++) {
+        for (let sig = 0; sig < N; sig++) {
+            for (let mu = 0; mu < N; mu++) {
+                for (let nu = 0; nu < N; nu++) {
+                    let sum = 0;
+                    for (let rho = 0; rho < N; rho++) {
+                        sum += gLowerFlat[a * N + rho] * R[rho][sig][mu][nu];
+                    }
+                    Rlow[a][sig][mu][nu] = sum;
+                }
+            }
+        }
+    }
+    return Rlow;
+}
+/**
+ * Sample the all-lower Riemann tensor R_{αβγδ}(x) — combines `riemannUpperAt`
+ * with `lowerFirstIndex` for the single-coordinate evaluation. Used both as
+ * the base sample and (with perturbed x) as the FD inputs for ∂_λ R_{αβγδ}.
+ */
+export function riemannLowerAt(x, gFn, gInverseFn, N, engine) {
+    const Rup = riemannUpperAt(x, gFn, gInverseFn, N, engine);
+    const gFlat = flattenNA(gFn(x));
+    if (gFlat.length !== N * N) {
+        throw new NumericalBackendError(`bianchi-residual: gFn returned wrong shape — expected ${N * N}, got ${gFlat.length}`);
+    }
+    return lowerFirstIndex(Rup, gFlat, N);
+}
+/**
+ * Compute ∂_λ R_{αβγδ}(x) via a 4th-order centered stencil on
+ * `riemannLowerAt`. Index order: `dR[λ][α][β][γ][δ] = ∂_λ R_{αβγδ}`.
+ *
+ * Uses the same outer step (`outerStep`) as Task 6's dGammaAt so the
+ * FD-noise compounding pattern matches. This is one extra layer of FD on top
+ * of the Christoffel-of-Christoffel double-FD already inside `riemannLowerAt`
+ * — total: ∂(∂(∂g)) at 4th order in each layer. Empirical noise floor reached
+ * by the per-component value is documented in the test report.
+ */
+export function dRiemannLowerAt(x, gFn, gInverseFn, N, engine) {
+    const dR = Array.from({ length: N }, () => Array.from({ length: N }, () => Array.from({ length: N }, () => Array.from({ length: N }, () => new Array(N).fill(0)))));
+    // v0.5.1 PD-7: 4th-order centered stencil via the shared
+    // `pderivNumericalFn(..., {order: 4})` — default h = 1e-4·max(|x|,1) matches
+    // the v0.5.0 `outerStep`. Each axis call returns a flat number[] of length
+    // N⁴ (riemannLowerAt produces [N,N,N,N]); unflatten into the
+    // `dR[λ][a][b][c][d]` axis-0 slab.
+    const RlowerFn = (xs) => riemannLowerAt(xs, gFn, gInverseFn, N, engine);
+    const N2 = N * N;
+    const N3 = N * N * N;
+    for (let lam = 0; lam < N; lam++) {
+        const flat = pderivNumericalFn(RlowerFn, x, lam, { order: 4 });
+        if (!Array.isArray(flat) || flat.length !== N2 * N2) {
+            throw new NumericalBackendError(`dRiemannLowerAt: derivative shape mismatch — expected ${N2 * N2}, got ${Array.isArray(flat) ? flat.length : 'scalar'}`);
+        }
+        for (let a = 0; a < N; a++) {
+            for (let b = 0; b < N; b++) {
+                for (let c = 0; c < N; c++) {
+                    for (let d = 0; d < N; d++) {
+                        dR[lam][a][b][c][d] = flat[a * N3 + b * N2 + c * N + d];
+                    }
+                }
+            }
+        }
+    }
+    return dR;
+}
+/**
+ * Compute ∇_λ R_{μνρσ} via partial + Christoffel-correction terms (Approach 1).
+ *
+ *   ∇_λ R_{μνρσ} = ∂_λ R_{μνρσ}
+ *                  − Γ^α_{λμ} R_{ανρσ}
+ *                  − Γ^α_{λν} R_{μαρσ}
+ *                  − Γ^α_{λρ} R_{μνασ}
+ *                  − Γ^α_{λσ} R_{μνρα}
+ *
+ * Index order: `covR[λ][μ][ν][ρ][σ] = ∇_λ R_{μνρσ}`.
+ *
+ * Approach choice (full ∇, not raw ∂): the second Bianchi identity
+ * ∇_{[λ} R_{μν]ρσ} = 0 is the canonical statement. Using raw ∂ would give a
+ * residual dominated by the (cyclic-non-cancelling) Christoffel-correction
+ * terms, masking the identity check with O(1) algebraic clutter. Full ∇
+ * makes the test a genuine self-consistency check of the lowered Riemann.
+ */
+export function covariantDerivRiemannLowerAt(x, gFn, gInverseFn, N, engine) {
+    const gamma = christoffelAt(x, gFn, gInverseFn, N, engine);
+    const R = riemannLowerAt(x, gFn, gInverseFn, N, engine);
+    const dR = dRiemannLowerAt(x, gFn, gInverseFn, N, engine);
+    const covR = Array.from({ length: N }, () => Array.from({ length: N }, () => Array.from({ length: N }, () => Array.from({ length: N }, () => new Array(N).fill(0)))));
+    for (let lam = 0; lam < N; lam++) {
+        for (let mu = 0; mu < N; mu++) {
+            for (let nu = 0; nu < N; nu++) {
+                for (let rho = 0; rho < N; rho++) {
+                    for (let sig = 0; sig < N; sig++) {
+                        let value = dR[lam][mu][nu][rho][sig];
+                        // Subtract Γ^α_{λμ} R_{ανρσ} (correction for the μ index)
+                        // Subtract Γ^α_{λν} R_{μαρσ} (correction for the ν index)
+                        // Subtract Γ^α_{λρ} R_{μνασ} (correction for the ρ index)
+                        // Subtract Γ^α_{λσ} R_{μνρα} (correction for the σ index)
+                        for (let a = 0; a < N; a++) {
+                            value -= gamma[a][lam][mu] * R[a][nu][rho][sig];
+                            value -= gamma[a][lam][nu] * R[mu][a][rho][sig];
+                            value -= gamma[a][lam][rho] * R[mu][nu][a][sig];
+                            value -= gamma[a][lam][sig] * R[mu][nu][rho][a];
+                        }
+                        covR[lam][mu][nu][rho][sig] = value;
+                    }
+                }
+            }
+        }
+    }
+    return covR;
+}
+/**
+ * Build the second-Bianchi-identity residual (cyclic over first three indices):
+ *
+ *   B_{λμνρσ} = ∇_λ R_{μνρσ} + ∇_μ R_{νλρσ} + ∇_ν R_{λμρσ}
+ *
+ * Carroll Eq. 3.95: B ≡ 0 in any (torsion-free) Lorentzian manifold. The
+ * residual measures the FD-truncation + cancellation noise on the lowered
+ * Riemann tensor through one extra ∂ layer.
+ *
+ * Index order: `B[λ][μ][ν][ρ][σ]` — all five lower.
+ */
+export function bianchiResidualAt(x, gFn, gInverseFn, N, engine) {
+    const covR = covariantDerivRiemannLowerAt(x, gFn, gInverseFn, N, engine);
+    const B = Array.from({ length: N }, () => Array.from({ length: N }, () => Array.from({ length: N }, () => Array.from({ length: N }, () => new Array(N).fill(0)))));
+    for (let lam = 0; lam < N; lam++) {
+        for (let mu = 0; mu < N; mu++) {
+            for (let nu = 0; nu < N; nu++) {
+                for (let rho = 0; rho < N; rho++) {
+                    for (let sig = 0; sig < N; sig++) {
+                        B[lam][mu][nu][rho][sig] =
+                            covR[lam][mu][nu][rho][sig] // ∇_λ R_{μνρσ}
+                                + covR[mu][nu][lam][rho][sig] // ∇_μ R_{νλρσ}
+                                + covR[nu][lam][mu][rho][sig]; // ∇_ν R_{λμρσ}
+                    }
+                }
+            }
+        }
+    }
+    return B;
+}
+/**
+ * JS-side contraction of a flat 4-axis tensor R[a][b][c][d] (length N⁴, row-major)
+ * down to a 2-tensor by summing one upper axis against one lower axis.
+ *
+ * Configuration:
+ *   - `upperAxis`: the axis index (0..3) that contracts WITH `lowerAxis` (dummy slot).
+ *   - `lowerAxis`: the partner axis (0..3) forming the dummy pair.
+ *   - `outAxes`:   the two surviving axes [i, j] that become the (free output)
+ *                  rows/cols of the returned N×N matrix.
+ *
+ * For Carroll Eq. 3.91 Ricci R_μν = R^λ_{μλν} on Riemann stored as R[ρ][σ][μ][ν]:
+ *   contractRiemannJS(flatR, N, { upperAxis: 0, lowerAxis: 2, outAxes: [1, 3] })
+ * sums `R[λ][μ_out][λ][ν_out]` over λ.
+ *
+ * Pre-flattened input convention: row-major `R[a][b][c][d] = flatR[a*N³ + b*N² + c*N + d]`.
+ * (Matches `flattenNestedArray` output for an [N,N,N,N] EngineTensor.)
+ *
+ * @internal — extracted from `lowering.ts` ricci-tensor arm (AS-1, v0.5.1). Kept
+ * deliberately narrow (4-axis source → 2-axis output) to match the v0.5.0 GR
+ * curvature pipeline; generalisation deferred to v0.6.0 if Weyl/Kretschmann
+ * arrives (see `docs/architecture/COMPONENTS.md` curvature-layer pattern note).
+ */
+export function contractRiemannJS(flatR, N, config) {
+    const { upperAxis, lowerAxis, outAxes } = config;
+    if (upperAxis === lowerAxis) {
+        throw new NumericalBackendError(`contractRiemannJS: upperAxis (${upperAxis}) and lowerAxis (${lowerAxis}) must differ`);
+    }
+    const [outA, outB] = outAxes;
+    if (outA === outB) {
+        throw new NumericalBackendError(`contractRiemannJS: outAxes [${outA},${outB}] must be distinct`);
+    }
+    if (outA === upperAxis || outA === lowerAxis || outB === upperAxis || outB === lowerAxis) {
+        throw new NumericalBackendError(`contractRiemannJS: outAxes [${outA},${outB}] must not overlap dummy pair [${upperAxis},${lowerAxis}]`);
+    }
+    if (flatR.length !== N * N * N * N) {
+        throw new NumericalBackendError(`contractRiemannJS: flatR length ${flatR.length} != N⁴ = ${N * N * N * N}`);
+    }
+    // Row-major strides for R[a][b][c][d].
+    const N2 = N * N;
+    const N3 = N * N * N;
+    const strides = [N3, N2, N, 1];
+    const strideDummy = strides[upperAxis] + strides[lowerAxis];
+    const strideOutA = strides[outA];
+    const strideOutB = strides[outB];
+    const out = Array.from({ length: N }, () => new Array(N).fill(0));
+    for (let i = 0; i < N; i++) {
+        const baseI = i * strideOutA;
+        for (let j = 0; j < N; j++) {
+            const baseIJ = baseI + j * strideOutB;
+            let sum = 0;
+            for (let lam = 0; lam < N; lam++) {
+                sum += flatR[baseIJ + lam * strideDummy];
+            }
+            out[i][j] = sum;
+        }
+    }
+    return out;
+}
+//# sourceMappingURL=curvature-lowering-helpers.js.map

package/dist/numerical/curvature-lowering-helpers.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"curvature-lowering-helpers.js","sourceRoot":"","sources":["../../src/numerical/curvature-lowering-helpers.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GA6BG;AAIH,OAAO,EAAE,qBAAqB,EAAE,MAAM,aAAa,CAAC;AACpD,OAAO,EAAE,wBAAwB,EAAE,SAAS,EAAE,MAAM,kCAAkC,CAAC;AACvF,OAAO,EAAE,iBAAiB,EAAE,MAAM,aAAa,CAAC;AAqBhD,8EAA8E;AAC9E,yBAAyB;AACzB,8EAA8E;AAE9E;;;;;;;;;;;;;;;;;;;;GAoBG;AACH,MAAM,UAAU,SAAS,CAAC,CAAS;IACjC,OAAO,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;AACzC,CAAC;AAED;;;;;;;;;;;;;;GAcG;AACH,SAAS,eAAe,CACtB,GAAa,EACb,EAAyB,EACzB,CAAS;IAET,4EAA4E;IAC5E,0EAA0E;IAC1E,qEAAqE;IACrE,2EAA2E;IAC3E,sEAAsE;IACtE,OAAO,CAAC,EAAU,EAAY,EAAE;QAC9B,MAAM,EAAE,GAAG,EAAE,CAAC,EAAE,CAAC,CAAC;QAClB,MAAM,CAAC,GAAG,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;QAC3C,MAAM,CAAC,GAAG,iBAAiB,CAAC,GAAG,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,KAAK,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;QAC1D,MAAM,IAAI,GAAG,KAAK,CAAC,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC,CAAE,CAAc,CAAC,CAAC,CAAC,CAAC,CAAW,CAAC,CAAC;QAChE,IAAI,IAAI,CAAC,MAAM,KAAK,CAAC,GAAG,CAAC,EAAE,CAAC;YAC1B,MAAM,IAAI,qBAAqB,CAC7B,sEAAsE,CAAC,GAAG,CAAC,IAAI;gBAC/E,OAAO,IAAI,CAAC,MAAM,EAAE,CACrB,CAAC;QACJ,CAAC;QACD,OAAO,IAAI,CAAC;IACd,CAAC,CAAC;AACJ,CAAC;AAED;;;;;;;;GAQG;AACH,MAAM,UAAU,aAAa,CAC3B,CAAwB,EACxB,GAAa,EACb,UAAoB,EACpB,CAAS,EACT,MAAoB;IAEpB,MAAM,QAAQ,GAAG,SAAS,CAAC,UAAU,CAAC,CAAC,CAAC,CAAC,CAAC;IAC1C,IAAI,QAAQ,CAAC,MAAM,KAAK,CAAC,GAAG,CAAC,EAAE,CAAC;QAC9B,MAAM,IAAI,qBAAqB,CAC7B,kEAAkE,CAAC,GAAG,CAAC,SAAS,QAAQ,CAAC,MAAM,EAAE,CAClG,CAAC;IACJ,CAAC;IACD,MAAM,MAAM,GAAG,eAAe,CAAC,GAAG,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;IAC1C,MAAM,MAAM,GAAG,wBAAwB,CAAC,QAAQ,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,CAAC,CAAC;IACrE,MAAM,MAAM,GAAG,MAAM,CAAC,QAAQ,CAAC,MAAM,CAAgB,CAAC;IACtD,MAAM,CAAC,OAAO,EAAE,CAAC,MAAM,CAAC,CAAC;IAEzB,uEAAuE;IACvE,oEAAoE;IACpE,oEAAoE;IACpE,uDAAuD;IACvD,OAAO,MAAsB,CAAC;AAChC,CAAC;AAED;;;;;;GAMG;AACH,MAAM,UAAU,QAAQ,CACtB,CAAwB,EACxB,GAAa,EACb,UAAoB,EACpB,CAAS,EACT,MAAoB;IAEpB,8BAA8B;IAC9B,MAAM,MAAM,GAAmB,KAAK,CAAC,IAAI,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,EAAE,GAAG,EAAE,CAC5D,KAAK,CAAC,IAAI,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,EAAE,GAAG,EAAE,CAC7B,KAAK,CAAC,IAAI,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,EAAE,GAAG,EAAE,CAAC,IAAI,KAAK,CAAS,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAC9D,CACF,CAAC;IAEF,yDAAyD;IACzD,4EAA4E;IAC5E,2EAA2E;IAC3E,uDAAuD;IACvD,MAAM,QAAQ,GAAG,CAAC,EAAyB,EAAe,EAAE,CAC1D,aAAa,CAAC,EAAE,EAAE,GAAG,EAAE,UAAU,EAAE,CAAC,EAAE,MAAM,CAAC,CAAC;IAChD,KAAK,IAAI,GAAG,GAAG,CAAC,EAAE,GAAG,GAAG,CAAC,EAAE,GAAG,EAAE,EAAE,CAAC;QACjC,MAAM,CAAC,GAAG,iBAAiB,CAAC,QAAQ,EAAE,CAAC,EAAE,GAAG,EAAE,EAAE,KAAK,EAAE,CAAC,EAAE,CAAC,CAAC;QAC5D,MAAM,IAAI,GAAG,KAAK,CAAC,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC,CAAE,CAAc,CAAC,CAAC,CAAC,CAAC,CAAW,CAAC,CAAC;QAChE,IAAI,IAAI,CAAC,MAAM,KAAK,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,CAAC;YAC9B,MAAM,IAAI,qBAAqB,CAC7B,gEAAgE,CAAC,GAAG,CAAC,GAAG,CAAC,SAAS,IAAI,CAAC,MAAM,EAAE,CAChG,CAAC;QACJ,CAAC;QACD,MAAM,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC;QACjB,KAAK,IAAI,GAAG,GAAG,CAAC,EAAE,GAAG,GAAG,CAAC,EAAE,GAAG,EAAE,EAAE,CAAC;YACjC,KAAK,IAAI,KAAK,GAAG,CAAC,EAAE,KAAK,GAAG,CAAC,EAAE,KAAK,EAAE,EAAE,CAAC;gBACvC,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,CAAC,EAAE,EAAE,EAAE,EAAE,CAAC;oBAC9B,MAAM,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,KAAK,CAAC,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,GAAG,GAAG,EAAE,GAAG,KAAK,GAAG,CAAC,GAAG,EAAE,CAAC,CAAC;gBAChE,CAAC;YACH,CAAC;QACH,CAAC;IACH,CAAC;IAED,0EAA0E;IAC1E,4EAA4E;IAC5E,0EAA0E;IAC1E,sBAAsB;IACtB,MAAM,KAAK,GAAG,MAAM,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;IACvC,IAAI,OAAO,KAAK,KAAK,QAAQ,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,KAAK,CAAC,EAAE,CAAC;QACzD,MAAM,IAAI,qBAAqB,CAC7B,yEAAyE;YACzE,GAAG,KAAK,+DAA+D,CACxE,CAAC;IACJ,CAAC;IAED,OAAO,MAAsB,CAAC;AAChC,CAAC;AAED;;;;;;;;GAQG;AACH,MAAM,UAAU,YAAY,CAC1B,KAAkB,EAClB,MAAoB,EACpB,CAAS;IAET,MAAM,CAAC,GAAmB,KAAK,CAAC,IAAI,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,EAAE,GAAG,EAAE,CACvD,KAAK,CAAC,IAAI,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,EAAE,GAAG,EAAE,CAC7B,KAAK,CAAC,IAAI,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,EAAE,GAAG,EAAE,CAAC,IAAI,KAAK,CAAS,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAC9D,CACF,CAAC;IAEF,KAAK,IAAI,GAAG,GAAG,CAAC,EAAE,GAAG,GAAG,CAAC,EAAE,GAAG,EAAE,EAAE,CAAC;QACjC,KAAK,IAAI,KAAK,GAAG,CAAC,EAAE,KAAK,GAAG,CAAC,EAAE,KAAK,EAAE,EAAE,CAAC;YACvC,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,CAAC,EAAE,EAAE,EAAE,EAAE,CAAC;gBAC9B,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,CAAC,EAAE,EAAE,EAAE,EAAE,CAAC;oBAC9B,8BAA8B;oBAC9B,IAAI,KAAK,GAAG,MAAM,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,KAAK,CAAC,CAAC,EAAE,CAAC,GAAG,MAAM,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC;oBACpE,8CAA8C;oBAC9C,KAAK,IAAI,GAAG,GAAG,CAAC,EAAE,GAAG,GAAG,CAAC,EAAE,GAAG,EAAE,EAAE,CAAC;wBACjC,KAAK,IAAI,KAAK,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,GAAG,KAAK,CAAC,GAAG,CAAC,CAAC,KAAK,CAAC,CAAC,EAAE,CAAC;8BAC3C,KAAK,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,GAAG,KAAK,CAAC,GAAG,CAAC,CAAC,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC;oBACvD,CAAC;oBACD,CAAC,CAAC,GAAG,CAAC,CAAC,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,GAAG,KAAK,CAAC;gBAChC,CAAC;YACH,CAAC;QACH,CAAC;IACH,CAAC;IAED,OAAO,CAAC,CAAC;AACX,CAAC;AAED,8EAA8E;AAC9E,+CAA+C;AAC/C,8EAA8E;AAE9E;;;;;;GAMG;AACH,MAAM,UAAU,cAAc,CAC5B,CAAwB,EACxB,GAAa,EACb,UAAoB,EACpB,CAAS,EACT,MAAoB;IAEpB,MAAM,KAAK,GAAG,aAAa,CAAC,CAAC,EAAE,GAAG,EAAE,UAAU,EAAE,CAAC,EAAE,MAAM,CAAC,CAAC;IAC3D,MAAM,MAAM,GAAG,QAAQ,CAAC,CAAC,EAAE,GAAG,EAAE,UAAU,EAAE,CAAC,EAAE,MAAM,CAAC,CAAC;IACvD,OAAO,YAAY,CAAC,KAAK,EAAE,MAAM,EAAE,CAAC,CAAC,CAAC;AACxC,CAAC;AAED;;;;;;;GAOG;AACH,MAAM,UAAU,eAAe,CAC7B,CAAiB,EACjB,UAAiC,EACjC,CAAS;IAET,MAAM,IAAI,GAAmB,KAAK,CAAC,IAAI,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,EAAE,GAAG,EAAE,CAC1D,KAAK,CAAC,IAAI,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,EAAE,GAAG,EAAE,CAC7B,KAAK,CAAC,IAAI,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,EAAE,GAAG,EAAE,CAAC,IAAI,KAAK,CAAS,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAC9D,CACF,CAAC;IACF,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;QAC3B,KAAK,IAAI,GAAG,GAAG,CAAC,EAAE,GAAG,GAAG,CAAC,EAAE,GAAG,EAAE,EAAE,CAAC;YACjC,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,CAAC,EAAE,EAAE,EAAE,EAAE,CAAC;gBAC9B,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,CAAC,EAAE,EAAE,EAAE,EAAE,CAAC;oBAC9B,IAAI,GAAG,GAAG,CAAC,CAAC;oBACZ,KAAK,IAAI,GAAG,GAAG,CAAC,EAAE,GAAG,GAAG,CAAC,EAAE,GAAG,EAAE,EAAE,CAAC;wBACjC,GAAG,IAAI,UAAU,CAAC,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC;oBACvD,CAAC;oBACD,IAAI,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,GAAG,GAAG,CAAC;gBAC7B,CAAC;YACH,CAAC;QACH,CAAC;IACH,CAAC;IACD,OAAO,IAAI,CAAC;AACd,CAAC;AAED;;;;GAIG;AACH,MAAM,UAAU,cAAc,CAC5B,CAAwB,EACxB,GAAa,EACb,UAAoB,EACpB,CAAS,EACT,MAAoB;IAEpB,MAAM,GAAG,GAAG,cAAc,CAAC,CAAC,EAAE,GAAG,EAAE,UAAU,EAAE,CAAC,EAAE,MAAM,CAAC,CAAC;IAC1D,MAAM,KAAK,GAAG,SAAS,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;IAChC,IAAI,KAAK,CAAC,MAAM,KAAK,CAAC,GAAG,CAAC,EAAE,CAAC;QAC3B,MAAM,IAAI,qBAAqB,CAC7B,yDAAyD,CAAC,GAAG,CAAC,SAAS,KAAK,CAAC,MAAM,EAAE,CACtF,CAAC;IACJ,CAAC;IACD,OAAO,eAAe,CAAC,GAAG,EAAE,KAAK,EAAE,CAAC,CAAC,CAAC;AACxC,CAAC;AAED;;;;;;;;;GASG;AACH,MAAM,UAAU,eAAe,CAC7B,CAAwB,EACxB,GAAa,EACb,UAAoB,EACpB,CAAS,EACT,MAAoB;IAEpB,MAAM,EAAE,GAAqB,KAAK,CAAC,IAAI,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,EAAE,GAAG,EAAE,CAC1D,KAAK,CAAC,IAAI,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,EAAE,GAAG,EAAE,CAC7B,KAAK,CAAC,IAAI,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,EAAE,GAAG,EAAE,CAC7B,KAAK,CAAC,IAAI,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,EAAE,GAAG,EAAE,CAAC,IAAI,KAAK,CAAS,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAC9D,CACF,CACF,CAAC;IAEF,yDAAyD;IACzD,6EAA6E;IAC7E,2EAA2E;IAC3E,6DAA6D;IAC7D,mCAAmC;IACnC,MAAM,QAAQ,GAAG,CAAC,EAAyB,EAAe,EAAE,CAC1D,cAAc,CAAC,EAAE,EAAE,GAAG,EAAE,UAAU,EAAE,CAAC,EAAE,MAAM,CAAC,CAAC;IACjD,MAAM,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC;IACjB,MAAM,EAAE,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;IACrB,KAAK,IAAI,GAAG,GAAG,CAAC,EAAE,GAAG,GAAG,CAAC,EAAE,GAAG,EAAE,EAAE,CAAC;QACjC,MAAM,IAAI,GAAG,iBAAiB,CAAC,QAAQ,EAAE,CAAC,EAAE,GAAG,EAAE,EAAE,KAAK,EAAE,CAAC,EAAE,CAAa,CAAC;QAC3E,IAAI,CAAC,KAAK,CAAC,OAAO,CAAC,IAAI,CAAC,IAAI,IAAI,CAAC,MAAM,KAAK,EAAE,GAAG,EAAE,EAAE,CAAC;YACpD,MAAM,IAAI,qBAAqB,CAC7B,yDAAyD,EAAE,GAAG,EAAE,SAC9D,KAAK,CAAC,OAAO,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,MAAM,CAAC,CAAC,CAAC,QACtC,EAAE,CACH,CAAC;QACJ,CAAC;QACD,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;YAC3B,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,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;wBAC3B,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC,GAAG,EAAE,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC;oBAC1D,CAAC;gBACH,CAAC;YACH,CAAC;QACH,CAAC;IACH,CAAC;IAED,OAAO,EAAE,CAAC;AACZ,CAAC;AAED;;;;;;;;;;;;;;;;GAgBG;AACH,MAAM,UAAU,4BAA4B,CAC1C,CAAwB,EACxB,GAAa,EACb,UAAoB,EACpB,CAAS,EACT,MAAoB;IAEpB,MAAM,KAAK,GAAG,aAAa,CAAC,CAAC,EAAE,GAAG,EAAE,UAAU,EAAE,CAAC,EAAE,MAAM,CAAC,CAAC;IAC3D,MAAM,CAAC,GAAG,cAAc,CAAC,CAAC,EAAE,GAAG,EAAE,UAAU,EAAE,CAAC,EAAE,MAAM,CAAC,CAAC;IACxD,MAAM,EAAE,GAAG,eAAe,CAAC,CAAC,EAAE,GAAG,EAAE,UAAU,EAAE,CAAC,EAAE,MAAM,CAAC,CAAC;IAE1D,MAAM,IAAI,GAAqB,KAAK,CAAC,IAAI,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,EAAE,GAAG,EAAE,CAC5D,KAAK,CAAC,IAAI,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,EAAE,GAAG,EAAE,CAC7B,KAAK,CAAC,IAAI,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,EAAE,GAAG,EAAE,CAC7B,KAAK,CAAC,IAAI,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,EAAE,GAAG,EAAE,CAAC,IAAI,KAAK,CAAS,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAC9D,CACF,CACF,CAAC;IAEF,KAAK,IAAI,GAAG,GAAG,CAAC,EAAE,GAAG,GAAG,CAAC,EAAE,GAAG,EAAE,EAAE,CAAC;QACjC,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,CAAC,EAAE,EAAE,EAAE,EAAE,CAAC;YAC9B,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,CAAC,EAAE,EAAE,EAAE,EAAE,CAAC;gBAC9B,KAAK,IAAI,GAAG,GAAG,CAAC,EAAE,GAAG,GAAG,CAAC,EAAE,GAAG,EAAE,EAAE,CAAC;oBACjC,KAAK,IAAI,GAAG,GAAG,CAAC,EAAE,GAAG,GAAG,CAAC,EAAE,GAAG,EAAE,EAAE,CAAC;wBACjC,IAAI,KAAK,GAAG,EAAE,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC;wBACtC,0DAA0D;wBAC1D,0DAA0D;wBAC1D,0DAA0D;wBAC1D,0DAA0D;wBAC1D,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;4BAC3B,KAAK,IAAI,KAAK,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC;4BAChD,KAAK,IAAI,KAAK,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC;4BAChD,KAAK,IAAI,KAAK,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;4BAChD,KAAK,IAAI,KAAK,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;wBAClD,CAAC;wBACD,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,KAAK,CAAC;oBACtC,CAAC;gBACH,CAAC;YACH,CAAC;QACH,CAAC;IACH,CAAC;IAED,OAAO,IAAI,CAAC;AACd,CAAC;AAED;;;;;;;;;;GAUG;AACH,MAAM,UAAU,iBAAiB,CAC/B,CAAwB,EACxB,GAAa,EACb,UAAoB,EACpB,CAAS,EACT,MAAoB;IAEpB,MAAM,IAAI,GAAG,4BAA4B,CAAC,CAAC,EAAE,GAAG,EAAE,UAAU,EAAE,CAAC,EAAE,MAAM,CAAC,CAAC;IAEzE,MAAM,CAAC,GAAqB,KAAK,CAAC,IAAI,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,EAAE,GAAG,EAAE,CACzD,KAAK,CAAC,IAAI,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,EAAE,GAAG,EAAE,CAC7B,KAAK,CAAC,IAAI,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,EAAE,GAAG,EAAE,CAC7B,KAAK,CAAC,IAAI,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,EAAE,GAAG,EAAE,CAAC,IAAI,KAAK,CAAS,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAC9D,CACF,CACF,CAAC;IAEF,KAAK,IAAI,GAAG,GAAG,CAAC,EAAE,GAAG,GAAG,CAAC,EAAE,GAAG,EAAE,EAAE,CAAC;QACjC,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,CAAC,EAAE,EAAE,EAAE,EAAE,CAAC;YAC9B,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,CAAC,EAAE,EAAE,EAAE,EAAE,CAAC;gBAC9B,KAAK,IAAI,GAAG,GAAG,CAAC,EAAE,GAAG,GAAG,CAAC,EAAE,GAAG,EAAE,EAAE,CAAC;oBACjC,KAAK,IAAI,GAAG,GAAG,CAAC,EAAE,GAAG,GAAG,CAAC,EAAE,GAAG,EAAE,EAAE,CAAC;wBACjC,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC;4BACtB,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAG,eAAe;kCAC3C,IAAI,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,eAAe;kCAC3C,IAAI,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,eAAe;oBAClD,CAAC;gBACH,CAAC;YACH,CAAC;QACH,CAAC;IACH,CAAC;IAED,OAAO,CAAC,CAAC;AACX,CAAC;AAED;;;;;;;;;;;;;;;;;;;;;GAqBG;AACH,MAAM,UAAU,iBAAiB,CAC/B,KAA4B,EAC5B,CAAS,EACT,MAIC;IAED,MAAM,EAAE,SAAS,EAAE,SAAS,EAAE,OAAO,EAAE,GAAG,MAAM,CAAC;IACjD,IAAI,SAAS,KAAK,SAAS,EAAE,CAAC;QAC5B,MAAM,IAAI,qBAAqB,CAC7B,iCAAiC,SAAS,oBAAoB,SAAS,eAAe,CACvF,CAAC;IACJ,CAAC;IACD,MAAM,CAAC,IAAI,EAAE,IAAI,CAAC,GAAG,OAAO,CAAC;IAC7B,IAAI,IAAI,KAAK,IAAI,EAAE,CAAC;QAClB,MAAM,IAAI,qBAAqB,CAC7B,+BAA+B,IAAI,IAAI,IAAI,oBAAoB,CAChE,CAAC;IACJ,CAAC;IACD,IAAI,IAAI,KAAK,SAAS,IAAI,IAAI,KAAK,SAAS,IAAI,IAAI,KAAK,SAAS,IAAI,IAAI,KAAK,SAAS,EAAE,CAAC;QACzF,MAAM,IAAI,qBAAqB,CAC7B,+BAA+B,IAAI,IAAI,IAAI,kCAAkC,SAAS,IAAI,SAAS,GAAG,CACvG,CAAC;IACJ,CAAC;IACD,IAAI,KAAK,CAAC,MAAM,KAAK,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,CAAC;QACnC,MAAM,IAAI,qBAAqB,CAC7B,mCAAmC,KAAK,CAAC,MAAM,YAAY,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,CAC3E,CAAC;IACJ,CAAC;IAED,uCAAuC;IACvC,MAAM,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC;IACjB,MAAM,EAAE,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;IACrB,MAAM,OAAO,GAAG,CAAC,EAAE,EAAE,EAAE,EAAE,CAAC,EAAE,CAAC,CAAU,CAAC;IACxC,MAAM,WAAW,GAAG,OAAO,CAAC,SAAS,CAAC,GAAG,OAAO,CAAC,SAAS,CAAC,CAAC;IAC5D,MAAM,UAAU,GAAG,OAAO,CAAC,IAAI,CAAC,CAAC;IACjC,MAAM,UAAU,GAAG,OAAO,CAAC,IAAI,CAAC,CAAC;IAEjC,MAAM,GAAG,GAAe,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;IACtF,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;QAC3B,MAAM,KAAK,GAAG,CAAC,GAAG,UAAU,CAAC;QAC7B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;YAC3B,MAAM,MAAM,GAAG,KAAK,GAAG,CAAC,GAAG,UAAU,CAAC;YACtC,IAAI,GAAG,GAAG,CAAC,CAAC;YACZ,KAAK,IAAI,GAAG,GAAG,CAAC,EAAE,GAAG,GAAG,CAAC,EAAE,GAAG,EAAE,EAAE,CAAC;gBACjC,GAAG,IAAI,KAAK,CAAC,MAAM,GAAG,GAAG,GAAG,WAAW,CAAC,CAAC;YAC3C,CAAC;YACD,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC;QAClB,CAAC;IACH,CAAC;IACD,OAAO,GAAG,CAAC;AACb,CAAC"}

package/dist/numerical/engine-registry.js

@@ -16,9 +16,9 @@ async function detectDefault() {
         return _resolvedDefault;
     try {
         await import('@danielsimonjr/mathts-tensor');
-        // eslint-disable-next-line @typescript-eslint/ban-ts-comment
-        // @ts-ignore — @danielsimonjr/mathts-autograd is an optional peer dep with
-        // no bundled type declarations; the try-catch handles the absent case.
+        // v0.5.1 TS-4: the prior `@ts-ignore` is replaced by the ambient module
+        // declaration at the top of this file; the try-catch still handles the
+        // runtime-absent case (Promise rejection from the bundler/loader).
         await import('@danielsimonjr/mathts-autograd');
         _resolvedDefault = 'mathts';
     }
@@ -26,7 +26,7 @@ async function detectDefault() {
         // I5 fix: `process` is undefined in browser bundlers; guard before reading.
         const quiet = typeof process !== 'undefined' && process.env?.UPT_QUIET_FALLBACK === '1';
         if (!quiet) {
-            console.warn('[universal-physics-tensor v0.4.0] @danielsimonjr/mathts-tensor and/or '
+            console.warn('[universal-physics-tensor] @danielsimonjr/mathts-tensor and/or '
                 + '@danielsimonjr/mathts-autograd not installed. '
                 + 'Falling back to Float64ReferenceEngine. '
                 + 'Install both peers for MathTS-default + autograd (AD) support. '

package/dist/numerical/engine-registry.js.map

@@ -1 +1 @@
-{"version":3,"file":"engine-registry.js","sourceRoot":"","sources":["../../src/numerical/engine-registry.ts"],"names":[],"mappings":"AAuBA,OAAO,EAAE,sBAAsB,EAAE,MAAM,qBAAqB,CAAC;AAE7D,2EAA2E;AAC3E,6EAA6E;AAC7E,uDAAuD;AACvD,IAAI,cAAiD,CAAC;AAEtD,8EAA8E;AAC9E,mEAAmE;AACnE,IAAI,gBAAkD,CAAC;AAEvD;;;;GAIG;AACH,KAAK,UAAU,aAAa;IAC1B,IAAI,gBAAgB,KAAK,SAAS;QAAE,OAAO,gBAAgB,CAAC;IAC5D,IAAI,CAAC;QACH,MAAM,MAAM,CAAC,8BAA8B,CAAC,CAAC;QAC7C,6DAA6D;QAC7D,2EAA2E;QAC3E,uEAAuE;QACvE,MAAM,MAAM,CAAC,gCAAgC,CAAC,CAAC;QAC/C,gBAAgB,GAAG,QAAQ,CAAC;IAC9B,CAAC;IAAC,MAAM,CAAC;QACP,4EAA4E;QAC5E,MAAM,KAAK,GACT,OAAO,OAAO,KAAK,WAAW,IAAI,OAAO,CAAC,GAAG,EAAE,kBAAkB,KAAK,GAAG,CAAC;QAC5E,IAAI,CAAC,KAAK,EAAE,CAAC;YACX,OAAO,CAAC,IAAI,CACV,wEAAwE;kBACtE,gDAAgD;kBAChD,0CAA0C;kBAC1C,iEAAiE;kBACjE,oDAAoD,CACvD,CAAC;QACJ,CAAC;QACD,gBAAgB,GAAG,SAAS,CAAC;IAC/B,CAAC;IACD,OAAO,gBAAgB,CAAC;AAC1B,CAAC;AAED,2EAA2E;AAC3E,gEAAgE;AAEhE;;;;;;;;;;;GAWG;AACH,MAAM,CAAC,KAAK,UAAU,eAAe;IACnC,cAAc,KAAK,CAAC,KAAK,IAAI,EAAE;QAC7B,MAAM,MAAM,GAAG,MAAM,aAAa,EAAE,CAAC;QACrC,IAAI,MAAM,KAAK,QAAQ,EAAE,CAAC;YACxB,MAAM,EAAE,YAAY,EAAE,GAAG,MAAM,MAAM,CAAC,oBAAoB,CAAC,CAAC;YAC5D,OAAO,IAAI,YAAY,EAAE,CAAC;QAC5B,CAAC;QACD,OAAO,IAAI,sBAAsB,EAAE,CAAC;IACtC,CAAC,CAAC,EAAE,CAAC;IACL,OAAO,cAAc,CAAC;AACxB,CAAC;AAED;;;;;GAKG;AACH,MAAM,UAAU,eAAe,CAAC,MAAoB;IAClD,cAAc,GAAG,OAAO,CAAC,OAAO,CAAC,MAAM,CAAC,CAAC;AAC3C,CAAC;AAED;;;;GAIG;AACH,MAAM,UAAU,qBAAqB;IACnC,cAAc,GAAG,SAAS,CAAC;IAC3B,gBAAgB,GAAG,SAAS,CAAC;AAC/B,CAAC"}
+{"version":3,"file":"engine-registry.js","sourceRoot":"","sources":["../../src/numerical/engine-registry.ts"],"names":[],"mappings":"AAuBA,OAAO,EAAE,sBAAsB,EAAE,MAAM,qBAAqB,CAAC;AAE7D,2EAA2E;AAC3E,6EAA6E;AAC7E,uDAAuD;AACvD,IAAI,cAAiD,CAAC;AAEtD,8EAA8E;AAC9E,mEAAmE;AACnE,IAAI,gBAAkD,CAAC;AAEvD;;;;GAIG;AACH,KAAK,UAAU,aAAa;IAC1B,IAAI,gBAAgB,KAAK,SAAS;QAAE,OAAO,gBAAgB,CAAC;IAC5D,IAAI,CAAC;QACH,MAAM,MAAM,CAAC,8BAA8B,CAAC,CAAC;QAC7C,wEAAwE;QACxE,uEAAuE;QACvE,mEAAmE;QACnE,MAAM,MAAM,CAAC,gCAAgC,CAAC,CAAC;QAC/C,gBAAgB,GAAG,QAAQ,CAAC;IAC9B,CAAC;IAAC,MAAM,CAAC;QACP,4EAA4E;QAC5E,MAAM,KAAK,GACT,OAAO,OAAO,KAAK,WAAW,IAAI,OAAO,CAAC,GAAG,EAAE,kBAAkB,KAAK,GAAG,CAAC;QAC5E,IAAI,CAAC,KAAK,EAAE,CAAC;YACX,OAAO,CAAC,IAAI,CACV,iEAAiE;kBAC/D,gDAAgD;kBAChD,0CAA0C;kBAC1C,iEAAiE;kBACjE,oDAAoD,CACvD,CAAC;QACJ,CAAC;QACD,gBAAgB,GAAG,SAAS,CAAC;IAC/B,CAAC;IACD,OAAO,gBAAgB,CAAC;AAC1B,CAAC;AAED,2EAA2E;AAC3E,gEAAgE;AAEhE;;;;;;;;;;;GAWG;AACH,MAAM,CAAC,KAAK,UAAU,eAAe;IACnC,cAAc,KAAK,CAAC,KAAK,IAAI,EAAE;QAC7B,MAAM,MAAM,GAAG,MAAM,aAAa,EAAE,CAAC;QACrC,IAAI,MAAM,KAAK,QAAQ,EAAE,CAAC;YACxB,MAAM,EAAE,YAAY,EAAE,GAAG,MAAM,MAAM,CAAC,oBAAoB,CAAC,CAAC;YAC5D,OAAO,IAAI,YAAY,EAAE,CAAC;QAC5B,CAAC;QACD,OAAO,IAAI,sBAAsB,EAAE,CAAC;IACtC,CAAC,CAAC,EAAE,CAAC;IACL,OAAO,cAAc,CAAC;AACxB,CAAC;AAED;;;;;GAKG;AACH,MAAM,UAAU,eAAe,CAAC,MAAoB;IAClD,cAAc,GAAG,OAAO,CAAC,OAAO,CAAC,MAAM,CAAC,CAAC;AAC3C,CAAC;AAED;;;;GAIG;AACH,MAAM,UAAU,qBAAqB;IACnC,cAAc,GAAG,SAAS,CAAC;IAC3B,gBAAgB,GAAG,SAAS,CAAC;AAC/B,CAAC"}

package/dist/numerical/errors.d.ts

@@ -25,4 +25,20 @@ export declare class EngineCapabilityError extends NumericalBackendError {
     readonly missingMethod: string;
     constructor(engineName: string, missingMethod: string);
 }
+/**
+ * Thrown when the GL4 implicit-stage Picard fixed-point solver fails to
+ * reach `picardTol` within `picardMaxIter` iterations. Carries a message
+ * matching `/Picard iteration did not converge/` for I7-style specific
+ * error-class assertions.
+ *
+ * Picard convergence is linear with contraction rate ≈ h·|∂f/∂x|; failure
+ * usually indicates the step size h is too large for the local curvature.
+ * Adaptive step-halving (Task 5) is the production response — direct
+ * callers of `solveGL4Stage` should adjust h or picardMaxIter.
+ *
+ * @public
+ */
+export declare class GL4ConvergenceError extends NumericalBackendError {
+    constructor(message: string);
+}
 //# sourceMappingURL=errors.d.ts.map

package/dist/numerical/errors.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"errors.d.ts","sourceRoot":"","sources":["../../src/numerical/errors.ts"],"names":[],"mappings":"AAAA;;;;;GAKG;AACH,OAAO,EAAE,QAAQ,EAAE,MAAM,0BAA0B,CAAC;AAEpD;;;;GAIG;AACH,qBAAa,qBAAsB,SAAQ,QAAQ;gBACrC,OAAO,EAAE,MAAM;CAK5B;AAED;;;;;;GAMG;AACH,qBAAa,qBAAsB,SAAQ,qBAAqB;aAE5C,UAAU,EAAE,MAAM;aAClB,aAAa,EAAE,MAAM;gBADrB,UAAU,EAAE,MAAM,EAClB,aAAa,EAAE,MAAM;CASxC"}
+{"version":3,"file":"errors.d.ts","sourceRoot":"","sources":["../../src/numerical/errors.ts"],"names":[],"mappings":"AAAA;;;;;GAKG;AACH,OAAO,EAAE,QAAQ,EAAE,MAAM,0BAA0B,CAAC;AAEpD;;;;GAIG;AACH,qBAAa,qBAAsB,SAAQ,QAAQ;gBACrC,OAAO,EAAE,MAAM;CAK5B;AAED;;;;;;GAMG;AACH,qBAAa,qBAAsB,SAAQ,qBAAqB;aAE5C,UAAU,EAAE,MAAM;aAClB,aAAa,EAAE,MAAM;gBADrB,UAAU,EAAE,MAAM,EAClB,aAAa,EAAE,MAAM;CASxC;AAED;;;;;;;;;;;;GAYG;AACH,qBAAa,mBAAoB,SAAQ,qBAAqB;gBAChD,OAAO,EAAE,MAAM;CAK5B"}

package/dist/numerical/errors.js

@@ -36,4 +36,24 @@ export class EngineCapabilityError extends NumericalBackendError {
         Object.setPrototypeOf(this, EngineCapabilityError.prototype);
     }
 }
+/**
+ * Thrown when the GL4 implicit-stage Picard fixed-point solver fails to
+ * reach `picardTol` within `picardMaxIter` iterations. Carries a message
+ * matching `/Picard iteration did not converge/` for I7-style specific
+ * error-class assertions.
+ *
+ * Picard convergence is linear with contraction rate ≈ h·|∂f/∂x|; failure
+ * usually indicates the step size h is too large for the local curvature.
+ * Adaptive step-halving (Task 5) is the production response — direct
+ * callers of `solveGL4Stage` should adjust h or picardMaxIter.
+ *
+ * @public
+ */
+export class GL4ConvergenceError extends NumericalBackendError {
+    constructor(message) {
+        super(message);
+        this.name = 'GL4ConvergenceError';
+        Object.setPrototypeOf(this, GL4ConvergenceError.prototype);
+    }
+}
 //# sourceMappingURL=errors.js.map

package/dist/numerical/errors.js.map

@@ -1 +1 @@
-{"version":3,"file":"errors.js","sourceRoot":"","sources":["../../src/numerical/errors.ts"],"names":[],"mappings":"AAAA;;;;;GAKG;AACH,OAAO,EAAE,QAAQ,EAAE,MAAM,0BAA0B,CAAC;AAEpD;;;;GAIG;AACH,MAAM,OAAO,qBAAsB,SAAQ,QAAQ;IACjD,YAAY,OAAe;QACzB,KAAK,CAAC,OAAO,CAAC,CAAC;QACf,IAAI,CAAC,IAAI,GAAG,uBAAuB,CAAC;QACpC,MAAM,CAAC,cAAc,CAAC,IAAI,EAAE,qBAAqB,CAAC,SAAS,CAAC,CAAC;IAC/D,CAAC;CACF;AAED;;;;;;GAMG;AACH,MAAM,OAAO,qBAAsB,SAAQ,qBAAqB;IAE5C;IACA;IAFlB,YACkB,UAAkB,EAClB,aAAqB;QAErC,KAAK,CACH,WAAW,UAAU,yBAAyB,aAAa,KAAK;YAChE,wEAAwE,CACzE,CAAC;QANc,eAAU,GAAV,UAAU,CAAQ;QAClB,kBAAa,GAAb,aAAa,CAAQ;QAMrC,IAAI,CAAC,IAAI,GAAG,uBAAuB,CAAC;QACpC,MAAM,CAAC,cAAc,CAAC,IAAI,EAAE,qBAAqB,CAAC,SAAS,CAAC,CAAC;IAC/D,CAAC;CACF"}
+{"version":3,"file":"errors.js","sourceRoot":"","sources":["../../src/numerical/errors.ts"],"names":[],"mappings":"AAAA;;;;;GAKG;AACH,OAAO,EAAE,QAAQ,EAAE,MAAM,0BAA0B,CAAC;AAEpD;;;;GAIG;AACH,MAAM,OAAO,qBAAsB,SAAQ,QAAQ;IACjD,YAAY,OAAe;QACzB,KAAK,CAAC,OAAO,CAAC,CAAC;QACf,IAAI,CAAC,IAAI,GAAG,uBAAuB,CAAC;QACpC,MAAM,CAAC,cAAc,CAAC,IAAI,EAAE,qBAAqB,CAAC,SAAS,CAAC,CAAC;IAC/D,CAAC;CACF;AAED;;;;;;GAMG;AACH,MAAM,OAAO,qBAAsB,SAAQ,qBAAqB;IAE5C;IACA;IAFlB,YACkB,UAAkB,EAClB,aAAqB;QAErC,KAAK,CACH,WAAW,UAAU,yBAAyB,aAAa,KAAK;YAChE,wEAAwE,CACzE,CAAC;QANc,eAAU,GAAV,UAAU,CAAQ;QAClB,kBAAa,GAAb,aAAa,CAAQ;QAMrC,IAAI,CAAC,IAAI,GAAG,uBAAuB,CAAC;QACpC,MAAM,CAAC,cAAc,CAAC,IAAI,EAAE,qBAAqB,CAAC,SAAS,CAAC,CAAC;IAC/D,CAAC;CACF;AAED;;;;;;;;;;;;GAYG;AACH,MAAM,OAAO,mBAAoB,SAAQ,qBAAqB;IAC5D,YAAY,OAAe;QACzB,KAAK,CAAC,OAAO,CAAC,CAAC;QACf,IAAI,CAAC,IAAI,GAAG,qBAAqB,CAAC;QAClC,MAAM,CAAC,cAAc,CAAC,IAAI,EAAE,mBAAmB,CAAC,SAAS,CAAC,CAAC;IAC7D,CAAC;CACF"}