universal-physics-tensor 8.0.0 → 10.0.0

package/dist/numerical/gl4-integrator.d.ts

@@ -26,8 +26,9 @@ export interface GL4State {
 }
 /**
  * Per-step snapshot recorded by the integrator. `v` (the contravariant
- * 4-velocity v^μ = g^{μν} p_ν) is optional — emitted when the caller asks
- * for it, since it requires an extra metric-inverse contraction.
+ * 4-velocity v^μ = g^{μν} p_ν) is declared optional and is never emitted:
+ * no option asks for it, and nothing in `src/` reads it. It stays so a caller's
+ * own snapshot type can carry one.
  *
  * @public
  */
@@ -90,34 +91,27 @@ export interface GL4Options {
     }) => void;
 }
 /**
- * Result of `solveGL4Stage` — the two converged stage values plus the
- * iteration count actually consumed. Consumed by the
- * `integrateGeodesicGL4` step driver.
- *
- * v0.6.1: dropped export — internal-only result shape (was already
- * @internal-tagged but had no external consumer).
+ * The state after one GL4 step, with the Picard iteration count that step
+ * consumed. MathTS `gaussLegendre4` does not expose its two internal stage
+ * values, so this is the advanced state and nothing else.
+ * @internal
  */
-interface StageSolveResult {
-    readonly stageX: readonly [readonly number[], readonly number[]];
-    readonly stageP: readonly [readonly number[], readonly number[]];
-    readonly stageDx: readonly [readonly number[], readonly number[]];
-    readonly stageDp: readonly [readonly number[], readonly number[]];
+export interface GL4StepResult {
+    readonly x: readonly number[];
+    readonly p: readonly number[];
     readonly iterations: number;
 }
 /**
- * One GL4 step, reported in the stage-result shape the stage tests read.
- *
- * The step itself is MathTS `gaussLegendre4`. In flat space `p` is constant,
- * so both stage momenta equal that advanced `p`. Throws `GL4ConvergenceError`
- * with a message matching `/Picard iteration did not converge/` when the
- * Picard cap is exhausted.
- *
+ * One GL4 step of the geodesic Hamiltonian, via MathTS `gaussLegendre4`.
+ * Throws `GL4ConvergenceError` with a message matching
+ * `/Picard iteration did not converge/` when the Picard cap is exhausted;
+ * any other error from the metric closures propagates unchanged.
  * @internal
  */
-export declare function solveGL4Stage(state: GL4State, h: number, gInverseFn: (x: readonly number[]) => Float64Array, dgInverseFn: (x: readonly number[]) => Float64Array, opts: {
+export declare function gl4Step(state: GL4State, h: number, gInverseFn: (x: readonly number[]) => Float64Array, dgInverseFn: (x: readonly number[]) => Float64Array, opts: {
     picardTol: number;
     picardMaxIter: number;
-}): StageSolveResult;
+}): GL4StepResult;
 /**
  * GL4 symplectic integrator on the canonical (x, p) geodesic Hamiltonian.
  *
@@ -163,10 +157,11 @@ export declare function solveGL4Stage(state: GL4State, h: number, gInverseFn: (x
  * @param initialState — canonical (x, p) at τ = 0. Units follow the
  *   `gInverseFn` convention (see above).
  * @param options — see {@link GL4Options}:
- *   - `steps` — integer step count (dimensionless).
- *   - `tauMax` — affine-parameter extent (seconds in canonical SI).
- *   - `gInverseFn(x)[μ][ν]` — inverse metric g^{μν}(x).
- *   - `dgInverseFn(x)[λ][μ][ν]` — ∂_λ g^{μν}(x).
+ *   - `steps` — positive integer step count (dimensionless).
+ *   - `tauMax` — positive finite affine-parameter extent (seconds in
+ *     canonical SI). The integrator runs forward from τ = 0 only.
+ *   - `gInverseFn(x)[μ*dim + ν]` — inverse metric g^{μν}(x), row-major flat.
+ *   - `dgInverseFn(x)[λ*dim² + μ*dim + ν]` — ∂_λ g^{μν}(x), flat.
  *   - `picardTol` — convergence tolerance (dimensionless, default 1e-12).
  *   - `picardMaxIter` — fixed-point iteration cap (dimensionless, default 50).
  *   - `hMin` — step-halving floor (same units as `tauMax / steps`).
@@ -174,8 +169,14 @@ export declare function solveGL4Stage(state: GL4State, h: number, gInverseFn: (x
  * @returns `steps + 1` snapshots: index `0` is the initial state, index `n`
  *   is the state after `n` steps (τ = n · h). Each snapshot carries `tau`,
  *   `x`, `p` (and optional `v` = g^{μν} p_ν) in the units chosen above.
- * @throws NumericalBackendError if `initialState.x[1] < domainMinRadius`.
+ * @throws NumericalBackendError if `steps` is not a positive integer, if
+ *   `tauMax` is not a positive finite number, or if
+ *   `initialState.x[1] < domainMinRadius`. The two option checks run before
+ *   any metric evaluation.
  * @throws GL4ConvergenceError if Picard fails even after step-halving to h_min.
+ *   Any other error thrown by `gInverseFn` or `dgInverseFn` (a domain
+ *   crossing the caller signals by throwing, or a programming error) is
+ *   rethrown unchanged; it is not halved away.
  *
  * @public
  */

package/dist/numerical/gl4-integrator.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"gl4-integrator.d.ts","sourceRoot":"","sources":["../../src/numerical/gl4-integrator.ts"],"names":[],"mappings":"AA4BA;;;;;;GAMG;AACH,OAAO,EAAE,KAAK,EAAE,KAAK,EAAE,KAAK,EAAE,MAAM,iCAAiC,CAAC;AAEtE;;;;;;;;;;;GAWG;AACH,MAAM,WAAW,QAAQ;IACvB,+BAA+B;IAC/B,QAAQ,CAAC,CAAC,EAAE,SAAS,MAAM,EAAE,CAAC;IAC9B,yCAAyC;IACzC,QAAQ,CAAC,CAAC,EAAE,SAAS,MAAM,EAAE,CAAC;CAC/B;AAED;;;;;;GAMG;AACH,MAAM,WAAW,WAAW;IAC1B,QAAQ,CAAC,GAAG,EAAE,MAAM,CAAC;IACrB,QAAQ,CAAC,CAAC,EAAE,SAAS,MAAM,EAAE,CAAC;IAC9B,QAAQ,CAAC,CAAC,EAAE,SAAS,MAAM,EAAE,CAAC;IAC9B,QAAQ,CAAC,CAAC,CAAC,EAAE,SAAS,MAAM,EAAE,CAAC;CAChC;AAED;;;;GAIG;AACH,MAAM,WAAW,UAAU;IACzB;qCACiC;IACjC,QAAQ,CAAC,KAAK,EAAE,MAAM,CAAC;IACvB,+CAA+C;IAC/C,QAAQ,CAAC,MAAM,EAAE,MAAM,CAAC;IACxB;gFAC4E;IAC5E,QAAQ,CAAC,UAAU,EAAE,CAAC,CAAC,EAAE,SAAS,MAAM,EAAE,KAAK,YAAY,CAAC;IAC5D;;;;;OAKG;IACH,QAAQ,CAAC,WAAW,EAAE,CAAC,CAAC,EAAE,SAAS,MAAM,EAAE,KAAK,YAAY,CAAC;IAC7D,+DAA+D;IAC/D,QAAQ,CAAC,SAAS,CAAC,EAAE,MAAM,CAAC;IAC5B,mEAAmE;IACnE,QAAQ,CAAC,aAAa,CAAC,EAAE,MAAM,CAAC;IAChC,+FAA+F;IAC/F,QAAQ,CAAC,IAAI,CAAC,EAAE,MAAM,CAAC;IACvB;;qCAEiC;IACjC,QAAQ,CAAC,eAAe,CAAC,EAAE,MAAM,CAAC;IAClC;;;;;;;;;;;;;;OAcG;IACH,QAAQ,CAAC,MAAM,CAAC,EAAE,CAAC,KAAK,EAAE;QAAE,IAAI,EAAE,MAAM,CAAC;QAAC,UAAU,EAAE,MAAM,CAAC;QAAC,QAAQ,EAAE,MAAM,CAAA;KAAE,KAAK,IAAI,CAAC;CAC3F;AAED;;;;;;;GAOG;AACH,UAAU,gBAAgB;IACxB,QAAQ,CAAC,MAAM,EAAE,SAAS,CAAC,SAAS,MAAM,EAAE,EAAE,SAAS,MAAM,EAAE,CAAC,CAAC;IACjE,QAAQ,CAAC,MAAM,EAAE,SAAS,CAAC,SAAS,MAAM,EAAE,EAAE,SAAS,MAAM,EAAE,CAAC,CAAC;IACjE,QAAQ,CAAC,OAAO,EAAE,SAAS,CAAC,SAAS,MAAM,EAAE,EAAE,SAAS,MAAM,EAAE,CAAC,CAAC;IAClE,QAAQ,CAAC,OAAO,EAAE,SAAS,CAAC,SAAS,MAAM,EAAE,EAAE,SAAS,MAAM,EAAE,CAAC,CAAC;IAClE,QAAQ,CAAC,UAAU,EAAE,MAAM,CAAC;CAC7B;AA4FD;;;;;;;;;GASG;AACH,wBAAgB,aAAa,CAC3B,KAAK,EAAE,QAAQ,EACf,CAAC,EAAE,MAAM,EACT,UAAU,EAAE,CAAC,CAAC,EAAE,SAAS,MAAM,EAAE,KAAK,YAAY,EAClD,WAAW,EAAE,CAAC,CAAC,EAAE,SAAS,MAAM,EAAE,KAAK,YAAY,EACnD,IAAI,EAAE;IAAE,SAAS,EAAE,MAAM,CAAC;IAAC,aAAa,EAAE,MAAM,CAAA;CAAE,GACjD,gBAAgB,CAYlB;AAMD;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GA4DG;AACH,wBAAgB,oBAAoB,CAClC,YAAY,EAAE,QAAQ,EACtB,OAAO,EAAE,UAAU,GAClB,SAAS,WAAW,EAAE,CAmGxB"}
+{"version":3,"file":"gl4-integrator.d.ts","sourceRoot":"","sources":["../../src/numerical/gl4-integrator.ts"],"names":[],"mappings":"AA4BA;;;;;;GAMG;AACH,OAAO,EAAE,KAAK,EAAE,KAAK,EAAE,KAAK,EAAE,MAAM,iCAAiC,CAAC;AAEtE;;;;;;;;;;;GAWG;AACH,MAAM,WAAW,QAAQ;IACvB,+BAA+B;IAC/B,QAAQ,CAAC,CAAC,EAAE,SAAS,MAAM,EAAE,CAAC;IAC9B,yCAAyC;IACzC,QAAQ,CAAC,CAAC,EAAE,SAAS,MAAM,EAAE,CAAC;CAC/B;AAED;;;;;;;GAOG;AACH,MAAM,WAAW,WAAW;IAC1B,QAAQ,CAAC,GAAG,EAAE,MAAM,CAAC;IACrB,QAAQ,CAAC,CAAC,EAAE,SAAS,MAAM,EAAE,CAAC;IAC9B,QAAQ,CAAC,CAAC,EAAE,SAAS,MAAM,EAAE,CAAC;IAC9B,QAAQ,CAAC,CAAC,CAAC,EAAE,SAAS,MAAM,EAAE,CAAC;CAChC;AAED;;;;GAIG;AACH,MAAM,WAAW,UAAU;IACzB;qCACiC;IACjC,QAAQ,CAAC,KAAK,EAAE,MAAM,CAAC;IACvB,+CAA+C;IAC/C,QAAQ,CAAC,MAAM,EAAE,MAAM,CAAC;IACxB;gFAC4E;IAC5E,QAAQ,CAAC,UAAU,EAAE,CAAC,CAAC,EAAE,SAAS,MAAM,EAAE,KAAK,YAAY,CAAC;IAC5D;;;;;OAKG;IACH,QAAQ,CAAC,WAAW,EAAE,CAAC,CAAC,EAAE,SAAS,MAAM,EAAE,KAAK,YAAY,CAAC;IAC7D,+DAA+D;IAC/D,QAAQ,CAAC,SAAS,CAAC,EAAE,MAAM,CAAC;IAC5B,mEAAmE;IACnE,QAAQ,CAAC,aAAa,CAAC,EAAE,MAAM,CAAC;IAChC,+FAA+F;IAC/F,QAAQ,CAAC,IAAI,CAAC,EAAE,MAAM,CAAC;IACvB;;qCAEiC;IACjC,QAAQ,CAAC,eAAe,CAAC,EAAE,MAAM,CAAC;IAClC;;;;;;;;;;;;;;OAcG;IACH,QAAQ,CAAC,MAAM,CAAC,EAAE,CAAC,KAAK,EAAE;QAAE,IAAI,EAAE,MAAM,CAAC;QAAC,UAAU,EAAE,MAAM,CAAC;QAAC,QAAQ,EAAE,MAAM,CAAA;KAAE,KAAK,IAAI,CAAC;CAC3F;AAgDD;;;;;GAKG;AACH,MAAM,WAAW,aAAa;IAC5B,QAAQ,CAAC,CAAC,EAAE,SAAS,MAAM,EAAE,CAAC;IAC9B,QAAQ,CAAC,CAAC,EAAE,SAAS,MAAM,EAAE,CAAC;IAC9B,QAAQ,CAAC,UAAU,EAAE,MAAM,CAAC;CAC7B;AAED;;;;;;GAMG;AACH,wBAAgB,OAAO,CACrB,KAAK,EAAE,QAAQ,EACf,CAAC,EAAE,MAAM,EACT,UAAU,EAAE,CAAC,CAAC,EAAE,SAAS,MAAM,EAAE,KAAK,YAAY,EAClD,WAAW,EAAE,CAAC,CAAC,EAAE,SAAS,MAAM,EAAE,KAAK,YAAY,EACnD,IAAI,EAAE;IAAE,SAAS,EAAE,MAAM,CAAC;IAAC,aAAa,EAAE,MAAM,CAAA;CAAE,GACjD,aAAa,CA6Bf;AAMD;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GAmEG;AACH,wBAAgB,oBAAoB,CAClC,YAAY,EAAE,QAAQ,EACtB,OAAO,EAAE,UAAU,GAClB,SAAS,WAAW,EAAE,CA0GxB"}

package/dist/numerical/gl4-integrator.js

@@ -17,13 +17,13 @@
  * p_μ = g_μν dx^ν/dτ — not (x, v). This is what makes the flow symplectic
  * on T*M.
  *
- * This module holds the types, the Butcher constants, the implicit Picard
- * stage solver (`solveGL4Stage`, internal), and the integrator entry point
- * `integrateGeodesicGL4`.
+ * This module holds the types, the Butcher constants, the one-step driver
+ * (`gl4Step`, internal; the implicit Picard stages are MathTS
+ * `gaussLegendre4`'s), and the integrator entry point `integrateGeodesicGL4`.
  *
  * @module numerical/gl4-integrator
  */
-import { gaussLegendre4 } from '@danielsimonjr/mathts-functions';
+import { gaussLegendre4, GL4ConvergenceError as MathTsGL4ConvergenceError } from '@danielsimonjr/mathts-functions';
 import { GL4ConvergenceError, NumericalBackendError } from './errors.js';
 /**
  * Gauss–Legendre order-4 Butcher tableau from MathTS. `GL4_A` is the 2×2
@@ -41,14 +41,15 @@ export { GL4_A, GL4_B, GL4_C } from '@danielsimonjr/mathts-functions';
  *
  * `dg[μ*dim² + ν*dim + ρ] = ∂_μ g^{νρ}`. A coefficient that is exactly 0 is
  * skipped, so `0 * NaN` does not poison a component the term does not enter.
+ * `xBuf` is a reused coordinate buffer of length `dim`; it is the array the
+ * metric closures receive, so it is a plain `number[]` and needs no cast.
  */
 function geodesicDeriv(y, gInverseFn, dgInverseFn, xBuf) {
     const dim = xBuf.length;
     for (let i = 0; i < dim; i++)
         xBuf[i] = y[i];
-    const coords = xBuf;
-    const gInv = gInverseFn(coords);
-    const dgInv = dgInverseFn(coords);
+    const gInv = gInverseFn(xBuf);
+    const dgInv = dgInverseFn(xBuf);
     const out = new Array(dim * 2);
     for (let mu = 0; mu < dim; mu++) {
         let dx = 0;
@@ -76,23 +77,29 @@ function geodesicDeriv(y, gInverseFn, dgInverseFn, xBuf) {
     }
     return out;
 }
-/** One GL4 step of the geodesic Hamiltonian, via MathTS `gaussLegendre4`. */
-function advanceGl4(state, h, gInverseFn, dgInverseFn, opts) {
+/**
+ * One GL4 step of the geodesic Hamiltonian, via MathTS `gaussLegendre4`.
+ * Throws `GL4ConvergenceError` with a message matching
+ * `/Picard iteration did not converge/` when the Picard cap is exhausted;
+ * any other error from the metric closures propagates unchanged.
+ * @internal
+ */
+export function gl4Step(state, h, gInverseFn, dgInverseFn, opts) {
     const dim = state.x.length;
     const y0 = new Array(dim * 2);
     for (let i = 0; i < dim; i++) {
         y0[i] = state.x[i];
         y0[dim + i] = state.p[i];
     }
-    const xBuf = new Float64Array(dim);
+    const xBuf = new Array(dim).fill(0);
     let solved;
     try {
         solved = gaussLegendre4((_t, y) => geodesicDeriv(y, gInverseFn, dgInverseFn, xBuf), y0, [0, h], { steps: 1, picardTol: opts.picardTol, picardMaxIter: opts.picardMaxIter });
     }
     catch (err) {
-        if (err instanceof Error && /Picard iteration did not converge/i.test(err.message)) {
+        // MathTS throws its own typed GL4ConvergenceError for a Picard step that does not converge.
+        if (err instanceof MathTsGL4ConvergenceError)
             throw new GL4ConvergenceError(err.message);
-        }
         throw err;
     }
     const y1 = solved.y[solved.y.length - 1];
@@ -102,29 +109,6 @@ function advanceGl4(state, h, gInverseFn, dgInverseFn, opts) {
     }
     return { x: y1.slice(0, dim), p: y1.slice(dim), iterations };
 }
-/**
- * One GL4 step, reported in the stage-result shape the stage tests read.
- *
- * The step itself is MathTS `gaussLegendre4`. In flat space `p` is constant,
- * so both stage momenta equal that advanced `p`. Throws `GL4ConvergenceError`
- * with a message matching `/Picard iteration did not converge/` when the
- * Picard cap is exhausted.
- *
- * @internal
- */
-export function solveGL4Stage(state, h, gInverseFn, dgInverseFn, opts) {
-    const step = advanceGl4(state, h, gInverseFn, dgInverseFn, opts);
-    const dim = step.p.length;
-    const zeros = new Array(dim).fill(0);
-    // In flat space p is constant, so both stage momenta equal the advanced p.
-    return {
-        stageX: [step.x.slice(), step.x.slice()],
-        stageP: [step.p.slice(), step.p.slice()],
-        stageDx: [zeros, zeros.slice()],
-        stageDp: [zeros.slice(), zeros.slice()],
-        iterations: step.iterations,
-    };
-}
 // ---------------------------------------------------------------------------
 // Integrator entry-point (Task 3, Phase 1a-iii)
 // ---------------------------------------------------------------------------
@@ -173,10 +157,11 @@ export function solveGL4Stage(state, h, gInverseFn, dgInverseFn, opts) {
  * @param initialState — canonical (x, p) at τ = 0. Units follow the
  *   `gInverseFn` convention (see above).
  * @param options — see {@link GL4Options}:
- *   - `steps` — integer step count (dimensionless).
- *   - `tauMax` — affine-parameter extent (seconds in canonical SI).
- *   - `gInverseFn(x)[μ][ν]` — inverse metric g^{μν}(x).
- *   - `dgInverseFn(x)[λ][μ][ν]` — ∂_λ g^{μν}(x).
+ *   - `steps` — positive integer step count (dimensionless).
+ *   - `tauMax` — positive finite affine-parameter extent (seconds in
+ *     canonical SI). The integrator runs forward from τ = 0 only.
+ *   - `gInverseFn(x)[μ*dim + ν]` — inverse metric g^{μν}(x), row-major flat.
+ *   - `dgInverseFn(x)[λ*dim² + μ*dim + ν]` — ∂_λ g^{μν}(x), flat.
  *   - `picardTol` — convergence tolerance (dimensionless, default 1e-12).
  *   - `picardMaxIter` — fixed-point iteration cap (dimensionless, default 50).
  *   - `hMin` — step-halving floor (same units as `tauMax / steps`).
@@ -184,21 +169,31 @@ export function solveGL4Stage(state, h, gInverseFn, dgInverseFn, opts) {
  * @returns `steps + 1` snapshots: index `0` is the initial state, index `n`
  *   is the state after `n` steps (τ = n · h). Each snapshot carries `tau`,
  *   `x`, `p` (and optional `v` = g^{μν} p_ν) in the units chosen above.
- * @throws NumericalBackendError if `initialState.x[1] < domainMinRadius`.
+ * @throws NumericalBackendError if `steps` is not a positive integer, if
+ *   `tauMax` is not a positive finite number, or if
+ *   `initialState.x[1] < domainMinRadius`. The two option checks run before
+ *   any metric evaluation.
  * @throws GL4ConvergenceError if Picard fails even after step-halving to h_min.
+ *   Any other error thrown by `gInverseFn` or `dgInverseFn` (a domain
+ *   crossing the caller signals by throwing, or a programming error) is
+ *   rethrown unchanged; it is not halved away.
  *
  * @public
  */
 export function integrateGeodesicGL4(initialState, options) {
     const { steps, tauMax, gInverseFn, dgInverseFn, picardTol = 1e-12, picardMaxIter = 50, hMin, domainMinRadius, onStep, } = options;
+    if (!Number.isInteger(steps) || steps <= 0) {
+        throw new NumericalBackendError(`GL4 integrator: steps must be a positive integer, got ${steps}`);
+    }
+    if (!Number.isFinite(tauMax) || tauMax <= 0) {
+        throw new NumericalBackendError(`GL4 integrator: tauMax must be a positive finite number, got ${tauMax}`);
+    }
     if (domainMinRadius !== undefined && initialState.x[1] < domainMinRadius) {
         throw new NumericalBackendError(`GL4 integrator: initial r=${initialState.x[1]} < domainMinRadius=${domainMinRadius} (domain violation)`);
     }
     const h = tauMax / steps;
     const hFloor = hMin ?? h * 1e-9;
-    const snapshots = [
-        { tau: 0, x: initialState.x.slice(), p: initialState.p.slice() },
-    ];
+    const snapshots = [{ tau: 0, x: initialState.x.slice(), p: initialState.p.slice() }];
     let x = initialState.x.slice();
     let p = initialState.p.slice();
     const stateDim = x.length;
@@ -225,14 +220,18 @@ export function integrateGeodesicGL4(initialState, options) {
             let subH = trialH;
             while (subH >= hFloor) {
                 try {
-                    advanced = advanceGl4({ x, p }, subH, gInverseFn, dgInverseFn, {
+                    advanced = gl4Step({ x, p }, subH, gInverseFn, dgInverseFn, {
                         picardTol,
                         picardMaxIter,
                     });
                     stepSucceeded = true;
                     break;
                 }
-                catch {
+                catch (err) {
+                    // Only a Picard failure is a reason to halve. A throw from the metric
+                    // closures (a domain crossing the caller signals, or a bug) is that error.
+                    if (!(err instanceof GL4ConvergenceError))
+                        throw err;
                     subH /= 2;
                     halvings++;
                 }

package/dist/numerical/gl4-integrator.js.map

@@ -1 +1 @@
-{"version":3,"file":"gl4-integrator.js","sourceRoot":"","sources":["../../src/numerical/gl4-integrator.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;;;;;;GAwBG;AACH,OAAO,EAAE,cAAc,EAAE,MAAM,iCAAiC,CAAC;AACjE,OAAO,EAAE,mBAAmB,EAAE,qBAAqB,EAAE,MAAM,aAAa,CAAC;AAEzE;;;;;;GAMG;AACH,OAAO,EAAE,KAAK,EAAE,KAAK,EAAE,KAAK,EAAE,MAAM,iCAAiC,CAAC;AAoGtE;;;;;;;;GAQG;AACH,SAAS,aAAa,CACpB,CAAoB,EACpB,UAAkD,EAClD,WAAmD,EACnD,IAAkB;IAElB,MAAM,GAAG,GAAG,IAAI,CAAC,MAAM,CAAC;IACxB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,GAAG,EAAE,CAAC,EAAE;QAAE,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAE,CAAC;IAC9C,MAAM,MAAM,GAAG,IAAoC,CAAC;IACpD,MAAM,IAAI,GAAG,UAAU,CAAC,MAAM,CAAC,CAAC;IAChC,MAAM,KAAK,GAAG,WAAW,CAAC,MAAM,CAAC,CAAC;IAClC,MAAM,GAAG,GAAG,IAAI,KAAK,CAAS,GAAG,GAAG,CAAC,CAAC,CAAC;IACvC,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,GAAG,EAAE,EAAE,EAAE,EAAE,CAAC;QAChC,IAAI,EAAE,GAAG,CAAC,CAAC;QACX,IAAI,EAAE,GAAG,CAAC,CAAC;QACX,MAAM,IAAI,GAAG,EAAE,GAAG,GAAG,CAAC;QACtB,MAAM,MAAM,GAAG,IAAI,GAAG,GAAG,CAAC;QAC1B,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,GAAG,EAAE,EAAE,EAAE,EAAE,CAAC;YAChC,MAAM,GAAG,GAAG,CAAC,CAAC,GAAG,GAAG,EAAE,CAAE,CAAC;YACzB,MAAM,CAAC,GAAG,IAAI,CAAC,IAAI,GAAG,EAAE,CAAE,CAAC;YAC3B,wFAAwF;YACxF,IAAI,CAAC,KAAK,CAAC;gBAAE,EAAE,IAAI,CAAC,GAAG,GAAG,CAAC;YAC3B,IAAI,IAAI,GAAG,CAAC,CAAC;YACb,MAAM,MAAM,GAAG,MAAM,GAAG,EAAE,GAAG,GAAG,CAAC;YACjC,KAAK,IAAI,GAAG,GAAG,CAAC,EAAE,GAAG,GAAG,GAAG,EAAE,GAAG,EAAE,EAAE,CAAC;gBACnC,MAAM,EAAE,GAAG,KAAK,CAAC,MAAM,GAAG,GAAG,CAAE,CAAC;gBAChC,IAAI,EAAE,KAAK,CAAC;oBAAE,IAAI,IAAI,EAAE,GAAG,CAAC,CAAC,GAAG,GAAG,GAAG,CAAE,CAAC;YAC3C,CAAC;YACD,IAAI,IAAI,KAAK,CAAC;gBAAE,EAAE,IAAI,IAAI,GAAG,GAAG,CAAC;QACnC,CAAC;QACD,GAAG,CAAC,EAAE,CAAC,GAAG,EAAE,CAAC;QACb,GAAG,CAAC,GAAG,GAAG,EAAE,CAAC,GAAG,EAAE,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,GAAG,EAAE,CAAC;IAC3C,CAAC;IACD,OAAO,GAAG,CAAC;AACb,CAAC;AAQD,6EAA6E;AAC7E,SAAS,UAAU,CACjB,KAAe,EACf,CAAS,EACT,UAAkD,EAClD,WAAmD,EACnD,IAAkD;IAElD,MAAM,GAAG,GAAG,KAAK,CAAC,CAAC,CAAC,MAAM,CAAC;IAC3B,MAAM,EAAE,GAAG,IAAI,KAAK,CAAS,GAAG,GAAG,CAAC,CAAC,CAAC;IACtC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,GAAG,EAAE,CAAC,EAAE,EAAE,CAAC;QAC7B,EAAE,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC,CAAE,CAAC;QACpB,EAAE,CAAC,GAAG,GAAG,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC,CAAE,CAAC;IAC5B,CAAC;IACD,MAAM,IAAI,GAAG,IAAI,YAAY,CAAC,GAAG,CAAC,CAAC;IACnC,IAAI,MAAyC,CAAC;IAC9C,IAAI,CAAC;QACH,MAAM,GAAG,cAAc,CACrB,CAAC,EAAE,EAAE,CAAC,EAAE,EAAE,CAAC,aAAa,CAAC,CAAC,EAAE,UAAU,EAAE,WAAW,EAAE,IAAI,CAAC,EAC1D,EAAE,EACF,CAAC,CAAC,EAAE,CAAC,CAAC,EACN,EAAE,KAAK,EAAE,CAAC,EAAE,SAAS,EAAE,IAAI,CAAC,SAAS,EAAE,aAAa,EAAE,IAAI,CAAC,aAAa,EAAE,CAC3E,CAAC;IACJ,CAAC;IAAC,OAAO,GAAG,EAAE,CAAC;QACb,IAAI,GAAG,YAAY,KAAK,IAAI,oCAAoC,CAAC,IAAI,CAAC,GAAG,CAAC,OAAO,CAAC,EAAE,CAAC;YACnF,MAAM,IAAI,mBAAmB,CAAC,GAAG,CAAC,OAAO,CAAC,CAAC;QAC7C,CAAC;QACD,MAAM,GAAG,CAAC;IACZ,CAAC;IACD,MAAM,EAAE,GAAG,MAAM,CAAC,CAAC,CAAC,MAAM,CAAC,CAAC,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC;IACzC,MAAM,UAAU,GAAG,MAAM,CAAC,UAAU,CAAC,CAAC,CAAC,CAAC;IACxC,IAAI,EAAE,KAAK,SAAS,IAAI,UAAU,KAAK,SAAS,EAAE,CAAC;QACjD,MAAM,IAAI,mBAAmB,CAC3B,wCAAwC,IAAI,CAAC,aAAa,yCAAyC,IAAI,CAAC,SAAS,GAAG,CACrH,CAAC;IACJ,CAAC;IACD,OAAO,EAAE,CAAC,EAAE,EAAE,CAAC,KAAK,CAAC,CAAC,EAAE,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,KAAK,CAAC,GAAG,CAAC,EAAE,UAAU,EAAE,CAAC;AAC/D,CAAC;AAED;;;;;;;;;GASG;AACH,MAAM,UAAU,aAAa,CAC3B,KAAe,EACf,CAAS,EACT,UAAkD,EAClD,WAAmD,EACnD,IAAkD;IAElD,MAAM,IAAI,GAAG,UAAU,CAAC,KAAK,EAAE,CAAC,EAAE,UAAU,EAAE,WAAW,EAAE,IAAI,CAAC,CAAC;IACjE,MAAM,GAAG,GAAG,IAAI,CAAC,CAAC,CAAC,MAAM,CAAC;IAC1B,MAAM,KAAK,GAAG,IAAI,KAAK,CAAS,GAAG,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;IAC7C,2EAA2E;IAC3E,OAAO;QACL,MAAM,EAAE,CAAC,IAAI,CAAC,CAAC,CAAC,KAAK,EAAE,EAAE,IAAI,CAAC,CAAC,CAAC,KAAK,EAAE,CAAC;QACxC,MAAM,EAAE,CAAC,IAAI,CAAC,CAAC,CAAC,KAAK,EAAE,EAAE,IAAI,CAAC,CAAC,CAAC,KAAK,EAAE,CAAC;QACxC,OAAO,EAAE,CAAC,KAAK,EAAE,KAAK,CAAC,KAAK,EAAE,CAAC;QAC/B,OAAO,EAAE,CAAC,KAAK,CAAC,KAAK,EAAE,EAAE,KAAK,CAAC,KAAK,EAAE,CAAC;QACvC,UAAU,EAAE,IAAI,CAAC,UAAU;KAC5B,CAAC;AACJ,CAAC;AAED,8EAA8E;AAC9E,gDAAgD;AAChD,8EAA8E;AAE9E;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GA4DG;AACH,MAAM,UAAU,oBAAoB,CAClC,YAAsB,EACtB,OAAmB;IAEnB,MAAM,EACJ,KAAK,EACL,MAAM,EACN,UAAU,EACV,WAAW,EACX,SAAS,GAAG,KAAK,EACjB,aAAa,GAAG,EAAE,EAClB,IAAI,EACJ,eAAe,EACf,MAAM,GACP,GAAG,OAAO,CAAC;IAEZ,IAAI,eAAe,KAAK,SAAS,IAAI,YAAY,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,eAAe,EAAE,CAAC;QACzE,MAAM,IAAI,qBAAqB,CAC7B,6BAA6B,YAAY,CAAC,CAAC,CAAC,CAAC,CAAC,sBAAsB,eAAe,qBAAqB,CACzG,CAAC;IACJ,CAAC;IAED,MAAM,CAAC,GAAG,MAAM,GAAG,KAAK,CAAC;IACzB,MAAM,MAAM,GAAG,IAAI,IAAI,CAAC,GAAG,IAAI,CAAC;IAChC,MAAM,SAAS,GAAkB;QAC/B,EAAE,GAAG,EAAE,CAAC,EAAE,CAAC,EAAE,YAAY,CAAC,CAAC,CAAC,KAAK,EAAc,EAAE,CAAC,EAAE,YAAY,CAAC,CAAC,CAAC,KAAK,EAAc,EAAE;KACzF,CAAC;IACF,IAAI,CAAC,GAAG,YAAY,CAAC,CAAC,CAAC,KAAK,EAAc,CAAC;IAC3C,IAAI,CAAC,GAAG,YAAY,CAAC,CAAC,CAAC,KAAK,EAAc,CAAC;IAE3C,MAAM,QAAQ,GAAG,CAAC,CAAC,MAAM,CAAC;IAC1B,IAAI,IAAI,GAAG,IAAI,KAAK,CAAS,QAAQ,CAAC,CAAC;IACvC,IAAI,IAAI,GAAG,IAAI,KAAK,CAAS,QAAQ,CAAC,CAAC;IAEvC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,KAAK,EAAE,CAAC,EAAE,EAAE,CAAC;QAC/B,8DAA8D;QAC9D,uEAAuE;QACvE,yEAAyE;QACzE,uEAAuE;QACvE,wEAAwE;QACxE,sEAAsE;QACtE,4DAA4D;QAC5D,IAAI,SAAS,GAAG,CAAC,CAAC;QAClB,IAAI,KAAK,GAAG,CAAC,CAAC;QACd,IAAI,QAAQ,GAAG,CAAC,CAAC;QACjB,IAAI,cAAc,GAAG,CAAC,CAAC;QACvB,+EAA+E;QAC/E,MAAM,SAAS,GAAG,CAAC,GAAG,KAAK,CAAC;QAC5B,OAAO,SAAS,GAAG,SAAS,EAAE,CAAC;YAC7B,MAAM,MAAM,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,EAAE,SAAS,CAAC,CAAC;YAC1C,IAAI,QAAgC,CAAC;YACrC,IAAI,aAAa,GAAG,KAAK,CAAC;YAC1B,IAAI,IAAI,GAAG,MAAM,CAAC;YAClB,OAAO,IAAI,IAAI,MAAM,EAAE,CAAC;gBACtB,IAAI,CAAC;oBACH,QAAQ,GAAG,UAAU,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,EAAE,IAAI,EAAE,UAAU,EAAE,WAAW,EAAE;wBAC7D,SAAS;wBACT,aAAa;qBACd,CAAC,CAAC;oBACH,aAAa,GAAG,IAAI,CAAC;oBACrB,MAAM;gBACR,CAAC;gBAAC,MAAM,CAAC;oBACP,IAAI,IAAI,CAAC,CAAC;oBACV,QAAQ,EAAE,CAAC;gBACb,CAAC;YACH,CAAC;YACD,IAAI,CAAC,aAAa,IAAI,QAAQ,KAAK,SAAS,EAAE,CAAC;gBAC7C,MAAM,IAAI,mBAAmB,CAC3B,mEAAmE,MAAM,UAAU,CAAC,8CAA8C,CACnI,CAAC;YACJ,CAAC;YACD,IAAI,QAAQ,KAAK,CAAC,EAAE,CAAC;gBACnB,IAAI,CAAC,CAAC,CAAC,GAAG,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAE,CAAC;gBAAC,IAAI,CAAC,CAAC,CAAC,GAAG,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAE,CAAC;gBACnD,IAAI,CAAC,CAAC,CAAC,GAAG,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAE,CAAC;gBAAC,IAAI,CAAC,CAAC,CAAC,GAAG,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAE,CAAC;gBACnD,IAAI,CAAC,CAAC,CAAC,GAAG,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAE,CAAC;gBAAC,IAAI,CAAC,CAAC,CAAC,GAAG,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAE,CAAC;gBACnD,IAAI,CAAC,CAAC,CAAC,GAAG,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAE,CAAC;gBAAC,IAAI,CAAC,CAAC,CAAC,GAAG,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAE,CAAC;YACrD,CAAC;iBAAM,CAAC;gBACN,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,QAAQ,EAAE,EAAE,EAAE,EAAE,CAAC;oBACrC,IAAI,CAAC,EAAE,CAAC,GAAG,QAAQ,CAAC,CAAC,CAAC,EAAE,CAAE,CAAC;oBAC3B,IAAI,CAAC,EAAE,CAAC,GAAG,QAAQ,CAAC,CAAC,CAAC,EAAE,CAAE,CAAC;gBAC7B,CAAC;YACH,CAAC;YAED,MAAM,IAAI,GAAG,CAAC,CAAC;YACf,CAAC,GAAG,IAAI,CAAC;YACT,IAAI,GAAG,IAAI,CAAC;YAEZ,MAAM,IAAI,GAAG,CAAC,CAAC;YACf,CAAC,GAAG,IAAI,CAAC;YACT,IAAI,GAAG,IAAI,CAAC;YAEZ,SAAS,IAAI,IAAI,CAAC;YAClB,KAAK,GAAG,IAAI,CAAC,CAAC,2DAA2D;YACzE,cAAc,GAAG,QAAQ,CAAC,UAAU,CAAC;QACvC,CAAC;QACD,SAAS,CAAC,IAAI,CAAC,EAAE,GAAG,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,KAAK,EAAE,EAAE,CAAC,EAAE,CAAC,CAAC,KAAK,EAAE,EAAE,CAAC,CAAC;QACjE,IAAI,MAAM,KAAK,SAAS,EAAE,CAAC;YACzB,MAAM,CAAC,EAAE,IAAI,EAAE,CAAC,EAAE,UAAU,EAAE,cAAc,EAAE,QAAQ,EAAE,CAAC,CAAC;QAC5D,CAAC;IACH,CAAC;IAED,OAAO,SAAS,CAAC;AACnB,CAAC"}
+{"version":3,"file":"gl4-integrator.js","sourceRoot":"","sources":["../../src/numerical/gl4-integrator.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;;;;;;GAwBG;AACH,OAAO,EAAE,cAAc,EAAE,mBAAmB,IAAI,yBAAyB,EAAE,MAAM,iCAAiC,CAAC;AACnH,OAAO,EAAE,mBAAmB,EAAE,qBAAqB,EAAE,MAAM,aAAa,CAAC;AAEzE;;;;;;GAMG;AACH,OAAO,EAAE,KAAK,EAAE,KAAK,EAAE,KAAK,EAAE,MAAM,iCAAiC,CAAC;AAqFtE;;;;;;;;;;GAUG;AACH,SAAS,aAAa,CACpB,CAAoB,EACpB,UAAkD,EAClD,WAAmD,EACnD,IAAc;IAEd,MAAM,GAAG,GAAG,IAAI,CAAC,MAAM,CAAC;IACxB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,GAAG,EAAE,CAAC,EAAE;QAAE,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAE,CAAC;IAC9C,MAAM,IAAI,GAAG,UAAU,CAAC,IAAI,CAAC,CAAC;IAC9B,MAAM,KAAK,GAAG,WAAW,CAAC,IAAI,CAAC,CAAC;IAChC,MAAM,GAAG,GAAG,IAAI,KAAK,CAAS,GAAG,GAAG,CAAC,CAAC,CAAC;IACvC,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,GAAG,EAAE,EAAE,EAAE,EAAE,CAAC;QAChC,IAAI,EAAE,GAAG,CAAC,CAAC;QACX,IAAI,EAAE,GAAG,CAAC,CAAC;QACX,MAAM,IAAI,GAAG,EAAE,GAAG,GAAG,CAAC;QACtB,MAAM,MAAM,GAAG,IAAI,GAAG,GAAG,CAAC;QAC1B,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,GAAG,EAAE,EAAE,EAAE,EAAE,CAAC;YAChC,MAAM,GAAG,GAAG,CAAC,CAAC,GAAG,GAAG,EAAE,CAAE,CAAC;YACzB,MAAM,CAAC,GAAG,IAAI,CAAC,IAAI,GAAG,EAAE,CAAE,CAAC;YAC3B,wFAAwF;YACxF,IAAI,CAAC,KAAK,CAAC;gBAAE,EAAE,IAAI,CAAC,GAAG,GAAG,CAAC;YAC3B,IAAI,IAAI,GAAG,CAAC,CAAC;YACb,MAAM,MAAM,GAAG,MAAM,GAAG,EAAE,GAAG,GAAG,CAAC;YACjC,KAAK,IAAI,GAAG,GAAG,CAAC,EAAE,GAAG,GAAG,GAAG,EAAE,GAAG,EAAE,EAAE,CAAC;gBACnC,MAAM,EAAE,GAAG,KAAK,CAAC,MAAM,GAAG,GAAG,CAAE,CAAC;gBAChC,IAAI,EAAE,KAAK,CAAC;oBAAE,IAAI,IAAI,EAAE,GAAG,CAAC,CAAC,GAAG,GAAG,GAAG,CAAE,CAAC;YAC3C,CAAC;YACD,IAAI,IAAI,KAAK,CAAC;gBAAE,EAAE,IAAI,IAAI,GAAG,GAAG,CAAC;QACnC,CAAC;QACD,GAAG,CAAC,EAAE,CAAC,GAAG,EAAE,CAAC;QACb,GAAG,CAAC,GAAG,GAAG,EAAE,CAAC,GAAG,EAAE,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,GAAG,EAAE,CAAC;IAC3C,CAAC;IACD,OAAO,GAAG,CAAC;AACb,CAAC;AAcD;;;;;;GAMG;AACH,MAAM,UAAU,OAAO,CACrB,KAAe,EACf,CAAS,EACT,UAAkD,EAClD,WAAmD,EACnD,IAAkD;IAElD,MAAM,GAAG,GAAG,KAAK,CAAC,CAAC,CAAC,MAAM,CAAC;IAC3B,MAAM,EAAE,GAAG,IAAI,KAAK,CAAS,GAAG,GAAG,CAAC,CAAC,CAAC;IACtC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,GAAG,EAAE,CAAC,EAAE,EAAE,CAAC;QAC7B,EAAE,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC,CAAE,CAAC;QACpB,EAAE,CAAC,GAAG,GAAG,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC,CAAE,CAAC;IAC5B,CAAC;IACD,MAAM,IAAI,GAAG,IAAI,KAAK,CAAS,GAAG,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;IAC5C,IAAI,MAAyC,CAAC;IAC9C,IAAI,CAAC;QACH,MAAM,GAAG,cAAc,CACrB,CAAC,EAAE,EAAE,CAAC,EAAE,EAAE,CAAC,aAAa,CAAC,CAAC,EAAE,UAAU,EAAE,WAAW,EAAE,IAAI,CAAC,EAC1D,EAAE,EACF,CAAC,CAAC,EAAE,CAAC,CAAC,EACN,EAAE,KAAK,EAAE,CAAC,EAAE,SAAS,EAAE,IAAI,CAAC,SAAS,EAAE,aAAa,EAAE,IAAI,CAAC,aAAa,EAAE,CAC3E,CAAC;IACJ,CAAC;IAAC,OAAO,GAAG,EAAE,CAAC;QACb,4FAA4F;QAC5F,IAAI,GAAG,YAAY,yBAAyB;YAAE,MAAM,IAAI,mBAAmB,CAAC,GAAG,CAAC,OAAO,CAAC,CAAC;QACzF,MAAM,GAAG,CAAC;IACZ,CAAC;IACD,MAAM,EAAE,GAAG,MAAM,CAAC,CAAC,CAAC,MAAM,CAAC,CAAC,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC;IACzC,MAAM,UAAU,GAAG,MAAM,CAAC,UAAU,CAAC,CAAC,CAAC,CAAC;IACxC,IAAI,EAAE,KAAK,SAAS,IAAI,UAAU,KAAK,SAAS,EAAE,CAAC;QACjD,MAAM,IAAI,mBAAmB,CAC3B,wCAAwC,IAAI,CAAC,aAAa,yCAAyC,IAAI,CAAC,SAAS,GAAG,CACrH,CAAC;IACJ,CAAC;IACD,OAAO,EAAE,CAAC,EAAE,EAAE,CAAC,KAAK,CAAC,CAAC,EAAE,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,KAAK,CAAC,GAAG,CAAC,EAAE,UAAU,EAAE,CAAC;AAC/D,CAAC;AAED,8EAA8E;AAC9E,gDAAgD;AAChD,8EAA8E;AAE9E;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GAmEG;AACH,MAAM,UAAU,oBAAoB,CAClC,YAAsB,EACtB,OAAmB;IAEnB,MAAM,EACJ,KAAK,EACL,MAAM,EACN,UAAU,EACV,WAAW,EACX,SAAS,GAAG,KAAK,EACjB,aAAa,GAAG,EAAE,EAClB,IAAI,EACJ,eAAe,EACf,MAAM,GACP,GAAG,OAAO,CAAC;IAEZ,IAAI,CAAC,MAAM,CAAC,SAAS,CAAC,KAAK,CAAC,IAAI,KAAK,IAAI,CAAC,EAAE,CAAC;QAC3C,MAAM,IAAI,qBAAqB,CAAC,yDAAyD,KAAK,EAAE,CAAC,CAAC;IACpG,CAAC;IACD,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,MAAM,CAAC,IAAI,MAAM,IAAI,CAAC,EAAE,CAAC;QAC5C,MAAM,IAAI,qBAAqB,CAAC,gEAAgE,MAAM,EAAE,CAAC,CAAC;IAC5G,CAAC;IACD,IAAI,eAAe,KAAK,SAAS,IAAI,YAAY,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,eAAe,EAAE,CAAC;QACzE,MAAM,IAAI,qBAAqB,CAC7B,6BAA6B,YAAY,CAAC,CAAC,CAAC,CAAC,CAAC,sBAAsB,eAAe,qBAAqB,CACzG,CAAC;IACJ,CAAC;IAED,MAAM,CAAC,GAAG,MAAM,GAAG,KAAK,CAAC;IACzB,MAAM,MAAM,GAAG,IAAI,IAAI,CAAC,GAAG,IAAI,CAAC;IAChC,MAAM,SAAS,GAAkB,CAAC,EAAE,GAAG,EAAE,CAAC,EAAE,CAAC,EAAE,YAAY,CAAC,CAAC,CAAC,KAAK,EAAE,EAAE,CAAC,EAAE,YAAY,CAAC,CAAC,CAAC,KAAK,EAAE,EAAE,CAAC,CAAC;IACpG,IAAI,CAAC,GAAG,YAAY,CAAC,CAAC,CAAC,KAAK,EAAE,CAAC;IAC/B,IAAI,CAAC,GAAG,YAAY,CAAC,CAAC,CAAC,KAAK,EAAE,CAAC;IAE/B,MAAM,QAAQ,GAAG,CAAC,CAAC,MAAM,CAAC;IAC1B,IAAI,IAAI,GAAG,IAAI,KAAK,CAAS,QAAQ,CAAC,CAAC;IACvC,IAAI,IAAI,GAAG,IAAI,KAAK,CAAS,QAAQ,CAAC,CAAC;IAEvC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,KAAK,EAAE,CAAC,EAAE,EAAE,CAAC;QAC/B,8DAA8D;QAC9D,uEAAuE;QACvE,yEAAyE;QACzE,uEAAuE;QACvE,wEAAwE;QACxE,sEAAsE;QACtE,4DAA4D;QAC5D,IAAI,SAAS,GAAG,CAAC,CAAC;QAClB,IAAI,KAAK,GAAG,CAAC,CAAC;QACd,IAAI,QAAQ,GAAG,CAAC,CAAC;QACjB,IAAI,cAAc,GAAG,CAAC,CAAC;QACvB,+EAA+E;QAC/E,MAAM,SAAS,GAAG,CAAC,GAAG,KAAK,CAAC;QAC5B,OAAO,SAAS,GAAG,SAAS,EAAE,CAAC;YAC7B,MAAM,MAAM,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,EAAE,SAAS,CAAC,CAAC;YAC1C,IAAI,QAAmC,CAAC;YACxC,IAAI,aAAa,GAAG,KAAK,CAAC;YAC1B,IAAI,IAAI,GAAG,MAAM,CAAC;YAClB,OAAO,IAAI,IAAI,MAAM,EAAE,CAAC;gBACtB,IAAI,CAAC;oBACH,QAAQ,GAAG,OAAO,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,EAAE,IAAI,EAAE,UAAU,EAAE,WAAW,EAAE;wBAC1D,SAAS;wBACT,aAAa;qBACd,CAAC,CAAC;oBACH,aAAa,GAAG,IAAI,CAAC;oBACrB,MAAM;gBACR,CAAC;gBAAC,OAAO,GAAG,EAAE,CAAC;oBACb,sEAAsE;oBACtE,2EAA2E;oBAC3E,IAAI,CAAC,CAAC,GAAG,YAAY,mBAAmB,CAAC;wBAAE,MAAM,GAAG,CAAC;oBACrD,IAAI,IAAI,CAAC,CAAC;oBACV,QAAQ,EAAE,CAAC;gBACb,CAAC;YACH,CAAC;YACD,IAAI,CAAC,aAAa,IAAI,QAAQ,KAAK,SAAS,EAAE,CAAC;gBAC7C,MAAM,IAAI,mBAAmB,CAC3B,mEAAmE,MAAM,UAAU,CAAC,8CAA8C,CACnI,CAAC;YACJ,CAAC;YACD,IAAI,QAAQ,KAAK,CAAC,EAAE,CAAC;gBACnB,IAAI,CAAC,CAAC,CAAC,GAAG,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAE,CAAC;gBAAC,IAAI,CAAC,CAAC,CAAC,GAAG,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAE,CAAC;gBACnD,IAAI,CAAC,CAAC,CAAC,GAAG,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAE,CAAC;gBAAC,IAAI,CAAC,CAAC,CAAC,GAAG,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAE,CAAC;gBACnD,IAAI,CAAC,CAAC,CAAC,GAAG,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAE,CAAC;gBAAC,IAAI,CAAC,CAAC,CAAC,GAAG,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAE,CAAC;gBACnD,IAAI,CAAC,CAAC,CAAC,GAAG,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAE,CAAC;gBAAC,IAAI,CAAC,CAAC,CAAC,GAAG,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAE,CAAC;YACrD,CAAC;iBAAM,CAAC;gBACN,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,QAAQ,EAAE,EAAE,EAAE,EAAE,CAAC;oBACrC,IAAI,CAAC,EAAE,CAAC,GAAG,QAAQ,CAAC,CAAC,CAAC,EAAE,CAAE,CAAC;oBAC3B,IAAI,CAAC,EAAE,CAAC,GAAG,QAAQ,CAAC,CAAC,CAAC,EAAE,CAAE,CAAC;gBAC7B,CAAC;YACH,CAAC;YAED,MAAM,IAAI,GAAG,CAAC,CAAC;YACf,CAAC,GAAG,IAAI,CAAC;YACT,IAAI,GAAG,IAAI,CAAC;YAEZ,MAAM,IAAI,GAAG,CAAC,CAAC;YACf,CAAC,GAAG,IAAI,CAAC;YACT,IAAI,GAAG,IAAI,CAAC;YAEZ,SAAS,IAAI,IAAI,CAAC;YAClB,KAAK,GAAG,IAAI,CAAC,CAAC,2DAA2D;YACzE,cAAc,GAAG,QAAQ,CAAC,UAAU,CAAC;QACvC,CAAC;QACD,SAAS,CAAC,IAAI,CAAC,EAAE,GAAG,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,KAAK,EAAE,EAAE,CAAC,EAAE,CAAC,CAAC,KAAK,EAAE,EAAE,CAAC,CAAC;QACjE,IAAI,MAAM,KAAK,SAAS,EAAE,CAAC;YACzB,MAAM,CAAC,EAAE,IAAI,EAAE,CAAC,EAAE,UAAU,EAAE,cAAc,EAAE,QAAQ,EAAE,CAAC,CAAC;QAC5D,CAAC;IACH,CAAC;IAED,OAAO,SAAS,CAAC;AACnB,CAAC","sourcesContent":["/**\n * Gauss-Legendre 4th-order (GL4) symplectic integrator — types + Butcher\n * tableau scaffold (v0.5.0 Task 1, Phase 1a-i).\n *\n * GL4 is a 2-stage implicit Runge-Kutta method built on the roots of the\n * shifted Legendre polynomial P₂, with order p = 4 and stage order s = 2.\n * It is symplectic for non-separable Hamiltonians — the property that\n * justifies its selection over Ruth-4 for the geodesic Hamiltonian\n *\n *   H(x, p) = ½ g^{μν}(x) p_μ p_ν,\n *\n * which is non-separable because g^{μν} depends on x. (v0.5.0 Decision #2,\n * post-adversarial reconciliation; Sanz-Serna 1988; Hairer/Lubich/Wanner\n * \"Geometric Numerical Integration\" §II.1.)\n *\n * The state is canonical (x, p) (Decision #3) — covariant momentum\n * p_μ = g_μν dx^ν/dτ — not (x, v). This is what makes the flow symplectic\n * on T*M.\n *\n * This module holds the types, the Butcher constants, the one-step driver\n * (`gl4Step`, internal; the implicit Picard stages are MathTS\n * `gaussLegendre4`'s), and the integrator entry point `integrateGeodesicGL4`.\n *\n * @module numerical/gl4-integrator\n */\nimport { gaussLegendre4, GL4ConvergenceError as MathTsGL4ConvergenceError } from '@danielsimonjr/mathts-functions';\nimport { GL4ConvergenceError, NumericalBackendError } from './errors.js';\n\n/**\n * Gauss–Legendre order-4 Butcher tableau from MathTS. `GL4_A` is the 2×2\n * stage matrix, `GL4_B` the weights (both ½), and `GL4_C` the two nodes.\n * `integrateGeodesicGL4` does not read these names; the tableau tests do.\n *\n * @internal\n */\nexport { GL4_A, GL4_B, GL4_C } from '@danielsimonjr/mathts-functions';\n\n/**\n * Canonical (x, p) phase-space state for the geodesic flow on T*M.\n *\n *   x^μ        — coordinate 4-vector\n *   p_μ        — covariant momentum, p_μ = g_μν dx^ν/dτ\n *\n * Decision #3 (v0.5.0): the canonical state is (x, p), not (x, v). The\n * symplectic 2-form ω = dp_μ ∧ dx^μ is preserved by the GL4 flow only on\n * this representation.\n *\n * @public\n */\nexport interface GL4State {\n  /** Coordinate 4-vector x^μ. */\n  readonly x: readonly number[];\n  /** Covariant momentum p_μ = g_μν v^ν. */\n  readonly p: readonly number[];\n}\n\n/**\n * Per-step snapshot recorded by the integrator. `v` (the contravariant\n * 4-velocity v^μ = g^{μν} p_ν) is declared optional and is never emitted:\n * no option asks for it, and nothing in `src/` reads it. It stays so a caller's\n * own snapshot type can carry one.\n *\n * @public\n */\nexport interface GL4Snapshot {\n  readonly tau: number;\n  readonly x: readonly number[];\n  readonly p: readonly number[];\n  readonly v?: readonly number[];\n}\n\n/**\n * Options for `integrateGeodesicGL4`.\n *\n * @public\n */\nexport interface GL4Options {\n  /** Number of integration steps (uniform-step baseline; may be subdivided\n   *  by adaptive step-halving). */\n  readonly steps: number;\n  /** Final proper time τ_max (initial τ = 0). */\n  readonly tauMax: number;\n  /** Inverse-metric closure (v0.9.0 O-1 flat layout):\n   *  `gInverseFn(x)[μ*dim + ν] = g^{μν}(x)`, row-major Float64Array(dim²). */\n  readonly gInverseFn: (x: readonly number[]) => Float64Array;\n  /**\n   * Partial derivatives of the inverse metric.\n   * Index order (v0.9.0 O-1 flat layout):\n   * `dgInverseFn(x)[lambda*dim² + mu*dim + nu] = ∂_lambda g^{mu nu}` at coords x.\n   * (I2: axis semantics pinned here to prevent silent transposition bugs.)\n   */\n  readonly dgInverseFn: (x: readonly number[]) => Float64Array;\n  /** Picard fixed-point tolerance (default chosen in Task 2). */\n  readonly picardTol?: number;\n  /** Picard fixed-point iteration cap (default chosen in Task 2). */\n  readonly picardMaxIter?: number;\n  /** Adaptive step-halving floor (I4). If step-halving reaches h_min, throws with diagnostic. */\n  readonly hMin?: number;\n  /** Minimum radial coordinate (or domain analog) — abort integration if\n   *  the trajectory crosses inside this radius (e.g., the Schwarzschild\n   *  event horizon at r = r_s). */\n  readonly domainMinRadius?: number;\n  /**\n   * v0.5.1 PD-4: opt-in per-step diagnostics callback. Fires once per\n   * successful step with the Picard iteration count consumed and whether\n   * adaptive step-halving had to subdivide (an \"exhaustion\" event from the\n   * caller's perspective: the original h failed Picard and was halved).\n   *\n   * Used by the gated `GL4_LONG=1` Mercury 100-orbit Picard-convergence\n   * test to measure the failure fraction across millions of steps without\n   * polluting the integrator's return shape for normal callers.\n   *\n   * - `iterations`: Picard iteration count actually consumed at the\n   *   successful step size (always 1..picardMaxIter).\n   * - `halvings`: number of step-halvings the step required before\n   *   succeeding (0 = first try; ≥1 means original h hit picardMaxIter).\n   */\n  readonly onStep?: (event: { step: number; iterations: number; halvings: number }) => void;\n}\n\n/**\n * Geodesic Hamiltonian derivative, packed as `y = [x..., p...]`.\n *\n *   dx^μ = g^{μν} p_ν\n *   dp_μ = −½ (∂_μ g^{νρ}) p_ν p_ρ\n *\n * `dg[μ*dim² + ν*dim + ρ] = ∂_μ g^{νρ}`. A coefficient that is exactly 0 is\n * skipped, so `0 * NaN` does not poison a component the term does not enter.\n * `xBuf` is a reused coordinate buffer of length `dim`; it is the array the\n * metric closures receive, so it is a plain `number[]` and needs no cast.\n */\nfunction geodesicDeriv(\n  y: readonly number[],\n  gInverseFn: (x: readonly number[]) => Float64Array,\n  dgInverseFn: (x: readonly number[]) => Float64Array,\n  xBuf: number[],\n): number[] {\n  const dim = xBuf.length;\n  for (let i = 0; i < dim; i++) xBuf[i] = y[i]!;\n  const gInv = gInverseFn(xBuf);\n  const dgInv = dgInverseFn(xBuf);\n  const out = new Array<number>(dim * 2);\n  for (let mu = 0; mu < dim; mu++) {\n    let dx = 0;\n    let dp = 0;\n    const gRow = mu * dim;\n    const dgBase = gRow * dim;\n    for (let nu = 0; nu < dim; nu++) {\n      const pNu = y[dim + nu]!;\n      const g = gInv[gRow + nu]!;\n      // Skip an exact zero so `0 * NaN` does not poison a component this term does not enter.\n      if (g !== 0) dx += g * pNu;\n      let pDot = 0;\n      const offset = dgBase + nu * dim;\n      for (let rho = 0; rho < dim; rho++) {\n        const dg = dgInv[offset + rho]!;\n        if (dg !== 0) pDot += dg * y[dim + rho]!;\n      }\n      if (pDot !== 0) dp += pDot * pNu;\n    }\n    out[mu] = dx;\n    out[dim + mu] = dp === 0 ? 0 : -0.5 * dp;\n  }\n  return out;\n}\n\n/**\n * The state after one GL4 step, with the Picard iteration count that step\n * consumed. MathTS `gaussLegendre4` does not expose its two internal stage\n * values, so this is the advanced state and nothing else.\n * @internal\n */\nexport interface GL4StepResult {\n  readonly x: readonly number[];\n  readonly p: readonly number[];\n  readonly iterations: number;\n}\n\n/**\n * One GL4 step of the geodesic Hamiltonian, via MathTS `gaussLegendre4`.\n * Throws `GL4ConvergenceError` with a message matching\n * `/Picard iteration did not converge/` when the Picard cap is exhausted;\n * any other error from the metric closures propagates unchanged.\n * @internal\n */\nexport function gl4Step(\n  state: GL4State,\n  h: number,\n  gInverseFn: (x: readonly number[]) => Float64Array,\n  dgInverseFn: (x: readonly number[]) => Float64Array,\n  opts: { picardTol: number; picardMaxIter: number },\n): GL4StepResult {\n  const dim = state.x.length;\n  const y0 = new Array<number>(dim * 2);\n  for (let i = 0; i < dim; i++) {\n    y0[i] = state.x[i]!;\n    y0[dim + i] = state.p[i]!;\n  }\n  const xBuf = new Array<number>(dim).fill(0);\n  let solved: ReturnType<typeof gaussLegendre4>;\n  try {\n    solved = gaussLegendre4(\n      (_t, y) => geodesicDeriv(y, gInverseFn, dgInverseFn, xBuf),\n      y0,\n      [0, h],\n      { steps: 1, picardTol: opts.picardTol, picardMaxIter: opts.picardMaxIter },\n    );\n  } catch (err) {\n    // MathTS throws its own typed GL4ConvergenceError for a Picard step that does not converge.\n    if (err instanceof MathTsGL4ConvergenceError) throw new GL4ConvergenceError(err.message);\n    throw err;\n  }\n  const y1 = solved.y[solved.y.length - 1];\n  const iterations = solved.iterations[0];\n  if (y1 === undefined || iterations === undefined) {\n    throw new GL4ConvergenceError(\n      `Picard iteration did not converge in ${opts.picardMaxIter} iterations (maxDelta above picardTol=${opts.picardTol})`,\n    );\n  }\n  return { x: y1.slice(0, dim), p: y1.slice(dim), iterations };\n}\n\n// ---------------------------------------------------------------------------\n// Integrator entry-point (Task 3, Phase 1a-iii)\n// ---------------------------------------------------------------------------\n\n/**\n * GL4 symplectic integrator on the canonical (x, p) geodesic Hamiltonian.\n *\n *   H(x, p) = ½ g^{μν}(x) p_μ p_ν\n *\n * Each accepted step is MathTS `gaussLegendre4` on the geodesic Hamiltonian.\n * Adaptive step-halving stays here. `onStep` reads that step's `iterations`.\n * with **adaptive step-halving on Picard non-convergence** (Adam+Eve I4,\n * replaces the single-retry R8): if Picard fails at step size h, retry at\n * h/2, h/4, … down to `hMin` (default `h · 1e-9`); throw\n * `GL4ConvergenceError` with a diagnostic message only when h_min is also\n * exhausted.\n *\n * Symplecticity (preservation of ω = dp_μ ∧ dx^μ) is a property of the\n * Butcher tableau, not of the inner Picard solver — see Sanz-Serna 1988,\n * Hairer/Lubich/Wanner §II.1. Hamiltonian drift over long integrations is\n * bounded; for non-resonant systems it remains O(h^p) over exponentially\n * long times (`p = 4` for GL4).\n *\n * **Domain guard.** If `domainMinRadius` is provided and `initialState.x[1]`\n * (radial coordinate) is below the bound, throws `NumericalBackendError`\n * synchronously with a `/domain/i`-matching message. The mid-trajectory\n * domain crossing is not checked here — callers needing that supply a\n * `gInverseFn` that throws on out-of-domain input.\n *\n * **Units.** The integrator is metric-agnostic — units follow the units of\n * the supplied `gInverseFn` and `initialState`. For UPT's canonical SI\n * Schwarzschild applications (the covariant-eikonal method used by\n * catalog relation be-37, and the perihelion method used by catalog\n * relation be-52):\n *   - `initialState.x` — `(t, r, θ, φ)` in **(s, m, rad, rad)** (SI).\n *   - `initialState.p` — covariant 4-momentum `p_μ = g_μν v^ν` in\n *     **(J·s, kg·m, kg·m², kg·m²)** under the affine normalization\n *     `p_t = −c²` used by `evaluateCovariantEikonalNumerical`.\n *   - `tauMax` — affine-parameter (proper-time for timelike, coordinate-\n *     time-like for the null normalization) extent in **seconds** under\n *     the covariant-eikonal method's convention (catalog relation be-37);\n *     **dimensionless** if the caller chose\n *     geometric units. The integrator does not enforce a choice.\n *   - `domainMinRadius` — radial coordinate lower bound in the same length\n *     units as `initialState.x[1]` (typically **metres** for SI).\n *\n * @param initialState — canonical (x, p) at τ = 0. Units follow the\n *   `gInverseFn` convention (see above).\n * @param options — see {@link GL4Options}:\n *   - `steps` — positive integer step count (dimensionless).\n *   - `tauMax` — positive finite affine-parameter extent (seconds in\n *     canonical SI). The integrator runs forward from τ = 0 only.\n *   - `gInverseFn(x)[μ*dim + ν]` — inverse metric g^{μν}(x), row-major flat.\n *   - `dgInverseFn(x)[λ*dim² + μ*dim + ν]` — ∂_λ g^{μν}(x), flat.\n *   - `picardTol` — convergence tolerance (dimensionless, default 1e-12).\n *   - `picardMaxIter` — fixed-point iteration cap (dimensionless, default 50).\n *   - `hMin` — step-halving floor (same units as `tauMax / steps`).\n *   - `domainMinRadius` — radial cutoff (same units as `initialState.x[1]`).\n * @returns `steps + 1` snapshots: index `0` is the initial state, index `n`\n *   is the state after `n` steps (τ = n · h). Each snapshot carries `tau`,\n *   `x`, `p` (and optional `v` = g^{μν} p_ν) in the units chosen above.\n * @throws NumericalBackendError if `steps` is not a positive integer, if\n *   `tauMax` is not a positive finite number, or if\n *   `initialState.x[1] < domainMinRadius`. The two option checks run before\n *   any metric evaluation.\n * @throws GL4ConvergenceError if Picard fails even after step-halving to h_min.\n *   Any other error thrown by `gInverseFn` or `dgInverseFn` (a domain\n *   crossing the caller signals by throwing, or a programming error) is\n *   rethrown unchanged; it is not halved away.\n *\n * @public\n */\nexport function integrateGeodesicGL4(\n  initialState: GL4State,\n  options: GL4Options,\n): readonly GL4Snapshot[] {\n  const {\n    steps,\n    tauMax,\n    gInverseFn,\n    dgInverseFn,\n    picardTol = 1e-12,\n    picardMaxIter = 50,\n    hMin,\n    domainMinRadius,\n    onStep,\n  } = options;\n\n  if (!Number.isInteger(steps) || steps <= 0) {\n    throw new NumericalBackendError(`GL4 integrator: steps must be a positive integer, got ${steps}`);\n  }\n  if (!Number.isFinite(tauMax) || tauMax <= 0) {\n    throw new NumericalBackendError(`GL4 integrator: tauMax must be a positive finite number, got ${tauMax}`);\n  }\n  if (domainMinRadius !== undefined && initialState.x[1] < domainMinRadius) {\n    throw new NumericalBackendError(\n      `GL4 integrator: initial r=${initialState.x[1]} < domainMinRadius=${domainMinRadius} (domain violation)`,\n    );\n  }\n\n  const h = tauMax / steps;\n  const hFloor = hMin ?? h * 1e-9;\n  const snapshots: GL4Snapshot[] = [{ tau: 0, x: initialState.x.slice(), p: initialState.p.slice() }];\n  let x = initialState.x.slice();\n  let p = initialState.p.slice();\n\n  const stateDim = x.length;\n  let newX = new Array<number>(stateDim);\n  let newP = new Array<number>(stateDim);\n\n  for (let n = 0; n < steps; n++) {\n    // I4: adaptive step-halving loop (not single-retry) on Picard\n    // non-convergence. When the full step fails, we advance by the SMALLER\n    // converging step and sub-step until the macro-step `h` is covered — the\n    // stages are solved at `stepH`, so the state/τ update MUST use `stepH`\n    // too (advancing by full `h` after a halved solve destroys accuracy and\n    // symplecticity). τ is reported on the fixed `(n+1)·h` grid; only the\n    // residual-to-cover accumulator tracks intra-step progress.\n    let remaining = h;\n    let stepH = h;\n    let halvings = 0;\n    let lastIterations = 0;\n    // Guard against FP residue: treat anything below a tiny fraction of h as done.\n    const remainEps = h * 1e-12;\n    while (remaining > remainEps) {\n      const trialH = Math.min(stepH, remaining);\n      let advanced: GL4StepResult | undefined;\n      let stepSucceeded = false;\n      let subH = trialH;\n      while (subH >= hFloor) {\n        try {\n          advanced = gl4Step({ x, p }, subH, gInverseFn, dgInverseFn, {\n            picardTol,\n            picardMaxIter,\n          });\n          stepSucceeded = true;\n          break;\n        } catch (err) {\n          // Only a Picard failure is a reason to halve. A throw from the metric\n          // closures (a domain crossing the caller signals, or a bug) is that error.\n          if (!(err instanceof GL4ConvergenceError)) throw err;\n          subH /= 2;\n          halvings++;\n        }\n      }\n      if (!stepSucceeded || advanced === undefined) {\n        throw new GL4ConvergenceError(\n          `GL4 integrator: Picard iteration did not converge even at h_min=${hFloor} (step ${n}). Diagnose step-size or metric singularity.`,\n        );\n      }\n      if (stateDim === 4) {\n        newX[0] = advanced.x[0]!; newP[0] = advanced.p[0]!;\n        newX[1] = advanced.x[1]!; newP[1] = advanced.p[1]!;\n        newX[2] = advanced.x[2]!; newP[2] = advanced.p[2]!;\n        newX[3] = advanced.x[3]!; newP[3] = advanced.p[3]!;\n      } else {\n        for (let mu = 0; mu < stateDim; mu++) {\n          newX[mu] = advanced.x[mu]!;\n          newP[mu] = advanced.p[mu]!;\n        }\n      }\n\n      const tmpX = x;\n      x = newX;\n      newX = tmpX;\n\n      const tmpP = p;\n      p = newP;\n      newP = tmpP;\n\n      remaining -= subH;\n      stepH = subH; // keep the converging size for the rest of this macro-step\n      lastIterations = advanced.iterations;\n    }\n    snapshots.push({ tau: (n + 1) * h, x: x.slice(), p: p.slice() });\n    if (onStep !== undefined) {\n      onStep({ step: n, iterations: lastIterations, halvings });\n    }\n  }\n\n  return snapshots;\n}\n"]}

package/dist/numerical/grid-field.js.map

@@ -1 +1 @@
-{"version":3,"file":"grid-field.js","sourceRoot":"","sources":["../../src/numerical/grid-field.ts"],"names":[],"mappings":""}
+{"version":3,"file":"grid-field.js","sourceRoot":"","sources":["../../src/numerical/grid-field.ts"],"names":[],"mappings":"","sourcesContent":["/**\n * GridField — a sampled field on a regular grid, for the 'grid' numericalForm\n * finite-difference path. See v0.3.5-Design.md §5.\n *\n * The interface now lives in the leaf `./types.js` module (it only depends on\n * `NestedArray`), which removes the former `types ↔ grid-field` type cycle. This\n * module re-exports it so the established `numerical/grid-field` import path —\n * used by `pderiv.ts`, `numerical/index.ts`, and the public `src/index.ts` — keeps\n * working unchanged.\n *\n * @module numerical/grid-field\n */\nexport type { GridField } from './types.js';\n"]}

package/dist/numerical/index.js.map

@@ -1 +1 @@
-{"version":3,"file":"index.js","sourceRoot":"","sources":["../../src/numerical/index.ts"],"names":[],"mappings":"AAQA,OAAO,EAAE,QAAQ,EAAE,MAAM,6BAA6B,CAAC;AAGvD,OAAO,EAAE,SAAS,EAAE,MAAM,eAAe,CAAC;AAC1C,OAAO,EAAE,eAAe,EAAE,MAAM,sBAAsB,CAAC;AACvD,OAAO,EAAE,qBAAqB,EAAE,MAAM,aAAa,CAAC;AACpD,OAAO,EAAE,qBAAqB,EAAE,iBAAiB,EAAE,MAAM,qBAAqB,CAAC;AAI/E,cAAc;AACd,OAAO,EAAE,kBAAkB,EAAE,qBAAqB,EAAE,MAAM,oBAAoB,CAAC;AAK/E,cAAc;AACd,OAAO,EAAE,eAAe,EAAE,eAAe,EAAE,MAAM,sBAAsB,CAAC;AACxE,cAAc;AACd,OAAO,EAAE,qBAAqB,EAAE,MAAM,aAAa,CAAC;AACpD;;0EAE0E;AAC1E,OAAO,EAAE,0BAA0B,EAAE,MAAM,0BAA0B,CAAC;AACtE;;;;GAIG;AACH,OAAO,EAAE,qBAAqB,EAAE,MAAM,qBAAqB,CAAC;AAC5D,mEAAmE;AACnE,OAAO,EACL,iCAAiC,GAClC,MAAM,wBAAwB,CAAC;AAOhC,kDAAkD;AAClD,OAAO,EAAE,oBAAoB,EAAE,MAAM,qBAAqB,CAAC;AAI3D,6DAA6D;AAC7D,OAAO,EAAE,cAAc,EAAE,MAAM,wBAAwB,CAAC;AAoCxD,wEAAwE;AACxE,KAAK,UAAU,OAAO,CAAC,IAAc,EAAE,OAAyB;IAC9D,MAAM,EAAE,GAAG,QAAQ,CAAC,IAAI,CAAC,CAAC;IAC1B,IAAI,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;QACX,MAAM,IAAI,qBAAqB,CAC7B,0DAA0D;cACxD,GAAG,EAAE,CAAC,UAAU,CAAC,MAAM,iBAAiB;cACxC,EAAE,CAAC,UAAU,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,IAAI,CAAC,CAAC,QAAQ,KAAK,CAAC,CAAC,IAAI,EAAE,CAAC,CAAC,IAAI,CAAC,IAAI,CAAC,CACnE,CAAC;IACJ,CAAC;IACD,MAAM,MAAM,GAAG,OAAO,EAAE,MAAM,IAAI,MAAM,eAAe,EAAE,CAAC;IAC1D,MAAM,QAAQ,GAAG,EAAE,CAAC,UAAU,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,QAAQ,KAAK,SAAS,CAAC,CAAC;IACvE,OAAO,EAAE,MAAM,EAAE,GAAG,EAAE,EAAE,CAAC,iBAA8B,EAAE,WAAW,EAAE,EAAE,CAAC,WAAW,EAAE,QAAQ,EAAE,CAAC;AACnG,CAAC;AAED;;mDAEmD;AACnD,KAAK,UAAU,4BAA4B,CACzC,IAAc,EACd,MAAuB,EACvB,MAAoB;IAEpB,MAAM,IAAI,GAAG,iBAAiB,CAAC,IAAI,CAAC,CAAC;IACrC,IAAI,CAAC,IAAI;QAAE,OAAO,EAAE,CAAC;IACrB,MAAM,EAAE,OAAO,EAAE,GAAG,MAAM,qBAAqB,CAC7C,IAAI,CAAC,MAAM,EAAE,IAAI,CAAC,MAAM,EAAE,MAAM,EAAE,SAAS,EAAE,EAAE,MAAM,EAAE,CACxD,CAAC;IACF,OAAO,OAAO,CAAC,CAAC,CAAC,CAAC,OAAO,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC;AAClC,CAAC;AAED;;;GAGG;AACH,MAAM,CAAC,KAAK,UAAU,iBAAiB,CACrC,IAAc,EACd,MAAuB,EACvB,OAAyB;IAEzB,MAAM,EAAE,MAAM,EAAE,GAAG,EAAE,WAAW,EAAE,QAAQ,EAAE,GAAG,MAAM,OAAO,CAAC,IAAI,EAAE,OAAO,CAAC,CAAC;IAC5E,MAAM,MAAM,GAAG,SAAS,CAAC,IAAI,EAAE,MAAM,EAAE,MAAM,CAAC,CAAC;IAC/C,MAAM,qBAAqB,GAAG,MAAM,4BAA4B,CAAC,IAAI,EAAE,MAAM,EAAE,MAAM,CAAC,CAAC;IACvF,OAAO;QACL,KAAK,EAAE,MAAM,CAAC,QAAQ,CAAC,MAAM,CAAC;QAC9B,GAAG;QACH,WAAW;QACX,QAAQ,EAAE,CAAC,GAAG,QAAQ,EAAE,GAAG,qBAAqB,CAAC;KAClD,CAAC;AACJ,CAAC;AAED;;;GAGG;AACH,MAAM,CAAC,KAAK,UAAU,oBAAoB,CACxC,IAAc,EACd,MAAuB,EACvB,OAAyB;IAEzB,MAAM,EAAE,MAAM,EAAE,GAAG,EAAE,WAAW,EAAE,QAAQ,EAAE,GAAG,MAAM,OAAO,CAAC,IAAI,EAAE,OAAO,CAAC,CAAC;IAC5E,MAAM,MAAM,GAAG,SAAS,CAAC,IAAI,EAAE,MAAM,EAAE,MAAM,CAAC,CAAC;IAC/C,MAAM,qBAAqB,GAAG,MAAM,4BAA4B,CAAC,IAAI,EAAE,MAAM,EAAE,MAAM,CAAC,CAAC;IACvF,OAAO;QACL,KAAK,EAAE,MAAM;QACb,GAAG;QACH,WAAW;QACX,QAAQ,EAAE,CAAC,GAAG,QAAQ,EAAE,GAAG,qBAAqB,CAAC;QACjD,OAAO,EAAE,GAAG,EAAE,GAAG,MAAM,CAAC,OAAO,EAAE,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC;KAC7C,CAAC;AACJ,CAAC"}
+{"version":3,"file":"index.js","sourceRoot":"","sources":["../../src/numerical/index.ts"],"names":[],"mappings":"AAQA,OAAO,EAAE,QAAQ,EAAE,MAAM,6BAA6B,CAAC;AAGvD,OAAO,EAAE,SAAS,EAAE,MAAM,eAAe,CAAC;AAC1C,OAAO,EAAE,eAAe,EAAE,MAAM,sBAAsB,CAAC;AACvD,OAAO,EAAE,qBAAqB,EAAE,MAAM,aAAa,CAAC;AACpD,OAAO,EAAE,qBAAqB,EAAE,iBAAiB,EAAE,MAAM,qBAAqB,CAAC;AAI/E,cAAc;AACd,OAAO,EAAE,kBAAkB,EAAE,qBAAqB,EAAE,MAAM,oBAAoB,CAAC;AAK/E,cAAc;AACd,OAAO,EAAE,eAAe,EAAE,eAAe,EAAE,MAAM,sBAAsB,CAAC;AACxE,cAAc;AACd,OAAO,EAAE,qBAAqB,EAAE,MAAM,aAAa,CAAC;AACpD;;0EAE0E;AAC1E,OAAO,EAAE,0BAA0B,EAAE,MAAM,0BAA0B,CAAC;AACtE;;;;GAIG;AACH,OAAO,EAAE,qBAAqB,EAAE,MAAM,qBAAqB,CAAC;AAC5D,mEAAmE;AACnE,OAAO,EACL,iCAAiC,GAClC,MAAM,wBAAwB,CAAC;AAOhC,kDAAkD;AAClD,OAAO,EAAE,oBAAoB,EAAE,MAAM,qBAAqB,CAAC;AAI3D,6DAA6D;AAC7D,OAAO,EAAE,cAAc,EAAE,MAAM,wBAAwB,CAAC;AAoCxD,wEAAwE;AACxE,KAAK,UAAU,OAAO,CAAC,IAAc,EAAE,OAAyB;IAC9D,MAAM,EAAE,GAAG,QAAQ,CAAC,IAAI,CAAC,CAAC;IAC1B,IAAI,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;QACX,MAAM,IAAI,qBAAqB,CAC7B,0DAA0D;cACxD,GAAG,EAAE,CAAC,UAAU,CAAC,MAAM,iBAAiB;cACxC,EAAE,CAAC,UAAU,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,IAAI,CAAC,CAAC,QAAQ,KAAK,CAAC,CAAC,IAAI,EAAE,CAAC,CAAC,IAAI,CAAC,IAAI,CAAC,CACnE,CAAC;IACJ,CAAC;IACD,MAAM,MAAM,GAAG,OAAO,EAAE,MAAM,IAAI,MAAM,eAAe,EAAE,CAAC;IAC1D,MAAM,QAAQ,GAAG,EAAE,CAAC,UAAU,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,QAAQ,KAAK,SAAS,CAAC,CAAC;IACvE,OAAO,EAAE,MAAM,EAAE,GAAG,EAAE,EAAE,CAAC,iBAA8B,EAAE,WAAW,EAAE,EAAE,CAAC,WAAW,EAAE,QAAQ,EAAE,CAAC;AACnG,CAAC;AAED;;mDAEmD;AACnD,KAAK,UAAU,4BAA4B,CACzC,IAAc,EACd,MAAuB,EACvB,MAAoB;IAEpB,MAAM,IAAI,GAAG,iBAAiB,CAAC,IAAI,CAAC,CAAC;IACrC,IAAI,CAAC,IAAI;QAAE,OAAO,EAAE,CAAC;IACrB,MAAM,EAAE,OAAO,EAAE,GAAG,MAAM,qBAAqB,CAC7C,IAAI,CAAC,MAAM,EAAE,IAAI,CAAC,MAAM,EAAE,MAAM,EAAE,SAAS,EAAE,EAAE,MAAM,EAAE,CACxD,CAAC;IACF,OAAO,OAAO,CAAC,CAAC,CAAC,CAAC,OAAO,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC;AAClC,CAAC;AAED;;;GAGG;AACH,MAAM,CAAC,KAAK,UAAU,iBAAiB,CACrC,IAAc,EACd,MAAuB,EACvB,OAAyB;IAEzB,MAAM,EAAE,MAAM,EAAE,GAAG,EAAE,WAAW,EAAE,QAAQ,EAAE,GAAG,MAAM,OAAO,CAAC,IAAI,EAAE,OAAO,CAAC,CAAC;IAC5E,MAAM,MAAM,GAAG,SAAS,CAAC,IAAI,EAAE,MAAM,EAAE,MAAM,CAAC,CAAC;IAC/C,MAAM,qBAAqB,GAAG,MAAM,4BAA4B,CAAC,IAAI,EAAE,MAAM,EAAE,MAAM,CAAC,CAAC;IACvF,OAAO;QACL,KAAK,EAAE,MAAM,CAAC,QAAQ,CAAC,MAAM,CAAC;QAC9B,GAAG;QACH,WAAW;QACX,QAAQ,EAAE,CAAC,GAAG,QAAQ,EAAE,GAAG,qBAAqB,CAAC;KAClD,CAAC;AACJ,CAAC;AAED;;;GAGG;AACH,MAAM,CAAC,KAAK,UAAU,oBAAoB,CACxC,IAAc,EACd,MAAuB,EACvB,OAAyB;IAEzB,MAAM,EAAE,MAAM,EAAE,GAAG,EAAE,WAAW,EAAE,QAAQ,EAAE,GAAG,MAAM,OAAO,CAAC,IAAI,EAAE,OAAO,CAAC,CAAC;IAC5E,MAAM,MAAM,GAAG,SAAS,CAAC,IAAI,EAAE,MAAM,EAAE,MAAM,CAAC,CAAC;IAC/C,MAAM,qBAAqB,GAAG,MAAM,4BAA4B,CAAC,IAAI,EAAE,MAAM,EAAE,MAAM,CAAC,CAAC;IACvF,OAAO;QACL,KAAK,EAAE,MAAM;QACb,GAAG;QACH,WAAW;QACX,QAAQ,EAAE,CAAC,GAAG,QAAQ,EAAE,GAAG,qBAAqB,CAAC;QACjD,OAAO,EAAE,GAAG,EAAE,GAAG,MAAM,CAAC,OAAO,EAAE,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC;KAC7C,CAAC;AACJ,CAAC","sourcesContent":["/**\n * Public surface of the UPT numerical-contraction backend (v0.3.5).\n * See docs/planning/v0.3.5-Design.md §5.\n *\n * @module numerical\n */\nimport type { Dimension } from '../dimensional/types.js';\nimport type { ExprNode, Violation } from '../dimensional/validator.js';\nimport { validate } from '../dimensional/validator.js';\nimport type { EngineTensor, TensorEngine } from './tensor-engine.js';\nimport type { NumericalInputs, NestedArray } from './types.js';\nimport { lowerNode } from './lowering.js';\nimport { getActiveEngine } from './engine-registry.js';\nimport { NumericalBackendError } from './errors.js';\nimport { evaluateMetricInverse, scanForMetricPair } from './metric-inverse.js';\n\n/** @public */\nexport type { TensorEngine, EngineTensor, EinsumSpec, ForwardGradResult, ReverseGradResult } from './tensor-engine.js';\n/** @public */\nexport { hasAutogradSupport, EngineCapabilityError } from './tensor-engine.js';\n/** @public */\nexport type { NumericalInputs, NestedArray } from './types.js';\n/** @public — part of the `NumericalInputs.grids` public contract. */\nexport type { GridField } from './grid-field.js';\n/** @public */\nexport { getActiveEngine, setActiveEngine } from './engine-registry.js';\n/** @public */\nexport { NumericalBackendError } from './errors.js';\n/** @public — re-exported from dimensional/errors to keep numerical/index.ts\n *  as the single public API surface, without creating a dimensional→numerical\n *  import cycle. See v0.4.0-Implementation-Plan Task 13 for rationale. */\nexport { DuplicateCoordinateWarning } from '../dimensional/errors.js';\n/**\n * Evaluate a metric inverse from a detected pair.\n *\n * @public\n */\nexport { evaluateMetricInverse } from './metric-inverse.js';\n/** Integrate the covariant eikonal on a numerical grid. @public */\nexport {\n  evaluateCovariantEikonalNumerical,\n} from './covariant-eikonal.js';\n/** Inputs and the result of the covariant eikonal integrator. */\nexport type {\n  CovariantEikonalInputs,\n  CovariantEikonalResult,\n} from './covariant-eikonal.js';\n\n/** @public — v0.5.0 GL4 symplectic integrator. */\nexport { integrateGeodesicGL4 } from './gl4-integrator.js';\n/** @public — v0.5.0 GL4 types. */\nexport type { GL4State, GL4Snapshot, GL4Options } from './gl4-integrator.js';\n\n/** @public — v0.5.0 perihelion finder (Task 4, Phase 1b). */\nexport { findPerihelion } from './perihelion-finder.js';\n/** @public — v0.5.0 perihelion-finder types. */\nexport type { PerihelionResult, FindPerihelionOptions } from './perihelion-finder.js';\n\n/**\n * Plain-JS result of `evaluateNumerical`.\n * @public\n */\nexport interface NumericalResult {\n  readonly value: NestedArray;\n  readonly dim: Dimension;\n  readonly freeIndices: ReadonlyMap<string, { upper: number; lower: number }>;\n  readonly warnings: ReadonlyArray<Violation>;\n}\n\n/**\n * Result of `evaluateNumericalRaw` — carries a live `EngineTensor` for\n * chaining workloads; the caller must `dispose()` it.\n * @public\n */\nexport interface NumericalRawResult {\n  readonly value: EngineTensor;\n  readonly dim: Dimension;\n  readonly freeIndices: ReadonlyMap<string, { upper: number; lower: number }>;\n  readonly warnings: ReadonlyArray<Violation>;\n  dispose(): void;\n}\n\n/**\n * Per-call options for the `evaluateNumerical*` entry points.\n * @public\n */\nexport interface EvaluateOptions {\n  readonly engine?: TensorEngine;\n}\n\n/** Validate the AST and return (engine, validation result) or throw. */\nasync function prepare(node: ExprNode, options?: EvaluateOptions) {\n  const vr = validate(node);\n  if (!vr.ok) {\n    throw new NumericalBackendError(\n      `evaluateNumerical: the AST is not dimensionally valid — `\n      + `${vr.violations.length} violation(s): `\n      + vr.violations.map((v) => `[${v.location}] ${v.note}`).join('; '),\n    );\n  }\n  const engine = options?.engine ?? await getActiveEngine();\n  const warnings = vr.violations.filter((v) => v.severity === 'warning');\n  return { engine, dim: vr.inferredDimension as Dimension, freeIndices: vr.freeIndices, warnings };\n}\n\n/** If the AST contains an identifiable lower/upper metric pair, run the\n *  numerical inverse-metric check and return any warning (deduplicated —\n *  scanForMetricPair returns at most one pair). */\nasync function collectInverseMetricWarnings(\n  node: ExprNode,\n  inputs: NumericalInputs,\n  engine: TensorEngine,\n): Promise<Violation[]> {\n  const pair = scanForMetricPair(node);\n  if (!pair) return [];\n  const { warning } = await evaluateMetricInverse(\n    pair.gUpper, pair.gLower, inputs, undefined, { engine },\n  );\n  return warning ? [warning] : [];\n}\n\n/**\n * Evaluate a validated AST to plain JS.\n * @public\n */\nexport async function evaluateNumerical(\n  node: ExprNode,\n  inputs: NumericalInputs,\n  options?: EvaluateOptions,\n): Promise<NumericalResult> {\n  const { engine, dim, freeIndices, warnings } = await prepare(node, options);\n  const tensor = lowerNode(node, inputs, engine);\n  const inverseMetricWarnings = await collectInverseMetricWarnings(node, inputs, engine);\n  return {\n    value: engine.toNested(tensor),\n    dim,\n    freeIndices,\n    warnings: [...warnings, ...inverseMetricWarnings],\n  };\n}\n\n/**\n * Evaluate to a live EngineTensor for chaining workloads.\n * @public\n */\nexport async function evaluateNumericalRaw(\n  node: ExprNode,\n  inputs: NumericalInputs,\n  options?: EvaluateOptions,\n): Promise<NumericalRawResult> {\n  const { engine, dim, freeIndices, warnings } = await prepare(node, options);\n  const tensor = lowerNode(node, inputs, engine);\n  const inverseMetricWarnings = await collectInverseMetricWarnings(node, inputs, engine);\n  return {\n    value: tensor,\n    dim,\n    freeIndices,\n    warnings: [...warnings, ...inverseMetricWarnings],\n    dispose: () => { engine.dispose?.(tensor); },\n  };\n}\n"]}

package/dist/numerical/input-validation.d.ts

@@ -1,12 +1,9 @@
 /**
  * Runtime input validation for numeric evaluators — `validateFiniteInputs`.
  *
- * A generic, dependency-free finite-and-in-range check shared by the numerical
- * evaluators (`klein-gordon.ts`) and the 42 bridge-equation evaluators
- * (`bridges/equations/be-*.ts`). It lived in `bridges/equations/_be-helpers.ts`,
- * which forced the lower numerical layer to import upward from `bridges`; moving
- * it here (a leaf module with no imports) removes that dependency while
- * `_be-helpers.ts` re-exports it so every bridge evaluator is unchanged.
+ * A generic, dependency-free finite-and-in-range check for a numeric
+ * evaluator's input record (`klein-gordon.ts` reads it). It is a leaf module
+ * with no imports, so any layer may read it without an upward edge.
  *
  * @module numerical/input-validation
  */

package/dist/numerical/input-validation.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"input-validation.d.ts","sourceRoot":"","sources":["../../src/numerical/input-validation.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;GAWG;AAEH;;;;;;;;;;GAUG;AACH,MAAM,WAAW,SAAS;IACxB,yDAAyD;IACzD,QAAQ,CAAC,IAAI,EAAE,MAAM,CAAC;IACtB,mDAAmD;IACnD,QAAQ,CAAC,GAAG,CAAC,EAAE,MAAM,CAAC;IACtB,kDAAkD;IAClD,QAAQ,CAAC,GAAG,CAAC,EAAE,MAAM,CAAC;IACtB,oDAAoD;IACpD,QAAQ,CAAC,UAAU,CAAC,EAAE,OAAO,CAAC;IAC9B,oDAAoD;IACpD,QAAQ,CAAC,UAAU,CAAC,EAAE,OAAO,CAAC;IAC9B,wEAAwE;IACxE,QAAQ,CAAC,SAAS,CAAC,EAAE,OAAO,CAAC;IAC7B;;;;;OAKG;IACH,QAAQ,CAAC,WAAW,CAAC,EAAE,MAAM,CAAC;CAC/B;AAgCD;;;;;;;;;GASG;AACH,wBAAgB,oBAAoB,CAAC,CAAC,SAAS,MAAM,EACnD,KAAK,EAAE,CAAC,EACR,UAAU,EAAE,aAAa,CAAC,SAAS,CAAC,EACpC,aAAa,EAAE,MAAM,GACpB,IAAI,CAgCN"}
+{"version":3,"file":"input-validation.d.ts","sourceRoot":"","sources":["../../src/numerical/input-validation.ts"],"names":[],"mappings":"AAAA;;;;;;;;GAQG;AAEH;;;;;;;;;;GAUG;AACH,MAAM,WAAW,SAAS;IACxB,yDAAyD;IACzD,QAAQ,CAAC,IAAI,EAAE,MAAM,CAAC;IACtB,mDAAmD;IACnD,QAAQ,CAAC,GAAG,CAAC,EAAE,MAAM,CAAC;IACtB,kDAAkD;IAClD,QAAQ,CAAC,GAAG,CAAC,EAAE,MAAM,CAAC;IACtB,oDAAoD;IACpD,QAAQ,CAAC,UAAU,CAAC,EAAE,OAAO,CAAC;IAC9B,oDAAoD;IACpD,QAAQ,CAAC,UAAU,CAAC,EAAE,OAAO,CAAC;IAC9B,wEAAwE;IACxE,QAAQ,CAAC,SAAS,CAAC,EAAE,OAAO,CAAC;IAC7B;;;;;OAKG;IACH,QAAQ,CAAC,WAAW,CAAC,EAAE,MAAM,CAAC;CAC/B;AAgCD;;;;;;;;;GASG;AACH,wBAAgB,oBAAoB,CAAC,CAAC,SAAS,MAAM,EACnD,KAAK,EAAE,CAAC,EACR,UAAU,EAAE,aAAa,CAAC,SAAS,CAAC,EACpC,aAAa,EAAE,MAAM,GACpB,IAAI,CA+BN"}

package/dist/numerical/input-validation.js

@@ -1,12 +1,9 @@
 /**
  * Runtime input validation for numeric evaluators — `validateFiniteInputs`.
  *
- * A generic, dependency-free finite-and-in-range check shared by the numerical
- * evaluators (`klein-gordon.ts`) and the 42 bridge-equation evaluators
- * (`bridges/equations/be-*.ts`). It lived in `bridges/equations/_be-helpers.ts`,
- * which forced the lower numerical layer to import upward from `bridges`; moving
- * it here (a leaf module with no imports) removes that dependency while
- * `_be-helpers.ts` re-exports it so every bridge evaluator is unchanged.
+ * A generic, dependency-free finite-and-in-range check for a numeric
+ * evaluator's input record (`klein-gordon.ts` reads it). It is a leaf module
+ * with no imports, so any layer may read it without an upward edge.
  *
  * @module numerical/input-validation
  */
@@ -50,8 +47,7 @@ function describeRange(spec) {
  */
 export function validateFiniteInputs(input, fieldSpecs, evaluatorName) {
     for (const spec of fieldSpecs) {
-        const raw = input[spec.name];
-        const value = raw;
+        const value = input[spec.name];
         if (typeof value !== 'number' || !Number.isFinite(value)) {
             const range = describeRange(spec);
             throw new RangeError(`${evaluatorName}: ${spec.name} must be a finite ${range}number, got ${String(value)}`);

package/dist/numerical/input-validation.js.map

@@ -1 +1 @@
-{"version":3,"file":"input-validation.js","sourceRoot":"","sources":["../../src/numerical/input-validation.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;GAWG;AAmCH;;;;;;;GAOG;AACH,SAAS,aAAa,CAAC,IAAe;IACpC,IAAI,IAAI,CAAC,WAAW,KAAK,SAAS;QAAE,OAAO,GAAG,IAAI,CAAC,WAAW,GAAG,CAAC;IAElE,MAAM,MAAM,GAAG,IAAI,CAAC,GAAG,KAAK,SAAS,CAAC;IACtC,MAAM,MAAM,GAAG,IAAI,CAAC,GAAG,KAAK,SAAS,CAAC;IACtC,MAAM,OAAO,GAAG,CAAC,CAAC,IAAI,CAAC,UAAU,CAAC;IAClC,MAAM,QAAQ,GAAG,CAAC,CAAC,IAAI,CAAC,UAAU,CAAC;IAEnC,IAAI,MAAM,IAAI,CAAC,MAAM,IAAI,IAAI,CAAC,GAAG,KAAK,CAAC,IAAI,CAAC,IAAI,CAAC,SAAS,EAAE,CAAC;QAC3D,IAAI,OAAO;YAAE,OAAO,WAAW,CAAC;QAChC,OAAO,eAAe,CAAC;IACzB,CAAC;IAED,IAAI,CAAC,MAAM,IAAI,CAAC,MAAM;QAAE,OAAO,EAAE,CAAC;IAElC,MAAM,EAAE,GAAG,MAAM,CAAC,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,QAAQ,CAAC;IACzC,MAAM,EAAE,GAAG,MAAM,CAAC,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,QAAQ,CAAC;IACxC,MAAM,EAAE,GAAG,OAAO,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC;IAC/B,MAAM,EAAE,GAAG,QAAQ,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC;IAChC,OAAO,YAAY,EAAE,GAAG,EAAE,KAAK,EAAE,GAAG,EAAE,GAAG,CAAC;AAC5C,CAAC;AAED;;;;;;;;;GASG;AACH,MAAM,UAAU,oBAAoB,CAClC,KAAQ,EACR,UAAoC,EACpC,aAAqB;IAErB,KAAK,MAAM,IAAI,IAAI,UAAU,EAAE,CAAC;QAC9B,MAAM,GAAG,GAAI,KAAiC,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC;QAC1D,MAAM,KAAK,GAAG,GAAa,CAAC;QAE5B,IAAI,OAAO,KAAK,KAAK,QAAQ,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,KAAK,CAAC,EAAE,CAAC;YACzD,MAAM,KAAK,GAAG,aAAa,CAAC,IAAI,CAAC,CAAC;YAClC,MAAM,IAAI,UAAU,CAClB,GAAG,aAAa,KAAK,IAAI,CAAC,IAAI,qBAAqB,KAAK,eAAe,MAAM,CAAC,KAAK,CAAC,EAAE,CACvF,CAAC;QACJ,CAAC;QAED,IAAI,IAAI,CAAC,GAAG,KAAK,SAAS,EAAE,CAAC;YAC3B,MAAM,WAAW,GAAG,IAAI,CAAC,UAAU,CAAC,CAAC,CAAC,KAAK,IAAI,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC;YAC3E,IAAI,WAAW,EAAE,CAAC;gBAChB,MAAM,KAAK,GAAG,aAAa,CAAC,IAAI,CAAC,CAAC;gBAClC,MAAM,IAAI,UAAU,CAClB,GAAG,aAAa,KAAK,IAAI,CAAC,IAAI,qBAAqB,KAAK,eAAe,KAAK,EAAE,CAC/E,CAAC;YACJ,CAAC;QACH,CAAC;QAED,IAAI,IAAI,CAAC,GAAG,KAAK,SAAS,EAAE,CAAC;YAC3B,MAAM,WAAW,GAAG,IAAI,CAAC,UAAU,CAAC,CAAC,CAAC,KAAK,IAAI,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC;YAC3E,IAAI,WAAW,EAAE,CAAC;gBAChB,MAAM,KAAK,GAAG,aAAa,CAAC,IAAI,CAAC,CAAC;gBAClC,MAAM,IAAI,UAAU,CAClB,GAAG,aAAa,KAAK,IAAI,CAAC,IAAI,qBAAqB,KAAK,eAAe,KAAK,EAAE,CAC/E,CAAC;YACJ,CAAC;QACH,CAAC;IACH,CAAC;AACH,CAAC"}
+{"version":3,"file":"input-validation.js","sourceRoot":"","sources":["../../src/numerical/input-validation.ts"],"names":[],"mappings":"AAAA;;;;;;;;GAQG;AAmCH;;;;;;;GAOG;AACH,SAAS,aAAa,CAAC,IAAe;IACpC,IAAI,IAAI,CAAC,WAAW,KAAK,SAAS;QAAE,OAAO,GAAG,IAAI,CAAC,WAAW,GAAG,CAAC;IAElE,MAAM,MAAM,GAAG,IAAI,CAAC,GAAG,KAAK,SAAS,CAAC;IACtC,MAAM,MAAM,GAAG,IAAI,CAAC,GAAG,KAAK,SAAS,CAAC;IACtC,MAAM,OAAO,GAAG,CAAC,CAAC,IAAI,CAAC,UAAU,CAAC;IAClC,MAAM,QAAQ,GAAG,CAAC,CAAC,IAAI,CAAC,UAAU,CAAC;IAEnC,IAAI,MAAM,IAAI,CAAC,MAAM,IAAI,IAAI,CAAC,GAAG,KAAK,CAAC,IAAI,CAAC,IAAI,CAAC,SAAS,EAAE,CAAC;QAC3D,IAAI,OAAO;YAAE,OAAO,WAAW,CAAC;QAChC,OAAO,eAAe,CAAC;IACzB,CAAC;IAED,IAAI,CAAC,MAAM,IAAI,CAAC,MAAM;QAAE,OAAO,EAAE,CAAC;IAElC,MAAM,EAAE,GAAG,MAAM,CAAC,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,QAAQ,CAAC;IACzC,MAAM,EAAE,GAAG,MAAM,CAAC,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,QAAQ,CAAC;IACxC,MAAM,EAAE,GAAG,OAAO,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC;IAC/B,MAAM,EAAE,GAAG,QAAQ,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC;IAChC,OAAO,YAAY,EAAE,GAAG,EAAE,KAAK,EAAE,GAAG,EAAE,GAAG,CAAC;AAC5C,CAAC;AAED;;;;;;;;;GASG;AACH,MAAM,UAAU,oBAAoB,CAClC,KAAQ,EACR,UAAoC,EACpC,aAAqB;IAErB,KAAK,MAAM,IAAI,IAAI,UAAU,EAAE,CAAC;QAC9B,MAAM,KAAK,GAAa,KAAiC,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC;QAErE,IAAI,OAAO,KAAK,KAAK,QAAQ,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,KAAK,CAAC,EAAE,CAAC;YACzD,MAAM,KAAK,GAAG,aAAa,CAAC,IAAI,CAAC,CAAC;YAClC,MAAM,IAAI,UAAU,CAClB,GAAG,aAAa,KAAK,IAAI,CAAC,IAAI,qBAAqB,KAAK,eAAe,MAAM,CAAC,KAAK,CAAC,EAAE,CACvF,CAAC;QACJ,CAAC;QAED,IAAI,IAAI,CAAC,GAAG,KAAK,SAAS,EAAE,CAAC;YAC3B,MAAM,WAAW,GAAG,IAAI,CAAC,UAAU,CAAC,CAAC,CAAC,KAAK,IAAI,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC;YAC3E,IAAI,WAAW,EAAE,CAAC;gBAChB,MAAM,KAAK,GAAG,aAAa,CAAC,IAAI,CAAC,CAAC;gBAClC,MAAM,IAAI,UAAU,CAClB,GAAG,aAAa,KAAK,IAAI,CAAC,IAAI,qBAAqB,KAAK,eAAe,KAAK,EAAE,CAC/E,CAAC;YACJ,CAAC;QACH,CAAC;QAED,IAAI,IAAI,CAAC,GAAG,KAAK,SAAS,EAAE,CAAC;YAC3B,MAAM,WAAW,GAAG,IAAI,CAAC,UAAU,CAAC,CAAC,CAAC,KAAK,IAAI,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC;YAC3E,IAAI,WAAW,EAAE,CAAC;gBAChB,MAAM,KAAK,GAAG,aAAa,CAAC,IAAI,CAAC,CAAC;gBAClC,MAAM,IAAI,UAAU,CAClB,GAAG,aAAa,KAAK,IAAI,CAAC,IAAI,qBAAqB,KAAK,eAAe,KAAK,EAAE,CAC/E,CAAC;YACJ,CAAC;QACH,CAAC;IACH,CAAC;AACH,CAAC","sourcesContent":["/**\n * Runtime input validation for numeric evaluators — `validateFiniteInputs`.\n *\n * A generic, dependency-free finite-and-in-range check for a numeric\n * evaluator's input record (`klein-gordon.ts` reads it). It is a leaf module\n * with no imports, so any layer may read it without an upward edge.\n *\n * @module numerical/input-validation\n */\n\n/**\n * Declarative spec for a single field of an evaluator's `input` record.\n *\n * The default behaviour (no `min`/`max`) only verifies finiteness\n * (`Number.isFinite`). The bound predicates compose with `Number.isFinite`\n * so the resulting error always reads \"finite ⟨adj⟩ number\" — drop-in\n * compatible with the per-evaluator message format\n * `\"<evaluator>: <field> must be a finite <adjective> number, got <v>\"`.\n *\n * @internal\n */\nexport interface FieldSpec {\n  /** Field name on the `input` record (e.g., `'M_kg'`). */\n  readonly name: string;\n  /** Inclusive lower bound (default `-Infinity`). */\n  readonly min?: number;\n  /** Inclusive upper bound (default `Infinity`). */\n  readonly max?: number;\n  /** Strict `>` instead of `≥` at the lower bound. */\n  readonly excludeMin?: boolean;\n  /** Strict `<` instead of `≤` at the upper bound. */\n  readonly excludeMax?: boolean;\n  /** Explicit override — when `min === 0` and you want to permit zero. */\n  readonly allowZero?: boolean;\n  /**\n   * Optional domain-specific noun phrase (e.g., 'probability', 'angle').\n   * When set, replaces the auto-generated \"in range [a, b]\" prose:\n   *   \"<evaluator>: <field> must be a finite <description> number, got <v>\"\n   * @internal\n   */\n  readonly description?: string;\n}\n\n/**\n * Compose the adjective phrase that appears in the error message, picked to\n * match the per-evaluator prose conventions:\n *   - `min=0, excludeMin=true`  → \"positive\"        (i.e., `> 0`)\n *   - `min=0`                   → \"non-negative\"    (i.e., `≥ 0`)\n *   - bounded but not a 0-pin   → \"in range [a, b]\" (closed/open as set)\n *   - unbounded                 → \"\"                (just \"finite\")\n */\nfunction describeRange(spec: FieldSpec): string {\n  if (spec.description !== undefined) return `${spec.description} `;\n\n  const hasMin = spec.min !== undefined;\n  const hasMax = spec.max !== undefined;\n  const lowOpen = !!spec.excludeMin;\n  const highOpen = !!spec.excludeMax;\n\n  if (hasMin && !hasMax && spec.min === 0 && !spec.allowZero) {\n    if (lowOpen) return 'positive ';\n    return 'non-negative ';\n  }\n\n  if (!hasMin && !hasMax) return '';\n\n  const lo = hasMin ? spec.min : -Infinity;\n  const hi = hasMax ? spec.max : Infinity;\n  const lb = lowOpen ? '(' : '[';\n  const rb = highOpen ? ')' : ']';\n  return `in range ${lb}${lo}, ${hi}${rb} `;\n}\n\n/**\n * Validate every field of `input` named by `fieldSpecs`. Throws a `RangeError`\n * on the first violation with a message containing the keywords (`\"finite\"`,\n * plus `\"positive\"` / `\"non-negative\"` / `\"in range\"`) the per-evaluator tests\n * pattern-match against. Format:\n *\n *   `<evaluatorName>: <fieldName> must be a finite <range> number, got <value>`\n *\n * @internal\n */\nexport function validateFiniteInputs<T extends object>(\n  input: T,\n  fieldSpecs: ReadonlyArray<FieldSpec>,\n  evaluatorName: string,\n): void {\n  for (const spec of fieldSpecs) {\n    const value: unknown = (input as Record<string, unknown>)[spec.name];\n\n    if (typeof value !== 'number' || !Number.isFinite(value)) {\n      const range = describeRange(spec);\n      throw new RangeError(\n        `${evaluatorName}: ${spec.name} must be a finite ${range}number, got ${String(value)}`,\n      );\n    }\n\n    if (spec.min !== undefined) {\n      const violatesMin = spec.excludeMin ? value <= spec.min : value < spec.min;\n      if (violatesMin) {\n        const range = describeRange(spec);\n        throw new RangeError(\n          `${evaluatorName}: ${spec.name} must be a finite ${range}number, got ${value}`,\n        );\n      }\n    }\n\n    if (spec.max !== undefined) {\n      const violatesMax = spec.excludeMax ? value >= spec.max : value > spec.max;\n      if (violatesMax) {\n        const range = describeRange(spec);\n        throw new RangeError(\n          `${evaluatorName}: ${spec.name} must be a finite ${range}number, got ${value}`,\n        );\n      }\n    }\n  }\n}\n"]}

package/dist/numerical/killing.d.ts

@@ -63,8 +63,10 @@ export interface KillingEquationOptions {
     readonly dMetricFn?: (x: ReadonlyArray<number>) => number[][][];
 }
 type Vec4 = [number, number, number, number];
-type KillingFn = () => (x: Vec4) => Vec4;
-type MetricFn = (x: ReadonlyArray<number>) => number[][];
+/** A Killing-field factory: returns the field ξ^μ(x) as a function of the 4-position. @public */
+export type KillingFn = () => (x: Vec4) => Vec4;
+/** A covariant metric closure for the Killing checks: `g(x)[μ][ν] = g_{μν}(x)`, nested. @public */
+export type KillingMetricFn = (x: ReadonlyArray<number>) => number[][];
 /**
  * Layout-agnostic Christoffel accessor.
  *
@@ -77,7 +79,8 @@ type MetricFn = (x: ReadonlyArray<number>) => number[][];
  * @public
  */
 export type ChristoffelAccess = (lambda: number, mu: number, nu: number) => number;
-type ChristoffelAtFn = (x: Vec4) => ChristoffelAccess;
+/** A Christoffel producer: the {@link ChristoffelAccess} at a 4-position. @public */
+export type ChristoffelAtFn = (x: Vec4) => ChristoffelAccess;
 /**
  * Compute `||∇_μ ξ_ν + ∇_ν ξ_μ||_∞` at point `x`, returning the maximum
  * absolute symmetrized residual over all 16 (μ, ν) pairs.
@@ -128,7 +131,7 @@ type ChristoffelAtFn = (x: Vec4) => ChristoffelAccess;
  *
  * @public
  */
-export declare function verifyKillingEquation(killingFn: KillingFn, metricFn: MetricFn, christoffelAt: ChristoffelAtFn, x: Vec4, opts?: KillingEquationOptions): number;
+export declare function verifyKillingEquation(killingFn: KillingFn, metricFn: KillingMetricFn, christoffelAt: ChristoffelAtFn, x: Vec4, opts?: KillingEquationOptions): number;
 /**
  * Result of {@link checkKillingEquation}.
  *
@@ -164,7 +167,7 @@ export interface KillingEquationCheck {
  *
  * @public
  */
-export declare function checkKillingEquation(killingFn: KillingFn, metricFn: MetricFn, christoffelAt: ChristoffelAtFn, x: Vec4, opts?: KillingEquationOptions): KillingEquationCheck;
+export declare function checkKillingEquation(killingFn: KillingFn, metricFn: KillingMetricFn, christoffelAt: ChristoffelAtFn, x: Vec4, opts?: KillingEquationOptions): KillingEquationCheck;
 /**
  * Evaluate the conserved charge Q = ξ^μ p_μ at point x along a geodesic.
  *

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

@@ -1 +1 @@
-{"version":3,"file":"killing.d.ts","sourceRoot":"","sources":["../../src/numerical/killing.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;;;;;;;GAyBG;AAKH;;cAEc;AACd,MAAM,WAAW,sBAAsB;IACrC;;;;;OAKG;IACH,QAAQ,CAAC,SAAS,CAAC,EAAE,MAAM,CAAC;IAC5B;;;;OAIG;IACH,QAAQ,CAAC,KAAK,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC;IACvB;;;;;;;;;;OAUG;IACH,QAAQ,CAAC,eAAe,CAAC,EAAE,OAAO,CAAC;IACnC;;;;;;OAMG;IACH,QAAQ,CAAC,SAAS,CAAC,EAAE,CAAC,CAAC,EAAE,aAAa,CAAC,MAAM,CAAC,KAAK,MAAM,EAAE,EAAE,EAAE,CAAC;CACjE;AAED,KAAK,IAAI,GAAG,CAAC,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,CAAC,CAAC;AAC7C,KAAK,SAAS,GAAG,MAAM,CAAC,CAAC,EAAE,IAAI,KAAK,IAAI,CAAC;AACzC,KAAK,QAAQ,GAAG,CAAC,CAAC,EAAE,aAAa,CAAC,MAAM,CAAC,KAAK,MAAM,EAAE,EAAE,CAAC;AAEzD;;;;;;;;;;GAUG;AACH,MAAM,MAAM,iBAAiB,GAAG,CAAC,MAAM,EAAE,MAAM,EAAE,EAAE,EAAE,MAAM,EAAE,EAAE,EAAE,MAAM,KAAK,MAAM,CAAC;AAEnF,KAAK,eAAe,GAAG,CAAC,CAAC,EAAE,IAAI,KAAK,iBAAiB,CAAC;AAEtD;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GAiDG;AACH,wBAAgB,qBAAqB,CACnC,SAAS,EAAE,SAAS,EACpB,QAAQ,EAAE,QAAQ,EAClB,aAAa,EAAE,eAAe,EAC9B,CAAC,EAAE,IAAI,EACP,IAAI,GAAE,sBAA2B,GAChC,MAAM,CAmIR;AAED;;;;GAIG;AACH,MAAM,WAAW,oBAAoB;IACnC;wDACoD;IACpD,QAAQ,CAAC,QAAQ,EAAE,MAAM,CAAC;IAC1B;qFACiF;IACjF,QAAQ,CAAC,gBAAgB,EAAE,MAAM,CAAC;IAClC,uDAAuD;IACvD,QAAQ,CAAC,eAAe,EAAE,OAAO,CAAC;CACnC;AAED;;;;;;;;;;;;;;;;;;;GAmBG;AACH,wBAAgB,oBAAoB,CAClC,SAAS,EAAE,SAAS,EACpB,QAAQ,EAAE,QAAQ,EAClB,aAAa,EAAE,eAAe,EAC9B,CAAC,EAAE,IAAI,EACP,IAAI,GAAE,sBAA2B,GAChC,oBAAoB,CAUtB;AAED;;;;;;;;;;;;;;;;;;;;;;;;;;;;GA4BG;AACH,wBAAgB,uBAAuB,CACrC,SAAS,EAAE,SAAS,EACpB,aAAa,EAAE,IAAI,EACnB,CAAC,EAAE,IAAI,GACN,MAAM,CAKR"}
+{"version":3,"file":"killing.d.ts","sourceRoot":"","sources":["../../src/numerical/killing.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;;;;;;;GAyBG;AAKH;;cAEc;AACd,MAAM,WAAW,sBAAsB;IACrC;;;;;OAKG;IACH,QAAQ,CAAC,SAAS,CAAC,EAAE,MAAM,CAAC;IAC5B;;;;OAIG;IACH,QAAQ,CAAC,KAAK,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC;IACvB;;;;;;;;;;OAUG;IACH,QAAQ,CAAC,eAAe,CAAC,EAAE,OAAO,CAAC;IACnC;;;;;;OAMG;IACH,QAAQ,CAAC,SAAS,CAAC,EAAE,CAAC,CAAC,EAAE,aAAa,CAAC,MAAM,CAAC,KAAK,MAAM,EAAE,EAAE,EAAE,CAAC;CACjE;AAED,KAAK,IAAI,GAAG,CAAC,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,CAAC,CAAC;AAC7C,iGAAiG;AACjG,MAAM,MAAM,SAAS,GAAG,MAAM,CAAC,CAAC,EAAE,IAAI,KAAK,IAAI,CAAC;AAChD,mGAAmG;AACnG,MAAM,MAAM,eAAe,GAAG,CAAC,CAAC,EAAE,aAAa,CAAC,MAAM,CAAC,KAAK,MAAM,EAAE,EAAE,CAAC;AAEvE;;;;;;;;;;GAUG;AACH,MAAM,MAAM,iBAAiB,GAAG,CAAC,MAAM,EAAE,MAAM,EAAE,EAAE,EAAE,MAAM,EAAE,EAAE,EAAE,MAAM,KAAK,MAAM,CAAC;AAEnF,qFAAqF;AACrF,MAAM,MAAM,eAAe,GAAG,CAAC,CAAC,EAAE,IAAI,KAAK,iBAAiB,CAAC;AAE7D;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GAiDG;AACH,wBAAgB,qBAAqB,CACnC,SAAS,EAAE,SAAS,EACpB,QAAQ,EAAE,eAAe,EACzB,aAAa,EAAE,eAAe,EAC9B,CAAC,EAAE,IAAI,EACP,IAAI,GAAE,sBAA2B,GAChC,MAAM,CAmIR;AAED;;;;GAIG;AACH,MAAM,WAAW,oBAAoB;IACnC;wDACoD;IACpD,QAAQ,CAAC,QAAQ,EAAE,MAAM,CAAC;IAC1B;qFACiF;IACjF,QAAQ,CAAC,gBAAgB,EAAE,MAAM,CAAC;IAClC,uDAAuD;IACvD,QAAQ,CAAC,eAAe,EAAE,OAAO,CAAC;CACnC;AAED;;;;;;;;;;;;;;;;;;;GAmBG;AACH,wBAAgB,oBAAoB,CAClC,SAAS,EAAE,SAAS,EACpB,QAAQ,EAAE,eAAe,EACzB,aAAa,EAAE,eAAe,EAC9B,CAAC,EAAE,IAAI,EACP,IAAI,GAAE,sBAA2B,GAChC,oBAAoB,CAUtB;AAED;;;;;;;;;;;;;;;;;;;;;;;;;;;;GA4BG;AACH,wBAAgB,uBAAuB,CACrC,SAAS,EAAE,SAAS,EACpB,aAAa,EAAE,IAAI,EACnB,CAAC,EAAE,IAAI,GACN,MAAM,CAKR"}

package/dist/numerical/killing.js.map

@@ -1 +1 @@
-{"version":3,"file":"killing.js","sourceRoot":"","sources":["../../src/numerical/killing.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;;;;;;;GAyBG;AAEH,OAAO,EAAE,iBAAiB,EAAE,MAAM,aAAa,CAAC;AAChD,OAAO,EAAE,oBAAoB,EAAE,MAAM,uBAAuB,CAAC;AA4D7D;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GAiDG;AACH,MAAM,UAAU,qBAAqB,CACnC,SAAoB,EACpB,QAAkB,EAClB,aAA8B,EAC9B,CAAO,EACP,IAAI,GAA2B,EAAE;IAEjC,MAAM,KAAK,GAAG,IAAI,CAAC,KAAK,IAAI,CAAC,CAAC;IAC9B,MAAM,eAAe,GAAG,IAAI,CAAC,eAAe,IAAI,IAAI,CAAC;IACrD,MAAM,EAAE,GAAG,SAAS,EAAE,CAAC;IACvB,MAAM,CAAC,GAAG,QAAQ,CAAC,CAAC,CAAC,CAAC;IACtB,MAAM,KAAK,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC;IACpB,MAAM,KAAK,GAAG,aAAa,CAAC,CAAC,CAAC,CAAC;IAE/B,0CAA0C;IAC1C,MAAM,OAAO,GAAS,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;IACnC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;QAC3B,IAAI,CAAC,GAAG,CAAC,CAAC;QACV,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE;YAAE,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC;QACpD,OAAO,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;IACjB,CAAC;IAED,4CAA4C;IAC5C,MAAM,GAAG,GAAe;QACtB,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;QACZ,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;QACZ,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;QACZ,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;KACb,CAAC;IAEF,IAAI,eAAe,EAAE,CAAC;QACpB,sDAAsD;QACtD,qCAAqC;QACrC,4DAA4D;QAC5D,yDAAyD;QACzD,0CAA0C;QAC1C,EAAE;QACF,+EAA+E;QAC/E,8DAA8D;QAC9D,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,CAAC,EAAE,EAAE,EAAE,EAAE,CAAC;YAC9B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;gBAC3B,IAAI,CAAC,GAAG,CAAC,CAAC;gBACV,KAAK,IAAI,GAAG,GAAG,CAAC,EAAE,GAAG,GAAG,CAAC,EAAE,GAAG,EAAE,EAAE,CAAC;oBACjC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;wBAC3B,CAAC,IAAI,KAAK,CAAC,GAAG,EAAE,EAAE,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC;oBAChD,CAAC;gBACH,CAAC;gBACD,GAAG,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;YACjB,CAAC;QACH,CAAC;IACH,CAAC;SAAM,IAAI,IAAI,CAAC,SAAS,KAAK,SAAS,EAAE,CAAC;QACxC,kEAAkE;QAClE,8CAA8C;QAC9C,MAAM,EAAE,GAAG,IAAI,CAAC,SAAS,CAAC,CAAC,CAAC,CAAC;QAE7B,oCAAoC;QACpC,MAAM,QAAQ,GAAe;YAC3B,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;YACZ,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;YACZ,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;YACZ,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;SACb,CAAC;QACF,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;YAC3B,MAAM,QAAQ,GAAG,CAAC,CAAwB,EAAU,EAAE,CAAC,EAAE,CAAC,CAAS,CAAC,CAAC,CAAC,CAAC,CAAC;YACxE,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,CAAC,EAAE,EAAE,EAAE,EAAE,CAAC;gBAC9B,MAAM,GAAG,GAAG,iBAAiB,CAAC,QAAQ,EAAE,CAAC,EAAE,EAAE,EAAE,EAAE,KAAK,EAAE,CAAC,CAAC;gBAC1D,QAAQ,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,OAAO,GAAG,KAAK,QAAQ,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAE,GAAgB,CAAC,CAAC,CAAC,CAAC;YACzE,CAAC;QACH,CAAC;QAED,MAAM,OAAO,GAAe;YAC1B,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;YACZ,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;YACZ,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;YACZ,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;SACb,CAAC;QACF,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,CAAC,EAAE,EAAE,EAAE,EAAE,CAAC;YAC9B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;gBAC3B,IAAI,CAAC,GAAG,CAAC,CAAC;gBACV,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;oBAC3B,CAAC,IAAI,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC;oBAC7B,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,QAAQ,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;gBACjC,CAAC;gBACD,OAAO,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;YACrB,CAAC;QACH,CAAC;QAED,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,CAAC,EAAE,EAAE,EAAE,EAAE,CAAC;YAC9B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;gBAC3B,IAAI,OAAO,GAAG,CAAC,CAAC;gBAChB,KAAK,IAAI,GAAG,GAAG,CAAC,EAAE,GAAG,GAAG,CAAC,EAAE,GAAG,EAAE;oBAAE,OAAO,IAAI,KAAK,CAAC,GAAG,EAAE,EAAE,EAAE,CAAC,CAAC,GAAG,OAAO,CAAC,GAAG,CAAC,CAAC;gBAC9E,GAAG,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,OAAO,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,OAAO,CAAC;YACxC,CAAC;QACH,CAAC;IACH,CAAC;SAAM,CAAC;QACN,yEAAyE;QACzE,wEAAwE;QACxE,kEAAkE;QAClE,yCAAyC;QACzC,MAAM,OAAO,GAAe;YAC1B,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;YACZ,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;YACZ,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;YACZ,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;SACb,CAAC;QACF,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;YAC3B,MAAM,QAAQ,GAAG,CAAC,CAAwB,EAAU,EAAE;gBACpD,MAAM,EAAE,GAAG,QAAQ,CAAC,CAAC,CAAC,CAAC;gBACvB,MAAM,EAAE,GAAG,EAAE,CAAC,CAAS,CAAC,CAAC;gBACzB,IAAI,CAAC,GAAG,CAAC,CAAC;gBACV,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE;oBAAE,CAAC,IAAI,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC;gBAClD,OAAO,CAAC,CAAC;YACX,CAAC,CAAC;YACF,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,CAAC,EAAE,EAAE,EAAE,EAAE,CAAC;gBAC9B,MAAM,GAAG,GAAG,iBAAiB,CAAC,QAAQ,EAAE,CAAC,EAAE,EAAE,EAAE,EAAE,KAAK,EAAE,CAAC,CAAC;gBAC1D,OAAO,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,OAAO,GAAG,KAAK,QAAQ,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAE,GAAgB,CAAC,CAAC,CAAC,CAAC;YACxE,CAAC;QACH,CAAC;QACD,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,CAAC,EAAE,EAAE,EAAE,EAAE,CAAC;YAC9B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;gBAC3B,IAAI,OAAO,GAAG,CAAC,CAAC;gBAChB,KAAK,IAAI,GAAG,GAAG,CAAC,EAAE,GAAG,GAAG,CAAC,EAAE,GAAG,EAAE;oBAAE,OAAO,IAAI,KAAK,CAAC,GAAG,EAAE,EAAE,EAAE,CAAC,CAAC,GAAG,OAAO,CAAC,GAAG,CAAC,CAAC;gBAC9E,GAAG,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,OAAO,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,OAAO,CAAC;YACxC,CAAC;QACH,CAAC;IACH,CAAC;IAED,8DAA8D;IAC9D,IAAI,QAAQ,GAAG,CAAC,CAAC;IACjB,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,CAAC,EAAE,EAAE,EAAE,EAAE,CAAC;QAC9B,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,CAAC,EAAE,EAAE,EAAE,EAAE,CAAC;YAC9B,MAAM,GAAG,GAAG,GAAG,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,GAAG,GAAG,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC;YACtC,MAAM,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;YACzB,IAAI,EAAE,GAAG,QAAQ;gBAAE,QAAQ,GAAG,EAAE,CAAC;QACnC,CAAC;IACH,CAAC;IACD,OAAO,QAAQ,CAAC;AAClB,CAAC;AAkBD;;;;;;;;;;;;;;;;;;;GAmBG;AACH,MAAM,UAAU,oBAAoB,CAClC,SAAoB,EACpB,QAAkB,EAClB,aAA8B,EAC9B,CAAO,EACP,IAAI,GAA2B,EAAE;IAEjC,IAAI,IAAI,CAAC,SAAS,KAAK,SAAS,EAAE,CAAC;QACjC,oBAAoB,CAAC,IAAI,EAAE,CAAC,EAAE,IAAI,EAAE,WAAW,EAAE,GAAG,EAAE,CAAC,EAAE,UAAU,EAAE,IAAI,EAAE,CAAC,EAAE,sBAAsB,CAAC,CAAC;IACxG,CAAC;IACD,MAAM,SAAS,GAAG,IAAI,CAAC,SAAS,IAAI,KAAK,CAAC;IAC1C,MAAM,QAAQ,GAAG,qBAAqB,CAAC,SAAS,EAAE,QAAQ,EAAE,aAAa,EAAE,CAAC,EAAE,IAAI,CAAC,CAAC;IACpF,IAAI,KAAK,GAAG,CAAC,CAAC;IACd,KAAK,MAAM,GAAG,IAAI,QAAQ,CAAC,CAAC,CAAC;QAAE,KAAK,MAAM,CAAC,IAAI,GAAG;YAAE,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,EAAE,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;IACzF,MAAM,gBAAgB,GAAG,QAAQ,GAAG,KAAK,CAAC;IAC1C,OAAO,EAAE,QAAQ,EAAE,gBAAgB,EAAE,eAAe,EAAE,gBAAgB,IAAI,SAAS,EAAE,CAAC;AACxF,CAAC;AAED;;;;;;;;;;;;;;;;;;;;;;;;;;;;GA4BG;AACH,MAAM,UAAU,uBAAuB,CACrC,SAAoB,EACpB,aAAmB,EACnB,CAAO;IAEP,MAAM,EAAE,GAAG,SAAS,EAAE,CAAC,CAAC,CAAC,CAAC;IAC1B,IAAI,CAAC,GAAG,CAAC,CAAC;IACV,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE;QAAE,CAAC,IAAI,EAAE,CAAC,CAAC,CAAC,GAAG,aAAa,CAAC,CAAC,CAAC,CAAC;IAC1D,OAAO,CAAC,CAAC;AACX,CAAC"}
+{"version":3,"file":"killing.js","sourceRoot":"","sources":["../../src/numerical/killing.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;;;;;;;GAyBG;AAEH,OAAO,EAAE,iBAAiB,EAAE,MAAM,aAAa,CAAC;AAChD,OAAO,EAAE,oBAAoB,EAAE,MAAM,uBAAuB,CAAC;AA+D7D;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GAiDG;AACH,MAAM,UAAU,qBAAqB,CACnC,SAAoB,EACpB,QAAyB,EACzB,aAA8B,EAC9B,CAAO,EACP,IAAI,GAA2B,EAAE;IAEjC,MAAM,KAAK,GAAG,IAAI,CAAC,KAAK,IAAI,CAAC,CAAC;IAC9B,MAAM,eAAe,GAAG,IAAI,CAAC,eAAe,IAAI,IAAI,CAAC;IACrD,MAAM,EAAE,GAAG,SAAS,EAAE,CAAC;IACvB,MAAM,CAAC,GAAG,QAAQ,CAAC,CAAC,CAAC,CAAC;IACtB,MAAM,KAAK,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC;IACpB,MAAM,KAAK,GAAG,aAAa,CAAC,CAAC,CAAC,CAAC;IAE/B,0CAA0C;IAC1C,MAAM,OAAO,GAAS,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;IACnC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;QAC3B,IAAI,CAAC,GAAG,CAAC,CAAC;QACV,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE;YAAE,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC;QACpD,OAAO,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;IACjB,CAAC;IAED,4CAA4C;IAC5C,MAAM,GAAG,GAAe;QACtB,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;QACZ,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;QACZ,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;QACZ,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;KACb,CAAC;IAEF,IAAI,eAAe,EAAE,CAAC;QACpB,sDAAsD;QACtD,qCAAqC;QACrC,4DAA4D;QAC5D,yDAAyD;QACzD,0CAA0C;QAC1C,EAAE;QACF,+EAA+E;QAC/E,8DAA8D;QAC9D,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,CAAC,EAAE,EAAE,EAAE,EAAE,CAAC;YAC9B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;gBAC3B,IAAI,CAAC,GAAG,CAAC,CAAC;gBACV,KAAK,IAAI,GAAG,GAAG,CAAC,EAAE,GAAG,GAAG,CAAC,EAAE,GAAG,EAAE,EAAE,CAAC;oBACjC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;wBAC3B,CAAC,IAAI,KAAK,CAAC,GAAG,EAAE,EAAE,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC;oBAChD,CAAC;gBACH,CAAC;gBACD,GAAG,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;YACjB,CAAC;QACH,CAAC;IACH,CAAC;SAAM,IAAI,IAAI,CAAC,SAAS,KAAK,SAAS,EAAE,CAAC;QACxC,kEAAkE;QAClE,8CAA8C;QAC9C,MAAM,EAAE,GAAG,IAAI,CAAC,SAAS,CAAC,CAAC,CAAC,CAAC;QAE7B,oCAAoC;QACpC,MAAM,QAAQ,GAAe;YAC3B,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;YACZ,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;YACZ,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;YACZ,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;SACb,CAAC;QACF,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;YAC3B,MAAM,QAAQ,GAAG,CAAC,CAAwB,EAAU,EAAE,CAAC,EAAE,CAAC,CAAS,CAAC,CAAC,CAAC,CAAC,CAAC;YACxE,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,CAAC,EAAE,EAAE,EAAE,EAAE,CAAC;gBAC9B,MAAM,GAAG,GAAG,iBAAiB,CAAC,QAAQ,EAAE,CAAC,EAAE,EAAE,EAAE,EAAE,KAAK,EAAE,CAAC,CAAC;gBAC1D,QAAQ,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,OAAO,GAAG,KAAK,QAAQ,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAE,GAAgB,CAAC,CAAC,CAAC,CAAC;YACzE,CAAC;QACH,CAAC;QAED,MAAM,OAAO,GAAe;YAC1B,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;YACZ,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;YACZ,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;YACZ,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;SACb,CAAC;QACF,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,CAAC,EAAE,EAAE,EAAE,EAAE,CAAC;YAC9B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;gBAC3B,IAAI,CAAC,GAAG,CAAC,CAAC;gBACV,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;oBAC3B,CAAC,IAAI,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC;oBAC7B,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,QAAQ,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;gBACjC,CAAC;gBACD,OAAO,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;YACrB,CAAC;QACH,CAAC;QAED,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,CAAC,EAAE,EAAE,EAAE,EAAE,CAAC;YAC9B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;gBAC3B,IAAI,OAAO,GAAG,CAAC,CAAC;gBAChB,KAAK,IAAI,GAAG,GAAG,CAAC,EAAE,GAAG,GAAG,CAAC,EAAE,GAAG,EAAE;oBAAE,OAAO,IAAI,KAAK,CAAC,GAAG,EAAE,EAAE,EAAE,CAAC,CAAC,GAAG,OAAO,CAAC,GAAG,CAAC,CAAC;gBAC9E,GAAG,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,OAAO,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,OAAO,CAAC;YACxC,CAAC;QACH,CAAC;IACH,CAAC;SAAM,CAAC;QACN,yEAAyE;QACzE,wEAAwE;QACxE,kEAAkE;QAClE,yCAAyC;QACzC,MAAM,OAAO,GAAe;YAC1B,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;YACZ,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;YACZ,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;YACZ,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;SACb,CAAC;QACF,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;YAC3B,MAAM,QAAQ,GAAG,CAAC,CAAwB,EAAU,EAAE;gBACpD,MAAM,EAAE,GAAG,QAAQ,CAAC,CAAC,CAAC,CAAC;gBACvB,MAAM,EAAE,GAAG,EAAE,CAAC,CAAS,CAAC,CAAC;gBACzB,IAAI,CAAC,GAAG,CAAC,CAAC;gBACV,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE;oBAAE,CAAC,IAAI,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC;gBAClD,OAAO,CAAC,CAAC;YACX,CAAC,CAAC;YACF,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,CAAC,EAAE,EAAE,EAAE,EAAE,CAAC;gBAC9B,MAAM,GAAG,GAAG,iBAAiB,CAAC,QAAQ,EAAE,CAAC,EAAE,EAAE,EAAE,EAAE,KAAK,EAAE,CAAC,CAAC;gBAC1D,OAAO,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,OAAO,GAAG,KAAK,QAAQ,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAE,GAAgB,CAAC,CAAC,CAAC,CAAC;YACxE,CAAC;QACH,CAAC;QACD,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,CAAC,EAAE,EAAE,EAAE,EAAE,CAAC;YAC9B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;gBAC3B,IAAI,OAAO,GAAG,CAAC,CAAC;gBAChB,KAAK,IAAI,GAAG,GAAG,CAAC,EAAE,GAAG,GAAG,CAAC,EAAE,GAAG,EAAE;oBAAE,OAAO,IAAI,KAAK,CAAC,GAAG,EAAE,EAAE,EAAE,CAAC,CAAC,GAAG,OAAO,CAAC,GAAG,CAAC,CAAC;gBAC9E,GAAG,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,OAAO,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,OAAO,CAAC;YACxC,CAAC;QACH,CAAC;IACH,CAAC;IAED,8DAA8D;IAC9D,IAAI,QAAQ,GAAG,CAAC,CAAC;IACjB,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,CAAC,EAAE,EAAE,EAAE,EAAE,CAAC;QAC9B,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,CAAC,EAAE,EAAE,EAAE,EAAE,CAAC;YAC9B,MAAM,GAAG,GAAG,GAAG,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,GAAG,GAAG,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC;YACtC,MAAM,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;YACzB,IAAI,EAAE,GAAG,QAAQ;gBAAE,QAAQ,GAAG,EAAE,CAAC;QACnC,CAAC;IACH,CAAC;IACD,OAAO,QAAQ,CAAC;AAClB,CAAC;AAkBD;;;;;;;;;;;;;;;;;;;GAmBG;AACH,MAAM,UAAU,oBAAoB,CAClC,SAAoB,EACpB,QAAyB,EACzB,aAA8B,EAC9B,CAAO,EACP,IAAI,GAA2B,EAAE;IAEjC,IAAI,IAAI,CAAC,SAAS,KAAK,SAAS,EAAE,CAAC;QACjC,oBAAoB,CAAC,IAAI,EAAE,CAAC,EAAE,IAAI,EAAE,WAAW,EAAE,GAAG,EAAE,CAAC,EAAE,UAAU,EAAE,IAAI,EAAE,CAAC,EAAE,sBAAsB,CAAC,CAAC;IACxG,CAAC;IACD,MAAM,SAAS,GAAG,IAAI,CAAC,SAAS,IAAI,KAAK,CAAC;IAC1C,MAAM,QAAQ,GAAG,qBAAqB,CAAC,SAAS,EAAE,QAAQ,EAAE,aAAa,EAAE,CAAC,EAAE,IAAI,CAAC,CAAC;IACpF,IAAI,KAAK,GAAG,CAAC,CAAC;IACd,KAAK,MAAM,GAAG,IAAI,QAAQ,CAAC,CAAC,CAAC;QAAE,KAAK,MAAM,CAAC,IAAI,GAAG;YAAE,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,EAAE,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;IACzF,MAAM,gBAAgB,GAAG,QAAQ,GAAG,KAAK,CAAC;IAC1C,OAAO,EAAE,QAAQ,EAAE,gBAAgB,EAAE,eAAe,EAAE,gBAAgB,IAAI,SAAS,EAAE,CAAC;AACxF,CAAC;AAED;;;;;;;;;;;;;;;;;;;;;;;;;;;;GA4BG;AACH,MAAM,UAAU,uBAAuB,CACrC,SAAoB,EACpB,aAAmB,EACnB,CAAO;IAEP,MAAM,EAAE,GAAG,SAAS,EAAE,CAAC,CAAC,CAAC,CAAC;IAC1B,IAAI,CAAC,GAAG,CAAC,CAAC;IACV,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE;QAAE,CAAC,IAAI,EAAE,CAAC,CAAC,CAAC,GAAG,aAAa,CAAC,CAAC,CAAC,CAAC;IAC1D,OAAO,CAAC,CAAC;AACX,CAAC","sourcesContent":["/**\n * Killing-equation numerical verification (v0.6.0 Phase 1, Task 1.3).\n *\n * F-3 reconciliation (hybrid impl):\n *   - Uses EXACT symbolic Christoffels for all Γ terms.\n *   - For the outer ∂_μ ξ_α = (∂_μ g_{αβ}) ξ^β + g_{αβ} ∂_μ ξ^β, two paths:\n *     1. constantKilling=true (default for constant ξ^β): uses metric\n *        compatibility ∂_μ g_{αβ} = Γ^λ_{μα} g_{λβ} + Γ^λ_{μβ} g_{αλ} to\n *        derive ∇_μ ξ_α = Γ^λ_{μβ} g_{αλ} ξ^β without any FD. Achieves\n *        machine precision RELATIVE to the metric scale for Schwarzschild.\n *        The residual itself is absolute, in the metric's units: with SI\n *        Schwarzschild (g_tt ≈ −c² ≈ −9e16) the exact time-translation field\n *        leaves 2.44e-4 at r = 3 r_s and 3.05e-5 at 10 r_s (≈ 3e-21 relative).\n *     2. constantKilling=false with dMetricFn supplied: uses exact analytic\n *        ∂_μ g_{αβ} + FD for ∂_μ ξ^β. Handles spatially-varying Killing fields.\n *     3. Fallback FD path: differentiates ξ_α(y) = g_{αβ}(y) ξ^β(y) via 4th-\n *        order FD. WARNING: suffers c²-cancellation noise for SI-scaled metrics.\n *\n * The Killing equation (Carroll Eq. 3.174):\n *   ∇_μ ξ_ν + ∇_ν ξ_μ = 0\n *\n * where ∇_μ ξ_ν = ∂_μ ξ_ν − Γ^λ_{μν} ξ_λ (covariant derivative of the\n * lowered Killing field ξ_ν = g_{νβ} ξ^β).\n *\n * @module numerical/killing\n */\n\nimport { pderivNumericalFn } from './pderiv.js';\nimport { validateFiniteInputs } from './input-validation.js';\n\n/** Options for {@link verifyKillingEquation}.\n *\n *  @public */\nexport interface KillingEquationOptions {\n  /**\n   * Maximum tolerated RELATIVE residual, read only by {@link checkKillingEquation}:\n   * ||∇_μ ξ_ν + ∇_ν ξ_μ||_∞ / max(max_{μν} |g_{μν}(x)|, 1) ≤ tolerance. Default 1e-10.\n   * Must be finite and > 0. {@link verifyKillingEquation} ignores it and returns the\n   * raw (absolute) residual.\n   */\n  readonly tolerance?: number;\n  /**\n   * Finite-difference stencil order for the outer ∂_μ ξ_α partials\n   * (used only when constantKilling=false and dMetricFn is not supplied).\n   * Default 4.\n   */\n  readonly order?: 2 | 4;\n  /**\n   * When true (default), assumes ξ^β is coordinate-independent (∂_μ ξ^β = 0)\n   * and uses the metric-compatibility identity to compute the covariant\n   * derivative purely algebraically:\n   *\n   *   ∇_μ ξ_α = Γ^λ_{μβ} g_{αλ} ξ^β  (for constant ξ^β)\n   *\n   * This avoids all FD and any c²-scale cancellation noise. Exactly correct for\n   * the Schwarzschild time-translation (ξ^t = 1) and axial (ξ^φ = 1) fields.\n   * Set to false for spatially-varying Killing fields.\n   */\n  readonly constantKilling?: boolean;\n  /**\n   * Optional exact partial derivatives of the covariant metric\n   * (used only when constantKilling=false).\n   * dMetricFn(x)[lambda][mu][nu] = ∂_lambda g_{mu nu}.\n   * When supplied with constantKilling=false, ∂_μ g_{αβ} is exact;\n   * ∂_μ ξ^β is computed via 4th-order FD.\n   */\n  readonly dMetricFn?: (x: ReadonlyArray<number>) => number[][][];\n}\n\ntype Vec4 = [number, number, number, number];\n/** A Killing-field factory: returns the field ξ^μ(x) as a function of the 4-position. @public */\nexport type KillingFn = () => (x: Vec4) => Vec4;\n/** A covariant metric closure for the Killing checks: `g(x)[μ][ν] = g_{μν}(x)`, nested. @public */\nexport type KillingMetricFn = (x: ReadonlyArray<number>) => number[][];\n\n/**\n * Layout-agnostic Christoffel accessor.\n *\n * P-3 fix: both nested (number[4][4][4]) and flat (Float64Array(64))\n * Christoffel producers can provide a ChristoffelAccess. This insulates\n * {@link verifyKillingEquation} from the Phase 2 BR-2 migration:\n * Task 2.9 changes how christoffelAt(x) is *constructed*, but consumers\n * only call `christoffelAt(x)(λ, μ, ν)`.\n *\n * @public\n */\nexport type ChristoffelAccess = (lambda: number, mu: number, nu: number) => number;\n\n/** A Christoffel producer: the {@link ChristoffelAccess} at a 4-position. @public */\nexport type ChristoffelAtFn = (x: Vec4) => ChristoffelAccess;\n\n/**\n * Compute `||∇_μ ξ_ν + ∇_ν ξ_μ||_∞` at point `x`, returning the maximum\n * absolute symmetrized residual over all 16 (μ, ν) pairs.\n *\n * For constant Killing fields (constantKilling=true, the default):\n *   Uses metric compatibility to derive ∇_μ ξ_α = Γ^λ_{μβ} g_{αλ} ξ^β\n *   algebraically — no FD at all. The residual is ABSOLUTE (metric units);\n *   for SI Schwarzschild it is at machine precision relative to |g_tt| ≈ c².\n *   `opts.tolerance` is not read here; use {@link checkKillingEquation} for a verdict.\n *\n * @param killingFn     - Factory returning ξ^μ(x) (upper index).\n * @param metricFn      - Covariant metric g_{μν}(x).\n * @param christoffelAt - Christoffel accessor factory (layout-agnostic).\n * @param x             - Evaluation point [t, r, θ, φ].\n * @param opts          - Optional overrides (see {@link KillingEquationOptions}).\n * @returns             Maximum absolute residual of the symmetrized Killing equation.\n *\n * @example\n * ```typescript\n * import { verifyKillingEquation } from 'universal-physics-tensor';\n * import {\n *   schwarzschildKillingT,\n *   schwarzschildGFn,\n *   schwarzschildChristoffelFn,\n *   schwarzschildRs,\n * } from '../tests/fixtures/schwarzschild.js'; // or your own metric closures\n *\n * const M = 1.989e30; // solar mass, kg\n * const r_s = schwarzschildRs(M);\n * const x: [number, number, number, number] = [0, 5 * r_s, Math.PI / 2, 0];\n *\n * // Christoffel accessor wrapping the flat Float64Array layout (BR-2)\n * const gammaFn = schwarzschildChristoffelFn(M);\n * const christoffelAccessAt = (pt: [number, number, number, number]) => {\n *   const arr = gammaFn(pt);\n *   return (lam: number, mu: number, nu: number) => arr[16 * lam + 4 * mu + nu];\n * };\n *\n * const residual = verifyKillingEquation(\n *   schwarzschildKillingT,\n *   schwarzschildGFn(M),\n *   christoffelAccessAt,\n *   x,\n *   { constantKilling: true },\n * );\n * // residual < 1e-12 for the exact time-translation Killing field\n * ```\n *\n * @public\n */\nexport function verifyKillingEquation(\n  killingFn: KillingFn,\n  metricFn: KillingMetricFn,\n  christoffelAt: ChristoffelAtFn,\n  x: Vec4,\n  opts: KillingEquationOptions = {},\n): number {\n  const order = opts.order ?? 4;\n  const constantKilling = opts.constantKilling ?? true;\n  const xi = killingFn();\n  const g = metricFn(x);\n  const xiAtX = xi(x);\n  const Gamma = christoffelAt(x);\n\n  // Lowered Killing field: ξ_α = g_{αβ} ξ^β\n  const xiLower: Vec4 = [0, 0, 0, 0];\n  for (let a = 0; a < 4; a++) {\n    let s = 0;\n    for (let b = 0; b < 4; b++) s += g[a][b] * xiAtX[b];\n    xiLower[a] = s;\n  }\n\n  // Covariant derivative ∇_μ ξ_α (4×4 matrix)\n  const cov: number[][] = [\n    [0, 0, 0, 0],\n    [0, 0, 0, 0],\n    [0, 0, 0, 0],\n    [0, 0, 0, 0],\n  ];\n\n  if (constantKilling) {\n    // Algebraic path (metric compatibility, ∂_μ ξ^β = 0):\n    //   ∇_μ ξ_α = ∂_μ ξ_α − Γ^λ_{μα} ξ_λ\n    //   where ∂_μ ξ_α = (Γ^λ_{μα} g_{λβ} + Γ^λ_{μβ} g_{αλ}) ξ^β\n    //   cancelling Γ^λ_{μα} g_{λβ} ξ^β = Γ^λ_{μα} ξ_λ gives:\n    //   ∇_μ ξ_α = Σ_β Σ_λ Γ^λ_{μβ} g_{αλ} ξ^β\n    //\n    // Symmetrized: ∇_μ ξ_α + ∇_α ξ_μ = Σ_β [Γ^λ_{μβ} g_{αλ} + Γ^λ_{αβ} g_{μλ}] ξ^β\n    // This is the Killing equation in terms of Christoffels only.\n    for (let mu = 0; mu < 4; mu++) {\n      for (let a = 0; a < 4; a++) {\n        let s = 0;\n        for (let lam = 0; lam < 4; lam++) {\n          for (let b = 0; b < 4; b++) {\n            s += Gamma(lam, mu, b) * g[a][lam] * xiAtX[b];\n          }\n        }\n        cov[mu][a] = s;\n      }\n    }\n  } else if (opts.dMetricFn !== undefined) {\n    // Exact dMetricFn + FD path for spatially-varying Killing fields.\n    // ∂_μ ξ_α = (∂_μ g_{αβ}) ξ^β + g_{αβ} ∂_μ ξ^β\n    const dg = opts.dMetricFn(x);\n\n    // FD for ∂_μ ξ^β (upper components)\n    const dxiUpper: number[][] = [\n      [0, 0, 0, 0],\n      [0, 0, 0, 0],\n      [0, 0, 0, 0],\n      [0, 0, 0, 0],\n    ];\n    for (let b = 0; b < 4; b++) {\n      const xiUpperB = (y: ReadonlyArray<number>): number => xi(y as Vec4)[b];\n      for (let mu = 0; mu < 4; mu++) {\n        const raw = pderivNumericalFn(xiUpperB, x, mu, { order });\n        dxiUpper[mu][b] = typeof raw === 'number' ? raw : (raw as number[])[0];\n      }\n    }\n\n    const partial: number[][] = [\n      [0, 0, 0, 0],\n      [0, 0, 0, 0],\n      [0, 0, 0, 0],\n      [0, 0, 0, 0],\n    ];\n    for (let mu = 0; mu < 4; mu++) {\n      for (let a = 0; a < 4; a++) {\n        let s = 0;\n        for (let b = 0; b < 4; b++) {\n          s += dg[mu][a][b] * xiAtX[b];\n          s += g[a][b] * dxiUpper[mu][b];\n        }\n        partial[mu][a] = s;\n      }\n    }\n\n    for (let mu = 0; mu < 4; mu++) {\n      for (let a = 0; a < 4; a++) {\n        let connSum = 0;\n        for (let lam = 0; lam < 4; lam++) connSum += Gamma(lam, mu, a) * xiLower[lam];\n        cov[mu][a] = partial[mu][a] - connSum;\n      }\n    }\n  } else {\n    // Pure FD fallback: differentiate ξ_α(y) = g_{αβ}(y) ξ^β(y) numerically.\n    // WARNING: large-scale SI metrics (c²·g_tt ~ 6e16) produce cancellation\n    // noise at the ~1-10 absolute level. Use constantKilling=true for\n    // constant Schwarzschild Killing fields.\n    const partial: number[][] = [\n      [0, 0, 0, 0],\n      [0, 0, 0, 0],\n      [0, 0, 0, 0],\n      [0, 0, 0, 0],\n    ];\n    for (let a = 0; a < 4; a++) {\n      const xiLowerA = (y: ReadonlyArray<number>): number => {\n        const gy = metricFn(y);\n        const xy = xi(y as Vec4);\n        let s = 0;\n        for (let b = 0; b < 4; b++) s += gy[a][b] * xy[b];\n        return s;\n      };\n      for (let mu = 0; mu < 4; mu++) {\n        const raw = pderivNumericalFn(xiLowerA, x, mu, { order });\n        partial[mu][a] = typeof raw === 'number' ? raw : (raw as number[])[0];\n      }\n    }\n    for (let mu = 0; mu < 4; mu++) {\n      for (let a = 0; a < 4; a++) {\n        let connSum = 0;\n        for (let lam = 0; lam < 4; lam++) connSum += Gamma(lam, mu, a) * xiLower[lam];\n        cov[mu][a] = partial[mu][a] - connSum;\n      }\n    }\n  }\n\n  // Symmetrize ∇_μ ξ_ν + ∇_ν ξ_μ and find max|·| over 16 pairs.\n  let maxResid = 0;\n  for (let mu = 0; mu < 4; mu++) {\n    for (let nu = 0; nu < 4; nu++) {\n      const sym = cov[mu][nu] + cov[nu][mu];\n      const ab = Math.abs(sym);\n      if (ab > maxResid) maxResid = ab;\n    }\n  }\n  return maxResid;\n}\n\n/**\n * Result of {@link checkKillingEquation}.\n *\n * @public\n */\nexport interface KillingEquationCheck {\n  /** Raw max residual ||∇_μ ξ_ν + ∇_ν ξ_μ||_∞ in the metric's units — identical to\n   *  {@link verifyKillingEquation}'s return value. */\n  readonly residual: number;\n  /** `residual / max(max_{μν} |g_{μν}(x)|, 1)`: dimensionless, the quantity the tolerance\n   *  is applied to (the normalization `evaluateEinsteinEquationResidual` uses). */\n  readonly relativeResidual: number;\n  /** `relativeResidual <= tolerance` (default 1e-10). */\n  readonly withinTolerance: boolean;\n}\n\n/**\n * Check the Killing equation ∇_μ ξ_ν + ∇_ν ξ_μ = 0 at `x` against a tolerance.\n *\n * Runs {@link verifyKillingEquation} with the same arguments and normalizes its absolute\n * residual by the metric scale at `x`: `relativeResidual = residual / max(max|g_μν(x)|, 1)`.\n * The floor of 1 keeps geometrized (order-1) metrics absolute. The verdict is\n * `relativeResidual <= opts.tolerance` (default 1e-10). A relative test matters for SI\n * metrics: an exact Killing field of SI Schwarzschild leaves an absolute residual up to\n * 2.44e-4 (g_tt ≈ −9e16), about 3e-21 relative.\n *\n * @param killingFn     - Factory returning ξ^μ(x) (upper index).\n * @param metricFn      - Covariant metric g_{μν}(x).\n * @param christoffelAt - Christoffel accessor factory (layout-agnostic).\n * @param x             - Evaluation point [t, r, θ, φ].\n * @param opts          - {@link KillingEquationOptions}; `tolerance` must be finite and > 0.\n * @returns `{ residual, relativeResidual, withinTolerance }`.\n * @throws RangeError when `opts.tolerance` is not a finite positive number.\n *\n * @public\n */\nexport function checkKillingEquation(\n  killingFn: KillingFn,\n  metricFn: KillingMetricFn,\n  christoffelAt: ChristoffelAtFn,\n  x: Vec4,\n  opts: KillingEquationOptions = {},\n): KillingEquationCheck {\n  if (opts.tolerance !== undefined) {\n    validateFiniteInputs(opts, [{ name: 'tolerance', min: 0, excludeMin: true }], 'checkKillingEquation');\n  }\n  const tolerance = opts.tolerance ?? 1e-10;\n  const residual = verifyKillingEquation(killingFn, metricFn, christoffelAt, x, opts);\n  let scale = 1;\n  for (const row of metricFn(x)) for (const v of row) scale = Math.max(scale, Math.abs(v));\n  const relativeResidual = residual / scale;\n  return { residual, relativeResidual, withinTolerance: relativeResidual <= tolerance };\n}\n\n/**\n * Evaluate the conserved charge Q = ξ^μ p_μ at point x along a geodesic.\n *\n * Sign convention: raw contraction ξ^μ p_μ. For the time-translation\n * Killing field in (−,+,+,+) signature, p_t < 0 for forward-time motion,\n * so the physical energy of common usage is E = −Q (Carroll Eq. 8.30).\n * Use `|ΔQ / Q|` for drift diagnostics to stay sign-agnostic.\n *\n * @param killingFn     - Factory returning ξ^μ(x) (upper index).\n * @param momentumLower - Co-momentum p_μ at the geodesic point (lower index).\n * @param x             - Geodesic coordinate [t, r, θ, φ].\n * @returns             Q = ξ^μ p_μ (dimensioned as energy if p is 4-momentum).\n *\n * @example\n * ```typescript\n * import { evaluateConservedCharge } from 'universal-physics-tensor';\n * import { schwarzschildKillingT } from '../tests/fixtures/schwarzschild.js';\n *\n * // Co-momentum p_μ at a geodesic snapshot (p_t < 0 for forward-time motion)\n * const pLower: [number, number, number, number] = [-8.988e16, 0, 0, 0];\n * const x: [number, number, number, number] = [0, 1e7, Math.PI / 2, 0];\n *\n * const Q = evaluateConservedCharge(schwarzschildKillingT, pLower, x);\n * // Q = ξ^t · p_t = 1 · p_t = p_t (energy with sign from metric)\n * // Physical energy E = −Q (Carroll sign convention for (−,+,+,+))\n * ```\n *\n * @public\n */\nexport function evaluateConservedCharge(\n  killingFn: KillingFn,\n  momentumLower: Vec4,\n  x: Vec4,\n): number {\n  const xi = killingFn()(x);\n  let q = 0;\n  for (let m = 0; m < 4; m++) q += xi[m] * momentumLower[m];\n  return q;\n}\n"]}