universal-physics-tensor 0.4.6 → 0.5.0

package/dist/numerical/curvature-lowering-helpers.d.ts

@@ -0,0 +1,162 @@
+/**
+ * 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 type { TensorEngine } from './tensor-engine.js';
+import type { NestedArray } from './types.js';
+/** Flat row-major N×N matrix (one g or g^{-1} sample). */
+export type FlatMatrix = ReadonlyArray<number>;
+/** A coordinate-dependent rank-2 closure: f(x) → N×N as a NestedArray
+ *  (number[][] or any nested form that `flattenNA` can flatten). */
+export type MetricFn = (x: ReadonlyArray<number>) => NestedArray;
+/**
+ * dGamma[λ][ρ][σ][ν] = ∂_λ Γ^ρ_{σν} at the given coordinates (I3).
+ *
+ * 4-deep nested number array — readonly to mark this as an output container
+ * the helper produces but does not mutate after construction.
+ */
+export type DGammaTensor = readonly (readonly (readonly (readonly number[])[])[])[];
+/** Gamma[ρ][σ][ν] = Γ^ρ_{σν} — same index order as `computeChristoffelTensor`
+ *  output (just renamed for the Riemann use). */
+export type GammaTensor = readonly (readonly (readonly number[])[])[];
+/**
+ * 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 declare function outerStep(x: number): number;
+/**
+ * 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 declare function christoffelAt(x: ReadonlyArray<number>, gFn: MetricFn, gInverseFn: MetricFn, N: number, engine: TensorEngine): number[][][];
+/**
+ * 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 declare function dGammaAt(x: ReadonlyArray<number>, gFn: MetricFn, gInverseFn: MetricFn, N: number, engine: TensorEngine): DGammaTensor;
+/**
+ * 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 declare function buildRiemann(gamma: GammaTensor, dGamma: DGammaTensor, N: number): number[][][][];
+/**
+ * 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 declare function riemannUpperAt(x: ReadonlyArray<number>, gFn: MetricFn, gInverseFn: MetricFn, N: number, engine: TensorEngine): number[][][][];
+/**
+ * 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 declare function lowerFirstIndex(R: number[][][][], gLowerFlat: ReadonlyArray<number>, N: number): number[][][][];
+/**
+ * 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 declare function riemannLowerAt(x: ReadonlyArray<number>, gFn: MetricFn, gInverseFn: MetricFn, N: number, engine: TensorEngine): number[][][][];
+/**
+ * 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 declare function dRiemannLowerAt(x: ReadonlyArray<number>, gFn: MetricFn, gInverseFn: MetricFn, N: number, engine: TensorEngine): number[][][][][];
+/**
+ * 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 declare function covariantDerivRiemannLowerAt(x: ReadonlyArray<number>, gFn: MetricFn, gInverseFn: MetricFn, N: number, engine: TensorEngine): number[][][][][];
+/**
+ * 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 declare function bianchiResidualAt(x: ReadonlyArray<number>, gFn: MetricFn, gInverseFn: MetricFn, N: number, engine: TensorEngine): number[][][][][];
+//# sourceMappingURL=curvature-lowering-helpers.d.ts.map

package/dist/numerical/curvature-lowering-helpers.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"curvature-lowering-helpers.d.ts","sourceRoot":"","sources":["../../src/numerical/curvature-lowering-helpers.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GA6BG;AAEH,OAAO,KAAK,EAAE,YAAY,EAAE,MAAM,oBAAoB,CAAC;AACvD,OAAO,KAAK,EAAE,WAAW,EAAE,MAAM,YAAY,CAAC;AAI9C,0DAA0D;AAC1D,MAAM,MAAM,UAAU,GAAG,aAAa,CAAC,MAAM,CAAC,CAAC;AAE/C;oEACoE;AACpE,MAAM,MAAM,QAAQ,GAAG,CAAC,CAAC,EAAE,aAAa,CAAC,MAAM,CAAC,KAAK,WAAW,CAAC;AAEjE;;;;;GAKG;AACH,MAAM,MAAM,YAAY,GAAG,SAAS,CAAC,SAAS,CAAC,SAAS,CAAC,SAAS,MAAM,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC;AAEpF;iDACiD;AACjD,MAAM,MAAM,WAAW,GAAG,SAAS,CAAC,SAAS,CAAC,SAAS,MAAM,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC;AAMtE;;;;;;;;;;;;;;;;;;;;GAoBG;AACH,wBAAgB,SAAS,CAAC,CAAC,EAAE,MAAM,GAAG,MAAM,CAE3C;AAmDD;;;;;;;;GAQG;AACH,wBAAgB,aAAa,CAC3B,CAAC,EAAE,aAAa,CAAC,MAAM,CAAC,EACxB,GAAG,EAAE,QAAQ,EACb,UAAU,EAAE,QAAQ,EACpB,CAAC,EAAE,MAAM,EACT,MAAM,EAAE,YAAY,GACnB,MAAM,EAAE,EAAE,EAAE,CAiBd;AAED;;;;;;GAMG;AACH,wBAAgB,QAAQ,CACtB,CAAC,EAAE,aAAa,CAAC,MAAM,CAAC,EACxB,GAAG,EAAE,QAAQ,EACb,UAAU,EAAE,QAAQ,EACpB,CAAC,EAAE,MAAM,EACT,MAAM,EAAE,YAAY,GACnB,YAAY,CAgDd;AAED;;;;;;;;GAQG;AACH,wBAAgB,YAAY,CAC1B,KAAK,EAAE,WAAW,EAClB,MAAM,EAAE,YAAY,EACpB,CAAC,EAAE,MAAM,GACR,MAAM,EAAE,EAAE,EAAE,EAAE,CAyBhB;AAMD;;;;;;GAMG;AACH,wBAAgB,cAAc,CAC5B,CAAC,EAAE,aAAa,CAAC,MAAM,CAAC,EACxB,GAAG,EAAE,QAAQ,EACb,UAAU,EAAE,QAAQ,EACpB,CAAC,EAAE,MAAM,EACT,MAAM,EAAE,YAAY,GACnB,MAAM,EAAE,EAAE,EAAE,EAAE,CAIhB;AAED;;;;;;;GAOG;AACH,wBAAgB,eAAe,CAC7B,CAAC,EAAE,MAAM,EAAE,EAAE,EAAE,EAAE,EACjB,UAAU,EAAE,aAAa,CAAC,MAAM,CAAC,EACjC,CAAC,EAAE,MAAM,GACR,MAAM,EAAE,EAAE,EAAE,EAAE,CAoBhB;AAED;;;;GAIG;AACH,wBAAgB,cAAc,CAC5B,CAAC,EAAE,aAAa,CAAC,MAAM,CAAC,EACxB,GAAG,EAAE,QAAQ,EACb,UAAU,EAAE,QAAQ,EACpB,CAAC,EAAE,MAAM,EACT,MAAM,EAAE,YAAY,GACnB,MAAM,EAAE,EAAE,EAAE,EAAE,CAShB;AAED;;;;;;;;;GASG;AACH,wBAAgB,eAAe,CAC7B,CAAC,EAAE,aAAa,CAAC,MAAM,CAAC,EACxB,GAAG,EAAE,QAAQ,EACb,UAAU,EAAE,QAAQ,EACpB,CAAC,EAAE,MAAM,EACT,MAAM,EAAE,YAAY,GACnB,MAAM,EAAE,EAAE,EAAE,EAAE,EAAE,CAwClB;AAED;;;;;;;;;;;;;;;;GAgBG;AACH,wBAAgB,4BAA4B,CAC1C,CAAC,EAAE,aAAa,CAAC,MAAM,CAAC,EACxB,GAAG,EAAE,QAAQ,EACb,UAAU,EAAE,QAAQ,EACpB,CAAC,EAAE,MAAM,EACT,MAAM,EAAE,YAAY,GACnB,MAAM,EAAE,EAAE,EAAE,EAAE,EAAE,CAqClB;AAED;;;;;;;;;;GAUG;AACH,wBAAgB,iBAAiB,CAC/B,CAAC,EAAE,aAAa,CAAC,MAAM,CAAC,EACxB,GAAG,EAAE,QAAQ,EACb,UAAU,EAAE,QAAQ,EACpB,CAAC,EAAE,MAAM,EACT,MAAM,EAAE,YAAY,GACnB,MAAM,EAAE,EAAE,EAAE,EAAE,EAAE,CA2BlB"}

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

@@ -0,0 +1,378 @@
+/**
+ * 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';
+// ---------------------------------------------------------------------------
+// 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) {
+    return (mu) => {
+        const xc = x0[mu];
+        const h = 1e-3 * Math.max(Math.abs(xc), 1);
+        const p1 = [...x0];
+        p1[mu] = xc + h;
+        const m1 = [...x0];
+        m1[mu] = xc - h;
+        const p2 = [...x0];
+        p2[mu] = xc + 2 * h;
+        const m2 = [...x0];
+        m2[mu] = xc - 2 * h;
+        const gp1 = flattenNA(gFn(p1));
+        const gm1 = flattenNA(gFn(m1));
+        const gp2 = flattenNA(gFn(p2));
+        const gm2 = flattenNA(gFn(m2));
+        if (gp1.length !== N * N || gm1.length !== N * N ||
+            gp2.length !== N * N || gm2.length !== N * N) {
+            throw new NumericalBackendError(`curvature-lowering: metric closure returned wrong shape — expected ${N * N}, ` +
+                `got [${gp1.length}, ${gm1.length}, ${gp2.length}, ${gm2.length}]`);
+        }
+        const inv12h = 1 / (12 * h);
+        const out = new Array(N * N);
+        for (let i = 0; i < N * N; i++) {
+            out[i] = (-gp2[i] + 8 * gp1[i] - 8 * gm1[i] + gm2[i]) * inv12h;
+        }
+        return out;
+    };
+}
+/**
+ * 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))));
+    // 4th-order centered stencil: f'(x) ≈ [−f₊₂ + 8 f₊₁ − 8 f₋₁ + f₋₂] / (12 h)
+    for (let lam = 0; lam < N; lam++) {
+        const xc = x[lam];
+        const h = outerStep(xc);
+        const xP1 = [...x];
+        xP1[lam] = xc + h;
+        const xM1 = [...x];
+        xM1[lam] = xc - h;
+        const xP2 = [...x];
+        xP2[lam] = xc + 2 * h;
+        const xM2 = [...x];
+        xM2[lam] = xc - 2 * h;
+        const Gp1 = christoffelAt(xP1, gFn, gInverseFn, N, engine);
+        const Gm1 = christoffelAt(xM1, gFn, gInverseFn, N, engine);
+        const Gp2 = christoffelAt(xP2, gFn, gInverseFn, N, engine);
+        const Gm2 = christoffelAt(xM2, gFn, gInverseFn, N, engine);
+        const inv12h = 1 / (12 * h);
+        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] = (-Gp2[rho][sigma][nu]
+                        + 8 * Gp1[rho][sigma][nu]
+                        - 8 * Gm1[rho][sigma][nu]
+                        + Gm2[rho][sigma][nu]) * inv12h;
+                }
+            }
+        }
+    }
+    // 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)))));
+    for (let lam = 0; lam < N; lam++) {
+        const xc = x[lam];
+        const h = outerStep(xc);
+        const xP1 = [...x];
+        xP1[lam] = xc + h;
+        const xM1 = [...x];
+        xM1[lam] = xc - h;
+        const xP2 = [...x];
+        xP2[lam] = xc + 2 * h;
+        const xM2 = [...x];
+        xM2[lam] = xc - 2 * h;
+        const Rp1 = riemannLowerAt(xP1, gFn, gInverseFn, N, engine);
+        const Rm1 = riemannLowerAt(xM1, gFn, gInverseFn, N, engine);
+        const Rp2 = riemannLowerAt(xP2, gFn, gInverseFn, N, engine);
+        const Rm2 = riemannLowerAt(xM2, gFn, gInverseFn, N, engine);
+        const inv12h = 1 / (12 * h);
+        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] = (-Rp2[a][b][c][d]
+                            + 8 * Rp1[a][b][c][d]
+                            - 8 * Rm1[a][b][c][d]
+                            + Rm2[a][b][c][d]) * inv12h;
+                    }
+                }
+            }
+        }
+    }
+    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;
+}
+//# 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;AAqBvF,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,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,EAAE,GAAG,CAAC,GAAG,EAAE,CAAC,CAAC;QAAC,EAAE,CAAC,EAAE,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC;QACpC,MAAM,EAAE,GAAG,CAAC,GAAG,EAAE,CAAC,CAAC;QAAC,EAAE,CAAC,EAAE,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC;QACpC,MAAM,EAAE,GAAG,CAAC,GAAG,EAAE,CAAC,CAAC;QAAC,EAAE,CAAC,EAAE,CAAC,GAAG,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC;QACxC,MAAM,EAAE,GAAG,CAAC,GAAG,EAAE,CAAC,CAAC;QAAC,EAAE,CAAC,EAAE,CAAC,GAAG,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC;QACxC,MAAM,GAAG,GAAG,SAAS,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,CAAC;QAC/B,MAAM,GAAG,GAAG,SAAS,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,CAAC;QAC/B,MAAM,GAAG,GAAG,SAAS,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,CAAC;QAC/B,MAAM,GAAG,GAAG,SAAS,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,CAAC;QAC/B,IACE,GAAG,CAAC,MAAM,KAAK,CAAC,GAAG,CAAC,IAAI,GAAG,CAAC,MAAM,KAAK,CAAC,GAAG,CAAC;YAC5C,GAAG,CAAC,MAAM,KAAK,CAAC,GAAG,CAAC,IAAI,GAAG,CAAC,MAAM,KAAK,CAAC,GAAG,CAAC,EAC5C,CAAC;YACD,MAAM,IAAI,qBAAqB,CAC7B,sEAAsE,CAAC,GAAG,CAAC,IAAI;gBAC/E,QAAQ,GAAG,CAAC,MAAM,KAAK,GAAG,CAAC,MAAM,KAAK,GAAG,CAAC,MAAM,KAAK,GAAG,CAAC,MAAM,GAAG,CACnE,CAAC;QACJ,CAAC;QACD,MAAM,MAAM,GAAG,CAAC,GAAG,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC;QAC5B,MAAM,GAAG,GAAG,IAAI,KAAK,CAAS,CAAC,GAAG,CAAC,CAAC,CAAC;QACrC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;YAC/B,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC;QACjE,CAAC;QACD,OAAO,GAAG,CAAC;IACb,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,4EAA4E;IAC5E,KAAK,IAAI,GAAG,GAAG,CAAC,EAAE,GAAG,GAAG,CAAC,EAAE,GAAG,EAAE,EAAE,CAAC;QACjC,MAAM,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC;QAClB,MAAM,CAAC,GAAG,SAAS,CAAC,EAAE,CAAC,CAAC;QACxB,MAAM,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC;QAAC,GAAG,CAAC,GAAG,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC;QACtC,MAAM,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC;QAAC,GAAG,CAAC,GAAG,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC;QACtC,MAAM,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC;QAAC,GAAG,CAAC,GAAG,CAAC,GAAG,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC;QAC1C,MAAM,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC;QAAC,GAAG,CAAC,GAAG,CAAC,GAAG,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC;QAC1C,MAAM,GAAG,GAAG,aAAa,CAAC,GAAG,EAAE,GAAG,EAAE,UAAU,EAAE,CAAC,EAAE,MAAM,CAAC,CAAC;QAC3D,MAAM,GAAG,GAAG,aAAa,CAAC,GAAG,EAAE,GAAG,EAAE,UAAU,EAAE,CAAC,EAAE,MAAM,CAAC,CAAC;QAC3D,MAAM,GAAG,GAAG,aAAa,CAAC,GAAG,EAAE,GAAG,EAAE,UAAU,EAAE,CAAC,EAAE,MAAM,CAAC,CAAC;QAC3D,MAAM,GAAG,GAAG,aAAa,CAAC,GAAG,EAAE,GAAG,EAAE,UAAU,EAAE,CAAC,EAAE,MAAM,CAAC,CAAC;QAC3D,MAAM,MAAM,GAAG,CAAC,GAAG,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC;QAC5B,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,CAC5B,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,KAAK,CAAC,CAAC,EAAE,CAAC;0BAClB,CAAC,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC,KAAK,CAAC,CAAC,EAAE,CAAC;0BACvB,CAAC,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC,KAAK,CAAC,CAAC,EAAE,CAAC;0BACvB,GAAG,CAAC,GAAG,CAAC,CAAC,KAAK,CAAC,CAAC,EAAE,CAAC,CACtB,GAAG,MAAM,CAAC;gBACb,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,KAAK,IAAI,GAAG,GAAG,CAAC,EAAE,GAAG,GAAG,CAAC,EAAE,GAAG,EAAE,EAAE,CAAC;QACjC,MAAM,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC;QAClB,MAAM,CAAC,GAAG,SAAS,CAAC,EAAE,CAAC,CAAC;QACxB,MAAM,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC;QAAC,GAAG,CAAC,GAAG,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC;QACtC,MAAM,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC;QAAC,GAAG,CAAC,GAAG,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC;QACtC,MAAM,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC;QAAC,GAAG,CAAC,GAAG,CAAC,GAAG,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC;QAC1C,MAAM,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC;QAAC,GAAG,CAAC,GAAG,CAAC,GAAG,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC;QAE1C,MAAM,GAAG,GAAG,cAAc,CAAC,GAAG,EAAE,GAAG,EAAE,UAAU,EAAE,CAAC,EAAE,MAAM,CAAC,CAAC;QAC5D,MAAM,GAAG,GAAG,cAAc,CAAC,GAAG,EAAE,GAAG,EAAE,UAAU,EAAE,CAAC,EAAE,MAAM,CAAC,CAAC;QAC5D,MAAM,GAAG,GAAG,cAAc,CAAC,GAAG,EAAE,GAAG,EAAE,UAAU,EAAE,CAAC,EAAE,MAAM,CAAC,CAAC;QAC5D,MAAM,GAAG,GAAG,cAAc,CAAC,GAAG,EAAE,GAAG,EAAE,UAAU,EAAE,CAAC,EAAE,MAAM,CAAC,CAAC;QAC5D,MAAM,MAAM,GAAG,CAAC,GAAG,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC;QAE5B,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,CACpB,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;8BACd,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;8BACnB,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;8BACnB,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAClB,GAAG,MAAM,CAAC;oBACb,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"}

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"}