universal-physics-tensor 0.5.0 → 0.7.0

package/dist/numerical/killing.js

@@ -0,0 +1,242 @@
+/**
+ * Killing-equation numerical verification (v0.6.0 Phase 1, Task 1.3).
+ *
+ * F-3 reconciliation (hybrid impl):
+ *   - Uses EXACT symbolic Christoffels for all Γ terms.
+ *   - For the outer ∂_μ ξ_α = (∂_μ g_{αβ}) ξ^β + g_{αβ} ∂_μ ξ^β, two paths:
+ *     1. constantKilling=true (default for constant ξ^β): uses metric
+ *        compatibility ∂_μ g_{αβ} = Γ^λ_{μα} g_{λβ} + Γ^λ_{μβ} g_{αλ} to
+ *        derive ∇_μ ξ_α = Γ^λ_{μβ} g_{αλ} ξ^β without any FD. Achieves
+ *        machine-precision (~1e-15 or exact 0) for Schwarzschild.
+ *     2. constantKilling=false with dMetricFn supplied: uses exact analytic
+ *        ∂_μ g_{αβ} + FD for ∂_μ ξ^β. Handles spatially-varying Killing fields.
+ *     3. Fallback FD path: differentiates ξ_α(y) = g_{αβ}(y) ξ^β(y) via 4th-
+ *        order FD. WARNING: suffers c²-cancellation noise for SI-scaled metrics.
+ *
+ * The Killing equation (Carroll Eq. 3.174):
+ *   ∇_μ ξ_ν + ∇_ν ξ_μ = 0
+ *
+ * where ∇_μ ξ_ν = ∂_μ ξ_ν − Γ^λ_{μν} ξ_λ (covariant derivative of the
+ * lowered Killing field ξ_ν = g_{νβ} ξ^β).
+ *
+ * @module numerical/killing
+ */
+import { pderivNumericalFn } from './pderiv.js';
+/**
+ * Compute `||∇_μ ξ_ν + ∇_ν ξ_μ||_∞` at point `x`, returning the maximum
+ * absolute symmetrized residual over all 16 (μ, ν) pairs.
+ *
+ * For constant Killing fields (constantKilling=true, the default):
+ *   Uses metric compatibility to derive ∇_μ ξ_α = Γ^λ_{μβ} g_{αλ} ξ^β
+ *   algebraically — no FD at all. Achieves machine-precision residuals
+ *   for Schwarzschild Killing fields in SI units.
+ *
+ * @param killingFn     - Factory returning ξ^μ(x) (upper index).
+ * @param metricFn      - Covariant metric g_{μν}(x).
+ * @param christoffelAt - Christoffel accessor factory (layout-agnostic).
+ * @param x             - Evaluation point [t, r, θ, φ].
+ * @param opts          - Optional overrides (see {@link KillingEquationOptions}).
+ * @returns             Maximum absolute residual of the symmetrized Killing equation.
+ *
+ * @example
+ * ```typescript
+ * import { verifyKillingEquation } from 'universal-physics-tensor';
+ * import {
+ *   schwarzschildKillingT,
+ *   schwarzschildGFn,
+ *   schwarzschildChristoffelFn,
+ *   schwarzschildRs,
+ * } from '../tests/fixtures/schwarzschild.js'; // or your own metric closures
+ *
+ * const M = 1.989e30; // solar mass, kg
+ * const r_s = schwarzschildRs(M);
+ * const x: [number, number, number, number] = [0, 5 * r_s, Math.PI / 2, 0];
+ *
+ * // Christoffel accessor wrapping the flat Float64Array layout (BR-2)
+ * const gammaFn = schwarzschildChristoffelFn(M);
+ * const christoffelAccessAt = (pt: [number, number, number, number]) => {
+ *   const arr = gammaFn(pt);
+ *   return (lam: number, mu: number, nu: number) => arr[16 * lam + 4 * mu + nu];
+ * };
+ *
+ * const residual = verifyKillingEquation(
+ *   schwarzschildKillingT,
+ *   schwarzschildGFn(M),
+ *   christoffelAccessAt,
+ *   x,
+ *   { constantKilling: true },
+ * );
+ * // residual < 1e-12 for the exact time-translation Killing field
+ * ```
+ *
+ * @public
+ */
+export function verifyKillingEquation(killingFn, metricFn, christoffelAt, x, opts = {}) {
+    const order = opts.order ?? 4;
+    const constantKilling = opts.constantKilling ?? true;
+    const xi = killingFn();
+    const g = metricFn(x);
+    const xiAtX = xi(x);
+    const Gamma = christoffelAt(x);
+    // Lowered Killing field: ξ_α = g_{αβ} ξ^β
+    const xiLower = [0, 0, 0, 0];
+    for (let a = 0; a < 4; a++) {
+        let s = 0;
+        for (let b = 0; b < 4; b++)
+            s += g[a][b] * xiAtX[b];
+        xiLower[a] = s;
+    }
+    // Covariant derivative ∇_μ ξ_α (4×4 matrix)
+    const cov = [
+        [0, 0, 0, 0],
+        [0, 0, 0, 0],
+        [0, 0, 0, 0],
+        [0, 0, 0, 0],
+    ];
+    if (constantKilling) {
+        // Algebraic path (metric compatibility, ∂_μ ξ^β = 0):
+        //   ∇_μ ξ_α = ∂_μ ξ_α − Γ^λ_{μα} ξ_λ
+        //   where ∂_μ ξ_α = (Γ^λ_{μα} g_{λβ} + Γ^λ_{μβ} g_{αλ}) ξ^β
+        //   cancelling Γ^λ_{μα} g_{λβ} ξ^β = Γ^λ_{μα} ξ_λ gives:
+        //   ∇_μ ξ_α = Σ_β Σ_λ Γ^λ_{μβ} g_{αλ} ξ^β
+        //
+        // Symmetrized: ∇_μ ξ_α + ∇_α ξ_μ = Σ_β [Γ^λ_{μβ} g_{αλ} + Γ^λ_{αβ} g_{μλ}] ξ^β
+        // This is the Killing equation in terms of Christoffels only.
+        for (let mu = 0; mu < 4; mu++) {
+            for (let a = 0; a < 4; a++) {
+                let s = 0;
+                for (let lam = 0; lam < 4; lam++) {
+                    for (let b = 0; b < 4; b++) {
+                        s += Gamma(lam, mu, b) * g[a][lam] * xiAtX[b];
+                    }
+                }
+                cov[mu][a] = s;
+            }
+        }
+    }
+    else if (opts.dMetricFn !== undefined) {
+        // Exact dMetricFn + FD path for spatially-varying Killing fields.
+        // ∂_μ ξ_α = (∂_μ g_{αβ}) ξ^β + g_{αβ} ∂_μ ξ^β
+        const dg = opts.dMetricFn(x);
+        // FD for ∂_μ ξ^β (upper components)
+        const dxiUpper = [
+            [0, 0, 0, 0],
+            [0, 0, 0, 0],
+            [0, 0, 0, 0],
+            [0, 0, 0, 0],
+        ];
+        for (let b = 0; b < 4; b++) {
+            const xiUpperB = (y) => xi(y)[b];
+            for (let mu = 0; mu < 4; mu++) {
+                const raw = pderivNumericalFn(xiUpperB, x, mu, { order });
+                dxiUpper[mu][b] = typeof raw === 'number' ? raw : raw[0];
+            }
+        }
+        const partial = [
+            [0, 0, 0, 0],
+            [0, 0, 0, 0],
+            [0, 0, 0, 0],
+            [0, 0, 0, 0],
+        ];
+        for (let mu = 0; mu < 4; mu++) {
+            for (let a = 0; a < 4; a++) {
+                let s = 0;
+                for (let b = 0; b < 4; b++) {
+                    s += dg[mu][a][b] * xiAtX[b];
+                    s += g[a][b] * dxiUpper[mu][b];
+                }
+                partial[mu][a] = s;
+            }
+        }
+        for (let mu = 0; mu < 4; mu++) {
+            for (let a = 0; a < 4; a++) {
+                let connSum = 0;
+                for (let lam = 0; lam < 4; lam++)
+                    connSum += Gamma(lam, mu, a) * xiLower[lam];
+                cov[mu][a] = partial[mu][a] - connSum;
+            }
+        }
+    }
+    else {
+        // Pure FD fallback: differentiate ξ_α(y) = g_{αβ}(y) ξ^β(y) numerically.
+        // WARNING: large-scale SI metrics (c²·g_tt ~ 6e16) produce cancellation
+        // noise at the ~1-10 absolute level. Use constantKilling=true for
+        // constant Schwarzschild Killing fields.
+        const partial = [
+            [0, 0, 0, 0],
+            [0, 0, 0, 0],
+            [0, 0, 0, 0],
+            [0, 0, 0, 0],
+        ];
+        for (let a = 0; a < 4; a++) {
+            const xiLowerA = (y) => {
+                const gy = metricFn(y);
+                const xy = xi(y);
+                let s = 0;
+                for (let b = 0; b < 4; b++)
+                    s += gy[a][b] * xy[b];
+                return s;
+            };
+            for (let mu = 0; mu < 4; mu++) {
+                const raw = pderivNumericalFn(xiLowerA, x, mu, { order });
+                partial[mu][a] = typeof raw === 'number' ? raw : raw[0];
+            }
+        }
+        for (let mu = 0; mu < 4; mu++) {
+            for (let a = 0; a < 4; a++) {
+                let connSum = 0;
+                for (let lam = 0; lam < 4; lam++)
+                    connSum += Gamma(lam, mu, a) * xiLower[lam];
+                cov[mu][a] = partial[mu][a] - connSum;
+            }
+        }
+    }
+    // Symmetrize ∇_μ ξ_ν + ∇_ν ξ_μ and find max|·| over 16 pairs.
+    let maxResid = 0;
+    for (let mu = 0; mu < 4; mu++) {
+        for (let nu = 0; nu < 4; nu++) {
+            const sym = cov[mu][nu] + cov[nu][mu];
+            const ab = Math.abs(sym);
+            if (ab > maxResid)
+                maxResid = ab;
+        }
+    }
+    return maxResid;
+}
+/**
+ * Evaluate the conserved charge Q = ξ^μ p_μ at point x along a geodesic.
+ *
+ * Sign convention: raw contraction ξ^μ p_μ. For the time-translation
+ * Killing field in (−,+,+,+) signature, p_t < 0 for forward-time motion,
+ * so the physical energy of common usage is E = −Q (Carroll Eq. 8.30).
+ * Use `|ΔQ / Q|` for drift diagnostics to stay sign-agnostic.
+ *
+ * @param killingFn     - Factory returning ξ^μ(x) (upper index).
+ * @param momentumLower - Co-momentum p_μ at the geodesic point (lower index).
+ * @param x             - Geodesic coordinate [t, r, θ, φ].
+ * @returns             Q = ξ^μ p_μ (dimensioned as energy if p is 4-momentum).
+ *
+ * @example
+ * ```typescript
+ * import { evaluateConservedCharge } from 'universal-physics-tensor';
+ * import { schwarzschildKillingT } from '../tests/fixtures/schwarzschild.js';
+ *
+ * // Co-momentum p_μ at a geodesic snapshot (p_t < 0 for forward-time motion)
+ * const pLower: [number, number, number, number] = [-8.988e16, 0, 0, 0];
+ * const x: [number, number, number, number] = [0, 1e7, Math.PI / 2, 0];
+ *
+ * const Q = evaluateConservedCharge(schwarzschildKillingT, pLower, x);
+ * // Q = ξ^t · p_t = 1 · p_t = p_t (energy with sign from metric)
+ * // Physical energy E = −Q (Carroll sign convention for (−,+,+,+))
+ * ```
+ *
+ * @public
+ */
+export function evaluateConservedCharge(killingFn, momentumLower, x) {
+    const xi = killingFn()(x);
+    let q = 0;
+    for (let m = 0; m < 4; m++)
+        q += xi[m] * momentumLower[m];
+    return q;
+}
+//# sourceMappingURL=killing.js.map

package/dist/numerical/killing.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"killing.js","sourceRoot":"","sources":["../../src/numerical/killing.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;;;;GAsBG;AAEH,OAAO,EAAE,iBAAiB,EAAE,MAAM,aAAa,CAAC;AA0DhD;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GAgDG;AACH,MAAM,UAAU,qBAAqB,CACnC,SAAoB,EACpB,QAAkB,EAClB,aAA8B,EAC9B,CAAO,EACP,OAA+B,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;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"}

package/dist/numerical/kretschmann.d.ts

@@ -0,0 +1,62 @@
+/**
+ * Kretschmann scalar numerical contraction (v0.6.0 Phase 3, Task 3.6).
+ *
+ * Computes K = R_{ρσμν} R^{ρσμν} from the fully-lowered Riemann tensor and
+ * the inverse metric by raising all four indices on the second Riemann factor.
+ *
+ * **PERFORMANCE WARNING**: O(4⁸) = 65536 multiplications per call. Acceptable
+ * for diagnostic / sample-point use, but DO NOT call inside tight integration
+ * loops (10⁶ steps × 65k ops = ~6.5e10 ops). Future optimization (v0.7.0
+ * candidate): exploit Riemann symmetries to reduce to ~20 independent
+ * components × symmetric weighting; estimated speedup ~12×. Tracked in
+ * v0.6.0 brainstorm carry-forward.
+ *
+ * @module numerical/kretschmann
+ */
+/**
+ * Kretschmann invariant K = R_{ρσμν} R^{ρσμν}.
+ *
+ * Algorithm:
+ *   1. Raise all four indices on the second Riemann factor:
+ *      R^{ρσμν} = g^{ρα} g^{σβ} g^{μγ} g^{νδ} R_{αβγδ}
+ *   2. Contract:
+ *      K = Σ_{ρσμν} R_{ρσμν} · R^{ρσμν}
+ *
+ * The outer sum has 4⁴ = 256 terms; each inner raise requires 4⁴ = 256
+ * metric products, giving the O(4⁸) = 65536 total multiplications.
+ *
+ * @param riemannLower - R_{ρσμν}: all-lower Riemann tensor, shape [4][4][4][4].
+ *   Obtain by lowering R^ρ_{σμν} via g_{ρα}: R_{αβγδ} = g_{αρ} R^ρ_{βγδ}.
+ * @param metricInverse - g^{αβ}: inverse metric, shape [4][4].
+ * @returns K, a scalar. Dimension [L⁻⁴] (Riemann is [L⁻²] per F8/I3 convention).
+ *
+ * @example
+ * ```typescript
+ * import { computeKretschmann, G_SI, C_SI } from 'universal-physics-tensor';
+ * import {
+ *   schwarzschildGFn,
+ *   schwarzschildGInverseFn,
+ *   schwarzschildRs,
+ * } from '../tests/fixtures/schwarzschild.js';
+ * import { riemannLowerAt } from '../src/numerical/curvature-lowering-helpers.js';
+ * import { Float64ReferenceEngine } from '../src/numerical/float64-engine.js';
+ *
+ * const M = 1.989e30;
+ * const r_s = schwarzschildRs(M);
+ * const r = 5 * r_s;
+ * const x: [number, number, number, number] = [0, r, Math.PI / 2, 0];
+ *
+ * const engine = new Float64ReferenceEngine();
+ * const rLower = riemannLowerAt(x, schwarzschildGFn(M), schwarzschildGInverseFn(M), 4, engine);
+ * const gInv = schwarzschildGInverseFn(M)(x);
+ * const K = computeKretschmann(rLower, gInv);
+ *
+ * // Schwarzschild closed-form: K = 48 G² M² / (c⁴ r⁶)
+ * const K_analytic = 48 * G_SI**2 * M**2 / (C_SI**4 * r**6);
+ * // |K - K_analytic| / K_analytic < 1e-4
+ * ```
+ *
+ * @public
+ */
+export declare function computeKretschmann(riemannLower: number[][][][], metricInverse: number[][]): number;
+//# sourceMappingURL=kretschmann.d.ts.map

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

@@ -0,0 +1 @@
+{"version":3,"file":"kretschmann.d.ts","sourceRoot":"","sources":["../../src/numerical/kretschmann.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;GAcG;AAEH;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GA4CG;AACH,wBAAgB,kBAAkB,CAChC,YAAY,EAAE,MAAM,EAAE,EAAE,EAAE,EAAE,EAC5B,aAAa,EAAE,MAAM,EAAE,EAAE,GACxB,MAAM,CA0BR"}

package/dist/numerical/kretschmann.js

@@ -0,0 +1,88 @@
+/**
+ * Kretschmann scalar numerical contraction (v0.6.0 Phase 3, Task 3.6).
+ *
+ * Computes K = R_{ρσμν} R^{ρσμν} from the fully-lowered Riemann tensor and
+ * the inverse metric by raising all four indices on the second Riemann factor.
+ *
+ * **PERFORMANCE WARNING**: O(4⁸) = 65536 multiplications per call. Acceptable
+ * for diagnostic / sample-point use, but DO NOT call inside tight integration
+ * loops (10⁶ steps × 65k ops = ~6.5e10 ops). Future optimization (v0.7.0
+ * candidate): exploit Riemann symmetries to reduce to ~20 independent
+ * components × symmetric weighting; estimated speedup ~12×. Tracked in
+ * v0.6.0 brainstorm carry-forward.
+ *
+ * @module numerical/kretschmann
+ */
+/**
+ * Kretschmann invariant K = R_{ρσμν} R^{ρσμν}.
+ *
+ * Algorithm:
+ *   1. Raise all four indices on the second Riemann factor:
+ *      R^{ρσμν} = g^{ρα} g^{σβ} g^{μγ} g^{νδ} R_{αβγδ}
+ *   2. Contract:
+ *      K = Σ_{ρσμν} R_{ρσμν} · R^{ρσμν}
+ *
+ * The outer sum has 4⁴ = 256 terms; each inner raise requires 4⁴ = 256
+ * metric products, giving the O(4⁸) = 65536 total multiplications.
+ *
+ * @param riemannLower - R_{ρσμν}: all-lower Riemann tensor, shape [4][4][4][4].
+ *   Obtain by lowering R^ρ_{σμν} via g_{ρα}: R_{αβγδ} = g_{αρ} R^ρ_{βγδ}.
+ * @param metricInverse - g^{αβ}: inverse metric, shape [4][4].
+ * @returns K, a scalar. Dimension [L⁻⁴] (Riemann is [L⁻²] per F8/I3 convention).
+ *
+ * @example
+ * ```typescript
+ * import { computeKretschmann, G_SI, C_SI } from 'universal-physics-tensor';
+ * import {
+ *   schwarzschildGFn,
+ *   schwarzschildGInverseFn,
+ *   schwarzschildRs,
+ * } from '../tests/fixtures/schwarzschild.js';
+ * import { riemannLowerAt } from '../src/numerical/curvature-lowering-helpers.js';
+ * import { Float64ReferenceEngine } from '../src/numerical/float64-engine.js';
+ *
+ * const M = 1.989e30;
+ * const r_s = schwarzschildRs(M);
+ * const r = 5 * r_s;
+ * const x: [number, number, number, number] = [0, r, Math.PI / 2, 0];
+ *
+ * const engine = new Float64ReferenceEngine();
+ * const rLower = riemannLowerAt(x, schwarzschildGFn(M), schwarzschildGInverseFn(M), 4, engine);
+ * const gInv = schwarzschildGInverseFn(M)(x);
+ * const K = computeKretschmann(rLower, gInv);
+ *
+ * // Schwarzschild closed-form: K = 48 G² M² / (c⁴ r⁶)
+ * const K_analytic = 48 * G_SI**2 * M**2 / (C_SI**4 * r**6);
+ * // |K - K_analytic| / K_analytic < 1e-4
+ * ```
+ *
+ * @public
+ */
+export function computeKretschmann(riemannLower, metricInverse) {
+    // Raise: R^{ρσμν} = g^{ρα} g^{σβ} g^{μγ} g^{νδ} R_{αβγδ}
+    // Then K = R_{ρσμν} R^{ρσμν} (sum over all 4⁴ = 256 components)
+    let K = 0;
+    for (let rho = 0; rho < 4; rho++) {
+        for (let sigma = 0; sigma < 4; sigma++) {
+            for (let mu = 0; mu < 4; mu++) {
+                for (let nu = 0; nu < 4; nu++) {
+                    let raised = 0;
+                    for (let a = 0; a < 4; a++) {
+                        for (let b = 0; b < 4; b++) {
+                            for (let c = 0; c < 4; c++) {
+                                for (let d = 0; d < 4; d++) {
+                                    raised += metricInverse[rho][a] * metricInverse[sigma][b]
+                                        * metricInverse[mu][c] * metricInverse[nu][d]
+                                        * riemannLower[a][b][c][d];
+                                }
+                            }
+                        }
+                    }
+                    K += riemannLower[rho][sigma][mu][nu] * raised;
+                }
+            }
+        }
+    }
+    return K;
+}
+//# sourceMappingURL=kretschmann.js.map

package/dist/numerical/kretschmann.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"kretschmann.js","sourceRoot":"","sources":["../../src/numerical/kretschmann.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;GAcG;AAEH;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GA4CG;AACH,MAAM,UAAU,kBAAkB,CAChC,YAA4B,EAC5B,aAAyB;IAEzB,yDAAyD;IACzD,gEAAgE;IAChE,IAAI,CAAC,GAAG,CAAC,CAAC;IACV,KAAK,IAAI,GAAG,GAAG,CAAC,EAAE,GAAG,GAAG,CAAC,EAAE,GAAG,EAAE,EAAE,CAAC;QACjC,KAAK,IAAI,KAAK,GAAG,CAAC,EAAE,KAAK,GAAG,CAAC,EAAE,KAAK,EAAE,EAAE,CAAC;YACvC,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,CAAC,EAAE,EAAE,EAAE,EAAE,CAAC;gBAC9B,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,CAAC,EAAE,EAAE,EAAE,EAAE,CAAC;oBAC9B,IAAI,MAAM,GAAG,CAAC,CAAC;oBACf,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;wBAC3B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;4BAC3B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;gCAC3B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;oCAC3B,MAAM,IAAI,aAAa,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,aAAa,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC;0CAC/C,aAAa,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,aAAa,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;0CAC3C,YAAY,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;gCACrC,CAAC;4BACH,CAAC;wBACH,CAAC;oBACH,CAAC;oBACD,CAAC,IAAI,YAAY,CAAC,GAAG,CAAC,CAAC,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,GAAG,MAAM,CAAC;gBACjD,CAAC;YACH,CAAC;QACH,CAAC;IACH,CAAC;IACD,OAAO,CAAC,CAAC;AACX,CAAC"}

package/dist/numerical/lowering-utils.d.ts

@@ -0,0 +1,47 @@
+/**
+ * Shared private utilities for the `numerical/lowering*.ts` modules.
+ *
+ * Carved out of `lowering.ts` in the v0.7 follow-up to its v0.6.1
+ * sprint-target miss (LOC 903 vs ≤890). Centralizes the four helpers
+ * that both `lowering.ts` and the v0.7 `derivative-lowering.ts`
+ * need, avoiding both a forward-import cycle and code duplication.
+ *
+ * No public surface — every export carries an `@internal` JSDoc tag.
+ *
+ * @module numerical/lowering-utils
+ */
+import type { NumericalInputs, NestedArray } from './types.js';
+import type { MetricTensorNode } from '../dimensional/metric-validators.js';
+/**
+ * Module-private type predicate for the `metric-tensor` AST kind.
+ * Replaces the duplicated `(n as { kind?: unknown }).kind ===
+ * 'metric-tensor'` cast pattern. TypeScript narrows the operand to
+ * `MetricTensorNode` after the predicate, removing follow-up `as
+ * MetricTensorNode` casts.
+ *
+ * @internal
+ */
+export declare function isMetricTensorNode(n: unknown): n is MetricTensorNode;
+/**
+ * Coordinate-space dimension supplied via `NumericalInputs.dimension`,
+ * defaulting to 4 (spacetime). Read at most case-arms in lowering.
+ *
+ * @internal
+ */
+export declare function dimensionOf(inputs: NumericalInputs): number;
+/**
+ * Look up a named tensor's concrete value from
+ * `NumericalInputs.tensors`, or throw with a clear message.
+ *
+ * @internal
+ */
+export declare function requireValue(name: string, inputs: NumericalInputs): NestedArray;
+/**
+ * Flatten a `NestedArray` to a plain `number[]` and check the result
+ * matches `expectedSize`. Wraps the canonical `flattenNA()` with a
+ * size-mismatch guard so consumers don't silently swallow bad input.
+ *
+ * @internal
+ */
+export declare function flattenNestedArray(data: NestedArray, expectedSize: number): number[];
+//# sourceMappingURL=lowering-utils.d.ts.map

package/dist/numerical/lowering-utils.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"lowering-utils.d.ts","sourceRoot":"","sources":["../../src/numerical/lowering-utils.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;GAWG;AAEH,OAAO,KAAK,EAAE,eAAe,EAAE,WAAW,EAAE,MAAM,YAAY,CAAC;AAC/D,OAAO,KAAK,EAAE,gBAAgB,EAAE,MAAM,qCAAqC,CAAC;AAI5E;;;;;;;;GAQG;AACH,wBAAgB,kBAAkB,CAAC,CAAC,EAAE,OAAO,GAAG,CAAC,IAAI,gBAAgB,CAGpE;AAED;;;;;GAKG;AACH,wBAAgB,WAAW,CAAC,MAAM,EAAE,eAAe,GAAG,MAAM,CAE3D;AAED;;;;;GAKG;AACH,wBAAgB,YAAY,CAAC,IAAI,EAAE,MAAM,EAAE,MAAM,EAAE,eAAe,GAAG,WAAW,CAM/E;AAED;;;;;;GAMG;AACH,wBAAgB,kBAAkB,CAAC,IAAI,EAAE,WAAW,EAAE,YAAY,EAAE,MAAM,GAAG,MAAM,EAAE,CAQpF"}

package/dist/numerical/lowering-utils.js

@@ -0,0 +1,64 @@
+/**
+ * Shared private utilities for the `numerical/lowering*.ts` modules.
+ *
+ * Carved out of `lowering.ts` in the v0.7 follow-up to its v0.6.1
+ * sprint-target miss (LOC 903 vs ≤890). Centralizes the four helpers
+ * that both `lowering.ts` and the v0.7 `derivative-lowering.ts`
+ * need, avoiding both a forward-import cycle and code duplication.
+ *
+ * No public surface — every export carries an `@internal` JSDoc tag.
+ *
+ * @module numerical/lowering-utils
+ */
+import { NumericalBackendError } from './errors.js';
+import { flattenNA } from './connection-lowering-helpers.js';
+/**
+ * Module-private type predicate for the `metric-tensor` AST kind.
+ * Replaces the duplicated `(n as { kind?: unknown }).kind ===
+ * 'metric-tensor'` cast pattern. TypeScript narrows the operand to
+ * `MetricTensorNode` after the predicate, removing follow-up `as
+ * MetricTensorNode` casts.
+ *
+ * @internal
+ */
+export function isMetricTensorNode(n) {
+    return typeof n === 'object' && n !== null
+        && n.kind === 'metric-tensor';
+}
+/**
+ * Coordinate-space dimension supplied via `NumericalInputs.dimension`,
+ * defaulting to 4 (spacetime). Read at most case-arms in lowering.
+ *
+ * @internal
+ */
+export function dimensionOf(inputs) {
+    return inputs.dimension ?? 4;
+}
+/**
+ * Look up a named tensor's concrete value from
+ * `NumericalInputs.tensors`, or throw with a clear message.
+ *
+ * @internal
+ */
+export function requireValue(name, inputs) {
+    const v = inputs.tensors.get(name);
+    if (v === undefined) {
+        throw new NumericalBackendError(`lowering: no value supplied for "${name}" in inputs.tensors`);
+    }
+    return v;
+}
+/**
+ * Flatten a `NestedArray` to a plain `number[]` and check the result
+ * matches `expectedSize`. Wraps the canonical `flattenNA()` with a
+ * size-mismatch guard so consumers don't silently swallow bad input.
+ *
+ * @internal
+ */
+export function flattenNestedArray(data, expectedSize) {
+    const out = flattenNA(data);
+    if (out.length !== expectedSize) {
+        throw new NumericalBackendError(`lowering: flattenNestedArray: got ${out.length} elements, expected ${expectedSize}`);
+    }
+    return out;
+}
+//# sourceMappingURL=lowering-utils.js.map

package/dist/numerical/lowering-utils.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"lowering-utils.js","sourceRoot":"","sources":["../../src/numerical/lowering-utils.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;GAWG;AAIH,OAAO,EAAE,qBAAqB,EAAE,MAAM,aAAa,CAAC;AACpD,OAAO,EAAE,SAAS,EAAE,MAAM,kCAAkC,CAAC;AAE7D;;;;;;;;GAQG;AACH,MAAM,UAAU,kBAAkB,CAAC,CAAU;IAC3C,OAAO,OAAO,CAAC,KAAK,QAAQ,IAAI,CAAC,KAAK,IAAI;WACpC,CAAwB,CAAC,IAAI,KAAK,eAAe,CAAC;AAC1D,CAAC;AAED;;;;;GAKG;AACH,MAAM,UAAU,WAAW,CAAC,MAAuB;IACjD,OAAO,MAAM,CAAC,SAAS,IAAI,CAAC,CAAC;AAC/B,CAAC;AAED;;;;;GAKG;AACH,MAAM,UAAU,YAAY,CAAC,IAAY,EAAE,MAAuB;IAChE,MAAM,CAAC,GAAG,MAAM,CAAC,OAAO,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC;IACnC,IAAI,CAAC,KAAK,SAAS,EAAE,CAAC;QACpB,MAAM,IAAI,qBAAqB,CAAC,oCAAoC,IAAI,qBAAqB,CAAC,CAAC;IACjG,CAAC;IACD,OAAO,CAAC,CAAC;AACX,CAAC;AAED;;;;;;GAMG;AACH,MAAM,UAAU,kBAAkB,CAAC,IAAiB,EAAE,YAAoB;IACxE,MAAM,GAAG,GAAG,SAAS,CAAC,IAAI,CAAC,CAAC;IAC5B,IAAI,GAAG,CAAC,MAAM,KAAK,YAAY,EAAE,CAAC;QAChC,MAAM,IAAI,qBAAqB,CAC7B,qCAAqC,GAAG,CAAC,MAAM,uBAAuB,YAAY,EAAE,CACrF,CAAC;IACJ,CAAC;IACD,OAAO,GAAG,CAAC;AACb,CAAC"}

package/dist/numerical/lowering.d.ts

@@ -11,37 +11,9 @@
  * @module numerical/lowering
  */
 import type { ExprNode } from '../dimensional/validator.js';
-import type { TensorIndex } from '../dimensional/tensor.js';
-import type { EngineTensor, TensorEngine, EinsumSpec } from './tensor-engine.js';
+import type { EngineTensor, TensorEngine } from './tensor-engine.js';
 import type { NumericalInputs } from './types.js';
-/** Operand kinds a flat tensor-product can lower in v0.3.5: the three
- *  index-carrying nodes plus tensor-partial-derivative (whose effective
- *  indices are `of`'s indices followed by its wrtIndex). */
-type ContractableNode = Extract<ExprNode, {
-    indices: ReadonlyArray<TensorIndex>;
-}> | Extract<ExprNode, {
-    kind: 'tensor-partial-derivative';
-}>;
-/**
- * Build the EinsumSpec for a flat tensor-product.
- *
- * NOTE: this function does NOT decide which indices contract. That authority
- * belongs to `computeContraction()` (v0.2.0 symbolic-layer) — variance-aware,
- * implicit-identity-metric rule — which already classified every label as
- * contracted or free. buildEinsumSpec only maps those already-classified
- * labels to their (operand, axis) sites. There is exactly one
- * contraction-decision implementation in the codebase.
- *
- * @internal — cross-module/test use only; not part of the consumer surface.
- */
-export declare function buildEinsumSpec(operands: ReadonlyArray<ContractableNode>, contractionPairs: ReadonlyArray<{
-    label: string;
-}>, freeIndices: ReadonlyMap<string, {
-    upper: number;
-    lower: number;
-}>): EinsumSpec;
 /** Lower a validated ExprNode to an EngineTensor.
  *  @internal — cross-module/test use only; not part of the consumer surface. */
 export declare function lowerNode(node: ExprNode, inputs: NumericalInputs, engine: TensorEngine): EngineTensor;
-export {};
 //# sourceMappingURL=lowering.d.ts.map

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

@@ -1 +1 @@
-{"version":3,"file":"lowering.d.ts","sourceRoot":"","sources":["../../src/numerical/lowering.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;GAWG;AAEH,OAAO,KAAK,EAAE,QAAQ,EAAE,MAAM,6BAA6B,CAAC;AAE5D,OAAO,KAAK,EAAE,WAAW,EAAoB,MAAM,0BAA0B,CAAC;AAY9E,OAAO,KAAK,EACV,YAAY,EAAE,YAAY,EAAE,UAAU,EACvC,MAAM,oBAAoB,CAAC;AAC5B,OAAO,KAAK,EAAE,eAAe,EAAe,MAAM,YAAY,CAAC;AAqB/D;;4DAE4D;AAC5D,KAAK,gBAAgB,GACjB,OAAO,CAAC,QAAQ,EAAE;IAAE,OAAO,EAAE,aAAa,CAAC,WAAW,CAAC,CAAA;CAAE,CAAC,GAC1D,OAAO,CAAC,QAAQ,EAAE;IAAE,IAAI,EAAE,2BAA2B,CAAA;CAAE,CAAC,CAAC;AAqD7D;;;;;;;;;;;GAWG;AACH,wBAAgB,eAAe,CAC7B,QAAQ,EAAE,aAAa,CAAC,gBAAgB,CAAC,EACzC,gBAAgB,EAAE,aAAa,CAAC;IAAE,KAAK,EAAE,MAAM,CAAA;CAAE,CAAC,EAClD,WAAW,EAAE,WAAW,CAAC,MAAM,EAAE;IAAE,KAAK,EAAE,MAAM,CAAC;IAAC,KAAK,EAAE,MAAM,CAAA;CAAE,CAAC,GACjE,UAAU,CA0CZ;AAiBD;gFACgF;AAChF,wBAAgB,SAAS,CACvB,IAAI,EAAE,QAAQ,EACd,MAAM,EAAE,eAAe,EACvB,MAAM,EAAE,YAAY,GACnB,YAAY,CA2lBd"}
+{"version":3,"file":"lowering.d.ts","sourceRoot":"","sources":["../../src/numerical/lowering.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;GAWG;AAEH,OAAO,KAAK,EAAE,QAAQ,EAAE,MAAM,6BAA6B,CAAC;AAiB5D,OAAO,KAAK,EACV,YAAY,EAAE,YAAY,EAC3B,MAAM,oBAAoB,CAAC;AAC5B,OAAO,KAAK,EAAE,eAAe,EAAe,MAAM,YAAY,CAAC;AAwW/D;gFACgF;AAChF,wBAAgB,SAAS,CACvB,IAAI,EAAE,QAAQ,EACd,MAAM,EAAE,eAAe,EACvB,MAAM,EAAE,YAAY,GACnB,YAAY,CAqMd"}