universal-physics-tensor 0.4.6 → 0.5.1

package/dist/numerical/be37-covariant-eikonal.js

@@ -1,51 +1,175 @@
 /**
- * v0.4.0 BE-37 covariant-eikonal STRUCTURAL PREVIEW.
+ * v0.5.0 BE-37 covariant-eikonal numerical evaluator.
  *
  * Provides `evaluateBE37CovariantEikonalNumerical`, which encodes the
  * covariant-eikonal equation
  *
  *   g^μν ∇_μ ∇_ν S = 0
  *
- * via the connection layer added in Tasks 2-3 + 12-14.  In v0.4.0 the
- * eikonal residual is **0 by construction**: the null wave-covector
+ * via the connection layer (Tasks 2-3 + 12-14 of v0.4.0). The eikonal
+ * residual is **0 by construction**: the null wave-covector
  *
  *   k_μ = ∂_μ S = (E, E, 0, 0)   (energy + radial component only)
  *
  * satisfies g^μν k_μ k_ν = 0 identically for the +,−,−,− Schwarzschild
- * metric.  No numerical integration is performed for this quantity.
+ * metric. No numerical integration is performed for this quantity.
  *
- * The `shapiroDelaySec` field is an **intentional v0.4.0 stub** (returns 0).
- * The full geodesic-integrated Shapiro cross-check — comparing
- * `shapiroDelaySec` against `evaluateShapiroDelay` — is the v0.5.0
- * deliverable; it requires wiring `integrateGeodesic` end-to-end through
- * the covariant-eikonal path, which is deferred.
+ * `shapiroDelaySec` is computed by **GL4-driven null-geodesic integration**
+ * on the canonical (x, p) state in the Schwarzschild spacetime (v0.5.0 Task
+ * 12, Phase 2c). The photon propagates from `r = R_far_m` down to
+ * `r = R_near_m`; coordinate time at the endpoint minus the flat-space
+ * straight-line transit time is returned as the Shapiro delay.
  *
- * @see tests/dimensional/covariant-derivative-preview.test.ts (Task 17 [U])
- * @see docs/planning/v0.4.0-Implementation-Plan.md Task 17
+ * For the default radial geometry (impact parameter `b_m = 0`) the
+ * integrator reproduces the closed-form `(2GM/c³) · ln(R_far/R_near)`
+ * (Task 11 cross-checks this to ±2×10⁻³).
+ *
+ * @see tests/numerical/be37-covariant-eikonal-real.test.ts (Task 12 [U])
+ * @see tests/dimensional/covariant-derivative-preview.test.ts (Task 17 [U]/Task 11 [U])
+ * @see docs/planning/v0.5.0-Implementation-Plan.md Task 12
  * @module numerical/be37-covariant-eikonal
  * @public
  */
+import { integrateGeodesicGL4 } from './gl4-integrator.js';
+import { C_SI, G_SI } from '../core/constants.js';
+// ─── Physical constants (canonical CODATA 2018 + exact-SI definitions) ────
+// Imported from src/core/constants.ts (v0.5.1 PC-1 canonicalization).
+const c_SI = C_SI;
+const c2_SI = c_SI * c_SI;
 /**
- * v0.4.0 structural preview — returns the eikonal residual (exactly 0 by
- * construction for a null ray) and a stub Shapiro delay (0 in v0.4.0).
+ * Schwarzschild radius `r_s = 2GM/c²`.
  *
- * Does NOT use the CovariantDerivativeNode or lowering infrastructure —
- * see the file-level module docstring for the structural rationale. The
- * domain guards (M_kg > 0, R_far_m > 0, R_near_m ∈ (0, R_far_m]) are
- * enforced for API consistency with evaluateShapiroDelay; the constants 0
- * do not depend on the inputs.
+ * Mirrors `tests/fixtures/schwarzschild.ts` (kept private here to avoid
+ * a src→tests dependency).
  *
- * The full geodesic-integrated Shapiro cross-check (wiring integrateGeodesic
- * through this path) is the v0.5.0 deliverable. Until then, both return
- * values are structural constants, not computed results.
+ * @internal
+ */
+function schwarzschildRs(M_kg) {
+    return (2 * G_SI * M_kg) / c2_SI;
+}
+/**
+ * Returns the contravariant Schwarzschild metric closure g^{μν}(x).
  *
- * @param inputs - Source mass + far/near radii (validated; unused in stub).
- * @returns `{eikonalResidual: 0, shapiroDelaySec: 0}`.
+ * SI convention with c² on `g_{tt}` — mirrors the test fixture so the
+ * coordinate time `t = x^0` carries seconds and `p_t = −c²` produces
+ * `dt/dτ ≈ 1` (τ ≈ coord time far from horizon).
+ *
+ * @internal
+ */
+function buildGInverseFn(M_kg) {
+    const r_s = schwarzschildRs(M_kg);
+    return function gInverse(x) {
+        const r = x[1];
+        const theta = x[2];
+        const f = 1 - r_s / r;
+        const sinT = Math.sin(theta);
+        const gInv = [
+            [0, 0, 0, 0],
+            [0, 0, 0, 0],
+            [0, 0, 0, 0],
+            [0, 0, 0, 0],
+        ];
+        gInv[0][0] = -1 / (f * c2_SI);
+        gInv[1][1] = f;
+        gInv[2][2] = 1 / (r * r);
+        gInv[3][3] = 1 / (r * r * sinT * sinT);
+        return gInv;
+    };
+}
+/**
+ * Returns ∂_λ g^{μν}(x) with index order `dg[λ][μ][ν]`.
+ *
+ * Only ∂_r and ∂_θ entries are non-zero (static, axisymmetric metric).
+ * Sign of `dg[1][0][0] = ∂_r g^{tt} = +r_s/(r²(1−r_s/r)²c²)` matches the
+ * Task 0 fixture fix (positive — incorrectly negative in the v0.5.0
+ * pre-Task-3 fixture; reversal would push the photon outward).
+ *
+ * @internal
+ */
+function buildDgInverseFn(M_kg) {
+    const r_s = schwarzschildRs(M_kg);
+    return function dgInverse(x) {
+        const r = x[1];
+        const theta = x[2];
+        const sinT = Math.sin(theta);
+        const cosT = Math.cos(theta);
+        const f = 1 - r_s / r;
+        const dg = Array.from({ length: 4 }, () => Array.from({ length: 4 }, () => [0, 0, 0, 0]));
+        // ∂_r entries (axis 1)
+        dg[1][0][0] = r_s / (r * r * f * f * c2_SI);
+        dg[1][1][1] = r_s / (r * r);
+        dg[1][2][2] = -2 / (r * r * r);
+        dg[1][3][3] = -2 / (r * r * r * sinT * sinT);
+        // ∂_θ entry (axis 2)
+        dg[2][3][3] = (-2 * cosT) / (r * r * sinT * sinT * sinT);
+        return dg;
+    };
+}
+/**
+ * v0.5.0 BE-37 covariant-eikonal numerical evaluator.
+ *
+ * Returns:
+ *  - `eikonalResidual = 0` by construction (null wave-covector identity).
+ *  - `shapiroDelaySec`: GL4-integrated coordinate-time delay relative to
+ *    the flat-space straight-line transit, in seconds.
+ *
+ * **Algorithm — null geodesic in Schwarzschild via GL4.**
+ *
+ * 1. Schwarzschild radius `r_s = 2GM/c²`. Domain guard `R_near_m > r_s`
+ *    (else the photon path crosses the event horizon).
+ * 2. Set initial state at `(t=0, r=R_far_m, θ=π/2, φ=0)`. Affine-parameter
+ *    normalization `p_t = −c²` → `dt/dτ = 1/(1−r_s/r) ≈ 1` far from the
+ *    horizon (τ ≈ coord time). Set `p_θ = 0` (motion confined to equatorial
+ *    plane), `p_φ = b_m · c` (Killing-vector angular momentum
+ *    corresponding to flat-limit impact parameter `b_m`).
+ * 3. Solve `p_r` from the null condition `g^μν p_μ p_ν = 0` at the
+ *    initial point, take the negative root (inward motion: `dr/dτ < 0`
+ *    since `dr/dτ = g^rr p_r = (1−r_s/r) p_r`).
+ * 4. Integrate via `integrateGeodesicGL4` for `tauMax` chosen to overshoot
+ *    `r = R_near_m` (safety factor 1.5 on the flat-space transit time).
+ * 5. Walk snapshots; find the first index where `r ≤ R_near_m`. Linearly
+ *    interpolate `t` between the bracketing snapshots to land exactly on
+ *    `r = R_near_m`. Take that `t` as the geodesic coord-time transit.
+ * 6. Flat-space reference: for a straight-line null ray with impact
+ *    parameter `b`, time to go from `r_far` to `r_near` (single leg, no
+ *    closest-approach crossing) is `(√(r_far²−b²) − √(r_near²−b²))/c`.
+ *    For `b = 0` (default) this reduces to `(R_far − R_near)/c`.
+ * 7. `shapiroDelaySec = t_geodesic − t_flat`.
+ *
+ * **Sign convention.** Shapiro delay is positive (light slowed by the
+ * gravitational potential). A negative result indicates an integration
+ * or reference-frame bug.
+ *
+ * **Impact parameter `b_m`.** Default `0` (radial). The default geometry
+ * reproduces the closed-form `evaluateShapiroDelay` to ≤1e-6 relative at
+ * 2048 steps for solar-scale Earth/Mars distances (Task 11 cross-check).
+ * Non-zero `b_m` requires `b_m ≤ R_near_m` (closest approach `r_min ≈ b`
+ * must not lie between R_far and R_near).
+ *
+ * @param inputs - Source mass + far/near radii (+ optional `b_m`, `steps`):
+ *   - `M_kg` — gravitational source mass in **kilograms** (SI). Must be a
+ *     finite positive number.
+ *   - `R_far_m` — origin radial coordinate in **metres** (SI). The signal
+ *     start point.
+ *   - `R_near_m` — destination radial coordinate in **metres** (SI). Must
+ *     satisfy `0 < R_near_m ≤ R_far_m` and `R_near_m > 1.01 · r_s` (stay
+ *     comfortably outside the Schwarzschild horizon `r_s = 2GM/c²`).
+ *   - `b_m` — optional impact parameter in **metres** (SI). Default `0`
+ *     (purely radial geodesic). Non-zero `b_m` requires `b_m ≤ R_near_m`.
+ *   - `steps` — optional GL4 integrator step count (dimensionless). Default
+ *     2048. Higher = tighter agreement with closed-form Shapiro at
+ *     additional CPU cost.
+ * @returns `{ eikonalResidual, shapiroDelaySec }`:
+ *   - `eikonalResidual` — **dimensionless** numerical residual of
+ *     `g^μν ∇_μ ∇_ν S = 0`; returns exactly `0` by construction.
+ *   - `shapiroDelaySec` — gravitational-time-delay in **seconds** (SI).
+ *     Positive by sign convention (light is slowed in the potential well).
+ * @throws RangeError on out-of-domain inputs.
  * @public
  */
 export async function evaluateBE37CovariantEikonalNumerical(inputs) {
-    const { M_kg, R_far_m, R_near_m } = inputs;
-    // Domain guards (mirrors evaluateShapiroDelay for API consistency).
+    const { M_kg, R_far_m, R_near_m, b_m = 0, steps = 2048 } = inputs;
+    // ─── Domain guards ──────────────────────────────────────────────────────
     if (!Number.isFinite(M_kg) || M_kg <= 0) {
         throw new RangeError(`evaluateBE37CovariantEikonalNumerical: M_kg must be a finite positive number, got ${M_kg}`);
     }
@@ -58,23 +182,78 @@ export async function evaluateBE37CovariantEikonalNumerical(inputs) {
     if (R_near_m > R_far_m) {
         throw new RangeError(`evaluateBE37CovariantEikonalNumerical: R_near_m (${R_near_m}) must be ≤ R_far_m (${R_far_m})`);
     }
+    if (!Number.isFinite(b_m) || b_m < 0) {
+        throw new RangeError(`evaluateBE37CovariantEikonalNumerical: b_m must be a finite non-negative number, got ${b_m}`);
+    }
+    if (b_m > R_near_m) {
+        throw new RangeError(`evaluateBE37CovariantEikonalNumerical: b_m (${b_m}) must be ≤ R_near_m (${R_near_m}); closest-approach geometry not supported`);
+    }
+    const r_s = schwarzschildRs(M_kg);
+    if (R_near_m <= r_s * 1.01) {
+        throw new RangeError(`evaluateBE37CovariantEikonalNumerical: R_near_m (${R_near_m}) must be safely outside the Schwarzschild horizon r_s=${r_s}`);
+    }
     // ─── Eikonal residual ────────────────────────────────────────────────────
-    //
-    // The covariant-eikonal equation g^μν ∇_μ ∇_ν S = 0 is evaluated at
-    // the null wave-covector k_μ = ∂_μ S = (E, E, 0, 0) in Schwarzschild
-    // coordinates.  For a null ray, by definition:
-    //
-    //   g^μν k_μ k_ν = 0
-    //
-    // This is a structural consequence of the null-ray construction — it does
-    // not depend on M_kg, R_far_m, or R_near_m.  The residual is exactly 0.
+    // g^μν k_μ k_ν = 0 by construction for the null wave-covector
+    // k_μ = (E, E, 0, 0). Structural — independent of M_kg, R_far, R_near.
     const eikonalResidual = 0;
-    // ─── Shapiro delay (v0.4.0 stub) ────────────────────────────────────────
-    //
-    // The geodesic-integrated Shapiro delay is deferred to v0.5.0.
-    // shapiroDelaySec returns 0 here; the v0.5.0 task will wire
-    // integrateGeodesic through this path and remove the stub.
-    const shapiroDelaySec = 0;
+    // ─── Shapiro delay via GL4 null-geodesic integration ────────────────────
+    const gInverseFn = buildGInverseFn(M_kg);
+    const dgInverseFn = buildDgInverseFn(M_kg);
+    // Initial state at (t=0, r=R_far, θ=π/2, φ=0). Equatorial plane (p_θ=0).
+    // Affine-parameter normalization: p_t = -c² makes dt/dτ ≈ 1 far from
+    // horizon, so τ ≈ coord time t and tauMax can be set in units of seconds.
+    const x0 = [0, R_far_m, Math.PI / 2, 0];
+    const p_t = -c2_SI;
+    const p_phi = b_m * c_SI; // Killing-vector conservation: classical b = p_φ/(|p_t|/c²·c) = p_φ/c.
+    const p_theta = 0;
+    // Solve null condition g^μν p_μ p_ν = 0 for p_r at the initial point:
+    //   g^tt p_t² + g^rr p_r² + g^φφ p_φ² = 0  (θ=π/2)
+    //   ⇒ p_r² = (−g^tt p_t² − g^φφ p_φ²) / g^rr
+    const gInv0 = gInverseFn(x0);
+    const numerator = -gInv0[0][0] * p_t * p_t - gInv0[3][3] * p_phi * p_phi;
+    if (numerator < 0) {
+        throw new RangeError(`evaluateBE37CovariantEikonalNumerical: null condition has no real p_r at R_far=${R_far_m} (impact parameter b=${b_m} too large for this geometry; numerator=${numerator})`);
+    }
+    const p_r_magnitude = Math.sqrt(numerator / gInv0[1][1]);
+    // Inward motion: dr/dτ = g^rr p_r = (1-r_s/r) p_r < 0 ⇒ p_r < 0.
+    const p_r = -p_r_magnitude;
+    const initialState = {
+        x: x0,
+        p: [p_t, p_r, p_theta, p_phi],
+    };
+    // Choose tauMax to overshoot R_near with safety margin. With p_t = -c²
+    // and far from horizon, dr/dτ ≈ -c for b=0 (radial); slightly slower for
+    // non-zero b. Flat-space single-leg time = (√(R_far² − b²) − √(R_near² − b²))/c.
+    const flatLegFar = Math.sqrt(R_far_m * R_far_m - b_m * b_m);
+    const flatLegNear = Math.sqrt(R_near_m * R_near_m - b_m * b_m);
+    const t_flat = (flatLegFar - flatLegNear) / c_SI;
+    const tauMax = t_flat * 1.5;
+    const snapshots = integrateGeodesicGL4(initialState, {
+        steps,
+        tauMax,
+        gInverseFn,
+        dgInverseFn,
+        domainMinRadius: r_s * 1.01,
+    });
+    // Find the first snapshot with r ≤ R_near_m and linearly interpolate t.
+    let t_geodesic;
+    for (let i = 1; i < snapshots.length; i++) {
+        const rPrev = snapshots[i - 1].x[1];
+        const rCurr = snapshots[i].x[1];
+        if (rPrev > R_near_m && rCurr <= R_near_m) {
+            const tPrev = snapshots[i - 1].x[0];
+            const tCurr = snapshots[i].x[0];
+            // Linear interpolation: alpha = (R_near - rPrev) / (rCurr - rPrev), in [0, 1].
+            const alpha = (R_near_m - rPrev) / (rCurr - rPrev);
+            t_geodesic = tPrev + alpha * (tCurr - tPrev);
+            break;
+        }
+    }
+    if (t_geodesic === undefined) {
+        const rFinal = snapshots[snapshots.length - 1].x[1];
+        throw new Error(`evaluateBE37CovariantEikonalNumerical: photon did not reach R_near=${R_near_m} within tauMax=${tauMax} s; final r=${rFinal}. Try increasing steps or tauMax.`);
+    }
+    const shapiroDelaySec = t_geodesic - t_flat;
     return { eikonalResidual, shapiroDelaySec };
 }
 //# sourceMappingURL=be37-covariant-eikonal.js.map

package/dist/numerical/be37-covariant-eikonal.js.map

@@ -1 +1 @@
-{"version":3,"file":"be37-covariant-eikonal.js","sourceRoot":"","sources":["../../src/numerical/be37-covariant-eikonal.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;;;;;;;;GA0BG;AAkDH;;;;;;;;;;;;;;;;;GAiBG;AACH,MAAM,CAAC,KAAK,UAAU,qCAAqC,CACzD,MAAkC;IAElC,MAAM,EAAE,IAAI,EAAE,OAAO,EAAE,QAAQ,EAAE,GAAG,MAAM,CAAC;IAE3C,oEAAoE;IACpE,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,IAAI,CAAC,IAAI,IAAI,IAAI,CAAC,EAAE,CAAC;QACxC,MAAM,IAAI,UAAU,CAClB,qFAAqF,IAAI,EAAE,CAC5F,CAAC;IACJ,CAAC;IACD,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,OAAO,CAAC,IAAI,OAAO,IAAI,CAAC,EAAE,CAAC;QAC9C,MAAM,IAAI,UAAU,CAClB,wFAAwF,OAAO,EAAE,CAClG,CAAC;IACJ,CAAC;IACD,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,QAAQ,CAAC,IAAI,QAAQ,IAAI,CAAC,EAAE,CAAC;QAChD,MAAM,IAAI,UAAU,CAClB,yFAAyF,QAAQ,EAAE,CACpG,CAAC;IACJ,CAAC;IACD,IAAI,QAAQ,GAAG,OAAO,EAAE,CAAC;QACvB,MAAM,IAAI,UAAU,CAClB,oDAAoD,QAAQ,wBAAwB,OAAO,GAAG,CAC/F,CAAC;IACJ,CAAC;IAED,4EAA4E;IAC5E,EAAE;IACF,oEAAoE;IACpE,qEAAqE;IACrE,+CAA+C;IAC/C,EAAE;IACF,qBAAqB;IACrB,EAAE;IACF,0EAA0E;IAC1E,wEAAwE;IACxE,MAAM,eAAe,GAAG,CAAC,CAAC;IAE1B,2EAA2E;IAC3E,EAAE;IACF,+DAA+D;IAC/D,4DAA4D;IAC5D,2DAA2D;IAC3D,MAAM,eAAe,GAAG,CAAC,CAAC;IAE1B,OAAO,EAAE,eAAe,EAAE,eAAe,EAAE,CAAC;AAC9C,CAAC"}
+{"version":3,"file":"be37-covariant-eikonal.js","sourceRoot":"","sources":["../../src/numerical/be37-covariant-eikonal.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GA+BG;AAEH,OAAO,EAAE,oBAAoB,EAAE,MAAM,qBAAqB,CAAC;AAC3D,OAAO,EAAE,IAAI,EAAE,IAAI,EAAE,MAAM,sBAAsB,CAAC;AAElD,6EAA6E;AAC7E,sEAAsE;AACtE,MAAM,IAAI,GAAG,IAAI,CAAC;AAClB,MAAM,KAAK,GAAG,IAAI,GAAG,IAAI,CAAC;AA0E1B;;;;;;;GAOG;AACH,SAAS,eAAe,CAAC,IAAY;IACnC,OAAO,CAAC,CAAC,GAAG,IAAI,GAAG,IAAI,CAAC,GAAG,KAAK,CAAC;AACnC,CAAC;AAED;;;;;;;;GAQG;AACH,SAAS,eAAe,CACtB,IAAY;IAEZ,MAAM,GAAG,GAAG,eAAe,CAAC,IAAI,CAAC,CAAC;IAClC,OAAO,SAAS,QAAQ,CAAC,CAAwB;QAC/C,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;QACf,MAAM,KAAK,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;QACnB,MAAM,CAAC,GAAG,CAAC,GAAG,GAAG,GAAG,CAAC,CAAC;QACtB,MAAM,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC;QAC7B,MAAM,IAAI,GAAe;YACvB,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,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC;QAC9B,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;QACf,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;QACzB,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,IAAI,GAAG,IAAI,CAAC,CAAC;QACvC,OAAO,IAAI,CAAC;IACd,CAAC,CAAC;AACJ,CAAC;AAED;;;;;;;;;GASG;AACH,SAAS,gBAAgB,CACvB,IAAY;IAEZ,MAAM,GAAG,GAAG,eAAe,CAAC,IAAI,CAAC,CAAC;IAClC,OAAO,SAAS,SAAS,CAAC,CAAwB;QAChD,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;QACf,MAAM,KAAK,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;QACnB,MAAM,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC;QAC7B,MAAM,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC;QAC7B,MAAM,CAAC,GAAG,CAAC,GAAG,GAAG,GAAG,CAAC,CAAC;QACtB,MAAM,EAAE,GAAiB,KAAK,CAAC,IAAI,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,EAAE,GAAG,EAAE,CACtD,KAAK,CAAC,IAAI,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,EAAE,GAAG,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,CAC9C,CAAC;QACF,uBAAuB;QACvB,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,GAAG,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,KAAK,CAAC,CAAC;QAC5C,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,GAAG,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;QAC5B,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC;QAC/B,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,IAAI,GAAG,IAAI,CAAC,CAAC;QAC7C,qBAAqB;QACrB,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,IAAI,GAAG,IAAI,GAAG,IAAI,CAAC,CAAC;QACzD,OAAO,EAAE,CAAC;IACZ,CAAC,CAAC;AACJ,CAAC;AAED;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GA6DG;AACH,MAAM,CAAC,KAAK,UAAU,qCAAqC,CACzD,MAAkC;IAElC,MAAM,EAAE,IAAI,EAAE,OAAO,EAAE,QAAQ,EAAE,GAAG,GAAG,CAAC,EAAE,KAAK,GAAG,IAAI,EAAE,GAAG,MAAM,CAAC;IAElE,2EAA2E;IAC3E,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,IAAI,CAAC,IAAI,IAAI,IAAI,CAAC,EAAE,CAAC;QACxC,MAAM,IAAI,UAAU,CAClB,qFAAqF,IAAI,EAAE,CAC5F,CAAC;IACJ,CAAC;IACD,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,OAAO,CAAC,IAAI,OAAO,IAAI,CAAC,EAAE,CAAC;QAC9C,MAAM,IAAI,UAAU,CAClB,wFAAwF,OAAO,EAAE,CAClG,CAAC;IACJ,CAAC;IACD,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,QAAQ,CAAC,IAAI,QAAQ,IAAI,CAAC,EAAE,CAAC;QAChD,MAAM,IAAI,UAAU,CAClB,yFAAyF,QAAQ,EAAE,CACpG,CAAC;IACJ,CAAC;IACD,IAAI,QAAQ,GAAG,OAAO,EAAE,CAAC;QACvB,MAAM,IAAI,UAAU,CAClB,oDAAoD,QAAQ,wBAAwB,OAAO,GAAG,CAC/F,CAAC;IACJ,CAAC;IACD,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,GAAG,CAAC,IAAI,GAAG,GAAG,CAAC,EAAE,CAAC;QACrC,MAAM,IAAI,UAAU,CAClB,wFAAwF,GAAG,EAAE,CAC9F,CAAC;IACJ,CAAC;IACD,IAAI,GAAG,GAAG,QAAQ,EAAE,CAAC;QACnB,MAAM,IAAI,UAAU,CAClB,+CAA+C,GAAG,yBAAyB,QAAQ,4CAA4C,CAChI,CAAC;IACJ,CAAC;IACD,MAAM,GAAG,GAAG,eAAe,CAAC,IAAI,CAAC,CAAC;IAClC,IAAI,QAAQ,IAAI,GAAG,GAAG,IAAI,EAAE,CAAC;QAC3B,MAAM,IAAI,UAAU,CAClB,oDAAoD,QAAQ,0DAA0D,GAAG,EAAE,CAC5H,CAAC;IACJ,CAAC;IAED,4EAA4E;IAC5E,8DAA8D;IAC9D,uEAAuE;IACvE,MAAM,eAAe,GAAG,CAAC,CAAC;IAE1B,2EAA2E;IAC3E,MAAM,UAAU,GAAG,eAAe,CAAC,IAAI,CAAC,CAAC;IACzC,MAAM,WAAW,GAAG,gBAAgB,CAAC,IAAI,CAAC,CAAC;IAE3C,yEAAyE;IACzE,qEAAqE;IACrE,0EAA0E;IAC1E,MAAM,EAAE,GAAa,CAAC,CAAC,EAAE,OAAO,EAAE,IAAI,CAAC,EAAE,GAAG,CAAC,EAAE,CAAC,CAAC,CAAC;IAClD,MAAM,GAAG,GAAG,CAAC,KAAK,CAAC;IACnB,MAAM,KAAK,GAAG,GAAG,GAAG,IAAI,CAAC,CAAC,uEAAuE;IACjG,MAAM,OAAO,GAAG,CAAC,CAAC;IAElB,sEAAsE;IACtE,mDAAmD;IACnD,6CAA6C;IAC7C,MAAM,KAAK,GAAG,UAAU,CAAC,EAAE,CAAC,CAAC;IAC7B,MAAM,SAAS,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,GAAG,GAAG,GAAG,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,GAAG,KAAK,CAAC;IACzE,IAAI,SAAS,GAAG,CAAC,EAAE,CAAC;QAClB,MAAM,IAAI,UAAU,CAClB,kFAAkF,OAAO,wBAAwB,GAAG,2CAA2C,SAAS,GAAG,CAC5K,CAAC;IACJ,CAAC;IACD,MAAM,aAAa,GAAG,IAAI,CAAC,IAAI,CAAC,SAAS,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IACzD,iEAAiE;IACjE,MAAM,GAAG,GAAG,CAAC,aAAa,CAAC;IAE3B,MAAM,YAAY,GAAG;QACnB,CAAC,EAAE,EAAE;QACL,CAAC,EAAE,CAAC,GAAG,EAAE,GAAG,EAAE,OAAO,EAAE,KAAK,CAAC;KAC9B,CAAC;IAEF,uEAAuE;IACvE,yEAAyE;IACzE,iFAAiF;IACjF,MAAM,UAAU,GAAG,IAAI,CAAC,IAAI,CAAC,OAAO,GAAG,OAAO,GAAG,GAAG,GAAG,GAAG,CAAC,CAAC;IAC5D,MAAM,WAAW,GAAG,IAAI,CAAC,IAAI,CAAC,QAAQ,GAAG,QAAQ,GAAG,GAAG,GAAG,GAAG,CAAC,CAAC;IAC/D,MAAM,MAAM,GAAG,CAAC,UAAU,GAAG,WAAW,CAAC,GAAG,IAAI,CAAC;IACjD,MAAM,MAAM,GAAG,MAAM,GAAG,GAAG,CAAC;IAE5B,MAAM,SAAS,GAAG,oBAAoB,CAAC,YAAY,EAAE;QACnD,KAAK;QACL,MAAM;QACN,UAAU;QACV,WAAW;QACX,eAAe,EAAE,GAAG,GAAG,IAAI;KAC5B,CAAC,CAAC;IAEH,wEAAwE;IACxE,IAAI,UAA8B,CAAC;IACnC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,SAAS,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;QAC1C,MAAM,KAAK,GAAG,SAAS,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QACpC,MAAM,KAAK,GAAG,SAAS,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QAChC,IAAI,KAAK,GAAG,QAAQ,IAAI,KAAK,IAAI,QAAQ,EAAE,CAAC;YAC1C,MAAM,KAAK,GAAG,SAAS,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;YACpC,MAAM,KAAK,GAAG,SAAS,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;YAChC,+EAA+E;YAC/E,MAAM,KAAK,GAAG,CAAC,QAAQ,GAAG,KAAK,CAAC,GAAG,CAAC,KAAK,GAAG,KAAK,CAAC,CAAC;YACnD,UAAU,GAAG,KAAK,GAAG,KAAK,GAAG,CAAC,KAAK,GAAG,KAAK,CAAC,CAAC;YAC7C,MAAM;QACR,CAAC;IACH,CAAC;IACD,IAAI,UAAU,KAAK,SAAS,EAAE,CAAC;QAC7B,MAAM,MAAM,GAAG,SAAS,CAAC,SAAS,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QACpD,MAAM,IAAI,KAAK,CACb,sEAAsE,QAAQ,kBAAkB,MAAM,eAAe,MAAM,mCAAmC,CAC/J,CAAC;IACJ,CAAC;IAED,MAAM,eAAe,GAAG,UAAU,GAAG,MAAM,CAAC;IAE5C,OAAO,EAAE,eAAe,EAAE,eAAe,EAAE,CAAC;AAC9C,CAAC"}

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

@@ -0,0 +1,189 @@
+/**
+ * Numerical helpers for Riemann-curvature lowering (Task 6 [U] / v0.5.0 1c-ii).
+ *
+ * Mirrors `connection-lowering-helpers.ts`: produces Γ and ∂Γ tensors from
+ * coordinate-dependent metric closures, then assembles the Riemann tensor via
+ * the corrected (Carroll-Ch.3 §3.4 / Adam+Eve F4-S3) formula:
+ *
+ *   R^ρ_{σμν} = ∂_μ Γ^ρ_{σν} − ∂_ν Γ^ρ_{σμ} + Σ_λ (Γ^ρ_{λμ} Γ^λ_{σν} − Γ^ρ_{λν} Γ^λ_{σμ})
+ *
+ * Index conventions (pinned via JSDoc + a runtime assert):
+ *   - Γ[ρ][σ][ν] = Γ^ρ_{σν}  (matches computeChristoffelTensor's [α][μ][ν] with
+ *     σ ↔ μ — i.e., the SECOND axis is σ, the THIRD is ν, per Carroll Ch.5).
+ *   - dGamma[λ][ρ][σ][ν] = ∂_λ Γ^ρ_{σν}     (I3)
+ *   - R[ρ][σ][μ][ν] = R^ρ_{σμν}              (Carroll)
+ *
+ * NOTE on the σ/μ swap vs `computeChristoffelTensor`:
+ *   `computeChristoffelTensor` was authored for v0.4.0 covariant-derivative
+ *   lowering, which sums Γ^α_{μλ} T^λ over the middle index `μ` (the wrt index
+ *   of the cov-deriv). v0.5.0 Riemann needs σ in the SECOND-lower slot of every
+ *   Γ. Because Γ is symmetric in its two lower indices (Christoffel symmetry),
+ *   the two storage conventions are numerically identical — we use the
+ *   existing [α][middle][last] storage and access it with σ in the middle slot.
+ *
+ * The fast path uses the v0.4.0 `computeChristoffelTensor` infrastructure
+ * (no new tensor-engine code) and the v0.4.0 `pderivNumericalFn` pattern
+ * (centered finite differences). M11: ∂Γ is computed by repeated centered
+ * finite-differences on `christoffel()`, not a new AST node type or AD pass.
+ *
+ * @module numerical/curvature-lowering-helpers
+ */
+import type { TensorEngine } from './tensor-engine.js';
+import type { NestedArray } from './types.js';
+/** Flat row-major N×N matrix (one g or g^{-1} sample). */
+export type FlatMatrix = ReadonlyArray<number>;
+/** A coordinate-dependent rank-2 closure: f(x) → N×N as a NestedArray
+ *  (number[][] or any nested form that `flattenNA` can flatten). */
+export type MetricFn = (x: ReadonlyArray<number>) => NestedArray;
+/**
+ * dGamma[λ][ρ][σ][ν] = ∂_λ Γ^ρ_{σν} at the given coordinates (I3).
+ *
+ * 4-deep nested number array — readonly to mark this as an output container
+ * the helper produces but does not mutate after construction.
+ */
+export type DGammaTensor = readonly (readonly (readonly (readonly number[])[])[])[];
+/** Gamma[ρ][σ][ν] = Γ^ρ_{σν} — same index order as `computeChristoffelTensor`
+ *  output (just renamed for the Riemann use). */
+export type GammaTensor = readonly (readonly (readonly number[])[])[];
+/**
+ * Step size for finite differences on the Christoffel function.
+ *
+ * Inner step (∂g for one Γ evaluation) uses the v0.4.0 default
+ * `1e-6·max(|x|,1)` via `pderivNumericalFn`-style centered FD. This gives ∂g
+ * with relative precision ≈ 1e-11 on Schwarzschild components (verified by
+ * inspection at r=3·r_s).
+ *
+ * Outer step (∂Γ for one R evaluation) is `1e-4·max(|x|,1)` and we use a
+ * **4th-order centered stencil**:
+ *
+ *   f'(x) ≈ [−f(x+2h) + 8 f(x+h) − 8 f(x-h) + f(x-2h)] / (12 h)
+ *
+ * Truncation O(h⁴) instead of O(h²). At h=0.886 m (r=3·r_s, M=M_sun),
+ * truncation ≈ h⁴·Γ⁽⁵⁾/Γ ~ 1e-13 relative, round-off ε·|Γ|/h ~ 1e-15 relative.
+ * Both sit well below the 1e-9 gate.
+ *
+ * The 2nd-order outer FD that v0.4.0 covariant-derivative lowering uses
+ * leaves ~3e-6 relative error on R^t_{rtr} due to noise propagation through
+ * the c²-scaled g_{tt} component — 4th-order outer FD recovers the precision.
+ */
+export declare function outerStep(x: number): number;
+/**
+ * Compute Γ^ρ_{σν} at coordinate x. Internally builds a centered-FD ∂g sampler
+ * around x and delegates to `computeChristoffelTensor` (the v0.4.0 helper).
+ *
+ * Returns a fully-materialised 3-deep nested array Γ[ρ][σ][ν]. (We do NOT keep
+ * the EngineTensor wrapper because we immediately rebuild dGamma and R as
+ * plain nested arrays — round-tripping through the engine on every FD sample
+ * would dominate the cost.)
+ */
+export declare function christoffelAt(x: ReadonlyArray<number>, gFn: MetricFn, gInverseFn: MetricFn, N: number, engine: TensorEngine): number[][][];
+/**
+ * dGamma[λ][ρ][σ][ν] = ∂_λ Γ^ρ_{σν} via centered FD on `christoffelAt`. (I3)
+ *
+ * Includes a runtime index-assert verifying `dGamma[1][1][1][1]` is finite
+ * before returning (catches silent shape-permutation bugs introduced by future
+ * edits to `computeChristoffelTensor`'s output ordering).
+ */
+export declare function dGammaAt(x: ReadonlyArray<number>, gFn: MetricFn, gInverseFn: MetricFn, N: number, engine: TensorEngine): DGammaTensor;
+/**
+ * Build the full 4×4×4×4 Riemann tensor R^ρ_{σμν} from Γ and ∂Γ per the
+ * Carroll formula (Adam+Eve F4-S3):
+ *
+ *   R[ρ][σ][μ][ν] = ∂_μ Γ[ρ][σ][ν] − ∂_ν Γ[ρ][σ][μ]
+ *                 + Σ_λ (Γ[ρ][λ][μ] Γ[λ][σ][ν] − Γ[ρ][λ][ν] Γ[λ][σ][μ])
+ *
+ * σ in the SECOND lower slot of each Γ (F4/S3 correction).
+ */
+export declare function buildRiemann(gamma: GammaTensor, dGamma: DGammaTensor, N: number): number[][][][];
+/**
+ * Compute the upper-Riemann tensor R^ρ_{σμν} at coordinate x using the same
+ * Γ + ∂Γ pipeline as the riemann-tensor lowering case. Encapsulates the
+ * full christoffelAt + dGammaAt + buildRiemann sequence so callers (Task 9
+ * Bianchi residual) can sample R at perturbed coordinates without
+ * re-implementing the FD machinery.
+ */
+export declare function riemannUpperAt(x: ReadonlyArray<number>, gFn: MetricFn, gInverseFn: MetricFn, N: number, engine: TensorEngine): number[][][][];
+/**
+ * Lower the upper-ρ index of R^ρ_{σμν} via the covariant metric g_{aρ}:
+ *
+ *   R_{aσμν} = Σ_ρ g_{aρ} R^ρ_{σμν}
+ *
+ * Output index order: `[a][σ][μ][ν]` — all four lower. (a is the freshly
+ * lowered index in the first slot.)
+ */
+export declare function lowerFirstIndex(R: number[][][][], gLowerFlat: ReadonlyArray<number>, N: number): number[][][][];
+/**
+ * Sample the all-lower Riemann tensor R_{αβγδ}(x) — combines `riemannUpperAt`
+ * with `lowerFirstIndex` for the single-coordinate evaluation. Used both as
+ * the base sample and (with perturbed x) as the FD inputs for ∂_λ R_{αβγδ}.
+ */
+export declare function riemannLowerAt(x: ReadonlyArray<number>, gFn: MetricFn, gInverseFn: MetricFn, N: number, engine: TensorEngine): number[][][][];
+/**
+ * Compute ∂_λ R_{αβγδ}(x) via a 4th-order centered stencil on
+ * `riemannLowerAt`. Index order: `dR[λ][α][β][γ][δ] = ∂_λ R_{αβγδ}`.
+ *
+ * Uses the same outer step (`outerStep`) as Task 6's dGammaAt so the
+ * FD-noise compounding pattern matches. This is one extra layer of FD on top
+ * of the Christoffel-of-Christoffel double-FD already inside `riemannLowerAt`
+ * — total: ∂(∂(∂g)) at 4th order in each layer. Empirical noise floor reached
+ * by the per-component value is documented in the test report.
+ */
+export declare function dRiemannLowerAt(x: ReadonlyArray<number>, gFn: MetricFn, gInverseFn: MetricFn, N: number, engine: TensorEngine): number[][][][][];
+/**
+ * Compute ∇_λ R_{μνρσ} via partial + Christoffel-correction terms (Approach 1).
+ *
+ *   ∇_λ R_{μνρσ} = ∂_λ R_{μνρσ}
+ *                  − Γ^α_{λμ} R_{ανρσ}
+ *                  − Γ^α_{λν} R_{μαρσ}
+ *                  − Γ^α_{λρ} R_{μνασ}
+ *                  − Γ^α_{λσ} R_{μνρα}
+ *
+ * Index order: `covR[λ][μ][ν][ρ][σ] = ∇_λ R_{μνρσ}`.
+ *
+ * Approach choice (full ∇, not raw ∂): the second Bianchi identity
+ * ∇_{[λ} R_{μν]ρσ} = 0 is the canonical statement. Using raw ∂ would give a
+ * residual dominated by the (cyclic-non-cancelling) Christoffel-correction
+ * terms, masking the identity check with O(1) algebraic clutter. Full ∇
+ * makes the test a genuine self-consistency check of the lowered Riemann.
+ */
+export declare function covariantDerivRiemannLowerAt(x: ReadonlyArray<number>, gFn: MetricFn, gInverseFn: MetricFn, N: number, engine: TensorEngine): number[][][][][];
+/**
+ * Build the second-Bianchi-identity residual (cyclic over first three indices):
+ *
+ *   B_{λμνρσ} = ∇_λ R_{μνρσ} + ∇_μ R_{νλρσ} + ∇_ν R_{λμρσ}
+ *
+ * Carroll Eq. 3.95: B ≡ 0 in any (torsion-free) Lorentzian manifold. The
+ * residual measures the FD-truncation + cancellation noise on the lowered
+ * Riemann tensor through one extra ∂ layer.
+ *
+ * Index order: `B[λ][μ][ν][ρ][σ]` — all five lower.
+ */
+export declare function bianchiResidualAt(x: ReadonlyArray<number>, gFn: MetricFn, gInverseFn: MetricFn, N: number, engine: TensorEngine): number[][][][][];
+/**
+ * JS-side contraction of a flat 4-axis tensor R[a][b][c][d] (length N⁴, row-major)
+ * down to a 2-tensor by summing one upper axis against one lower axis.
+ *
+ * Configuration:
+ *   - `upperAxis`: the axis index (0..3) that contracts WITH `lowerAxis` (dummy slot).
+ *   - `lowerAxis`: the partner axis (0..3) forming the dummy pair.
+ *   - `outAxes`:   the two surviving axes [i, j] that become the (free output)
+ *                  rows/cols of the returned N×N matrix.
+ *
+ * For Carroll Eq. 3.91 Ricci R_μν = R^λ_{μλν} on Riemann stored as R[ρ][σ][μ][ν]:
+ *   contractRiemannJS(flatR, N, { upperAxis: 0, lowerAxis: 2, outAxes: [1, 3] })
+ * sums `R[λ][μ_out][λ][ν_out]` over λ.
+ *
+ * Pre-flattened input convention: row-major `R[a][b][c][d] = flatR[a*N³ + b*N² + c*N + d]`.
+ * (Matches `flattenNestedArray` output for an [N,N,N,N] EngineTensor.)
+ *
+ * @internal — extracted from `lowering.ts` ricci-tensor arm (AS-1, v0.5.1). Kept
+ * deliberately narrow (4-axis source → 2-axis output) to match the v0.5.0 GR
+ * curvature pipeline; generalisation deferred to v0.6.0 if Weyl/Kretschmann
+ * arrives (see `docs/architecture/COMPONENTS.md` curvature-layer pattern note).
+ */
+export declare function contractRiemannJS(flatR: ReadonlyArray<number>, N: number, config: {
+    upperAxis: 0 | 1 | 2 | 3;
+    lowerAxis: 0 | 1 | 2 | 3;
+    outAxes: readonly [number, number];
+}): number[][];
+//# sourceMappingURL=curvature-lowering-helpers.d.ts.map

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

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