universal-physics-tensor 0.4.5

package/dist/numerical/be37-covariant-eikonal.d.ts

@@ -0,0 +1,90 @@
+/**
+ * v0.4.0 BE-37 covariant-eikonal STRUCTURAL PREVIEW.
+ *
+ * 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
+ *
+ *   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.
+ *
+ * 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.
+ *
+ * @see tests/dimensional/covariant-derivative-preview.test.ts (Task 17 [U])
+ * @see docs/planning/v0.4.0-Implementation-Plan.md Task 17
+ * @module numerical/be37-covariant-eikonal
+ * @public
+ */
+/**
+ * Input parameters for the BE-37 covariant-eikonal preview evaluator.
+ * @public
+ */
+export interface BE37CovariantEikonalInputs {
+    /** Gravitational source mass in kilograms (must be > 0). */
+    readonly M_kg: number;
+    /** Far-point radius in metres, i.e. the signal origin/destination (must be > 0). */
+    readonly R_far_m: number;
+    /** Near-point radius in metres, i.e. closest approach to the source (must be > 0, ≤ R_far_m). */
+    readonly R_near_m: number;
+}
+/**
+ * Result of `evaluateBE37CovariantEikonalNumerical`.
+ * @public
+ */
+export interface BE37CovariantEikonalResult {
+    /**
+     * Numerical residual of the covariant-eikonal equation
+     * g^μν ∇_μ ∇_ν S = 0 evaluated at the null wave-covector.
+     *
+     * **v0.4.0 structural preview**: returns exactly 0 by construction —
+     * the null wave-covector k_μ = ∂_μ S = (E, E, 0, 0) satisfies
+     * g^μν k_μ k_ν = 0 identically for the Schwarzschild +,−,−,− metric.
+     * No numerical integration is required to establish this; the result
+     * is a structural consequence of the null-ray construction.
+     */
+    readonly eikonalResidual: number;
+    /**
+     * Shapiro time delay in seconds via geodesic integration.
+     *
+     * **v0.4.0 stub — returns 0.**  The full geodesic-integrated Shapiro
+     * cross-check (comparing this field against `evaluateShapiroDelay`) is
+     * the v0.5.0 deliverable; see `it.skip` in
+     * `tests/dimensional/covariant-derivative-preview.test.ts`.
+     */
+    readonly shapiroDelaySec: number;
+    /**
+     * Optional closed-form cross-check value (seconds).
+     * Populated when a closed-form comparator is available.
+     * Not populated in the v0.4.0 structural preview.
+     */
+    readonly closedFormDelaySec?: number;
+}
+/**
+ * v0.4.0 BE-37 covariant-eikonal STRUCTURAL PREVIEW.
+ *
+ * Encodes the eikonal as ∇_μ ∇^μ S using Task 2's `CovariantDerivativeNode`
+ * and Task 12's lowering path.  The numerical eikonal residual is exactly 0
+ * by construction of the null wave-covector (same structural ground as v0.3.5
+ * `evaluateBE37EikonalNumerical`, but now expressed through the covariant
+ * rather than the partial-derivative form).
+ *
+ * **shapiroDelaySec = 0** in v0.4.0 (intentional stub).  The
+ * geodesic-integrated Shapiro cross-check via Task 14's `integrateGeodesic`
+ * is the v0.5.0 deliverable.
+ *
+ * @param inputs - Source mass + far/near radii.
+ * @returns Structural eikonal residual (0 by construction) + stub Shapiro delay.
+ * @public
+ */
+export declare function evaluateBE37CovariantEikonalNumerical(inputs: BE37CovariantEikonalInputs): Promise<BE37CovariantEikonalResult>;
+//# sourceMappingURL=be37-covariant-eikonal.d.ts.map

package/dist/numerical/be37-covariant-eikonal.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"be37-covariant-eikonal.d.ts","sourceRoot":"","sources":["../../src/numerical/be37-covariant-eikonal.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;;;;;;;;GA0BG;AAEH;;;GAGG;AACH,MAAM,WAAW,0BAA0B;IACzC,4DAA4D;IAC5D,QAAQ,CAAC,IAAI,EAAE,MAAM,CAAC;IACtB,oFAAoF;IACpF,QAAQ,CAAC,OAAO,EAAE,MAAM,CAAC;IACzB,iGAAiG;IACjG,QAAQ,CAAC,QAAQ,EAAE,MAAM,CAAC;CAC3B;AAED;;;GAGG;AACH,MAAM,WAAW,0BAA0B;IACzC;;;;;;;;;OASG;IACH,QAAQ,CAAC,eAAe,EAAE,MAAM,CAAC;IAEjC;;;;;;;OAOG;IACH,QAAQ,CAAC,eAAe,EAAE,MAAM,CAAC;IAEjC;;;;OAIG;IACH,QAAQ,CAAC,kBAAkB,CAAC,EAAE,MAAM,CAAC;CACtC;AAED;;;;;;;;;;;;;;;;GAgBG;AACH,wBAAsB,qCAAqC,CACzD,MAAM,EAAE,0BAA0B,GACjC,OAAO,CAAC,0BAA0B,CAAC,CA6CrC"}

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

@@ -0,0 +1,79 @@
+/**
+ * v0.4.0 BE-37 covariant-eikonal STRUCTURAL PREVIEW.
+ *
+ * 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
+ *
+ *   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.
+ *
+ * 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.
+ *
+ * @see tests/dimensional/covariant-derivative-preview.test.ts (Task 17 [U])
+ * @see docs/planning/v0.4.0-Implementation-Plan.md Task 17
+ * @module numerical/be37-covariant-eikonal
+ * @public
+ */
+/**
+ * v0.4.0 BE-37 covariant-eikonal STRUCTURAL PREVIEW.
+ *
+ * Encodes the eikonal as ∇_μ ∇^μ S using Task 2's `CovariantDerivativeNode`
+ * and Task 12's lowering path.  The numerical eikonal residual is exactly 0
+ * by construction of the null wave-covector (same structural ground as v0.3.5
+ * `evaluateBE37EikonalNumerical`, but now expressed through the covariant
+ * rather than the partial-derivative form).
+ *
+ * **shapiroDelaySec = 0** in v0.4.0 (intentional stub).  The
+ * geodesic-integrated Shapiro cross-check via Task 14's `integrateGeodesic`
+ * is the v0.5.0 deliverable.
+ *
+ * @param inputs - Source mass + far/near radii.
+ * @returns Structural eikonal residual (0 by construction) + stub Shapiro delay.
+ * @public
+ */
+export async function evaluateBE37CovariantEikonalNumerical(inputs) {
+    const { M_kg, R_far_m, R_near_m } = inputs;
+    // Domain guards (mirrors evaluateShapiroDelay for API consistency).
+    if (!Number.isFinite(M_kg) || M_kg <= 0) {
+        throw new RangeError(`evaluateBE37CovariantEikonalNumerical: M_kg must be a finite positive number, got ${M_kg}`);
+    }
+    if (!Number.isFinite(R_far_m) || R_far_m <= 0) {
+        throw new RangeError(`evaluateBE37CovariantEikonalNumerical: R_far_m must be a finite positive number, got ${R_far_m}`);
+    }
+    if (!Number.isFinite(R_near_m) || R_near_m <= 0) {
+        throw new RangeError(`evaluateBE37CovariantEikonalNumerical: R_near_m must be a finite positive number, got ${R_near_m}`);
+    }
+    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})`);
+    }
+    // ─── 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.
+    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;
+    return { eikonalResidual, shapiroDelaySec };
+}
+//# sourceMappingURL=be37-covariant-eikonal.js.map

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

@@ -0,0 +1 @@
+{"version":3,"file":"be37-covariant-eikonal.js","sourceRoot":"","sources":["../../src/numerical/be37-covariant-eikonal.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;;;;;;;;GA0BG;AAkDH;;;;;;;;;;;;;;;;GAgBG;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"}

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

@@ -0,0 +1,107 @@
+/**
+ * Numerical helpers for covariant-derivative lowering (Task 12 [U]).
+ *
+ * These functions compute Christoffel symbols and tensor corrections from raw
+ * metric data WITHOUT going through the christoffel() composite AST. The
+ * composite AST is designed for symbolic/dimensional use; numerical lowering
+ * builds Γ directly from metric components and metric-derivative arrays.
+ *
+ * All helpers are engine-agnostic: they operate on flat number[] arrays and
+ * produce EngineTensor results via the engine.fromNested() interface.
+ *
+ * @module numerical/connection-lowering-helpers
+ */
+import type { EngineTensor, TensorEngine } from './tensor-engine.js';
+import type { NestedArray } from './types.js';
+/**
+ * Flatten a NestedArray to a plain number[] in row-major order.
+ * Canonical implementation. `flattenNestedArray` in lowering.ts is NOT a
+ * duplicate — it delegates here and adds an expectedSize invariant check
+ * (throws NumericalBackendError if the flattened count mismatches). This
+ * function deliberately omits that check for callers that do not know the
+ * expected count ahead of time.
+ */
+export declare function flattenNA(data: NestedArray): number[];
+/**
+ * Return a same-shape EngineTensor whose data is all zeros.
+ * S2(c) helper spec: returns a fresh zero tensor of the same shape as t.
+ */
+export declare function zeroTensorLike(t: EngineTensor, engine: TensorEngine): EngineTensor;
+/**
+ * Return an all-zero EngineTensor of the given shape.
+ */
+export declare function zeroTensor(shape: ReadonlyArray<number>, engine: TensorEngine): EngineTensor;
+/**
+ * Convert a flat number[] to a properly-nested NestedArray matching `shape`.
+ * Exported for use in lowering.ts where the engine requires nested structure.
+ */
+export declare function flatToNested(flat: number[], shape: ReadonlyArray<number>): NestedArray;
+/**
+ * Elementwise a + b; throws on shape mismatch.
+ * S2(c) helper spec.
+ */
+export declare function tensorAdd(a: EngineTensor, b: EngineTensor, engine: TensorEngine): EngineTensor;
+/**
+ * Elementwise a[i] + sign * b[i]; throws on shape mismatch.
+ * S2(c) helper spec.
+ */
+export declare function tensorAddScaled(a: EngineTensor, b: EngineTensor, sign: 1 | -1, engine: TensorEngine): EngineTensor;
+/**
+ * Compute the full Christoffel symbol Γ^α_{μν} from raw metric arrays.
+ *
+ *   Γ^α_{μν} = (1/2) g^{αρ} (∂_μ g_{ρν} + ∂_ν g_{ρμ} − ∂_ρ g_{μν})
+ *
+ * Returns an EngineTensor of shape [N, N, N] indexed [α][μ][ν] (row-major).
+ *
+ * @param gInverseFlat  Flat [N*N] array for g^{αρ} (row-major, α outer)
+ * @param getMetricDeriv  Function (mu: number) => flat [N*N] ∂_{mu} g (row-major)
+ * @param N  Dimension
+ * @param engine  Active TensorEngine
+ */
+export declare function computeChristoffelTensor(gInverseFlat: number[], getMetricDeriv: (mu: number) => number[], N: number, engine: TensorEngine): EngineTensor;
+/**
+ * Contract Γ^α_{μλ} or Γ^λ_{μα} with the operand tensor along the axis
+ * corresponding to `freeIdxPos` (the position of the contracted free index
+ * in the operand's flat axis list).
+ *
+ * S2(c) spec — explicit loop body:
+ *
+ * Upper case (∇_μ T^α_β):
+ *   out[α, μ, β, ...] += Γ[α, μ, λ] * T[λ_at_pos, β, ...]
+ *   i.e., contract over λ (axis `freeIdxPos` of T, axis 2 of Γ)
+ *
+ * Lower case (∇_μ T_α_β):
+ *   out[α, μ, β, ...] += Γ[λ, μ, α] * T[λ_at_pos, β, ...]
+ *   i.e., contract over λ (axis `freeIdxPos` of T, axis 0 of Γ)
+ *   and α indexes the fixed axis of Γ at position 2.
+ *
+ * Shape conventions (all row-major):
+ *   Gamma: shape [N, N, N] — [α or λ, μ, ν or α]
+ *   ofFlat: flat [totalSize] with strides computed from `ofShape`
+ *   outputShape: same as the partial-derivative output: [N, ...ofShape]
+ *     (wrt axis prepended, the freeIdxPos axis contracted out,
+ *      all other axes kept — but since ∇_μ output equals ∂_μ output shape,
+ *      we match that shape exactly)
+ *
+ * @param GammaFlat  Flat [N^3] Γ^α_{μν} — [α][μ][ν] row-major
+ * @param ofFlat     Flat operand data, shape `ofShape`
+ * @param ofShape    Shape of the operand tensor
+ * @param freeIdxPos Position (0-based) of the free index being contracted in `ofShape`
+ * @param variance   'upper' → Γ^α_{μλ} contracts λ; 'lower' → Γ^λ_{μα} contracts λ
+ * @param N          Dimension
+ * @param engine     Active TensorEngine
+ *
+ * Returns an EngineTensor of shape [N, ...ofShape_without_freeIdxPos]
+ * i.e., the same shape as the partial-derivative output (wrt axis first, others follow).
+ */
+export declare function contractChristoffelWithOperand(GammaFlat: number[], ofFlat: number[], ofShape: ReadonlyArray<number>, freeIdxPos: number, variance: 'upper' | 'lower', N: number, engine: TensorEngine): EngineTensor;
+/**
+ * Look up ∂_{mu} g from inputs.metricDerivatives.
+ * Key format: `${metricName}/μ_${mu}` (e.g. 'g/μ_0' for ∂_0 g).
+ * Returns a flat [N*N] array.
+ *
+ * Strategy 'zero': returns [N*N] zeros.
+ * Strategy 'supplied': looks up the key; throws if absent.
+ */
+export declare function getMetricDerivFlat(metricName: string, coordLabel: string, mu: number, strategy: 'zero' | 'supplied', N: number, metricDerivatives: ReadonlyMap<string, NestedArray> | undefined): number[];
+//# sourceMappingURL=connection-lowering-helpers.d.ts.map

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

@@ -0,0 +1 @@
+{"version":3,"file":"connection-lowering-helpers.d.ts","sourceRoot":"","sources":["../../src/numerical/connection-lowering-helpers.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;GAYG;AAEH,OAAO,KAAK,EAAE,YAAY,EAAE,YAAY,EAAE,MAAM,oBAAoB,CAAC;AACrE,OAAO,KAAK,EAAE,WAAW,EAAE,MAAM,YAAY,CAAC;AAO9C;;;;;;;GAOG;AACH,wBAAgB,SAAS,CAAC,IAAI,EAAE,WAAW,GAAG,MAAM,EAAE,CASrD;AAsBD;;;GAGG;AACH,wBAAgB,cAAc,CAAC,CAAC,EAAE,YAAY,EAAE,MAAM,EAAE,YAAY,GAAG,YAAY,CAGlF;AAED;;GAEG;AACH,wBAAgB,UAAU,CAAC,KAAK,EAAE,aAAa,CAAC,MAAM,CAAC,EAAE,MAAM,EAAE,YAAY,GAAG,YAAY,CAE3F;AAED;;;GAGG;AACH,wBAAgB,YAAY,CAAC,IAAI,EAAE,MAAM,EAAE,EAAE,KAAK,EAAE,aAAa,CAAC,MAAM,CAAC,GAAG,WAAW,CAEtF;AAaD;;;GAGG;AACH,wBAAgB,SAAS,CAAC,CAAC,EAAE,YAAY,EAAE,CAAC,EAAE,YAAY,EAAE,MAAM,EAAE,YAAY,GAAG,YAAY,CAO9F;AAED;;;GAGG;AACH,wBAAgB,eAAe,CAC7B,CAAC,EAAE,YAAY,EACf,CAAC,EAAE,YAAY,EACf,IAAI,EAAE,CAAC,GAAG,CAAC,CAAC,EACZ,MAAM,EAAE,YAAY,GACnB,YAAY,CAQd;AAaD;;;;;;;;;;;GAWG;AACH,wBAAgB,wBAAwB,CACtC,YAAY,EAAE,MAAM,EAAE,EACtB,cAAc,EAAE,CAAC,EAAE,EAAE,MAAM,KAAK,MAAM,EAAE,EACxC,CAAC,EAAE,MAAM,EACT,MAAM,EAAE,YAAY,GACnB,YAAY,CA4Bd;AAMD;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GAkCG;AACH,wBAAgB,8BAA8B,CAC5C,SAAS,EAAE,MAAM,EAAE,EACnB,MAAM,EAAE,MAAM,EAAE,EAChB,OAAO,EAAE,aAAa,CAAC,MAAM,CAAC,EAC9B,UAAU,EAAE,MAAM,EAClB,QAAQ,EAAE,OAAO,GAAG,OAAO,EAC3B,CAAC,EAAE,MAAM,EACT,MAAM,EAAE,YAAY,GACnB,YAAY,CA0Dd;AA0CD;;;;;;;GAOG;AACH,wBAAgB,kBAAkB,CAChC,UAAU,EAAE,MAAM,EAClB,UAAU,EAAE,MAAM,EAClB,EAAE,EAAE,MAAM,EACV,QAAQ,EAAE,MAAM,GAAG,UAAU,EAC7B,CAAC,EAAE,MAAM,EACT,iBAAiB,EAAE,WAAW,CAAC,MAAM,EAAE,WAAW,CAAC,GAAG,SAAS,GAC9D,MAAM,EAAE,CAsBV"}

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

@@ -0,0 +1,315 @@
+/**
+ * Numerical helpers for covariant-derivative lowering (Task 12 [U]).
+ *
+ * These functions compute Christoffel symbols and tensor corrections from raw
+ * metric data WITHOUT going through the christoffel() composite AST. The
+ * composite AST is designed for symbolic/dimensional use; numerical lowering
+ * builds Γ directly from metric components and metric-derivative arrays.
+ *
+ * All helpers are engine-agnostic: they operate on flat number[] arrays and
+ * produce EngineTensor results via the engine.fromNested() interface.
+ *
+ * @module numerical/connection-lowering-helpers
+ */
+import { NumericalBackendError } from './errors.js';
+// ---------------------------------------------------------------------------
+// Flat-array utilities
+// ---------------------------------------------------------------------------
+/**
+ * Flatten a NestedArray to a plain number[] in row-major order.
+ * Canonical implementation. `flattenNestedArray` in lowering.ts is NOT a
+ * duplicate — it delegates here and adds an expectedSize invariant check
+ * (throws NumericalBackendError if the flattened count mismatches). This
+ * function deliberately omits that check for callers that do not know the
+ * expected count ahead of time.
+ */
+export function flattenNA(data) {
+    if (typeof data === 'number')
+        return [data];
+    const out = [];
+    const walk = (n) => {
+        if (typeof n === 'number')
+            out.push(n);
+        else
+            for (const c of n)
+                walk(c);
+    };
+    walk(data);
+    return out;
+}
+// ---------------------------------------------------------------------------
+// S2(c) helpers — specified by the Task 12 plan
+// ---------------------------------------------------------------------------
+/**
+ * Build a properly-nested zero array matching the given shape.
+ * engine.fromNested() expects nested arrays, not flat arrays, for rank > 1.
+ */
+function buildNestedZeros(shape) {
+    if (shape.length === 0)
+        return 0;
+    if (shape.length === 1)
+        return new Array(shape[0]).fill(0);
+    const inner = buildNestedZeros(shape.slice(1));
+    // Deep-clone each inner node so they are independent arrays.
+    const cloneNested = (n) => {
+        if (typeof n === 'number')
+            return n;
+        return n.map(cloneNested);
+    };
+    return Array.from({ length: shape[0] }, () => cloneNested(inner));
+}
+/**
+ * Return a same-shape EngineTensor whose data is all zeros.
+ * S2(c) helper spec: returns a fresh zero tensor of the same shape as t.
+ */
+export function zeroTensorLike(t, engine) {
+    const shape = [...t.shape];
+    return engine.fromNested(buildNestedZeros(shape), shape);
+}
+/**
+ * Return an all-zero EngineTensor of the given shape.
+ */
+export function zeroTensor(shape, engine) {
+    return engine.fromNested(buildNestedZeros([...shape]), [...shape]);
+}
+/**
+ * Convert a flat number[] to a properly-nested NestedArray matching `shape`.
+ * Exported for use in lowering.ts where the engine requires nested structure.
+ */
+export function flatToNested(flat, shape) {
+    return flatToNestedImpl(flat, shape);
+}
+function flatToNestedImpl(flat, shape) {
+    if (shape.length === 0)
+        return flat[0] ?? 0;
+    if (shape.length === 1)
+        return flat.slice();
+    const stride = flat.length / shape[0];
+    const result = [];
+    for (let i = 0; i < shape[0]; i++) {
+        result.push(flatToNestedImpl(flat.slice(i * stride, (i + 1) * stride), shape.slice(1)));
+    }
+    return result;
+}
+/**
+ * Elementwise a + b; throws on shape mismatch.
+ * S2(c) helper spec.
+ */
+export function tensorAdd(a, b, engine) {
+    if (!sameShape(a.shape, b.shape)) {
+        throw new NumericalBackendError(`tensorAdd: shape mismatch [${a.shape}] vs [${b.shape}]`);
+    }
+    return engine.add(a, b);
+}
+/**
+ * Elementwise a[i] + sign * b[i]; throws on shape mismatch.
+ * S2(c) helper spec.
+ */
+export function tensorAddScaled(a, b, sign, engine) {
+    if (!sameShape(a.shape, b.shape)) {
+        throw new NumericalBackendError(`tensorAddScaled: shape mismatch [${a.shape}] vs [${b.shape}]`);
+    }
+    if (sign === 1)
+        return engine.add(a, b);
+    return engine.sub(a, b);
+}
+function sameShape(a, b) {
+    return a.length === b.length && a.every((v, i) => v === b[i]);
+}
+// ---------------------------------------------------------------------------
+// Christoffel numerical computation
+// ---------------------------------------------------------------------------
+/**
+ * Compute the full Christoffel symbol Γ^α_{μν} from raw metric arrays.
+ *
+ *   Γ^α_{μν} = (1/2) g^{αρ} (∂_μ g_{ρν} + ∂_ν g_{ρμ} − ∂_ρ g_{μν})
+ *
+ * Returns an EngineTensor of shape [N, N, N] indexed [α][μ][ν] (row-major).
+ *
+ * @param gInverseFlat  Flat [N*N] array for g^{αρ} (row-major, α outer)
+ * @param getMetricDeriv  Function (mu: number) => flat [N*N] ∂_{mu} g (row-major)
+ * @param N  Dimension
+ * @param engine  Active TensorEngine
+ */
+export function computeChristoffelTensor(gInverseFlat, getMetricDeriv, N, engine) {
+    // Output: Gamma[alpha][mu][nu] stored in row-major order
+    const size = N * N * N;
+    const Gamma = new Array(size).fill(0);
+    for (let alpha = 0; alpha < N; alpha++) {
+        for (let mu = 0; mu < N; mu++) {
+            for (let nu = 0; nu < N; nu++) {
+                let sum = 0;
+                const dmu = getMetricDeriv(mu); // ∂_μ g_{ρν} — shape [N,N], ρ outer
+                const dnu = getMetricDeriv(nu); // ∂_ν g_{ρμ}
+                for (let rho = 0; rho < N; rho++) {
+                    const gInvAlphaRho = gInverseFlat[alpha * N + rho]; // g^{α ρ}
+                    const drho = getMetricDeriv(rho); // ∂_ρ g_{μν}
+                    // ∂_μ g_{ρν} → index [ρ,ν] in flat ∂_μ g
+                    const term1 = dmu[rho * N + nu];
+                    // ∂_ν g_{ρμ} → index [ρ,μ] in flat ∂_ν g
+                    const term2 = dnu[rho * N + mu];
+                    // ∂_ρ g_{μν} → index [μ,ν] in flat ∂_ρ g
+                    const term3 = drho[mu * N + nu];
+                    sum += gInvAlphaRho * (term1 + term2 - term3);
+                }
+                Gamma[alpha * N * N + mu * N + nu] = 0.5 * sum;
+            }
+        }
+    }
+    return engine.fromNested(flatToNestedImpl(Gamma, [N, N, N]), [N, N, N]);
+}
+// ---------------------------------------------------------------------------
+// Contraction of Γ with the operand tensor
+// ---------------------------------------------------------------------------
+/**
+ * Contract Γ^α_{μλ} or Γ^λ_{μα} with the operand tensor along the axis
+ * corresponding to `freeIdxPos` (the position of the contracted free index
+ * in the operand's flat axis list).
+ *
+ * S2(c) spec — explicit loop body:
+ *
+ * Upper case (∇_μ T^α_β):
+ *   out[α, μ, β, ...] += Γ[α, μ, λ] * T[λ_at_pos, β, ...]
+ *   i.e., contract over λ (axis `freeIdxPos` of T, axis 2 of Γ)
+ *
+ * Lower case (∇_μ T_α_β):
+ *   out[α, μ, β, ...] += Γ[λ, μ, α] * T[λ_at_pos, β, ...]
+ *   i.e., contract over λ (axis `freeIdxPos` of T, axis 0 of Γ)
+ *   and α indexes the fixed axis of Γ at position 2.
+ *
+ * Shape conventions (all row-major):
+ *   Gamma: shape [N, N, N] — [α or λ, μ, ν or α]
+ *   ofFlat: flat [totalSize] with strides computed from `ofShape`
+ *   outputShape: same as the partial-derivative output: [N, ...ofShape]
+ *     (wrt axis prepended, the freeIdxPos axis contracted out,
+ *      all other axes kept — but since ∇_μ output equals ∂_μ output shape,
+ *      we match that shape exactly)
+ *
+ * @param GammaFlat  Flat [N^3] Γ^α_{μν} — [α][μ][ν] row-major
+ * @param ofFlat     Flat operand data, shape `ofShape`
+ * @param ofShape    Shape of the operand tensor
+ * @param freeIdxPos Position (0-based) of the free index being contracted in `ofShape`
+ * @param variance   'upper' → Γ^α_{μλ} contracts λ; 'lower' → Γ^λ_{μα} contracts λ
+ * @param N          Dimension
+ * @param engine     Active TensorEngine
+ *
+ * Returns an EngineTensor of shape [N, ...ofShape_without_freeIdxPos]
+ * i.e., the same shape as the partial-derivative output (wrt axis first, others follow).
+ */
+export function contractChristoffelWithOperand(GammaFlat, ofFlat, ofShape, freeIdxPos, variance, N, engine) {
+    // Output shape: [N (wrt/mu axis), ...ofShape with freeIdxPos axis removed and
+    //               replaced by the alpha axis of the Christoffel]
+    // Concretely: output has the same shape as ∂_μ T which is [...ofShape, N]
+    // with wrtIndex appended. But freeIdx is now contracted (replaced by α):
+    // upper: T had axis freeIdxPos as α (upper) → output keeps α at same position, prepends μ
+    // lower: T had axis freeIdxPos as α (lower) → output keeps α at same position, prepends μ
+    //
+    // The output axis order matches ∂_μ T: [...ofShape, N] (wrt axis last? or first?)
+    // Actually ∂_μ appends wrt as the LAST axis per v0.3.5 design.
+    // So output shape of ∇_μ T = [...ofShape, N] where N is the wrt axis.
+    // But the correction term must have the same shape.
+    //
+    // For the correction: we keep all of T's axes, with freeIdxPos contracting,
+    // and add a new μ axis (last, like ∂_μ). The α index replaces the contracted
+    // position (same position in the output as freeIdxPos was in T).
+    //
+    // Output shape = [...ofShape, N] where ofShape[freeIdxPos] stays (α replaces λ)
+    // and the last axis N is the μ axis.
+    const rank = ofShape.length;
+    const outShape = [...ofShape, N]; // [...T_shape, mu]
+    const outSize = outShape.reduce((a, b) => a * b, 1);
+    const out = new Array(outSize).fill(0);
+    // Build strides for ofFlat and out
+    const ofStrides = buildStrides(ofShape);
+    const outStrides = buildStrides(outShape);
+    // Iterate over all output indices
+    forEachMultiIndex(outShape, (outIdx) => {
+        const mu = outIdx[rank]; // last axis of output is μ
+        // The other axes map to T's indices, with freeIdxPos being the α/result axis
+        const alpha = outIdx[freeIdxPos]; // the "free" index value (α for upper, α for lower)
+        let sum = 0;
+        for (let lam = 0; lam < N; lam++) {
+            // Build T's index: same as outIdx[0..rank-1] but with freeIdxPos = lam
+            const tIdx = outIdx.slice(0, rank);
+            tIdx[freeIdxPos] = lam;
+            const tFlat = flatIndex(tIdx, ofStrides);
+            let gammaFlat;
+            if (variance === 'upper') {
+                // Γ^α_{μλ} — shape [α][μ][λ]
+                gammaFlat = alpha * N * N + mu * N + lam;
+            }
+            else {
+                // Γ^λ_{μα} — shape [λ][μ][α]
+                gammaFlat = lam * N * N + mu * N + alpha;
+            }
+            sum += GammaFlat[gammaFlat] * ofFlat[tFlat];
+        }
+        out[flatIndex(outIdx, outStrides)] = sum;
+    });
+    return engine.fromNested(flatToNestedImpl(out, outShape), outShape);
+}
+// ---------------------------------------------------------------------------
+// Internal utilities
+// ---------------------------------------------------------------------------
+function buildStrides(shape) {
+    const strides = new Array(shape.length);
+    let s = 1;
+    for (let k = shape.length - 1; k >= 0; k--) {
+        strides[k] = s;
+        s *= shape[k];
+    }
+    return strides;
+}
+function flatIndex(idx, strides) {
+    let f = 0;
+    for (let k = 0; k < idx.length; k++)
+        f += idx[k] * strides[k];
+    return f;
+}
+function forEachMultiIndex(shape, visit) {
+    if (shape.length === 0) {
+        visit([]);
+        return;
+    }
+    const idx = new Array(shape.length).fill(0);
+    const total = shape.reduce((a, b) => a * b, 1);
+    for (let n = 0; n < total; n++) {
+        visit([...idx]);
+        for (let k = shape.length - 1; k >= 0; k--) {
+            if (++idx[k] < shape[k])
+                break;
+            idx[k] = 0;
+        }
+    }
+}
+// ---------------------------------------------------------------------------
+// Metric-partial-derivative lookup (for 'supplied' strategy)
+// ---------------------------------------------------------------------------
+/**
+ * Look up ∂_{mu} g from inputs.metricDerivatives.
+ * Key format: `${metricName}/μ_${mu}` (e.g. 'g/μ_0' for ∂_0 g).
+ * Returns a flat [N*N] array.
+ *
+ * Strategy 'zero': returns [N*N] zeros.
+ * Strategy 'supplied': looks up the key; throws if absent.
+ */
+export function getMetricDerivFlat(metricName, coordLabel, mu, strategy, N, metricDerivatives) {
+    if (strategy === 'zero') {
+        return new Array(N * N).fill(0);
+    }
+    // strategy === 'supplied'
+    const key = `${metricName}/${coordLabel}_${mu}`;
+    const d = metricDerivatives?.get(key);
+    if (d === undefined) {
+        throw new NumericalBackendError(`connection-lowering: no metric derivative supplied for "${key}" — ` +
+            `a covariant-derivative with derivativeStrategy='supplied' requires ` +
+            `inputs.metricDerivatives to contain "${key}" (flat [N×N] of ∂_{${mu}} g)`);
+    }
+    const flat = flattenNA(d);
+    if (flat.length !== N * N) {
+        throw new NumericalBackendError(`connection-lowering: metricDerivatives["${key}"] has ${flat.length} elements, ` +
+            `expected ${N * N} (N=${N})`);
+    }
+    return flat;
+}
+//# sourceMappingURL=connection-lowering-helpers.js.map

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

@@ -0,0 +1 @@
+{"version":3,"file":"connection-lowering-helpers.js","sourceRoot":"","sources":["../../src/numerical/connection-lowering-helpers.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;GAYG;AAIH,OAAO,EAAE,qBAAqB,EAAE,MAAM,aAAa,CAAC;AAEpD,8EAA8E;AAC9E,uBAAuB;AACvB,8EAA8E;AAE9E;;;;;;;GAOG;AACH,MAAM,UAAU,SAAS,CAAC,IAAiB;IACzC,IAAI,OAAO,IAAI,KAAK,QAAQ;QAAE,OAAO,CAAC,IAAI,CAAC,CAAC;IAC5C,MAAM,GAAG,GAAa,EAAE,CAAC;IACzB,MAAM,IAAI,GAAG,CAAC,CAAc,EAAQ,EAAE;QACpC,IAAI,OAAO,CAAC,KAAK,QAAQ;YAAE,GAAG,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;;YAClC,KAAK,MAAM,CAAC,IAAI,CAAC;gBAAE,IAAI,CAAC,CAAC,CAAC,CAAC;IAClC,CAAC,CAAC;IACF,IAAI,CAAC,IAAI,CAAC,CAAC;IACX,OAAO,GAAG,CAAC;AACb,CAAC;AAED,8EAA8E;AAC9E,gDAAgD;AAChD,8EAA8E;AAE9E;;;GAGG;AACH,SAAS,gBAAgB,CAAC,KAA4B;IACpD,IAAI,KAAK,CAAC,MAAM,KAAK,CAAC;QAAE,OAAO,CAAC,CAAC;IACjC,IAAI,KAAK,CAAC,MAAM,KAAK,CAAC;QAAE,OAAO,IAAI,KAAK,CAAS,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;IACnE,MAAM,KAAK,GAAG,gBAAgB,CAAC,KAAK,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC;IAC/C,6DAA6D;IAC7D,MAAM,WAAW,GAAG,CAAC,CAAc,EAAe,EAAE;QAClD,IAAI,OAAO,CAAC,KAAK,QAAQ;YAAE,OAAO,CAAC,CAAC;QACpC,OAAQ,CAAmB,CAAC,GAAG,CAAC,WAAW,CAAC,CAAC;IAC/C,CAAC,CAAC;IACF,OAAO,KAAK,CAAC,IAAI,CAAC,EAAE,MAAM,EAAE,KAAK,CAAC,CAAC,CAAC,EAAE,EAAE,GAAG,EAAE,CAAC,WAAW,CAAC,KAAK,CAAC,CAAC,CAAC;AACpE,CAAC;AAED;;;GAGG;AACH,MAAM,UAAU,cAAc,CAAC,CAAe,EAAE,MAAoB;IAClE,MAAM,KAAK,GAAG,CAAC,GAAG,CAAC,CAAC,KAAK,CAAC,CAAC;IAC3B,OAAO,MAAM,CAAC,UAAU,CAAC,gBAAgB,CAAC,KAAK,CAAC,EAAE,KAAK,CAAC,CAAC;AAC3D,CAAC;AAED;;GAEG;AACH,MAAM,UAAU,UAAU,CAAC,KAA4B,EAAE,MAAoB;IAC3E,OAAO,MAAM,CAAC,UAAU,CAAC,gBAAgB,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,EAAE,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC;AACrE,CAAC;AAED;;;GAGG;AACH,MAAM,UAAU,YAAY,CAAC,IAAc,EAAE,KAA4B;IACvE,OAAO,gBAAgB,CAAC,IAAI,EAAE,KAAK,CAAC,CAAC;AACvC,CAAC;AAED,SAAS,gBAAgB,CAAC,IAAc,EAAE,KAA4B;IACpE,IAAI,KAAK,CAAC,MAAM,KAAK,CAAC;QAAE,OAAO,IAAI,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC;IAC5C,IAAI,KAAK,CAAC,MAAM,KAAK,CAAC;QAAE,OAAO,IAAI,CAAC,KAAK,EAAE,CAAC;IAC5C,MAAM,MAAM,GAAG,IAAI,CAAC,MAAM,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC;IACtC,MAAM,MAAM,GAAkB,EAAE,CAAC;IACjC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;QAClC,MAAM,CAAC,IAAI,CAAC,gBAAgB,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC,GAAG,MAAM,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,MAAM,CAAC,EAAE,KAAK,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IAC1F,CAAC;IACD,OAAO,MAAM,CAAC;AAChB,CAAC;AAED;;;GAGG;AACH,MAAM,UAAU,SAAS,CAAC,CAAe,EAAE,CAAe,EAAE,MAAoB;IAC9E,IAAI,CAAC,SAAS,CAAC,CAAC,CAAC,KAAK,EAAE,CAAC,CAAC,KAAK,CAAC,EAAE,CAAC;QACjC,MAAM,IAAI,qBAAqB,CAC7B,8BAA8B,CAAC,CAAC,KAAK,SAAS,CAAC,CAAC,KAAK,GAAG,CACzD,CAAC;IACJ,CAAC;IACD,OAAO,MAAM,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;AAC1B,CAAC;AAED;;;GAGG;AACH,MAAM,UAAU,eAAe,CAC7B,CAAe,EACf,CAAe,EACf,IAAY,EACZ,MAAoB;IAEpB,IAAI,CAAC,SAAS,CAAC,CAAC,CAAC,KAAK,EAAE,CAAC,CAAC,KAAK,CAAC,EAAE,CAAC;QACjC,MAAM,IAAI,qBAAqB,CAC7B,oCAAoC,CAAC,CAAC,KAAK,SAAS,CAAC,CAAC,KAAK,GAAG,CAC/D,CAAC;IACJ,CAAC;IACD,IAAI,IAAI,KAAK,CAAC;QAAE,OAAO,MAAM,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;IACxC,OAAO,MAAM,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;AAC1B,CAAC;AAED,SAAS,SAAS,CAChB,CAAwB,EACxB,CAAwB;IAExB,OAAO,CAAC,CAAC,MAAM,KAAK,CAAC,CAAC,MAAM,IAAI,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;AAChE,CAAC;AAED,8EAA8E;AAC9E,oCAAoC;AACpC,8EAA8E;AAE9E;;;;;;;;;;;GAWG;AACH,MAAM,UAAU,wBAAwB,CACtC,YAAsB,EACtB,cAAwC,EACxC,CAAS,EACT,MAAoB;IAEpB,yDAAyD;IACzD,MAAM,IAAI,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;IACvB,MAAM,KAAK,GAAG,IAAI,KAAK,CAAS,IAAI,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;IAE9C,KAAK,IAAI,KAAK,GAAG,CAAC,EAAE,KAAK,GAAG,CAAC,EAAE,KAAK,EAAE,EAAE,CAAC;QACvC,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,CAAC,EAAE,EAAE,EAAE,EAAE,CAAC;YAC9B,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,CAAC,EAAE,EAAE,EAAE,EAAE,CAAC;gBAC9B,IAAI,GAAG,GAAG,CAAC,CAAC;gBACZ,MAAM,GAAG,GAAG,cAAc,CAAC,EAAE,CAAC,CAAC,CAAE,oCAAoC;gBACrE,MAAM,GAAG,GAAG,cAAc,CAAC,EAAE,CAAC,CAAC,CAAE,aAAa;gBAC9C,KAAK,IAAI,GAAG,GAAG,CAAC,EAAE,GAAG,GAAG,CAAC,EAAE,GAAG,EAAE,EAAE,CAAC;oBACjC,MAAM,YAAY,GAAG,YAAY,CAAC,KAAK,GAAG,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,UAAU;oBAC9D,MAAM,IAAI,GAAG,cAAc,CAAC,GAAG,CAAC,CAAC,CAAmB,aAAa;oBACjE,yCAAyC;oBACzC,MAAM,KAAK,GAAG,GAAG,CAAC,GAAG,GAAG,CAAC,GAAG,EAAE,CAAC,CAAC;oBAChC,yCAAyC;oBACzC,MAAM,KAAK,GAAG,GAAG,CAAC,GAAG,GAAG,CAAC,GAAG,EAAE,CAAC,CAAC;oBAChC,yCAAyC;oBACzC,MAAM,KAAK,GAAG,IAAI,CAAC,EAAE,GAAG,CAAC,GAAG,EAAE,CAAC,CAAC;oBAChC,GAAG,IAAI,YAAY,GAAG,CAAC,KAAK,GAAG,KAAK,GAAG,KAAK,CAAC,CAAC;gBAChD,CAAC;gBACD,KAAK,CAAC,KAAK,GAAG,CAAC,GAAG,CAAC,GAAG,EAAE,GAAG,CAAC,GAAG,EAAE,CAAC,GAAG,GAAG,GAAG,GAAG,CAAC;YACjD,CAAC;QACH,CAAC;IACH,CAAC;IAED,OAAO,MAAM,CAAC,UAAU,CAAC,gBAAgB,CAAC,KAAK,EAAE,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;AAC1E,CAAC;AAED,8EAA8E;AAC9E,2CAA2C;AAC3C,8EAA8E;AAE9E;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GAkCG;AACH,MAAM,UAAU,8BAA8B,CAC5C,SAAmB,EACnB,MAAgB,EAChB,OAA8B,EAC9B,UAAkB,EAClB,QAA2B,EAC3B,CAAS,EACT,MAAoB;IAEpB,8EAA8E;IAC9E,+DAA+D;IAC/D,0EAA0E;IAC1E,yEAAyE;IACzE,0FAA0F;IAC1F,0FAA0F;IAC1F,EAAE;IACF,kFAAkF;IAClF,+DAA+D;IAC/D,sEAAsE;IACtE,oDAAoD;IACpD,EAAE;IACF,4EAA4E;IAC5E,6EAA6E;IAC7E,iEAAiE;IACjE,EAAE;IACF,gFAAgF;IAChF,qCAAqC;IAErC,MAAM,IAAI,GAAG,OAAO,CAAC,MAAM,CAAC;IAC5B,MAAM,QAAQ,GAAG,CAAC,GAAG,OAAO,EAAE,CAAC,CAAC,CAAC,CAAC,mBAAmB;IACrD,MAAM,OAAO,GAAG,QAAQ,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,CAAC;IACpD,MAAM,GAAG,GAAG,IAAI,KAAK,CAAS,OAAO,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;IAE/C,mCAAmC;IACnC,MAAM,SAAS,GAAG,YAAY,CAAC,OAAO,CAAC,CAAC;IACxC,MAAM,UAAU,GAAG,YAAY,CAAC,QAAQ,CAAC,CAAC;IAE1C,kCAAkC;IAClC,iBAAiB,CAAC,QAAQ,EAAE,CAAC,MAAM,EAAE,EAAE;QACrC,MAAM,EAAE,GAAG,MAAM,CAAC,IAAI,CAAC,CAAC,CAAC,2BAA2B;QACpD,6EAA6E;QAC7E,MAAM,KAAK,GAAG,MAAM,CAAC,UAAU,CAAC,CAAC,CAAC,oDAAoD;QAEtF,IAAI,GAAG,GAAG,CAAC,CAAC;QACZ,KAAK,IAAI,GAAG,GAAG,CAAC,EAAE,GAAG,GAAG,CAAC,EAAE,GAAG,EAAE,EAAE,CAAC;YACjC,uEAAuE;YACvE,MAAM,IAAI,GAAG,MAAM,CAAC,KAAK,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC;YACnC,IAAI,CAAC,UAAU,CAAC,GAAG,GAAG,CAAC;YACvB,MAAM,KAAK,GAAG,SAAS,CAAC,IAAI,EAAE,SAAS,CAAC,CAAC;YAEzC,IAAI,SAAiB,CAAC;YACtB,IAAI,QAAQ,KAAK,OAAO,EAAE,CAAC;gBACzB,6BAA6B;gBAC7B,SAAS,GAAG,KAAK,GAAG,CAAC,GAAG,CAAC,GAAG,EAAE,GAAG,CAAC,GAAG,GAAG,CAAC;YAC3C,CAAC;iBAAM,CAAC;gBACN,6BAA6B;gBAC7B,SAAS,GAAG,GAAG,GAAG,CAAC,GAAG,CAAC,GAAG,EAAE,GAAG,CAAC,GAAG,KAAK,CAAC;YAC3C,CAAC;YAED,GAAG,IAAI,SAAS,CAAC,SAAS,CAAC,GAAG,MAAM,CAAC,KAAK,CAAC,CAAC;QAC9C,CAAC;QAED,GAAG,CAAC,SAAS,CAAC,MAAM,EAAE,UAAU,CAAC,CAAC,GAAG,GAAG,CAAC;IAC3C,CAAC,CAAC,CAAC;IAEH,OAAO,MAAM,CAAC,UAAU,CAAC,gBAAgB,CAAC,GAAG,EAAE,QAAQ,CAAC,EAAE,QAAQ,CAAC,CAAC;AACtE,CAAC;AAED,8EAA8E;AAC9E,qBAAqB;AACrB,8EAA8E;AAE9E,SAAS,YAAY,CAAC,KAA4B;IAChD,MAAM,OAAO,GAAG,IAAI,KAAK,CAAS,KAAK,CAAC,MAAM,CAAC,CAAC;IAChD,IAAI,CAAC,GAAG,CAAC,CAAC;IACV,KAAK,IAAI,CAAC,GAAG,KAAK,CAAC,MAAM,GAAG,CAAC,EAAE,CAAC,IAAI,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;QAC3C,OAAO,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;QACf,CAAC,IAAI,KAAK,CAAC,CAAC,CAAC,CAAC;IAChB,CAAC;IACD,OAAO,OAAO,CAAC;AACjB,CAAC;AAED,SAAS,SAAS,CAAC,GAA0B,EAAE,OAA8B;IAC3E,IAAI,CAAC,GAAG,CAAC,CAAC;IACV,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,GAAG,CAAC,MAAM,EAAE,CAAC,EAAE;QAAE,CAAC,IAAI,GAAG,CAAC,CAAC,CAAC,GAAG,OAAO,CAAC,CAAC,CAAC,CAAC;IAC9D,OAAO,CAAC,CAAC;AACX,CAAC;AAED,SAAS,iBAAiB,CACxB,KAA4B,EAC5B,KAA8B;IAE9B,IAAI,KAAK,CAAC,MAAM,KAAK,CAAC,EAAE,CAAC;QAAC,KAAK,CAAC,EAAE,CAAC,CAAC;QAAC,OAAO;IAAC,CAAC;IAC9C,MAAM,GAAG,GAAG,IAAI,KAAK,CAAS,KAAK,CAAC,MAAM,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;IACpD,MAAM,KAAK,GAAG,KAAK,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,CAAC;IAC/C,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,KAAK,EAAE,CAAC,EAAE,EAAE,CAAC;QAC/B,KAAK,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC;QAChB,KAAK,IAAI,CAAC,GAAG,KAAK,CAAC,MAAM,GAAG,CAAC,EAAE,CAAC,IAAI,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;YAC3C,IAAI,EAAE,GAAG,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC;gBAAE,MAAM;YAC/B,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;QACb,CAAC;IACH,CAAC;AACH,CAAC;AAED,8EAA8E;AAC9E,6DAA6D;AAC7D,8EAA8E;AAE9E;;;;;;;GAOG;AACH,MAAM,UAAU,kBAAkB,CAChC,UAAkB,EAClB,UAAkB,EAClB,EAAU,EACV,QAA6B,EAC7B,CAAS,EACT,iBAA+D;IAE/D,IAAI,QAAQ,KAAK,MAAM,EAAE,CAAC;QACxB,OAAO,IAAI,KAAK,CAAS,CAAC,GAAG,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;IAC1C,CAAC;IACD,0BAA0B;IAC1B,MAAM,GAAG,GAAG,GAAG,UAAU,IAAI,UAAU,IAAI,EAAE,EAAE,CAAC;IAChD,MAAM,CAAC,GAAG,iBAAiB,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC;IACtC,IAAI,CAAC,KAAK,SAAS,EAAE,CAAC;QACpB,MAAM,IAAI,qBAAqB,CAC7B,2DAA2D,GAAG,MAAM;YACpE,qEAAqE;YACrE,wCAAwC,GAAG,uBAAuB,EAAE,MAAM,CAC3E,CAAC;IACJ,CAAC;IACD,MAAM,IAAI,GAAG,SAAS,CAAC,CAAC,CAAC,CAAC;IAC1B,IAAI,IAAI,CAAC,MAAM,KAAK,CAAC,GAAG,CAAC,EAAE,CAAC;QAC1B,MAAM,IAAI,qBAAqB,CAC7B,2CAA2C,GAAG,UAAU,IAAI,CAAC,MAAM,aAAa;YAChF,YAAY,CAAC,GAAG,CAAC,OAAO,CAAC,GAAG,CAC7B,CAAC;IACJ,CAAC;IACD,OAAO,IAAI,CAAC;AACd,CAAC"}