universal-physics-tensor 0.4.5 → 0.5.0

package/dist/index.d.ts

@@ -13,6 +13,12 @@ export type { BridgeEquationEntry, BridgeEquationStatus, BridgeIssueSeverity, Br
 export { evaluateGravitationalLensing, type GravitationalLensingInputs, type GravitationalLensingResult, evaluatePerihelionPrecession, type PerihelionPrecessionInputs, type PerihelionPrecessionResult, } from './bridges/index.js';
 export { christoffel } from './dimensional/connection.js';
 export type { CovariantDerivativeNode } from './dimensional/validator.js';
+export { ricci } from './dimensional/curvature.js';
+export type { RicciTensorNode } from './dimensional/validator.js';
+export { einstein } from './dimensional/curvature.js';
+export type { EinsteinTensorNode } from './dimensional/validator.js';
+export { bianchiResidual } from './dimensional/curvature.js';
+export type { BianchiResidualNode } from './dimensional/validator.js';
 export { integrateGeodesic, type GeodesicIntegratorInputs, type GeodesicIntegratorResult, } from './numerical/geodesic-integrator.js';
 export type { Dimension } from './dimensional/types.js';
 export { DIMENSIONLESS, LENGTH, AREA, TIME, FREQUENCY, MASS, VELOCITY, ACCELERATION, FORCE, ENERGY, POWER, ACTION, TEMPERATURE, ENTROPY, CHARGE, } from './dimensional/types.js';
@@ -20,6 +26,6 @@ export { multiply, divide, power, add, subtract, equals, format, DimensionMismat
 export type { ExprNode, ValidationResult, Violation } from './dimensional/validator.js';
 export { validate, validateEquation, validateInverseMetricPair } from './dimensional/validator.js';
 export { inferDimensionForBridge } from './dimensional/bridge-check.js';
-export { evaluateNumerical, evaluateNumericalRaw, evaluateMetricInverse, Float64ReferenceEngine, getActiveEngine, setActiveEngine, NumericalBackendError, DuplicateCoordinateWarning, EngineCapabilityError, hasAutogradSupport, evaluateBE37CovariantEikonalNumerical, } from './numerical/index.js';
-export type { NumericalResult, NumericalRawResult, EvaluateOptions, NumericalInputs, TensorEngine, EngineTensor, EinsumSpec, NestedArray, GridField, ForwardGradResult, ReverseGradResult, } from './numerical/index.js';
+export { evaluateNumerical, evaluateNumericalRaw, evaluateMetricInverse, Float64ReferenceEngine, getActiveEngine, setActiveEngine, NumericalBackendError, DuplicateCoordinateWarning, EngineCapabilityError, hasAutogradSupport, evaluateBE37CovariantEikonalNumerical, integrateGeodesicGL4, findPerihelion, } from './numerical/index.js';
+export type { NumericalResult, NumericalRawResult, EvaluateOptions, NumericalInputs, TensorEngine, EngineTensor, EinsumSpec, NestedArray, GridField, ForwardGradResult, ReverseGradResult, GL4State, GL4Snapshot, GL4Options, PerihelionResult, FindPerihelionOptions, } from './numerical/index.js';
 //# sourceMappingURL=index.d.ts.map

package/dist/index.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"index.d.ts","sourceRoot":"","sources":["../src/index.ts"],"names":[],"mappings":"AAAA;;;;;;GAMG;AAEH,OAAO,EAAE,eAAe,EAAE,MAAM,kBAAkB,CAAC;AACnD,YAAY,EACV,YAAY,EACZ,aAAa,EACb,WAAW,EACX,cAAc,EACd,kBAAkB,EAClB,aAAa,EACb,KAAK,EACL,QAAQ,EACR,kBAAkB,GACnB,MAAM,iBAAiB,CAAC;AACzB,OAAO,EAAE,iBAAiB,EAAE,MAAM,iBAAiB,CAAC;AAOpD,OAAO,EAAE,gBAAgB,EAAE,MAAM,oBAAoB,CAAC;AACtD,YAAY,EACV,mBAAmB,EACnB,oBAAoB,EACpB,mBAAmB,EACnB,kBAAkB,EAClB,UAAU,GACX,MAAM,oBAAoB,CAAC;AAG5B,OAAO,EACL,4BAA4B,EAC5B,KAAK,0BAA0B,EAC/B,KAAK,0BAA0B,EAC/B,4BAA4B,EAC5B,KAAK,0BAA0B,EAC/B,KAAK,0BAA0B,GAChC,MAAM,oBAAoB,CAAC;AAM5B,OAAO,EAAE,WAAW,EAAE,MAAM,6BAA6B,CAAC;AAG1D,YAAY,EAAE,uBAAuB,EAAE,MAAM,4BAA4B,CAAC;AAG1E,OAAO,EACL,iBAAiB,EACjB,KAAK,wBAAwB,EAC7B,KAAK,wBAAwB,GAC9B,MAAM,oCAAoC,CAAC;AAM5C,YAAY,EAAE,SAAS,EAAE,MAAM,wBAAwB,CAAC;AACxD,OAAO,EACL,aAAa,EACb,MAAM,EACN,IAAI,EACJ,IAAI,EACJ,SAAS,EACT,IAAI,EACJ,QAAQ,EACR,YAAY,EACZ,KAAK,EACL,MAAM,EACN,KAAK,EACL,MAAM,EACN,WAAW,EACX,OAAO,EACP,MAAM,GACP,MAAM,wBAAwB,CAAC;AAChC,OAAO,EACL,QAAQ,EACR,MAAM,EACN,KAAK,EACL,GAAG,EACH,QAAQ,EACR,MAAM,EACN,MAAM,EACN,sBAAsB,GACvB,MAAM,0BAA0B,CAAC;AAClC,YAAY,EAAE,QAAQ,EAAE,gBAAgB,EAAE,SAAS,EAAE,MAAM,4BAA4B,CAAC;AACxF,OAAO,EAAE,QAAQ,EAAE,gBAAgB,EAAE,yBAAyB,EAAE,MAAM,4BAA4B,CAAC;AACnG,OAAO,EAAE,uBAAuB,EAAE,MAAM,+BAA+B,CAAC;AAOxE,OAAO,EACL,iBAAiB,EACjB,oBAAoB,EACpB,qBAAqB,EACrB,sBAAsB,EACtB,eAAe,EACf,eAAe,EACf,qBAAqB,EAErB,0BAA0B,EAC1B,qBAAqB,EACrB,kBAAkB,EAClB,qCAAqC,GACtC,MAAM,sBAAsB,CAAC;AAC9B,YAAY,EACV,eAAe,EACf,kBAAkB,EAClB,eAAe,EACf,eAAe,EACf,YAAY,EACZ,YAAY,EACZ,UAAU,EACV,WAAW,EACX,SAAS,EAET,iBAAiB,EACjB,iBAAiB,GAClB,MAAM,sBAAsB,CAAC"}
+{"version":3,"file":"index.d.ts","sourceRoot":"","sources":["../src/index.ts"],"names":[],"mappings":"AAAA;;;;;;GAMG;AAEH,OAAO,EAAE,eAAe,EAAE,MAAM,kBAAkB,CAAC;AACnD,YAAY,EACV,YAAY,EACZ,aAAa,EACb,WAAW,EACX,cAAc,EACd,kBAAkB,EAClB,aAAa,EACb,KAAK,EACL,QAAQ,EACR,kBAAkB,GACnB,MAAM,iBAAiB,CAAC;AACzB,OAAO,EAAE,iBAAiB,EAAE,MAAM,iBAAiB,CAAC;AAOpD,OAAO,EAAE,gBAAgB,EAAE,MAAM,oBAAoB,CAAC;AACtD,YAAY,EACV,mBAAmB,EACnB,oBAAoB,EACpB,mBAAmB,EACnB,kBAAkB,EAClB,UAAU,GACX,MAAM,oBAAoB,CAAC;AAG5B,OAAO,EACL,4BAA4B,EAC5B,KAAK,0BAA0B,EAC/B,KAAK,0BAA0B,EAC/B,4BAA4B,EAC5B,KAAK,0BAA0B,EAC/B,KAAK,0BAA0B,GAChC,MAAM,oBAAoB,CAAC;AAM5B,OAAO,EAAE,WAAW,EAAE,MAAM,6BAA6B,CAAC;AAG1D,YAAY,EAAE,uBAAuB,EAAE,MAAM,4BAA4B,CAAC;AAK1E,OAAO,EAAE,KAAK,EAAE,MAAM,4BAA4B,CAAC;AACnD,YAAY,EAAE,eAAe,EAAE,MAAM,4BAA4B,CAAC;AAOlE,OAAO,EAAE,QAAQ,EAAE,MAAM,4BAA4B,CAAC;AACtD,YAAY,EAAE,kBAAkB,EAAE,MAAM,4BAA4B,CAAC;AAKrE,OAAO,EAAE,eAAe,EAAE,MAAM,4BAA4B,CAAC;AAC7D,YAAY,EAAE,mBAAmB,EAAE,MAAM,4BAA4B,CAAC;AAGtE,OAAO,EACL,iBAAiB,EACjB,KAAK,wBAAwB,EAC7B,KAAK,wBAAwB,GAC9B,MAAM,oCAAoC,CAAC;AAM5C,YAAY,EAAE,SAAS,EAAE,MAAM,wBAAwB,CAAC;AACxD,OAAO,EACL,aAAa,EACb,MAAM,EACN,IAAI,EACJ,IAAI,EACJ,SAAS,EACT,IAAI,EACJ,QAAQ,EACR,YAAY,EACZ,KAAK,EACL,MAAM,EACN,KAAK,EACL,MAAM,EACN,WAAW,EACX,OAAO,EACP,MAAM,GACP,MAAM,wBAAwB,CAAC;AAChC,OAAO,EACL,QAAQ,EACR,MAAM,EACN,KAAK,EACL,GAAG,EACH,QAAQ,EACR,MAAM,EACN,MAAM,EACN,sBAAsB,GACvB,MAAM,0BAA0B,CAAC;AAClC,YAAY,EAAE,QAAQ,EAAE,gBAAgB,EAAE,SAAS,EAAE,MAAM,4BAA4B,CAAC;AACxF,OAAO,EAAE,QAAQ,EAAE,gBAAgB,EAAE,yBAAyB,EAAE,MAAM,4BAA4B,CAAC;AACnG,OAAO,EAAE,uBAAuB,EAAE,MAAM,+BAA+B,CAAC;AAQxE,OAAO,EACL,iBAAiB,EACjB,oBAAoB,EACpB,qBAAqB,EACrB,sBAAsB,EACtB,eAAe,EACf,eAAe,EACf,qBAAqB,EAErB,0BAA0B,EAC1B,qBAAqB,EACrB,kBAAkB,EAClB,qCAAqC,EAErC,oBAAoB,EAEpB,cAAc,GACf,MAAM,sBAAsB,CAAC;AAC9B,YAAY,EACV,eAAe,EACf,kBAAkB,EAClB,eAAe,EACf,eAAe,EACf,YAAY,EACZ,YAAY,EACZ,UAAU,EACV,WAAW,EACX,SAAS,EAET,iBAAiB,EACjB,iBAAiB,EAEjB,QAAQ,EACR,WAAW,EACX,UAAU,EAEV,gBAAgB,EAChB,qBAAqB,GACtB,MAAM,sBAAsB,CAAC"}

package/dist/index.js

@@ -20,18 +20,37 @@ export { evaluateGravitationalLensing, evaluatePerihelionPrecession, } from './b
 // callers compose Γ trees directly; `CovariantDerivativeNode` is the structural
 // type for the new ∇_μ AST kind.
 export { christoffel } from './dimensional/connection.js';
+// v0.5.0 curvature layer — Ricci tensor helper (Task 7). `ricci(R)` wraps a
+// RiemannTensorNode and produces the contracted R_μν = R^λ_{λμν} as a
+// composite ExprNode (own validator + lowering arms).
+export { ricci } from './dimensional/curvature.js';
+// v0.5.0 curvature layer — Einstein tensor helper (Task 8). `einstein(R, g,
+// gInverse)` wraps a RiemannTensorNode + metric pair and produces the
+// composite G_μν = R_μν − ½ R g_μν as an ExprNode (own validator +
+// lowering arms — no AST rewrite). Vacuum-Einstein scope; matter-coupled
+// `G_μν = κ T_μν` is deferred to v0.6.0+.
+export { einstein } from './dimensional/curvature.js';
+// v0.5.0 curvature layer — Bianchi residual helper (Task 9). `bianchiResidual(R)`
+// returns {residual, evaluate, evaluateMax} for the cyclic second-Bianchi-identity
+// check ∇_{[λ} R_{μν]ρσ} = 0. Closes Phase 1 — Foundations.
+export { bianchiResidual } from './dimensional/curvature.js';
 // v0.4.0 geodesic integrator (RK4 solver — headline feature of v0.4.0)
 export { integrateGeodesic, } from './numerical/geodesic-integrator.js';
 export { DIMENSIONLESS, LENGTH, AREA, TIME, FREQUENCY, MASS, VELOCITY, ACCELERATION, FORCE, ENERGY, POWER, ACTION, TEMPERATURE, ENTROPY, CHARGE, } from './dimensional/types.js';
 export { multiply, divide, power, add, subtract, equals, format, DimensionMismatchError, } from './dimensional/algebra.js';
 export { validate, validateEquation, validateInverseMetricPair } from './dimensional/validator.js';
 export { inferDimensionForBridge } from './dimensional/bridge-check.js';
-// v0.3.5 numerical-contraction backend. See docs/planning/v0.3.5-Design.md.
+// Numerical-contraction backend (v0.3.5+). See docs/planning/v0.3.5-Design.md
+// and docs/planning/v0.4.0-Implementation-Plan.md (v0.4.0 additions below).
 // Every symbol below is `@public` — the consumer-facing (TensorJS) surface.
 // `MathTSEngine` is intentionally NOT re-exported here: it lives behind the
 // `@danielsimonjr/mathts-tensor` optionalDependency and is reachable only via
 // the `universal-physics-tensor/numerical/mathts-engine` exports subpath.
 export { evaluateNumerical, evaluateNumericalRaw, evaluateMetricInverse, Float64ReferenceEngine, getActiveEngine, setActiveEngine, NumericalBackendError, 
 // v0.4.0 additions to the numerical surface
-DuplicateCoordinateWarning, EngineCapabilityError, hasAutogradSupport, evaluateBE37CovariantEikonalNumerical, } from './numerical/index.js';
+DuplicateCoordinateWarning, EngineCapabilityError, hasAutogradSupport, evaluateBE37CovariantEikonalNumerical, 
+// v0.5.0 GL4 symplectic integrator
+integrateGeodesicGL4, 
+// v0.5.0 perihelion finder (Task 4)
+findPerihelion, } from './numerical/index.js';
 //# sourceMappingURL=index.js.map

package/dist/index.js.map

@@ -1 +1 @@
-{"version":3,"file":"index.js","sourceRoot":"","sources":["../src/index.ts"],"names":[],"mappings":"AAAA;;;;;;GAMG;AAEH,OAAO,EAAE,eAAe,EAAE,MAAM,kBAAkB,CAAC;AAYnD,OAAO,EAAE,iBAAiB,EAAE,MAAM,iBAAiB,CAAC;AAEpD,yEAAyE;AACzE,4EAA4E;AAC5E,2EAA2E;AAC3E,2EAA2E;AAC3E,yDAAyD;AACzD,OAAO,EAAE,gBAAgB,EAAE,MAAM,oBAAoB,CAAC;AAStD,8EAA8E;AAC9E,OAAO,EACL,4BAA4B,EAG5B,4BAA4B,GAG7B,MAAM,oBAAoB,CAAC;AAE5B,iFAAiF;AACjF,+EAA+E;AAC/E,gFAAgF;AAChF,iCAAiC;AACjC,OAAO,EAAE,WAAW,EAAE,MAAM,6BAA6B,CAAC;AAK1D,uEAAuE;AACvE,OAAO,EACL,iBAAiB,GAGlB,MAAM,oCAAoC,CAAC;AAO5C,OAAO,EACL,aAAa,EACb,MAAM,EACN,IAAI,EACJ,IAAI,EACJ,SAAS,EACT,IAAI,EACJ,QAAQ,EACR,YAAY,EACZ,KAAK,EACL,MAAM,EACN,KAAK,EACL,MAAM,EACN,WAAW,EACX,OAAO,EACP,MAAM,GACP,MAAM,wBAAwB,CAAC;AAChC,OAAO,EACL,QAAQ,EACR,MAAM,EACN,KAAK,EACL,GAAG,EACH,QAAQ,EACR,MAAM,EACN,MAAM,EACN,sBAAsB,GACvB,MAAM,0BAA0B,CAAC;AAElC,OAAO,EAAE,QAAQ,EAAE,gBAAgB,EAAE,yBAAyB,EAAE,MAAM,4BAA4B,CAAC;AACnG,OAAO,EAAE,uBAAuB,EAAE,MAAM,+BAA+B,CAAC;AAExE,4EAA4E;AAC5E,4EAA4E;AAC5E,4EAA4E;AAC5E,8EAA8E;AAC9E,0EAA0E;AAC1E,OAAO,EACL,iBAAiB,EACjB,oBAAoB,EACpB,qBAAqB,EACrB,sBAAsB,EACtB,eAAe,EACf,eAAe,EACf,qBAAqB;AACrB,4CAA4C;AAC5C,0BAA0B,EAC1B,qBAAqB,EACrB,kBAAkB,EAClB,qCAAqC,GACtC,MAAM,sBAAsB,CAAC"}
+{"version":3,"file":"index.js","sourceRoot":"","sources":["../src/index.ts"],"names":[],"mappings":"AAAA;;;;;;GAMG;AAEH,OAAO,EAAE,eAAe,EAAE,MAAM,kBAAkB,CAAC;AAYnD,OAAO,EAAE,iBAAiB,EAAE,MAAM,iBAAiB,CAAC;AAEpD,yEAAyE;AACzE,4EAA4E;AAC5E,2EAA2E;AAC3E,2EAA2E;AAC3E,yDAAyD;AACzD,OAAO,EAAE,gBAAgB,EAAE,MAAM,oBAAoB,CAAC;AAStD,8EAA8E;AAC9E,OAAO,EACL,4BAA4B,EAG5B,4BAA4B,GAG7B,MAAM,oBAAoB,CAAC;AAE5B,iFAAiF;AACjF,+EAA+E;AAC/E,gFAAgF;AAChF,iCAAiC;AACjC,OAAO,EAAE,WAAW,EAAE,MAAM,6BAA6B,CAAC;AAK1D,4EAA4E;AAC5E,sEAAsE;AACtE,sDAAsD;AACtD,OAAO,EAAE,KAAK,EAAE,MAAM,4BAA4B,CAAC;AAGnD,4EAA4E;AAC5E,sEAAsE;AACtE,mEAAmE;AACnE,yEAAyE;AACzE,0CAA0C;AAC1C,OAAO,EAAE,QAAQ,EAAE,MAAM,4BAA4B,CAAC;AAGtD,kFAAkF;AAClF,mFAAmF;AACnF,4DAA4D;AAC5D,OAAO,EAAE,eAAe,EAAE,MAAM,4BAA4B,CAAC;AAG7D,uEAAuE;AACvE,OAAO,EACL,iBAAiB,GAGlB,MAAM,oCAAoC,CAAC;AAO5C,OAAO,EACL,aAAa,EACb,MAAM,EACN,IAAI,EACJ,IAAI,EACJ,SAAS,EACT,IAAI,EACJ,QAAQ,EACR,YAAY,EACZ,KAAK,EACL,MAAM,EACN,KAAK,EACL,MAAM,EACN,WAAW,EACX,OAAO,EACP,MAAM,GACP,MAAM,wBAAwB,CAAC;AAChC,OAAO,EACL,QAAQ,EACR,MAAM,EACN,KAAK,EACL,GAAG,EACH,QAAQ,EACR,MAAM,EACN,MAAM,EACN,sBAAsB,GACvB,MAAM,0BAA0B,CAAC;AAElC,OAAO,EAAE,QAAQ,EAAE,gBAAgB,EAAE,yBAAyB,EAAE,MAAM,4BAA4B,CAAC;AACnG,OAAO,EAAE,uBAAuB,EAAE,MAAM,+BAA+B,CAAC;AAExE,8EAA8E;AAC9E,4EAA4E;AAC5E,4EAA4E;AAC5E,4EAA4E;AAC5E,8EAA8E;AAC9E,0EAA0E;AAC1E,OAAO,EACL,iBAAiB,EACjB,oBAAoB,EACpB,qBAAqB,EACrB,sBAAsB,EACtB,eAAe,EACf,eAAe,EACf,qBAAqB;AACrB,4CAA4C;AAC5C,0BAA0B,EAC1B,qBAAqB,EACrB,kBAAkB,EAClB,qCAAqC;AACrC,mCAAmC;AACnC,oBAAoB;AACpB,oCAAoC;AACpC,cAAc,GACf,MAAM,sBAAsB,CAAC"}

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

@@ -1,41 +1,69 @@
 /**
- * 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
  */
 /**
- * Input parameters for the BE-37 covariant-eikonal preview evaluator.
+ * Input parameters for the BE-37 covariant-eikonal 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). */
+    /** Far-point radius in metres, i.e. the signal origin (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). */
+    /** Near-point radius in metres, i.e. the signal destination
+     *  (must be > 0, ≤ R_far_m, and outside the Schwarzschild horizon). */
     readonly R_near_m: number;
+    /**
+     * Optional impact parameter `b` in metres for the null geodesic. The
+     * photon's angular momentum is set such that the minimum r (closest
+     * approach) is `b` in the flat-space limit; `p_φ = b · c` in the
+     * affine-parameter normalization used here (`p_t = −c²`).
+     *
+     * Default `0` — purely radial null geodesic, which reproduces the
+     * closed-form `(2GM/c³) · ln(R_far/R_near)` Shapiro delay. Non-zero
+     * `b` modifies the geometry and the closed-form comparator changes
+     * to `(2GM/c³) · ln((R_far + √(R_far²−b²)) / (R_near + √(R_near²−b²)))`
+     * for the single-leg path (closest-approach not crossed).
+     *
+     * Must satisfy `b ≤ R_near` if non-zero (otherwise the photon cannot
+     * reach `R_near` without crossing closest approach).
+     */
+    readonly b_m?: number;
+    /**
+     * Optional GL4 step count override. Default `2048` — sufficient for
+     * the default radial geometry to reach ≤1e-6 relative agreement with
+     * `evaluateShapiroDelay` on Earth/Mars/Sun scales (Task 11).
+     */
+    readonly steps?: number;
 }
 /**
  * Result of `evaluateBE37CovariantEikonalNumerical`.
@@ -46,44 +74,74 @@ 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.
+     * 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.
+     * Shapiro coordinate-time delay in seconds, computed via GL4
+     * null-geodesic integration on the canonical (x, p) state.
      *
-     * **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`.
+     * **v0.5.0 Task 12:** real value (was a documented stub returning 0 in
+     * v0.4.0). For the default radial geometry (`b_m = 0`) this agrees with
+     * the closed-form `(2GM/c³) · ln(R_far/R_near)` to ≤1e-6 relative at
+     * the default step count (Task 11 [U]).
      */
     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.
+     * Optional closed-form cross-check value (seconds). Not populated by
+     * this evaluator — callers cross-check externally against
+     * `evaluateShapiroDelay`.
      */
     readonly closedFormDelaySec?: number;
 }
 /**
- * v0.4.0 BE-37 covariant-eikonal STRUCTURAL PREVIEW.
+ * 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`.
  *
- * 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).
+ * **Sign convention.** Shapiro delay is positive (light slowed by the
+ * gravitational potential). A negative result indicates an integration
+ * or reference-frame bug.
  *
- * **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.
+ * **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.
- * @returns Structural eikonal residual (0 by construction) + stub Shapiro delay.
+ * @param inputs - Source mass + far/near radii (+ optional `b_m`, `steps`).
+ * @returns `{eikonalResidual: 0, shapiroDelaySec: <real, positive>}`.
+ * @throws RangeError on out-of-domain inputs.
  * @public
  */
 export declare function evaluateBE37CovariantEikonalNumerical(inputs: BE37CovariantEikonalInputs): Promise<BE37CovariantEikonalResult>;

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

@@ -1 +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"}
+{"version":3,"file":"be37-covariant-eikonal.d.ts","sourceRoot":"","sources":["../../src/numerical/be37-covariant-eikonal.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GA+BG;AASH;;;GAGG;AACH,MAAM,WAAW,0BAA0B;IACzC,4DAA4D;IAC5D,QAAQ,CAAC,IAAI,EAAE,MAAM,CAAC;IACtB,wEAAwE;IACxE,QAAQ,CAAC,OAAO,EAAE,MAAM,CAAC;IACzB;2EACuE;IACvE,QAAQ,CAAC,QAAQ,EAAE,MAAM,CAAC;IAC1B;;;;;;;;;;;;;;OAcG;IACH,QAAQ,CAAC,GAAG,CAAC,EAAE,MAAM,CAAC;IACtB;;;;OAIG;IACH,QAAQ,CAAC,KAAK,CAAC,EAAE,MAAM,CAAC;CACzB;AAED;;;GAGG;AACH,MAAM,WAAW,0BAA0B;IACzC;;;;;;;;;OASG;IACH,QAAQ,CAAC,eAAe,EAAE,MAAM,CAAC;IAEjC;;;;;;;;OAQG;IACH,QAAQ,CAAC,eAAe,EAAE,MAAM,CAAC;IAEjC;;;;OAIG;IACH,QAAQ,CAAC,kBAAkB,CAAC,EAAE,MAAM,CAAC;CACtC;AAgFD;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GA6CG;AACH,wBAAsB,qCAAqC,CACzD,MAAM,EAAE,0BAA0B,GACjC,OAAO,CAAC,0BAA0B,CAAC,CAqHrC"}

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

@@ -1,50 +1,158 @@
 /**
- * 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';
+// ─── Physical constants (CODATA 2018 + SI definitions) ────────────────────
+const G_SI = 6.6743e-11; // m³ kg⁻¹ s⁻²  (matches tests/fixtures/schwarzschild.ts)
+const c_SI = 2.998e8; // m s⁻¹          (matches tests/fixtures/schwarzschild.ts)
+const c2_SI = c_SI * c_SI;
 /**
- * v0.4.0 BE-37 covariant-eikonal STRUCTURAL PREVIEW.
+ * Schwarzschild radius `r_s = 2GM/c²`.
  *
- * 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).
+ * Mirrors `tests/fixtures/schwarzschild.ts` (kept private here to avoid
+ * a src→tests dependency).
  *
- * **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.
+ * @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.
- * @returns Structural eikonal residual (0 by construction) + stub Shapiro delay.
+ * 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`).
+ * @returns `{eikonalResidual: 0, shapiroDelaySec: <real, positive>}`.
+ * @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}`);
     }
@@ -57,23 +165,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;;;;;;;;;;;;;;;;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"}
+{"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;AAE3D,6EAA6E;AAC7E,MAAM,IAAI,GAAG,UAAU,CAAC,CAAC,yDAAyD;AAClF,MAAM,IAAI,GAAG,OAAO,CAAC,CAAC,2DAA2D;AACjF,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;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GA6CG;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/connection-lowering-helpers.d.ts

@@ -15,11 +15,10 @@ 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.
+ * Canonical implementation for this module. `flattenNestedArray` in
+ * lowering.ts wraps this and adds an expectedSize invariant check.
+ * `pderiv.ts` previously had its own identical `flattenToNumbers` —
+ * consolidated to this function in v0.4.6 (see Task 9 of the v0.4.6 plan).
  */
 export declare function flattenNA(data: NestedArray): number[];
 /**
@@ -97,7 +96,10 @@ export declare function computeChristoffelTensor(gInverseFlat: number[], getMetr
 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).
+ * Key format: `${metricName}/${coordLabel}_${mu}` (e.g. 'g/x_0' where
+ * 'x' is the coordLabel string and '0' is the mu index). The coordLabel
+ * is whatever string the caller passes — typically a coordinate label
+ * like 'x', 't', 'r'; the literal Greek letter μ is NOT used as a key.
  * Returns a flat [N*N] array.
  *
  * Strategy 'zero': returns [N*N] zeros.

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

@@ -1 +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"}
+{"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;AAQ9C;;;;;;GAMG;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;AAMD;;;;;;;;;;;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,CAqCd;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;AA8BD;;;;;;;;;;GAUG;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"}