universal-physics-tensor 10.0.0 → 12.0.0

package/dist/numerical/finite-difference.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"finite-difference.d.ts","sourceRoot":"","sources":["../../src/numerical/finite-difference.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;GAYG;AAEH;;;;;;GAMG;AACH,wBAAgB,cAAc,CAC5B,CAAC,EAAE,CAAC,KAAK,EAAE,MAAM,CAAC,MAAM,EAAE,MAAM,CAAC,KAAK,MAAM,EAC5C,KAAK,EAAE,QAAQ,CAAC,MAAM,CAAC,MAAM,EAAE,MAAM,CAAC,CAAC,EACvC,GAAG,EAAE,MAAM,EACX,IAAI,EAAE,MAAM,GACX,MAAM,CAKR"}

package/dist/numerical/finite-difference.js

@@ -0,0 +1,27 @@
+/**
+ * One central-difference rule.
+ *
+ * Three readers differentiated a closed form by central differences, each
+ * with its own step rule and its own denominator (`2h` at one, the
+ * representable step at another). The derivative is here, once: the caller
+ * passes the step its rule chose, and the denominator is the step the
+ * floating-point point actually took, `dx = (x + h) − (x − h)`, so a
+ * perturbation that rounds at an astrophysical `x` is divided by what was
+ * applied, not by what was asked.
+ *
+ * @module numerical/finite-difference
+ */
+/**
+ * `∂f/∂x_key` at `point` by the central difference with half-step `step`:
+ * `(f(x + h) − f(x − h)) / ((x + h) − (x − h))`. The other keys of `point` are
+ * held. The step rule is the caller's; a zero or non-finite `step` makes the
+ * denominator 0 or NaN, which the caller refuses before calling.
+ * @internal
+ */
+export function centralPartial(f, point, key, step) {
+    const x = point[key];
+    const up = x + step;
+    const down = x - step;
+    return (f({ ...point, [key]: up }) - f({ ...point, [key]: down })) / (up - down);
+}
+//# sourceMappingURL=finite-difference.js.map

package/dist/numerical/finite-difference.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"finite-difference.js","sourceRoot":"","sources":["../../src/numerical/finite-difference.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;GAYG;AAEH;;;;;;GAMG;AACH,MAAM,UAAU,cAAc,CAC5B,CAA4C,EAC5C,KAAuC,EACvC,GAAW,EACX,IAAY;IAEZ,MAAM,CAAC,GAAG,KAAK,CAAC,GAAG,CAAE,CAAC;IACtB,MAAM,EAAE,GAAG,CAAC,GAAG,IAAI,CAAC;IACpB,MAAM,IAAI,GAAG,CAAC,GAAG,IAAI,CAAC;IACtB,OAAO,CAAC,CAAC,CAAC,EAAE,GAAG,KAAK,EAAE,CAAC,GAAG,CAAC,EAAE,EAAE,EAAE,CAAC,GAAG,CAAC,CAAC,EAAE,GAAG,KAAK,EAAE,CAAC,GAAG,CAAC,EAAE,IAAI,EAAE,CAAC,CAAC,GAAG,CAAC,EAAE,GAAG,IAAI,CAAC,CAAC;AACnF,CAAC","sourcesContent":["/**\n * One central-difference rule.\n *\n * Three readers differentiated a closed form by central differences, each\n * with its own step rule and its own denominator (`2h` at one, the\n * representable step at another). The derivative is here, once: the caller\n * passes the step its rule chose, and the denominator is the step the\n * floating-point point actually took, `dx = (x + h) − (x − h)`, so a\n * perturbation that rounds at an astrophysical `x` is divided by what was\n * applied, not by what was asked.\n *\n * @module numerical/finite-difference\n */\n\n/**\n * `∂f/∂x_key` at `point` by the central difference with half-step `step`:\n * `(f(x + h) − f(x − h)) / ((x + h) − (x − h))`. The other keys of `point` are\n * held. The step rule is the caller's; a zero or non-finite `step` makes the\n * denominator 0 or NaN, which the caller refuses before calling.\n * @internal\n */\nexport function centralPartial(\n  f: (point: Record<string, number>) => number,\n  point: Readonly<Record<string, number>>,\n  key: string,\n  step: number,\n): number {\n  const x = point[key]!;\n  const up = x + step;\n  const down = x - step;\n  return (f({ ...point, [key]: up }) - f({ ...point, [key]: down })) / (up - down);\n}\n"]}

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

@@ -1 +1 @@
-{"version":3,"file":"geodesic-integrator.d.ts","sourceRoot":"","sources":["../../src/numerical/geodesic-integrator.ts"],"names":[],"mappings":"AA2BA;;;GAGG;AACH,MAAM,WAAW,wBAAwB;IACvC;;;;;;;;;;;;OAYG;IACH,aAAa,EAAE,CAAC,CAAC,EAAE,aAAa,CAAC,MAAM,CAAC,EAAE,GAAG,CAAC,EAAE,YAAY,KAAK,YAAY,CAAC;IAE9E,iDAAiD;IACjD,EAAE,EAAE,SAAS,CAAC,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,CAAC,CAAC;IAE9C;;;;;;;;OAQG;IACH,EAAE,EAAE,SAAS,CAAC,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,CAAC,CAAC;IAE9C,yDAAyD;IACzD,QAAQ,EAAE,MAAM,CAAC;IAEjB,uDAAuD;IACvD,MAAM,EAAE,MAAM,CAAC;IAEf,8DAA8D;IAC9D,KAAK,EAAE,MAAM,CAAC;IAEd;;;;;;OAMG;IACH,eAAe,CAAC,EAAE,MAAM,CAAC;CAC1B;AAED;;;GAGG;AACH,MAAM,WAAW,wBAAwB;IACvC,kCAAkC;IAClC,QAAQ,CAAC,MAAM,EAAE,SAAS,CAAC,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,CAAC,CAAC;IAE3D,kCAAkC;IAClC,QAAQ,CAAC,MAAM,EAAE,SAAS,CAAC,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,CAAC,CAAC;IAE3D;;;OAGG;IACH,QAAQ,CAAC,UAAU,EAAE,aAAa,CAAC,SAAS,CAAC,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,CAAC,CAAC,CAAC;CAC/E;AAuFD;;;;;;GAMG;AACH,wBAAgB,iBAAiB,CAC/B,MAAM,EAAE,wBAAwB,GAC/B,wBAAwB,CAuD1B"}
+{"version":3,"file":"geodesic-integrator.d.ts","sourceRoot":"","sources":["../../src/numerical/geodesic-integrator.ts"],"names":[],"mappings":"AA4BA;;;GAGG;AACH,MAAM,WAAW,wBAAwB;IACvC;;;;;;;;;;;;OAYG;IACH,aAAa,EAAE,CAAC,CAAC,EAAE,aAAa,CAAC,MAAM,CAAC,EAAE,GAAG,CAAC,EAAE,YAAY,KAAK,YAAY,CAAC;IAE9E,iDAAiD;IACjD,EAAE,EAAE,SAAS,CAAC,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,CAAC,CAAC;IAE9C;;;;;;;;OAQG;IACH,EAAE,EAAE,SAAS,CAAC,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,CAAC,CAAC;IAE9C,yDAAyD;IACzD,QAAQ,EAAE,MAAM,CAAC;IAEjB,uDAAuD;IACvD,MAAM,EAAE,MAAM,CAAC;IAEf,8DAA8D;IAC9D,KAAK,EAAE,MAAM,CAAC;IAEd;;;;;;OAMG;IACH,eAAe,CAAC,EAAE,MAAM,CAAC;CAC1B;AAED;;;GAGG;AACH,MAAM,WAAW,wBAAwB;IACvC,kCAAkC;IAClC,QAAQ,CAAC,MAAM,EAAE,SAAS,CAAC,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,CAAC,CAAC;IAE3D,kCAAkC;IAClC,QAAQ,CAAC,MAAM,EAAE,SAAS,CAAC,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,CAAC,CAAC;IAE3D;;;OAGG;IACH,QAAQ,CAAC,UAAU,EAAE,aAAa,CAAC,SAAS,CAAC,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,CAAC,CAAC,CAAC;CAC/E;AAuFD;;;;;;GAMG;AACH,wBAAgB,iBAAiB,CAC/B,MAAM,EAAE,wBAAwB,GAC/B,wBAAwB,CA8D1B"}

package/dist/numerical/geodesic-integrator.js

@@ -14,16 +14,17 @@
  * The Christoffel symbol Γ^μ_{νρ} is supplied as a plain-JS closure that
  * maps a 4-coordinate array to a [4][4][4] tensor.  No TensorEngine
  * dependency — the integrator is self-contained. The fixed steps are
- * MathTS `solveODESystem` with `dt` set, which is that library's classical
- * RK4. Adaptive `solveODE` is not this step.
+ * `rk4Trajectory` (MathTS `solveODESystem` with `dt` set, that library's
+ * classical RK4), which takes exactly `steps` steps. Adaptive `solveODE` is
+ * not this step.
  *
  * Validated against the cycloid parametric form for radial Schwarzschild
  * infall to ±1e-6 (Task 14 test, η = 0.5).
  *
  * @module numerical/geodesic-integrator
  */
-import { solveODESystem } from '@danielsimonjr/mathts-functions';
 import { NumericalBackendError } from './errors.js';
+import { rk4Trajectory } from './rk4-trajectory.js';
 /**
  * Evaluate the Christoffel acceleration −Γ^μ_{νρ} v^ν v^ρ at position x
  * with 4-velocity v.  Returns (dx/dτ, dv/dτ) = (v, accel).
@@ -149,28 +150,35 @@ export function integrateGeodesic(inputs) {
     if (!Number.isInteger(steps) || steps <= 0) {
         throw new NumericalBackendError(`integrateGeodesic: steps must be a positive integer, got ${steps}`);
     }
-    const h = (tauEnd - tauStart) / steps;
+    // A NaN or infinite end made every step NaN and the result the start state; an
+    // end equal to the start is a zero step, which is not an integration either.
+    if (!Number.isFinite(tauStart)) {
+        throw new NumericalBackendError(`integrateGeodesic: tauStart must be finite, got ${tauStart}`);
+    }
+    if (!Number.isFinite(tauEnd)) {
+        throw new NumericalBackendError(`integrateGeodesic: tauEnd must be finite, got ${tauEnd}`);
+    }
+    if (tauEnd === tauStart) {
+        throw new NumericalBackendError(`integrateGeodesic: tauEnd must differ from tauStart, both are ${tauStart}`);
+    }
     const y0 = [x0[0], x0[1], x0[2], x0[3], v0[0], v0[1], v0[2], v0[3]];
     // One scratch buffer reused across all Christoffel evaluations (4 per step).
     // Each geodesicRHS consumes its Γ before the next call, so reuse is safe and
     // avoids ~4·steps Float64Array(64) allocations per integration.
     const scratch = new Float64Array(64);
-    const sol = solveODESystem((_tau, y) => {
+    const states = rk4Trajectory((_tau, y) => {
         const x = [y[0], y[1], y[2], y[3]];
         const v = [y[4], y[5], y[6], y[7]];
         const s = geodesicRHS(christoffelFn, x, v, scratch);
         return [s.dx[0], s.dx[1], s.dx[2], s.dx[3], s.dv[0], s.dv[1], s.dv[2], s.dv[3]];
-    }, y0, [tauStart, tauEnd], { dt: h });
+    }, y0, tauStart, tauEnd, steps);
     const asVec4 = (row) => [row[0], row[1], row[2], row[3]];
     const sampleEvery = Math.max(1, Math.floor(steps / 100));
-    const trajectory = [asVec4(sol.y[0] ?? y0)];
-    for (let i = 0; i < steps; i++) {
-        if ((i + 1) % sampleEvery === 0) {
-            const row = sol.y[i + 1] ?? sol.y[sol.y.length - 1] ?? y0;
-            trajectory.push(asVec4(row));
-        }
-    }
-    const last = sol.y[sol.y.length - 1] ?? y0;
+    // `rk4Trajectory` returns exactly steps + 1 states, so state i is the one after i steps.
+    const trajectory = [asVec4(states[0] ?? y0)];
+    for (let i = sampleEvery; i <= steps; i += sampleEvery)
+        trajectory.push(asVec4(states[i] ?? y0));
+    const last = states[steps] ?? y0;
     return {
         xFinal: asVec4(last),
         vFinal: [last[4], last[5], last[6], last[7]],

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

@@ -1 +1 @@
-{"version":3,"file":"geodesic-integrator.js","sourceRoot":"","sources":["../../src/numerical/geodesic-integrator.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;;;;;GAuBG;AACH,OAAO,EAAE,cAAc,EAAE,MAAM,iCAAiC,CAAC;AACjE,OAAO,EAAE,qBAAqB,EAAE,MAAM,aAAa,CAAC;AA+EpD;;;;;;GAMG;AACH,SAAS,WAAW,CAClB,aAA6E,EAC7E,CAAO,EACP,CAAO;AACP,yEAAyE;AACzE,4EAA4E;AAC5E,8DAA8D;AAC9D,OAAsB;IAEtB,MAAM,CAAC,GAAG,aAAa,CAAC,CAAC,EAAE,OAAO,CAAC,CAAC;IACpC,MAAM,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,EAAE,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,EAAE,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,EAAE,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;IACjD,IAAI,GAAG,GAAG,CAAC,EAAE,GAAG,GAAG,CAAC,EAAE,GAAG,GAAG,CAAC,EAAE,GAAG,GAAG,CAAC,CAAC;IAEvC,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,CAAC,EAAE,EAAE,EAAE,EAAE,CAAC;QAC9B,6EAA6E;QAC7E,MAAM,IAAI,GAAG,EAAE,GAAG,EAAE,CAAC;QACrB,IAAI,GAAG,GAAG,CAAC,CAAC;QACZ,IAAI,OAAO,GAAG,CAAC,CAAC;QAChB,IAAI,CAAC,GAAG,CAAC,CAAC;QAEV,6EAA6E;QAC7E,sFAAsF;QAEtF,SAAS;QACT,IAAI,MAAM,GAAG,IAAI,CAAC;QAClB,OAAO,GAAG,CAAC,CAAC;QACZ,CAAC,GAAG,CAAC,CAAC,MAAM,CAAC,CAAC;QAAC,IAAI,CAAC,KAAK,CAAC;YAAE,OAAO,IAAI,CAAC,GAAG,EAAE,CAAC;QAC9C,CAAC,GAAG,CAAC,CAAC,MAAM,GAAC,CAAC,CAAC,CAAC;QAAC,IAAI,CAAC,KAAK,CAAC;YAAE,OAAO,IAAI,CAAC,GAAG,EAAE,CAAC;QAChD,CAAC,GAAG,CAAC,CAAC,MAAM,GAAC,CAAC,CAAC,CAAC;QAAC,IAAI,CAAC,KAAK,CAAC;YAAE,OAAO,IAAI,CAAC,GAAG,EAAE,CAAC;QAChD,CAAC,GAAG,CAAC,CAAC,MAAM,GAAC,CAAC,CAAC,CAAC;QAAC,IAAI,CAAC,KAAK,CAAC;YAAE,OAAO,IAAI,CAAC,GAAG,EAAE,CAAC;QAChD,IAAI,OAAO,KAAK,CAAC;YAAE,GAAG,IAAI,OAAO,GAAG,EAAE,CAAC;QAEvC,SAAS;QACT,MAAM,IAAI,CAAC,CAAC;QACZ,OAAO,GAAG,CAAC,CAAC;QACZ,CAAC,GAAG,CAAC,CAAC,MAAM,CAAC,CAAC;QAAC,IAAI,CAAC,KAAK,CAAC;YAAE,OAAO,IAAI,CAAC,GAAG,EAAE,CAAC;QAC9C,CAAC,GAAG,CAAC,CAAC,MAAM,GAAC,CAAC,CAAC,CAAC;QAAC,IAAI,CAAC,KAAK,CAAC;YAAE,OAAO,IAAI,CAAC,GAAG,EAAE,CAAC;QAChD,CAAC,GAAG,CAAC,CAAC,MAAM,GAAC,CAAC,CAAC,CAAC;QAAC,IAAI,CAAC,KAAK,CAAC;YAAE,OAAO,IAAI,CAAC,GAAG,EAAE,CAAC;QAChD,CAAC,GAAG,CAAC,CAAC,MAAM,GAAC,CAAC,CAAC,CAAC;QAAC,IAAI,CAAC,KAAK,CAAC;YAAE,OAAO,IAAI,CAAC,GAAG,EAAE,CAAC;QAChD,IAAI,OAAO,KAAK,CAAC;YAAE,GAAG,IAAI,OAAO,GAAG,EAAE,CAAC;QAEvC,SAAS;QACT,MAAM,IAAI,CAAC,CAAC;QACZ,OAAO,GAAG,CAAC,CAAC;QACZ,CAAC,GAAG,CAAC,CAAC,MAAM,CAAC,CAAC;QAAC,IAAI,CAAC,KAAK,CAAC;YAAE,OAAO,IAAI,CAAC,GAAG,EAAE,CAAC;QAC9C,CAAC,GAAG,CAAC,CAAC,MAAM,GAAC,CAAC,CAAC,CAAC;QAAC,IAAI,CAAC,KAAK,CAAC;YAAE,OAAO,IAAI,CAAC,GAAG,EAAE,CAAC;QAChD,CAAC,GAAG,CAAC,CAAC,MAAM,GAAC,CAAC,CAAC,CAAC;QAAC,IAAI,CAAC,KAAK,CAAC;YAAE,OAAO,IAAI,CAAC,GAAG,EAAE,CAAC;QAChD,CAAC,GAAG,CAAC,CAAC,MAAM,GAAC,CAAC,CAAC,CAAC;QAAC,IAAI,CAAC,KAAK,CAAC;YAAE,OAAO,IAAI,CAAC,GAAG,EAAE,CAAC;QAChD,IAAI,OAAO,KAAK,CAAC;YAAE,GAAG,IAAI,OAAO,GAAG,EAAE,CAAC;QAEvC,SAAS;QACT,MAAM,IAAI,CAAC,CAAC;QACZ,OAAO,GAAG,CAAC,CAAC;QACZ,CAAC,GAAG,CAAC,CAAC,MAAM,CAAC,CAAC;QAAC,IAAI,CAAC,KAAK,CAAC;YAAE,OAAO,IAAI,CAAC,GAAG,EAAE,CAAC;QAC9C,CAAC,GAAG,CAAC,CAAC,MAAM,GAAC,CAAC,CAAC,CAAC;QAAC,IAAI,CAAC,KAAK,CAAC;YAAE,OAAO,IAAI,CAAC,GAAG,EAAE,CAAC;QAChD,CAAC,GAAG,CAAC,CAAC,MAAM,GAAC,CAAC,CAAC,CAAC;QAAC,IAAI,CAAC,KAAK,CAAC;YAAE,OAAO,IAAI,CAAC,GAAG,EAAE,CAAC;QAChD,CAAC,GAAG,CAAC,CAAC,MAAM,GAAC,CAAC,CAAC,CAAC;QAAC,IAAI,CAAC,KAAK,CAAC;YAAE,OAAO,IAAI,CAAC,GAAG,EAAE,CAAC;QAChD,IAAI,OAAO,KAAK,CAAC;YAAE,GAAG,IAAI,OAAO,GAAG,EAAE,CAAC;QAEvC,IAAI,EAAE,KAAK,CAAC;YAAE,GAAG,GAAG,CAAC,GAAG,CAAC;aACpB,IAAI,EAAE,KAAK,CAAC;YAAE,GAAG,GAAG,CAAC,GAAG,CAAC;aACzB,IAAI,EAAE,KAAK,CAAC;YAAE,GAAG,GAAG,CAAC,GAAG,CAAC;;YACzB,GAAG,GAAG,CAAC,GAAG,CAAC;IAClB,CAAC;IAED,OAAO,EAAE,EAAE,EAAE,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,EAAE,EAAE,EAAE,CAAC,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,CAAC,EAAE,CAAC;AAC5D,CAAC;AAED,8EAA8E;AAC9E,aAAa;AACb,8EAA8E;AAE9E;;;;;;GAMG;AACH,MAAM,UAAU,iBAAiB,CAC/B,MAAgC;IAEhC,MAAM,EAAE,aAAa,EAAE,EAAE,EAAE,EAAE,EAAE,QAAQ,EAAE,MAAM,EAAE,KAAK,EAAE,eAAe,EAAE,GAAG,MAAM,CAAC;IACnF,mFAAmF;IAEnF,2EAA2E;IAC3E,yEAAyE;IACzE,IAAI,OAAO,eAAe,KAAK,QAAQ,IAAI,EAAE,CAAC,CAAC,CAAC,GAAG,eAAe,EAAE,CAAC;QACnE,MAAM,IAAI,qBAAqB,CAC7B,wCAAwC,EAAE,CAAC,CAAC,CAAC,sBAAsB,eAAe,IAAI;YACtF,oDAAoD;YACpD,2EAA2E,CAC5E,CAAC;IACJ,CAAC;IAED,IAAI,CAAC,MAAM,CAAC,SAAS,CAAC,KAAK,CAAC,IAAI,KAAK,IAAI,CAAC,EAAE,CAAC;QAC3C,MAAM,IAAI,qBAAqB,CAC7B,4DAA4D,KAAK,EAAE,CACpE,CAAC;IACJ,CAAC;IAED,MAAM,CAAC,GAAG,CAAC,MAAM,GAAG,QAAQ,CAAC,GAAG,KAAK,CAAC;IACtC,MAAM,EAAE,GAAG,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;IAEpE,6EAA6E;IAC7E,6EAA6E;IAC7E,gEAAgE;IAChE,MAAM,OAAO,GAAG,IAAI,YAAY,CAAC,EAAE,CAAC,CAAC;IACrC,MAAM,GAAG,GAAG,cAAc,CACxB,CAAC,IAAI,EAAE,CAAC,EAAE,EAAE;QACV,MAAM,CAAC,GAAS,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QACzC,MAAM,CAAC,GAAS,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QACzC,MAAM,CAAC,GAAG,WAAW,CAAC,aAAa,EAAE,CAAC,EAAE,CAAC,EAAE,OAAO,CAAC,CAAC;QACpD,OAAO,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;IAClF,CAAC,EACD,EAAE,EACF,CAAC,QAAQ,EAAE,MAAM,CAAC,EAClB,EAAE,EAAE,EAAE,CAAC,EAAE,CACV,CAAC;IAEF,MAAM,MAAM,GAAG,CAAC,GAAsB,EAAQ,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;IAClF,MAAM,WAAW,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,IAAI,CAAC,KAAK,CAAC,KAAK,GAAG,GAAG,CAAC,CAAC,CAAC;IACzD,MAAM,UAAU,GAAW,CAAC,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI,EAAE,CAAC,CAAC,CAAC;IACpD,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,KAAK,EAAE,CAAC,EAAE,EAAE,CAAC;QAC/B,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,WAAW,KAAK,CAAC,EAAE,CAAC;YAChC,MAAM,GAAG,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,IAAI,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,MAAM,GAAG,CAAC,CAAC,IAAI,EAAE,CAAC;YAC1D,UAAU,CAAC,IAAI,CAAC,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC;QAC/B,CAAC;IACH,CAAC;IACD,MAAM,IAAI,GAAG,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,MAAM,GAAG,CAAC,CAAC,IAAI,EAAE,CAAC;IAE3C,OAAO;QACL,MAAM,EAAE,MAAM,CAAC,IAAI,CAAC;QACpB,MAAM,EAAE,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,CAAC,CAAC;QAC5C,UAAU;KACX,CAAC;AACJ,CAAC","sourcesContent":["/**\n * Fixed-step RK4 integrator for the geodesic equation in an arbitrary\n * Lorentzian manifold (v0.4.0, Task 14 [U]).\n *\n * Solves the second-order ODE system\n *\n *   d²x^μ/dτ² + Γ^μ_{νρ}(x) (dx^ν/dτ)(dx^ρ/dτ) = 0\n *\n * by casting it as a first-order system on the (x, v) phase space:\n *\n *   dx^μ/dτ  = v^μ\n *   dv^μ/dτ  = −Γ^μ_{νρ}(x) v^ν v^ρ    (Einstein summation over ν, ρ)\n *\n * The Christoffel symbol Γ^μ_{νρ} is supplied as a plain-JS closure that\n * maps a 4-coordinate array to a [4][4][4] tensor.  No TensorEngine\n * dependency — the integrator is self-contained. The fixed steps are\n * MathTS `solveODESystem` with `dt` set, which is that library's classical\n * RK4. Adaptive `solveODE` is not this step.\n *\n * Validated against the cycloid parametric form for radial Schwarzschild\n * infall to ±1e-6 (Task 14 test, η = 0.5).\n *\n * @module numerical/geodesic-integrator\n */\nimport { solveODESystem } from '@danielsimonjr/mathts-functions';\nimport { NumericalBackendError } from './errors.js';\n\n/**\n * Input bundle for `integrateGeodesic`.\n * @public\n */\nexport interface GeodesicIntegratorInputs {\n  /**\n   * Christoffel-symbol closure.  Maps a 4-coordinate array x^μ to a\n   * Float64Array(64) containing Γ^μ_{νρ}(x) in λ-major layout:\n   * index (λ, μ, ν) → 16·λ + 4·μ + ν.\n   *\n   * BR-2 (v0.6.0 Task 2.9): migrated from nested number[][][] to flat\n   * Float64Array(64). Consumers call `arr[16*mu + 4*nu + rho]`.\n   *\n   * The optional `out` buffer lets the integrator pass a reused scratch\n   * Float64Array(64) so the hot loop avoids per-call allocation. A function\n   * that ignores `out` (e.g. allocates fresh each call) still satisfies this\n   * type — the optimization is opt-in.\n   */\n  christoffelFn: (x: ReadonlyArray<number>, out?: Float64Array) => Float64Array;\n\n  /** Initial 4-position x^μ(τ₀) = [t, r, θ, φ]. */\n  x0: readonly [number, number, number, number];\n\n  /**\n   * Initial 4-velocity dx^μ/dτ at τ₀.\n   *\n   * I2 fix: for a Schwarzschild radial timelike geodesic starting at\n   * coordinate radius r₀, use\n   *   v0 = [1/√(1 − r_s/r₀), 0, 0, 0]\n   * so that g_{μν} u^μ u^ν = −c² holds.  At r₀ = 100·r_s this gives\n   * dt/dτ ≈ 1.005, not 1.\n   */\n  v0: readonly [number, number, number, number];\n\n  /** Proper-time parameter at the start of integration. */\n  tauStart: number;\n\n  /** Proper-time parameter at the end of integration. */\n  tauEnd: number;\n\n  /** Number of fixed RK4 steps.  Must be a positive integer. */\n  steps: number;\n\n  /**\n   * E11 fix: explicit domain-validity option in place of the\n   * `(christoffelFn as any).r_s` monkey-patch.  If supplied, the\n   * integrator throws `NumericalBackendError` when `x0[1] < domainMinRadius`\n   * (e.g. for Schwarzschild: caller passes `3 * r_s` to enforce the\n   * v0.4.0 domain restriction r ≥ 3·r_Schwarz).\n   */\n  domainMinRadius?: number;\n}\n\n/**\n * Result of `integrateGeodesic`.\n * @public\n */\nexport interface GeodesicIntegratorResult {\n  /** Final 4-position x^μ(τEnd). */\n  readonly xFinal: readonly [number, number, number, number];\n\n  /** Final 4-velocity v^μ(τEnd). */\n  readonly vFinal: readonly [number, number, number, number];\n\n  /**\n   * Sampled trajectory positions (approximately 100 evenly-spaced snapshots\n   * plus the initial point).  Useful for visualisation and diagnostics.\n   */\n  readonly trajectory: ReadonlyArray<readonly [number, number, number, number]>;\n}\n\n// ---------------------------------------------------------------------------\n// Internal helpers\n// ---------------------------------------------------------------------------\n\ntype Vec4 = [number, number, number, number];\n\n/**\n * Evaluate the Christoffel acceleration −Γ^μ_{νρ} v^ν v^ρ at position x\n * with 4-velocity v.  Returns (dx/dτ, dv/dτ) = (v, accel).\n *\n * BR-2 (v0.6.0 Task 2.9): christoffelFn now returns Float64Array(64).\n * Index access: G[16*mu + 4*nu + rho] (λ-major layout).\n */\nfunction geodesicRHS(\n  christoffelFn: (x: ReadonlyArray<number>, out?: Float64Array) => Float64Array,\n  x: Vec4,\n  v: Vec4,\n  // Reused scratch buffer for the Christoffel evaluation. Safe: G is fully\n  // consumed into `dv` below before geodesicRHS returns, so the next call may\n  // overwrite it. `dx`/`dv` are fresh arrays and never alias G.\n  scratch?: Float64Array,\n): { dx: Vec4; dv: Vec4 } {\n  const G = christoffelFn(x, scratch);\n  const v0 = v[0], v1 = v[1], v2 = v[2], v3 = v[3];\n  let dv0 = 0, dv1 = 0, dv2 = 0, dv3 = 0;\n\n  for (let mu = 0; mu < 4; mu++) {\n    // Bolt: Hoist the invariant 16 * mu multiplication outside the nu, rho loops\n    const mu16 = mu * 16;\n    let acc = 0;\n    let accTerm = 0;\n    let g = 0;\n\n    // Bolt: Fully unroll the inner `nu` and `rho` loops for a fixed 4D spacetime\n    // to eliminate loop overhead and index arithmetic, caching `v` into scalar variables.\n\n    // nu = 0\n    let offset = mu16;\n    accTerm = 0;\n    g = G[offset]; if (g !== 0) accTerm += g * v0;\n    g = G[offset+1]; if (g !== 0) accTerm += g * v1;\n    g = G[offset+2]; if (g !== 0) accTerm += g * v2;\n    g = G[offset+3]; if (g !== 0) accTerm += g * v3;\n    if (accTerm !== 0) acc += accTerm * v0;\n\n    // nu = 1\n    offset += 4;\n    accTerm = 0;\n    g = G[offset]; if (g !== 0) accTerm += g * v0;\n    g = G[offset+1]; if (g !== 0) accTerm += g * v1;\n    g = G[offset+2]; if (g !== 0) accTerm += g * v2;\n    g = G[offset+3]; if (g !== 0) accTerm += g * v3;\n    if (accTerm !== 0) acc += accTerm * v1;\n\n    // nu = 2\n    offset += 4;\n    accTerm = 0;\n    g = G[offset]; if (g !== 0) accTerm += g * v0;\n    g = G[offset+1]; if (g !== 0) accTerm += g * v1;\n    g = G[offset+2]; if (g !== 0) accTerm += g * v2;\n    g = G[offset+3]; if (g !== 0) accTerm += g * v3;\n    if (accTerm !== 0) acc += accTerm * v2;\n\n    // nu = 3\n    offset += 4;\n    accTerm = 0;\n    g = G[offset]; if (g !== 0) accTerm += g * v0;\n    g = G[offset+1]; if (g !== 0) accTerm += g * v1;\n    g = G[offset+2]; if (g !== 0) accTerm += g * v2;\n    g = G[offset+3]; if (g !== 0) accTerm += g * v3;\n    if (accTerm !== 0) acc += accTerm * v3;\n\n    if (mu === 0) dv0 = -acc;\n    else if (mu === 1) dv1 = -acc;\n    else if (mu === 2) dv2 = -acc;\n    else dv3 = -acc;\n  }\n\n  return { dx: [v0, v1, v2, v3], dv: [dv0, dv1, dv2, dv3] };\n}\n\n// ---------------------------------------------------------------------------\n// Public API\n// ---------------------------------------------------------------------------\n\n/**\n * Numerically integrate the geodesic equation for a timelike (or null)\n * worldline in an arbitrary curved spacetime using a fixed-step classical\n * fourth-order Runge-Kutta scheme.\n *\n * @public\n */\nexport function integrateGeodesic(\n  inputs: GeodesicIntegratorInputs,\n): GeodesicIntegratorResult {\n  const { christoffelFn, x0, v0, tauStart, tauEnd, steps, domainMinRadius } = inputs;\n  // BR-2: christoffelFn now returns Float64Array(64); geodesicRHS reads flat layout.\n\n  // E11 fix: explicit option in place of monkey-patched .r_s on the closure.\n  // Callers (e.g. the Schwarzschild test fixture) pass `3 * r_s` directly.\n  if (typeof domainMinRadius === 'number' && x0[1] < domainMinRadius) {\n    throw new NumericalBackendError(\n      `integrateGeodesic: initial radius r0=${x0[1]} < domainMinRadius=${domainMinRadius}. ` +\n      `Position is outside the valid integration domain. ` +\n      `For Schwarzschild geodesics this restriction is r ≥ 3·r_Schwarz (v0.4.0).`,\n    );\n  }\n\n  if (!Number.isInteger(steps) || steps <= 0) {\n    throw new NumericalBackendError(\n      `integrateGeodesic: steps must be a positive integer, got ${steps}`,\n    );\n  }\n\n  const h = (tauEnd - tauStart) / steps;\n  const y0 = [x0[0], x0[1], x0[2], x0[3], v0[0], v0[1], v0[2], v0[3]];\n\n  // One scratch buffer reused across all Christoffel evaluations (4 per step).\n  // Each geodesicRHS consumes its Γ before the next call, so reuse is safe and\n  // avoids ~4·steps Float64Array(64) allocations per integration.\n  const scratch = new Float64Array(64);\n  const sol = solveODESystem(\n    (_tau, y) => {\n      const x: Vec4 = [y[0], y[1], y[2], y[3]];\n      const v: Vec4 = [y[4], y[5], y[6], y[7]];\n      const s = geodesicRHS(christoffelFn, x, v, scratch);\n      return [s.dx[0], s.dx[1], s.dx[2], s.dx[3], s.dv[0], s.dv[1], s.dv[2], s.dv[3]];\n    },\n    y0,\n    [tauStart, tauEnd],\n    { dt: h },\n  );\n\n  const asVec4 = (row: readonly number[]): Vec4 => [row[0], row[1], row[2], row[3]];\n  const sampleEvery = Math.max(1, Math.floor(steps / 100));\n  const trajectory: Vec4[] = [asVec4(sol.y[0] ?? y0)];\n  for (let i = 0; i < steps; i++) {\n    if ((i + 1) % sampleEvery === 0) {\n      const row = sol.y[i + 1] ?? sol.y[sol.y.length - 1] ?? y0;\n      trajectory.push(asVec4(row));\n    }\n  }\n  const last = sol.y[sol.y.length - 1] ?? y0;\n\n  return {\n    xFinal: asVec4(last),\n    vFinal: [last[4], last[5], last[6], last[7]],\n    trajectory,\n  };\n}\n"]}
+{"version":3,"file":"geodesic-integrator.js","sourceRoot":"","sources":["../../src/numerical/geodesic-integrator.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;;;;;;GAwBG;AACH,OAAO,EAAE,qBAAqB,EAAE,MAAM,aAAa,CAAC;AACpD,OAAO,EAAE,aAAa,EAAE,MAAM,qBAAqB,CAAC;AA+EpD;;;;;;GAMG;AACH,SAAS,WAAW,CAClB,aAA6E,EAC7E,CAAO,EACP,CAAO;AACP,yEAAyE;AACzE,4EAA4E;AAC5E,8DAA8D;AAC9D,OAAsB;IAEtB,MAAM,CAAC,GAAG,aAAa,CAAC,CAAC,EAAE,OAAO,CAAC,CAAC;IACpC,MAAM,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,EAAE,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,EAAE,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,EAAE,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;IACjD,IAAI,GAAG,GAAG,CAAC,EAAE,GAAG,GAAG,CAAC,EAAE,GAAG,GAAG,CAAC,EAAE,GAAG,GAAG,CAAC,CAAC;IAEvC,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,CAAC,EAAE,EAAE,EAAE,EAAE,CAAC;QAC9B,6EAA6E;QAC7E,MAAM,IAAI,GAAG,EAAE,GAAG,EAAE,CAAC;QACrB,IAAI,GAAG,GAAG,CAAC,CAAC;QACZ,IAAI,OAAO,GAAG,CAAC,CAAC;QAChB,IAAI,CAAC,GAAG,CAAC,CAAC;QAEV,6EAA6E;QAC7E,sFAAsF;QAEtF,SAAS;QACT,IAAI,MAAM,GAAG,IAAI,CAAC;QAClB,OAAO,GAAG,CAAC,CAAC;QACZ,CAAC,GAAG,CAAC,CAAC,MAAM,CAAC,CAAC;QAAC,IAAI,CAAC,KAAK,CAAC;YAAE,OAAO,IAAI,CAAC,GAAG,EAAE,CAAC;QAC9C,CAAC,GAAG,CAAC,CAAC,MAAM,GAAC,CAAC,CAAC,CAAC;QAAC,IAAI,CAAC,KAAK,CAAC;YAAE,OAAO,IAAI,CAAC,GAAG,EAAE,CAAC;QAChD,CAAC,GAAG,CAAC,CAAC,MAAM,GAAC,CAAC,CAAC,CAAC;QAAC,IAAI,CAAC,KAAK,CAAC;YAAE,OAAO,IAAI,CAAC,GAAG,EAAE,CAAC;QAChD,CAAC,GAAG,CAAC,CAAC,MAAM,GAAC,CAAC,CAAC,CAAC;QAAC,IAAI,CAAC,KAAK,CAAC;YAAE,OAAO,IAAI,CAAC,GAAG,EAAE,CAAC;QAChD,IAAI,OAAO,KAAK,CAAC;YAAE,GAAG,IAAI,OAAO,GAAG,EAAE,CAAC;QAEvC,SAAS;QACT,MAAM,IAAI,CAAC,CAAC;QACZ,OAAO,GAAG,CAAC,CAAC;QACZ,CAAC,GAAG,CAAC,CAAC,MAAM,CAAC,CAAC;QAAC,IAAI,CAAC,KAAK,CAAC;YAAE,OAAO,IAAI,CAAC,GAAG,EAAE,CAAC;QAC9C,CAAC,GAAG,CAAC,CAAC,MAAM,GAAC,CAAC,CAAC,CAAC;QAAC,IAAI,CAAC,KAAK,CAAC;YAAE,OAAO,IAAI,CAAC,GAAG,EAAE,CAAC;QAChD,CAAC,GAAG,CAAC,CAAC,MAAM,GAAC,CAAC,CAAC,CAAC;QAAC,IAAI,CAAC,KAAK,CAAC;YAAE,OAAO,IAAI,CAAC,GAAG,EAAE,CAAC;QAChD,CAAC,GAAG,CAAC,CAAC,MAAM,GAAC,CAAC,CAAC,CAAC;QAAC,IAAI,CAAC,KAAK,CAAC;YAAE,OAAO,IAAI,CAAC,GAAG,EAAE,CAAC;QAChD,IAAI,OAAO,KAAK,CAAC;YAAE,GAAG,IAAI,OAAO,GAAG,EAAE,CAAC;QAEvC,SAAS;QACT,MAAM,IAAI,CAAC,CAAC;QACZ,OAAO,GAAG,CAAC,CAAC;QACZ,CAAC,GAAG,CAAC,CAAC,MAAM,CAAC,CAAC;QAAC,IAAI,CAAC,KAAK,CAAC;YAAE,OAAO,IAAI,CAAC,GAAG,EAAE,CAAC;QAC9C,CAAC,GAAG,CAAC,CAAC,MAAM,GAAC,CAAC,CAAC,CAAC;QAAC,IAAI,CAAC,KAAK,CAAC;YAAE,OAAO,IAAI,CAAC,GAAG,EAAE,CAAC;QAChD,CAAC,GAAG,CAAC,CAAC,MAAM,GAAC,CAAC,CAAC,CAAC;QAAC,IAAI,CAAC,KAAK,CAAC;YAAE,OAAO,IAAI,CAAC,GAAG,EAAE,CAAC;QAChD,CAAC,GAAG,CAAC,CAAC,MAAM,GAAC,CAAC,CAAC,CAAC;QAAC,IAAI,CAAC,KAAK,CAAC;YAAE,OAAO,IAAI,CAAC,GAAG,EAAE,CAAC;QAChD,IAAI,OAAO,KAAK,CAAC;YAAE,GAAG,IAAI,OAAO,GAAG,EAAE,CAAC;QAEvC,SAAS;QACT,MAAM,IAAI,CAAC,CAAC;QACZ,OAAO,GAAG,CAAC,CAAC;QACZ,CAAC,GAAG,CAAC,CAAC,MAAM,CAAC,CAAC;QAAC,IAAI,CAAC,KAAK,CAAC;YAAE,OAAO,IAAI,CAAC,GAAG,EAAE,CAAC;QAC9C,CAAC,GAAG,CAAC,CAAC,MAAM,GAAC,CAAC,CAAC,CAAC;QAAC,IAAI,CAAC,KAAK,CAAC;YAAE,OAAO,IAAI,CAAC,GAAG,EAAE,CAAC;QAChD,CAAC,GAAG,CAAC,CAAC,MAAM,GAAC,CAAC,CAAC,CAAC;QAAC,IAAI,CAAC,KAAK,CAAC;YAAE,OAAO,IAAI,CAAC,GAAG,EAAE,CAAC;QAChD,CAAC,GAAG,CAAC,CAAC,MAAM,GAAC,CAAC,CAAC,CAAC;QAAC,IAAI,CAAC,KAAK,CAAC;YAAE,OAAO,IAAI,CAAC,GAAG,EAAE,CAAC;QAChD,IAAI,OAAO,KAAK,CAAC;YAAE,GAAG,IAAI,OAAO,GAAG,EAAE,CAAC;QAEvC,IAAI,EAAE,KAAK,CAAC;YAAE,GAAG,GAAG,CAAC,GAAG,CAAC;aACpB,IAAI,EAAE,KAAK,CAAC;YAAE,GAAG,GAAG,CAAC,GAAG,CAAC;aACzB,IAAI,EAAE,KAAK,CAAC;YAAE,GAAG,GAAG,CAAC,GAAG,CAAC;;YACzB,GAAG,GAAG,CAAC,GAAG,CAAC;IAClB,CAAC;IAED,OAAO,EAAE,EAAE,EAAE,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,EAAE,EAAE,EAAE,CAAC,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,CAAC,EAAE,CAAC;AAC5D,CAAC;AAED,8EAA8E;AAC9E,aAAa;AACb,8EAA8E;AAE9E;;;;;;GAMG;AACH,MAAM,UAAU,iBAAiB,CAC/B,MAAgC;IAEhC,MAAM,EAAE,aAAa,EAAE,EAAE,EAAE,EAAE,EAAE,QAAQ,EAAE,MAAM,EAAE,KAAK,EAAE,eAAe,EAAE,GAAG,MAAM,CAAC;IACnF,mFAAmF;IAEnF,2EAA2E;IAC3E,yEAAyE;IACzE,IAAI,OAAO,eAAe,KAAK,QAAQ,IAAI,EAAE,CAAC,CAAC,CAAC,GAAG,eAAe,EAAE,CAAC;QACnE,MAAM,IAAI,qBAAqB,CAC7B,wCAAwC,EAAE,CAAC,CAAC,CAAC,sBAAsB,eAAe,IAAI;YACtF,oDAAoD;YACpD,2EAA2E,CAC5E,CAAC;IACJ,CAAC;IAED,IAAI,CAAC,MAAM,CAAC,SAAS,CAAC,KAAK,CAAC,IAAI,KAAK,IAAI,CAAC,EAAE,CAAC;QAC3C,MAAM,IAAI,qBAAqB,CAC7B,4DAA4D,KAAK,EAAE,CACpE,CAAC;IACJ,CAAC;IACD,+EAA+E;IAC/E,6EAA6E;IAC7E,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,QAAQ,CAAC,EAAE,CAAC;QAC/B,MAAM,IAAI,qBAAqB,CAAC,mDAAmD,QAAQ,EAAE,CAAC,CAAC;IACjG,CAAC;IACD,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,MAAM,CAAC,EAAE,CAAC;QAC7B,MAAM,IAAI,qBAAqB,CAAC,iDAAiD,MAAM,EAAE,CAAC,CAAC;IAC7F,CAAC;IACD,IAAI,MAAM,KAAK,QAAQ,EAAE,CAAC;QACxB,MAAM,IAAI,qBAAqB,CAAC,iEAAiE,QAAQ,EAAE,CAAC,CAAC;IAC/G,CAAC;IAED,MAAM,EAAE,GAAG,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;IAEpE,6EAA6E;IAC7E,6EAA6E;IAC7E,gEAAgE;IAChE,MAAM,OAAO,GAAG,IAAI,YAAY,CAAC,EAAE,CAAC,CAAC;IACrC,MAAM,MAAM,GAAG,aAAa,CAC1B,CAAC,IAAI,EAAE,CAAC,EAAE,EAAE;QACV,MAAM,CAAC,GAAS,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QACzC,MAAM,CAAC,GAAS,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QACzC,MAAM,CAAC,GAAG,WAAW,CAAC,aAAa,EAAE,CAAC,EAAE,CAAC,EAAE,OAAO,CAAC,CAAC;QACpD,OAAO,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;IAClF,CAAC,EACD,EAAE,EACF,QAAQ,EACR,MAAM,EACN,KAAK,CACN,CAAC;IAEF,MAAM,MAAM,GAAG,CAAC,GAAsB,EAAQ,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;IAClF,MAAM,WAAW,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,IAAI,CAAC,KAAK,CAAC,KAAK,GAAG,GAAG,CAAC,CAAC,CAAC;IACzD,yFAAyF;IACzF,MAAM,UAAU,GAAW,CAAC,MAAM,CAAC,MAAM,CAAC,CAAC,CAAC,IAAI,EAAE,CAAC,CAAC,CAAC;IACrD,KAAK,IAAI,CAAC,GAAG,WAAW,EAAE,CAAC,IAAI,KAAK,EAAE,CAAC,IAAI,WAAW;QAAE,UAAU,CAAC,IAAI,CAAC,MAAM,CAAC,MAAM,CAAC,CAAC,CAAC,IAAI,EAAE,CAAC,CAAC,CAAC;IACjG,MAAM,IAAI,GAAG,MAAM,CAAC,KAAK,CAAC,IAAI,EAAE,CAAC;IAEjC,OAAO;QACL,MAAM,EAAE,MAAM,CAAC,IAAI,CAAC;QACpB,MAAM,EAAE,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,CAAC,CAAC;QAC5C,UAAU;KACX,CAAC;AACJ,CAAC","sourcesContent":["/**\n * Fixed-step RK4 integrator for the geodesic equation in an arbitrary\n * Lorentzian manifold (v0.4.0, Task 14 [U]).\n *\n * Solves the second-order ODE system\n *\n *   d²x^μ/dτ² + Γ^μ_{νρ}(x) (dx^ν/dτ)(dx^ρ/dτ) = 0\n *\n * by casting it as a first-order system on the (x, v) phase space:\n *\n *   dx^μ/dτ  = v^μ\n *   dv^μ/dτ  = −Γ^μ_{νρ}(x) v^ν v^ρ    (Einstein summation over ν, ρ)\n *\n * The Christoffel symbol Γ^μ_{νρ} is supplied as a plain-JS closure that\n * maps a 4-coordinate array to a [4][4][4] tensor.  No TensorEngine\n * dependency — the integrator is self-contained. The fixed steps are\n * `rk4Trajectory` (MathTS `solveODESystem` with `dt` set, that library's\n * classical RK4), which takes exactly `steps` steps. Adaptive `solveODE` is\n * not this step.\n *\n * Validated against the cycloid parametric form for radial Schwarzschild\n * infall to ±1e-6 (Task 14 test, η = 0.5).\n *\n * @module numerical/geodesic-integrator\n */\nimport { NumericalBackendError } from './errors.js';\nimport { rk4Trajectory } from './rk4-trajectory.js';\n\n/**\n * Input bundle for `integrateGeodesic`.\n * @public\n */\nexport interface GeodesicIntegratorInputs {\n  /**\n   * Christoffel-symbol closure.  Maps a 4-coordinate array x^μ to a\n   * Float64Array(64) containing Γ^μ_{νρ}(x) in λ-major layout:\n   * index (λ, μ, ν) → 16·λ + 4·μ + ν.\n   *\n   * BR-2 (v0.6.0 Task 2.9): migrated from nested number[][][] to flat\n   * Float64Array(64). Consumers call `arr[16*mu + 4*nu + rho]`.\n   *\n   * The optional `out` buffer lets the integrator pass a reused scratch\n   * Float64Array(64) so the hot loop avoids per-call allocation. A function\n   * that ignores `out` (e.g. allocates fresh each call) still satisfies this\n   * type — the optimization is opt-in.\n   */\n  christoffelFn: (x: ReadonlyArray<number>, out?: Float64Array) => Float64Array;\n\n  /** Initial 4-position x^μ(τ₀) = [t, r, θ, φ]. */\n  x0: readonly [number, number, number, number];\n\n  /**\n   * Initial 4-velocity dx^μ/dτ at τ₀.\n   *\n   * I2 fix: for a Schwarzschild radial timelike geodesic starting at\n   * coordinate radius r₀, use\n   *   v0 = [1/√(1 − r_s/r₀), 0, 0, 0]\n   * so that g_{μν} u^μ u^ν = −c² holds.  At r₀ = 100·r_s this gives\n   * dt/dτ ≈ 1.005, not 1.\n   */\n  v0: readonly [number, number, number, number];\n\n  /** Proper-time parameter at the start of integration. */\n  tauStart: number;\n\n  /** Proper-time parameter at the end of integration. */\n  tauEnd: number;\n\n  /** Number of fixed RK4 steps.  Must be a positive integer. */\n  steps: number;\n\n  /**\n   * E11 fix: explicit domain-validity option in place of the\n   * `(christoffelFn as any).r_s` monkey-patch.  If supplied, the\n   * integrator throws `NumericalBackendError` when `x0[1] < domainMinRadius`\n   * (e.g. for Schwarzschild: caller passes `3 * r_s` to enforce the\n   * v0.4.0 domain restriction r ≥ 3·r_Schwarz).\n   */\n  domainMinRadius?: number;\n}\n\n/**\n * Result of `integrateGeodesic`.\n * @public\n */\nexport interface GeodesicIntegratorResult {\n  /** Final 4-position x^μ(τEnd). */\n  readonly xFinal: readonly [number, number, number, number];\n\n  /** Final 4-velocity v^μ(τEnd). */\n  readonly vFinal: readonly [number, number, number, number];\n\n  /**\n   * Sampled trajectory positions (approximately 100 evenly-spaced snapshots\n   * plus the initial point).  Useful for visualisation and diagnostics.\n   */\n  readonly trajectory: ReadonlyArray<readonly [number, number, number, number]>;\n}\n\n// ---------------------------------------------------------------------------\n// Internal helpers\n// ---------------------------------------------------------------------------\n\ntype Vec4 = [number, number, number, number];\n\n/**\n * Evaluate the Christoffel acceleration −Γ^μ_{νρ} v^ν v^ρ at position x\n * with 4-velocity v.  Returns (dx/dτ, dv/dτ) = (v, accel).\n *\n * BR-2 (v0.6.0 Task 2.9): christoffelFn now returns Float64Array(64).\n * Index access: G[16*mu + 4*nu + rho] (λ-major layout).\n */\nfunction geodesicRHS(\n  christoffelFn: (x: ReadonlyArray<number>, out?: Float64Array) => Float64Array,\n  x: Vec4,\n  v: Vec4,\n  // Reused scratch buffer for the Christoffel evaluation. Safe: G is fully\n  // consumed into `dv` below before geodesicRHS returns, so the next call may\n  // overwrite it. `dx`/`dv` are fresh arrays and never alias G.\n  scratch?: Float64Array,\n): { dx: Vec4; dv: Vec4 } {\n  const G = christoffelFn(x, scratch);\n  const v0 = v[0], v1 = v[1], v2 = v[2], v3 = v[3];\n  let dv0 = 0, dv1 = 0, dv2 = 0, dv3 = 0;\n\n  for (let mu = 0; mu < 4; mu++) {\n    // Bolt: Hoist the invariant 16 * mu multiplication outside the nu, rho loops\n    const mu16 = mu * 16;\n    let acc = 0;\n    let accTerm = 0;\n    let g = 0;\n\n    // Bolt: Fully unroll the inner `nu` and `rho` loops for a fixed 4D spacetime\n    // to eliminate loop overhead and index arithmetic, caching `v` into scalar variables.\n\n    // nu = 0\n    let offset = mu16;\n    accTerm = 0;\n    g = G[offset]; if (g !== 0) accTerm += g * v0;\n    g = G[offset+1]; if (g !== 0) accTerm += g * v1;\n    g = G[offset+2]; if (g !== 0) accTerm += g * v2;\n    g = G[offset+3]; if (g !== 0) accTerm += g * v3;\n    if (accTerm !== 0) acc += accTerm * v0;\n\n    // nu = 1\n    offset += 4;\n    accTerm = 0;\n    g = G[offset]; if (g !== 0) accTerm += g * v0;\n    g = G[offset+1]; if (g !== 0) accTerm += g * v1;\n    g = G[offset+2]; if (g !== 0) accTerm += g * v2;\n    g = G[offset+3]; if (g !== 0) accTerm += g * v3;\n    if (accTerm !== 0) acc += accTerm * v1;\n\n    // nu = 2\n    offset += 4;\n    accTerm = 0;\n    g = G[offset]; if (g !== 0) accTerm += g * v0;\n    g = G[offset+1]; if (g !== 0) accTerm += g * v1;\n    g = G[offset+2]; if (g !== 0) accTerm += g * v2;\n    g = G[offset+3]; if (g !== 0) accTerm += g * v3;\n    if (accTerm !== 0) acc += accTerm * v2;\n\n    // nu = 3\n    offset += 4;\n    accTerm = 0;\n    g = G[offset]; if (g !== 0) accTerm += g * v0;\n    g = G[offset+1]; if (g !== 0) accTerm += g * v1;\n    g = G[offset+2]; if (g !== 0) accTerm += g * v2;\n    g = G[offset+3]; if (g !== 0) accTerm += g * v3;\n    if (accTerm !== 0) acc += accTerm * v3;\n\n    if (mu === 0) dv0 = -acc;\n    else if (mu === 1) dv1 = -acc;\n    else if (mu === 2) dv2 = -acc;\n    else dv3 = -acc;\n  }\n\n  return { dx: [v0, v1, v2, v3], dv: [dv0, dv1, dv2, dv3] };\n}\n\n// ---------------------------------------------------------------------------\n// Public API\n// ---------------------------------------------------------------------------\n\n/**\n * Numerically integrate the geodesic equation for a timelike (or null)\n * worldline in an arbitrary curved spacetime using a fixed-step classical\n * fourth-order Runge-Kutta scheme.\n *\n * @public\n */\nexport function integrateGeodesic(\n  inputs: GeodesicIntegratorInputs,\n): GeodesicIntegratorResult {\n  const { christoffelFn, x0, v0, tauStart, tauEnd, steps, domainMinRadius } = inputs;\n  // BR-2: christoffelFn now returns Float64Array(64); geodesicRHS reads flat layout.\n\n  // E11 fix: explicit option in place of monkey-patched .r_s on the closure.\n  // Callers (e.g. the Schwarzschild test fixture) pass `3 * r_s` directly.\n  if (typeof domainMinRadius === 'number' && x0[1] < domainMinRadius) {\n    throw new NumericalBackendError(\n      `integrateGeodesic: initial radius r0=${x0[1]} < domainMinRadius=${domainMinRadius}. ` +\n      `Position is outside the valid integration domain. ` +\n      `For Schwarzschild geodesics this restriction is r ≥ 3·r_Schwarz (v0.4.0).`,\n    );\n  }\n\n  if (!Number.isInteger(steps) || steps <= 0) {\n    throw new NumericalBackendError(\n      `integrateGeodesic: steps must be a positive integer, got ${steps}`,\n    );\n  }\n  // A NaN or infinite end made every step NaN and the result the start state; an\n  // end equal to the start is a zero step, which is not an integration either.\n  if (!Number.isFinite(tauStart)) {\n    throw new NumericalBackendError(`integrateGeodesic: tauStart must be finite, got ${tauStart}`);\n  }\n  if (!Number.isFinite(tauEnd)) {\n    throw new NumericalBackendError(`integrateGeodesic: tauEnd must be finite, got ${tauEnd}`);\n  }\n  if (tauEnd === tauStart) {\n    throw new NumericalBackendError(`integrateGeodesic: tauEnd must differ from tauStart, both are ${tauStart}`);\n  }\n\n  const y0 = [x0[0], x0[1], x0[2], x0[3], v0[0], v0[1], v0[2], v0[3]];\n\n  // One scratch buffer reused across all Christoffel evaluations (4 per step).\n  // Each geodesicRHS consumes its Γ before the next call, so reuse is safe and\n  // avoids ~4·steps Float64Array(64) allocations per integration.\n  const scratch = new Float64Array(64);\n  const states = rk4Trajectory(\n    (_tau, y) => {\n      const x: Vec4 = [y[0], y[1], y[2], y[3]];\n      const v: Vec4 = [y[4], y[5], y[6], y[7]];\n      const s = geodesicRHS(christoffelFn, x, v, scratch);\n      return [s.dx[0], s.dx[1], s.dx[2], s.dx[3], s.dv[0], s.dv[1], s.dv[2], s.dv[3]];\n    },\n    y0,\n    tauStart,\n    tauEnd,\n    steps,\n  );\n\n  const asVec4 = (row: readonly number[]): Vec4 => [row[0], row[1], row[2], row[3]];\n  const sampleEvery = Math.max(1, Math.floor(steps / 100));\n  // `rk4Trajectory` returns exactly steps + 1 states, so state i is the one after i steps.\n  const trajectory: Vec4[] = [asVec4(states[0] ?? y0)];\n  for (let i = sampleEvery; i <= steps; i += sampleEvery) trajectory.push(asVec4(states[i] ?? y0));\n  const last = states[steps] ?? y0;\n\n  return {\n    xFinal: asVec4(last),\n    vFinal: [last[4], last[5], last[6], last[7]],\n    trajectory,\n  };\n}\n"]}

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

@@ -1 +1 @@
-{"version":3,"file":"geometrized.d.ts","sourceRoot":"","sources":["../../src/numerical/geometrized.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;;;;;;GAwBG;AAEH,OAAO,KAAK,EAAE,SAAS,EAAE,MAAM,yBAAyB,CAAC;AACzD,OAAO,EAAE,QAAQ,EAAE,MAAM,0BAA0B,CAAC;AAGpD;;;;;GAKG;AACH,qBAAa,8BAA+B,SAAQ,QAAQ;IAC1D,YAAY,OAAO,EAAE,MAAM,EAI1B;CACF;AAED;;;;;;;;GAQG;AACH,wBAAgB,iBAAiB,CAAC,GAAG,EAAE,SAAS,GAAG,MAAM,CASxD;AAED,4CAA4C;AAC5C,wBAAgB,aAAa,CAAC,OAAO,EAAE,MAAM,EAAE,GAAG,EAAE,SAAS,GAAG,MAAM,CAErE;AAED,4CAA4C;AAC5C,wBAAgB,eAAe,CAAC,SAAS,EAAE,MAAM,EAAE,GAAG,EAAE,SAAS,GAAG,MAAM,CAEzE"}
+{"version":3,"file":"geometrized.d.ts","sourceRoot":"","sources":["../../src/numerical/geometrized.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;;;;;;GAwBG;AAEH,OAAO,KAAK,EAAE,SAAS,EAAE,MAAM,yBAAyB,CAAC;AACzD,OAAO,EAAE,QAAQ,EAAE,MAAM,0BAA0B,CAAC;AAGpD;;;;;GAKG;AACH,qBAAa,8BAA+B,SAAQ,QAAQ;IAC1D,YAAY,OAAO,EAAE,MAAM,EAG1B;CACF;AAED;;;;;;;;GAQG;AACH,wBAAgB,iBAAiB,CAAC,GAAG,EAAE,SAAS,GAAG,MAAM,CASxD;AAED,4CAA4C;AAC5C,wBAAgB,aAAa,CAAC,OAAO,EAAE,MAAM,EAAE,GAAG,EAAE,SAAS,GAAG,MAAM,CAErE;AAED,4CAA4C;AAC5C,wBAAgB,eAAe,CAAC,SAAS,EAAE,MAAM,EAAE,GAAG,EAAE,SAAS,GAAG,MAAM,CAEzE"}

package/dist/numerical/geometrized.js

@@ -35,7 +35,6 @@ export class NonGeometrizableDimensionError extends UPTError {
     constructor(message) {
         super(message);
         this.name = 'NonGeometrizableDimensionError';
-        Object.setPrototypeOf(this, NonGeometrizableDimensionError.prototype);
     }
 }
 /**

package/dist/numerical/geometrized.js.map

@@ -1 +1 @@
-{"version":3,"file":"geometrized.js","sourceRoot":"","sources":["../../src/numerical/geometrized.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;;;;;;GAwBG;AAGH,OAAO,EAAE,QAAQ,EAAE,MAAM,0BAA0B,CAAC;AACpD,OAAO,EAAE,IAAI,EAAE,IAAI,EAAE,MAAM,sBAAsB,CAAC;AAElD;;;;;GAKG;AACH,MAAM,OAAO,8BAA+B,SAAQ,QAAQ;IAC1D,YAAY,OAAe;QACzB,KAAK,CAAC,OAAO,CAAC,CAAC;QACf,IAAI,CAAC,IAAI,GAAG,gCAAgC,CAAC;QAC7C,MAAM,CAAC,cAAc,CAAC,IAAI,EAAE,8BAA8B,CAAC,SAAS,CAAC,CAAC;IACxE,CAAC;CACF;AAED;;;;;;;;GAQG;AACH,MAAM,UAAU,iBAAiB,CAAC,GAAc;IAC9C,IAAI,GAAG,CAAC,CAAC,KAAK,CAAC,IAAI,GAAG,CAAC,KAAK,KAAK,CAAC,IAAI,GAAG,CAAC,CAAC,KAAK,CAAC,IAAI,GAAG,CAAC,CAAC,KAAK,CAAC,EAAE,CAAC;QACjE,MAAM,IAAI,8BAA8B,CACtC,qEAAqE;YACnE,wCAAwC,GAAG,CAAC,CAAC,OAAO,GAAG,CAAC,KAAK,IAAI;YACjE,KAAK,GAAG,CAAC,CAAC,OAAO,GAAG,CAAC,CAAC,iDAAiD,CAC1E,CAAC;IACJ,CAAC;IACD,OAAO,IAAI,CAAC,GAAG,CAAC,IAAI,EAAE,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC;AACnE,CAAC;AAED,4CAA4C;AAC5C,MAAM,UAAU,aAAa,CAAC,OAAe,EAAE,GAAc;IAC3D,OAAO,OAAO,GAAG,iBAAiB,CAAC,GAAG,CAAC,CAAC;AAC1C,CAAC;AAED,4CAA4C;AAC5C,MAAM,UAAU,eAAe,CAAC,SAAiB,EAAE,GAAc;IAC/D,OAAO,SAAS,GAAG,iBAAiB,CAAC,GAAG,CAAC,CAAC;AAC5C,CAAC","sourcesContent":["/**\n * Geometrized-units boundary adapters (G-9 increment 1).\n *\n * Converts a scalar between SI and geometrized (G = c = 1) units, driven\n * MECHANICALLY by the `Dimension` exponent vector — the dimension functor\n * earning its keep. For an SI quantity of dimension L^l M^m T^t the single\n * conversion factor is `G^m · c^(t − 2m)` (each kg → G/c² metres, each second\n * → c metres); the geometrized quantity is a pure length L^(l+m+t).\n *\n * This is the self-contained FOUNDATION of the units-normalization layer\n * (design: docs/planning/v0.10.0-Units-Normalization-Design-Note.md, Adam-\n * vetted r2; plan: docs/planning/v0.13-G9-Adapters-Plan.md, Eve-vetted). Only c\n * and G are in scope: a nonzero electromagnetic / thermal / molar / luminous\n * exponent (I/Θ/N/J) throws.\n *\n * PUBLIC boundary API (G-9 increment 2): a consumer running the GR pipeline in\n * geometrized (G = c = 1) units converts scalars at the boundary with these\n * adapters. The consumer-wide DEFAULT-pipeline migration (increment 3) was\n * DECLINED after adversarial review — it is measured to give no precision win\n * (and far worse far-field; see docs/planning/v0.14-G9-Increment3-Disposition.md).\n * The SI pipeline stays the default; this geometrized layer is the optional\n * boundary API.\n *\n * @module numerical/geometrized\n */\n\nimport type { Dimension } from '../dimensional/types.js';\nimport { UPTError } from '../dimensional/errors.js';\nimport { C_SI, G_SI } from '../core/constants.js';\n\n/**\n * A dimension carries a non-c/G base (I/Θ/N/J ≠ 0) and cannot be geometrized\n * by c and G alone (k_B / ε₀ extensions are out of scope for this increment).\n *\n * @public\n */\nexport class NonGeometrizableDimensionError extends UPTError {\n  constructor(message: string) {\n    super(message);\n    this.name = 'NonGeometrizableDimensionError';\n    Object.setPrototypeOf(this, NonGeometrizableDimensionError.prototype);\n  }\n}\n\n/**\n * The single SI→geometrized conversion factor for `dim`: `G^M · c^(T − 2M)`\n * (only the mass and time exponents enter). Shared by `toGeometrized` /\n * `fromGeometrized` so the round-trip is a clean inverse. A nonzero I/Θ/N/J\n * exponent throws `NonGeometrizableDimensionError` BEFORE any factor is\n * computed.\n *\n * @public\n */\nexport function geometrizedFactor(dim: Dimension): number {\n  if (dim.I !== 0 || dim.Theta !== 0 || dim.N !== 0 || dim.J !== 0) {\n    throw new NonGeometrizableDimensionError(\n      `geometrizedFactor: only c and G are in scope; this dimension has a ` +\n        `nonzero non-mechanizable exponent (I=${dim.I}, Θ=${dim.Theta}, ` +\n        `N=${dim.N}, J=${dim.J}). Geometrization handles length / mass / time.`,\n    );\n  }\n  return Math.pow(G_SI, dim.M) * Math.pow(C_SI, dim.T - 2 * dim.M);\n}\n\n/** SI → geometrized (G = c = 1). @public */\nexport function toGeometrized(valueSI: number, dim: Dimension): number {\n  return valueSI * geometrizedFactor(dim);\n}\n\n/** Geometrized (G = c = 1) → SI. @public */\nexport function fromGeometrized(valueGeom: number, dim: Dimension): number {\n  return valueGeom / geometrizedFactor(dim);\n}\n"]}
+{"version":3,"file":"geometrized.js","sourceRoot":"","sources":["../../src/numerical/geometrized.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;;;;;;GAwBG;AAGH,OAAO,EAAE,QAAQ,EAAE,MAAM,0BAA0B,CAAC;AACpD,OAAO,EAAE,IAAI,EAAE,IAAI,EAAE,MAAM,sBAAsB,CAAC;AAElD;;;;;GAKG;AACH,MAAM,OAAO,8BAA+B,SAAQ,QAAQ;IAC1D,YAAY,OAAe;QACzB,KAAK,CAAC,OAAO,CAAC,CAAC;QACf,IAAI,CAAC,IAAI,GAAG,gCAAgC,CAAC;IAC/C,CAAC;CACF;AAED;;;;;;;;GAQG;AACH,MAAM,UAAU,iBAAiB,CAAC,GAAc;IAC9C,IAAI,GAAG,CAAC,CAAC,KAAK,CAAC,IAAI,GAAG,CAAC,KAAK,KAAK,CAAC,IAAI,GAAG,CAAC,CAAC,KAAK,CAAC,IAAI,GAAG,CAAC,CAAC,KAAK,CAAC,EAAE,CAAC;QACjE,MAAM,IAAI,8BAA8B,CACtC,qEAAqE;YACnE,wCAAwC,GAAG,CAAC,CAAC,OAAO,GAAG,CAAC,KAAK,IAAI;YACjE,KAAK,GAAG,CAAC,CAAC,OAAO,GAAG,CAAC,CAAC,iDAAiD,CAC1E,CAAC;IACJ,CAAC;IACD,OAAO,IAAI,CAAC,GAAG,CAAC,IAAI,EAAE,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC;AACnE,CAAC;AAED,4CAA4C;AAC5C,MAAM,UAAU,aAAa,CAAC,OAAe,EAAE,GAAc;IAC3D,OAAO,OAAO,GAAG,iBAAiB,CAAC,GAAG,CAAC,CAAC;AAC1C,CAAC;AAED,4CAA4C;AAC5C,MAAM,UAAU,eAAe,CAAC,SAAiB,EAAE,GAAc;IAC/D,OAAO,SAAS,GAAG,iBAAiB,CAAC,GAAG,CAAC,CAAC;AAC5C,CAAC","sourcesContent":["/**\n * Geometrized-units boundary adapters (G-9 increment 1).\n *\n * Converts a scalar between SI and geometrized (G = c = 1) units, driven\n * MECHANICALLY by the `Dimension` exponent vector — the dimension functor\n * earning its keep. For an SI quantity of dimension L^l M^m T^t the single\n * conversion factor is `G^m · c^(t − 2m)` (each kg → G/c² metres, each second\n * → c metres); the geometrized quantity is a pure length L^(l+m+t).\n *\n * This is the self-contained FOUNDATION of the units-normalization layer\n * (design: docs/planning/v0.10.0-Units-Normalization-Design-Note.md, Adam-\n * vetted r2; plan: docs/planning/v0.13-G9-Adapters-Plan.md, Eve-vetted). Only c\n * and G are in scope: a nonzero electromagnetic / thermal / molar / luminous\n * exponent (I/Θ/N/J) throws.\n *\n * PUBLIC boundary API (G-9 increment 2): a consumer running the GR pipeline in\n * geometrized (G = c = 1) units converts scalars at the boundary with these\n * adapters. The consumer-wide DEFAULT-pipeline migration (increment 3) was\n * DECLINED after adversarial review — it is measured to give no precision win\n * (and far worse far-field; see docs/planning/v0.14-G9-Increment3-Disposition.md).\n * The SI pipeline stays the default; this geometrized layer is the optional\n * boundary API.\n *\n * @module numerical/geometrized\n */\n\nimport type { Dimension } from '../dimensional/types.js';\nimport { UPTError } from '../dimensional/errors.js';\nimport { C_SI, G_SI } from '../core/constants.js';\n\n/**\n * A dimension carries a non-c/G base (I/Θ/N/J ≠ 0) and cannot be geometrized\n * by c and G alone (k_B / ε₀ extensions are out of scope for this increment).\n *\n * @public\n */\nexport class NonGeometrizableDimensionError extends UPTError {\n  constructor(message: string) {\n    super(message);\n    this.name = 'NonGeometrizableDimensionError';\n  }\n}\n\n/**\n * The single SI→geometrized conversion factor for `dim`: `G^M · c^(T − 2M)`\n * (only the mass and time exponents enter). Shared by `toGeometrized` /\n * `fromGeometrized` so the round-trip is a clean inverse. A nonzero I/Θ/N/J\n * exponent throws `NonGeometrizableDimensionError` BEFORE any factor is\n * computed.\n *\n * @public\n */\nexport function geometrizedFactor(dim: Dimension): number {\n  if (dim.I !== 0 || dim.Theta !== 0 || dim.N !== 0 || dim.J !== 0) {\n    throw new NonGeometrizableDimensionError(\n      `geometrizedFactor: only c and G are in scope; this dimension has a ` +\n        `nonzero non-mechanizable exponent (I=${dim.I}, Θ=${dim.Theta}, ` +\n        `N=${dim.N}, J=${dim.J}). Geometrization handles length / mass / time.`,\n    );\n  }\n  return Math.pow(G_SI, dim.M) * Math.pow(C_SI, dim.T - 2 * dim.M);\n}\n\n/** SI → geometrized (G = c = 1). @public */\nexport function toGeometrized(valueSI: number, dim: Dimension): number {\n  return valueSI * geometrizedFactor(dim);\n}\n\n/** Geometrized (G = c = 1) → SI. @public */\nexport function fromGeometrized(valueGeom: number, dim: Dimension): number {\n  return valueGeom / geometrizedFactor(dim);\n}\n"]}

package/dist/numerical/mathts-engine.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"mathts-engine.d.ts","sourceRoot":"","sources":["../../src/numerical/mathts-engine.ts"],"names":[],"mappings":"AAYA,OAAO,KAAK,EAAE,YAAY,EAAE,YAAY,EAAE,UAAU,EAAE,iBAAiB,EAAE,iBAAiB,EAAE,MAAM,oBAAoB,CAAC;AACvH,OAAO,KAAK,EAAE,WAAW,EAAE,MAAM,YAAY,CAAC;AAmC9C;;;;;;;GAOG;AACH,qBAAa,YAAa,YAAW,YAAY;IAC/C,QAAQ,CAAC,IAAI,kBAAkB;IAE/B,UAAU,CAAC,IAAI,EAAE,WAAW,EAAE,KAAK,EAAE,aAAa,CAAC,MAAM,CAAC,GAAG,YAAY,CAGxE;IACD,QAAQ,CAAC,CAAC,EAAE,YAAY,GAAG,WAAW,CAErC;IAED,MAAM,CAAC,IAAI,EAAE,UAAU,EAAE,GAAG,QAAQ,EAAE,YAAY,EAAE,GAAG,YAAY,CAGlE;IACD,MAAM,CAAC,CAAC,EAAE,YAAY,EAAE,CAAC,EAAE,YAAY,GAAG,YAAY,CAErD;IACD,SAAS,CAAC,CAAC,EAAE,YAAY,EAAE,IAAI,CAAC,EAAE,aAAa,CAAC,MAAM,CAAC,GAAG,YAAY,CAErE;IACD,OAAO,CAAC,CAAC,EAAE,YAAY,EAAE,KAAK,EAAE,aAAa,CAAC,MAAM,CAAC,GAAG,YAAY,CAEnE;IASD,GAAG,CAAC,CAAC,EAAE,YAAY,EAAE,CAAC,EAAE,YAAY,GAAG,YAAY,CASlD;IACD,GAAG,CAAC,CAAC,EAAE,YAAY,EAAE,CAAC,EAAE,YAAY,GAAG,YAAY,CAQlD;IACD,GAAG,CAAC,CAAC,EAAE,YAAY,EAAE,CAAC,EAAE,YAAY,GAAG,YAAY,CAQlD;IACD,KAAK,CAAC,CAAC,EAAE,YAAY,EAAE,CAAC,EAAE,MAAM,GAAG,YAAY,CAK9C;IAED,QAAQ,CAAC,CAAC,EAAE,MAAM,GAAG,YAAY,CAEhC;IACD,OAAO,CAAC,CAAC,EAAE,YAAY,GAAG,MAAM,CAE/B;IAMD,oFAAoF;IACpF,OAAO,CAAC,cAAc;IAItB,2EAA2E;IAC3E,OAAO,CAAC,gBAAgB;IAQxB;;;;;;;;;;;OAWG;IACG,WAAW,CACf,EAAE,EAAE,CAAC,CAAC,EAAE,YAAY,KAAK,YAAY,EACrC,CAAC,EAAE,YAAY,GACd,OAAO,CAAC,iBAAiB,CAAC,CAkB5B;IAMD;;;;;;;;;OASG;IACG,WAAW,CACf,EAAE,EAAE,CAAC,CAAC,EAAE,YAAY,KAAK,YAAY,EACrC,CAAC,EAAE,YAAY,EACf,SAAS,CAAC,EAAE,YAAY,GACvB,OAAO,CAAC,iBAAiB,CAAC,CAiB5B;CACF"}
+{"version":3,"file":"mathts-engine.d.ts","sourceRoot":"","sources":["../../src/numerical/mathts-engine.ts"],"names":[],"mappings":"AAYA,OAAO,KAAK,EAAE,YAAY,EAAE,YAAY,EAAE,UAAU,EAAE,iBAAiB,EAAE,iBAAiB,EAAE,MAAM,oBAAoB,CAAC;AACvH,OAAO,KAAK,EAAE,WAAW,EAAE,MAAM,YAAY,CAAC;AAmC9C;;;;;;;GAOG;AACH,qBAAa,YAAa,YAAW,YAAY;IAC/C,QAAQ,CAAC,IAAI,kBAAkB;IAE/B,UAAU,CAAC,IAAI,EAAE,WAAW,EAAE,KAAK,EAAE,aAAa,CAAC,MAAM,CAAC,GAAG,YAAY,CAGxE;IACD,QAAQ,CAAC,CAAC,EAAE,YAAY,GAAG,WAAW,CAErC;IAED,MAAM,CAAC,IAAI,EAAE,UAAU,EAAE,GAAG,QAAQ,EAAE,YAAY,EAAE,GAAG,YAAY,CAGlE;IACD,MAAM,CAAC,CAAC,EAAE,YAAY,EAAE,CAAC,EAAE,YAAY,GAAG,YAAY,CAErD;IACD,SAAS,CAAC,CAAC,EAAE,YAAY,EAAE,IAAI,CAAC,EAAE,aAAa,CAAC,MAAM,CAAC,GAAG,YAAY,CAErE;IACD,OAAO,CAAC,CAAC,EAAE,YAAY,EAAE,KAAK,EAAE,aAAa,CAAC,MAAM,CAAC,GAAG,YAAY,CAEnE;IASD,GAAG,CAAC,CAAC,EAAE,YAAY,EAAE,CAAC,EAAE,YAAY,GAAG,YAAY,CAGlD;IACD,GAAG,CAAC,CAAC,EAAE,YAAY,EAAE,CAAC,EAAE,YAAY,GAAG,YAAY,CAGlD;IACD,GAAG,CAAC,CAAC,EAAE,YAAY,EAAE,CAAC,EAAE,YAAY,GAAG,YAAY,CAGlD;IACD,KAAK,CAAC,CAAC,EAAE,YAAY,EAAE,CAAC,EAAE,MAAM,GAAG,YAAY,CAG9C;IAED,QAAQ,CAAC,CAAC,EAAE,MAAM,GAAG,YAAY,CAEhC;IACD,OAAO,CAAC,CAAC,EAAE,YAAY,GAAG,MAAM,CAE/B;IAMD,oFAAoF;IACpF,OAAO,CAAC,cAAc;IAItB,2EAA2E;IAC3E,OAAO,CAAC,gBAAgB;IAQxB;;;;;;;;;;;OAWG;IACG,WAAW,CACf,EAAE,EAAE,CAAC,CAAC,EAAE,YAAY,KAAK,YAAY,EACrC,CAAC,EAAE,YAAY,GACd,OAAO,CAAC,iBAAiB,CAAC,CAkB5B;IAMD;;;;;;;;;OASG;IACG,WAAW,CACf,EAAE,EAAE,CAAC,CAAC,EAAE,YAAY,KAAK,YAAY,EACrC,CAAC,EAAE,YAAY,EACf,SAAS,CAAC,EAAE,YAAY,GACvB,OAAO,CAAC,iBAAiB,CAAC,CAiB5B;CACF"}

package/dist/numerical/mathts-engine.js

@@ -11,7 +11,7 @@
  */
 import { Tensor } from '@danielsimonjr/mathts-tensor';
 import { NumericalBackendError } from './errors.js';
-import { EngineCapabilityError } from './tensor-engine.js';
+import { EngineCapabilityError, isTracedTensor, tracedMode } from './tensor-engine.js';
 /** EngineTensor handle wrapping a MathTS Tensor. */
 class MathTSEngineTensor {
     inner;
@@ -58,43 +58,29 @@ export class MathTSEngine {
     }
     // AD dispatch. `@danielsimonjr/mathts-autograd` wraps the traced input as a
     // DualTensor (forward mode, carries `tangent`) or a TapedTensor (reverse
-    // mode, carries `tape`) and calls these methods with it. The autograd
-    // package is loaded at the AD call, not at module load, so its classes
-    // cannot be named here; the structural casts below state the invariant
-    // the package keeps: a wrapped tensor implements add/sub/mul/scale over
-    // another wrapped tensor and returns a wrapped tensor.
+    // mode, carries `tape`) and calls these methods with it. The operations a
+    // wrapped tensor implements are declared on `TracedEngineTensor`; the guard
+    // `isTracedTensor` is the dispatch. Two operands trace together only when
+    // both are wrapped in the same mode; a mixed or half-wrapped pair falls
+    // through to the plain path, whose `unwrap` refuses it.
     add(a, b) {
-        if ('tangent' in a && 'tangent' in b) {
+        if (isTracedTensor(a) && isTracedTensor(b) && tracedMode(a) === tracedMode(b))
             return a.add(b);
-        }
-        // AD dispatch: TapedTensor (reverse-mode)
-        if ('tape' in a && 'tape' in b) {
-            return a.add(b);
-        }
         return new MathTSEngineTensor(unwrap(a, 'add').add(unwrap(b, 'add')));
     }
     sub(a, b) {
-        if ('tangent' in a && 'tangent' in b) {
-            return a.sub(b);
-        }
-        if ('tape' in a && 'tape' in b) {
+        if (isTracedTensor(a) && isTracedTensor(b) && tracedMode(a) === tracedMode(b))
             return a.sub(b);
-        }
         return new MathTSEngineTensor(unwrap(a, 'sub').sub(unwrap(b, 'sub')));
     }
     mul(a, b) {
-        if ('tangent' in a && 'tangent' in b) {
-            return a.mul(b);
-        }
-        if ('tape' in a && 'tape' in b) {
+        if (isTracedTensor(a) && isTracedTensor(b) && tracedMode(a) === tracedMode(b))
             return a.mul(b);
-        }
         return new MathTSEngineTensor(unwrap(a, 'mul').mul(unwrap(b, 'mul')));
     }
     scale(t, k) {
-        if ('tangent' in t || 'tape' in t) {
+        if (isTracedTensor(t))
             return t.scale(k);
-        }
         return new MathTSEngineTensor(unwrap(t, 'scale').scale(k));
     }
     identity(n) {

package/dist/numerical/mathts-engine.js.map

@@ -1 +1 @@
-{"version":3,"file":"mathts-engine.js","sourceRoot":"","sources":["../../src/numerical/mathts-engine.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;GAUG;AACH,OAAO,EAAE,MAAM,EAAE,MAAM,8BAA8B,CAAC;AAGtD,OAAO,EAAE,qBAAqB,EAAE,MAAM,aAAa,CAAC;AACpD,OAAO,EAAE,qBAAqB,EAAE,MAAM,oBAAoB,CAAC;AAE3D,oDAAoD;AACpD,MAAM,kBAAkB;IACD,KAAK;IAA1B,YAAqB,KAAa;qBAAb,KAAK;IAAW,CAAC;IACtC,IAAI,KAAK,KAA4B,OAAO,IAAI,CAAC,KAAK,CAAC,KAAK,CAAC,CAAC,CAAC;CAChE;AAED,SAAS,MAAM,CAAC,CAAe,EAAE,EAAU;IACzC,IAAI,CAAC,CAAC,CAAC,YAAY,kBAAkB,CAAC,EAAE,CAAC;QACvC,MAAM,IAAI,qBAAqB,CAAC,gBAAgB,EAAE,uCAAuC,CAAC,CAAC;IAC7F,CAAC;IACD,OAAO,CAAC,CAAC,KAAK,CAAC;AACjB,CAAC;AAoBD;;;;;;;GAOG;AACH,MAAM,OAAO,YAAY;IACd,IAAI,GAAG,cAAc,CAAC;IAE/B,UAAU,CAAC,IAAiB,EAAE,KAA4B;QACxD,iFAAiF;QACjF,OAAO,IAAI,kBAAkB,CAAC,MAAM,CAAC,UAAU,CAAC,IAAa,EAAE,KAAK,CAAC,CAAC,CAAC;IACzE,CAAC;IACD,QAAQ,CAAC,CAAe;QACtB,OAAO,MAAM,CAAC,CAAC,EAAE,UAAU,CAAC,CAAC,QAAQ,EAAE,CAAC;IAC1C,CAAC;IAED,MAAM,CAAC,IAAgB,EAAE,GAAG,QAAwB;QAClD,MAAM,KAAK,GAAG,QAAQ,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,MAAM,CAAC,CAAC,EAAE,mBAAmB,CAAC,GAAG,CAAC,CAAC,CAAC;QACzE,OAAO,IAAI,kBAAkB,CAAC,MAAM,CAAC,MAAM,CAAC,IAAI,EAAE,GAAG,KAAK,CAAC,CAAC,CAAC;IAC/D,CAAC;IACD,MAAM,CAAC,CAAe,EAAE,CAAe;QACrC,OAAO,IAAI,kBAAkB,CAAC,MAAM,CAAC,CAAC,EAAE,QAAQ,CAAC,CAAC,MAAM,CAAC,MAAM,CAAC,CAAC,EAAE,QAAQ,CAAC,CAAC,CAAC,CAAC;IACjF,CAAC;IACD,SAAS,CAAC,CAAe,EAAE,IAA4B;QACrD,OAAO,IAAI,kBAAkB,CAAC,MAAM,CAAC,CAAC,EAAE,WAAW,CAAC,CAAC,SAAS,CAAC,IAAI,CAAC,CAAC,CAAC;IACxE,CAAC;IACD,OAAO,CAAC,CAAe,EAAE,KAA4B;QACnD,OAAO,IAAI,kBAAkB,CAAC,MAAM,CAAC,CAAC,EAAE,SAAS,CAAC,CAAC,OAAO,CAAC,KAAK,CAAC,CAAC,CAAC;IACrE,CAAC;IAED,4EAA4E;IAC5E,yEAAyE;IACzE,sEAAsE;IACtE,uEAAuE;IACvE,uEAAuE;IACvE,wEAAwE;IACxE,uDAAuD;IACvD,GAAG,CAAC,CAAe,EAAE,CAAe;QAClC,IAAI,SAAS,IAAI,CAAC,IAAI,SAAS,IAAI,CAAC,EAAE,CAAC;YACrC,OAAQ,CAAkD,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;QACpE,CAAC;QACD,0CAA0C;QAC1C,IAAI,MAAM,IAAI,CAAC,IAAI,MAAM,IAAI,CAAC,EAAE,CAAC;YAC/B,OAAQ,CAAkD,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;QACpE,CAAC;QACD,OAAO,IAAI,kBAAkB,CAAC,MAAM,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,GAAG,CAAC,MAAM,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC,CAAC;IACxE,CAAC;IACD,GAAG,CAAC,CAAe,EAAE,CAAe;QAClC,IAAI,SAAS,IAAI,CAAC,IAAI,SAAS,IAAI,CAAC,EAAE,CAAC;YACrC,OAAQ,CAAkD,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;QACpE,CAAC;QACD,IAAI,MAAM,IAAI,CAAC,IAAI,MAAM,IAAI,CAAC,EAAE,CAAC;YAC/B,OAAQ,CAAkD,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;QACpE,CAAC;QACD,OAAO,IAAI,kBAAkB,CAAC,MAAM,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,GAAG,CAAC,MAAM,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC,CAAC;IACxE,CAAC;IACD,GAAG,CAAC,CAAe,EAAE,CAAe;QAClC,IAAI,SAAS,IAAI,CAAC,IAAI,SAAS,IAAI,CAAC,EAAE,CAAC;YACrC,OAAQ,CAAkD,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;QACpE,CAAC;QACD,IAAI,MAAM,IAAI,CAAC,IAAI,MAAM,IAAI,CAAC,EAAE,CAAC;YAC/B,OAAQ,CAAkD,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;QACpE,CAAC;QACD,OAAO,IAAI,kBAAkB,CAAC,MAAM,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,GAAG,CAAC,MAAM,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC,CAAC;IACxE,CAAC;IACD,KAAK,CAAC,CAAe,EAAE,CAAS;QAC9B,IAAI,SAAS,IAAI,CAAC,IAAI,MAAM,IAAI,CAAC,EAAE,CAAC;YAClC,OAAQ,CAAmD,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC;QACvE,CAAC;QACD,OAAO,IAAI,kBAAkB,CAAC,MAAM,CAAC,CAAC,EAAE,OAAO,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC;IAC7D,CAAC;IAED,QAAQ,CAAC,CAAS;QAChB,OAAO,IAAI,kBAAkB,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC;IACpD,CAAC;IACD,OAAO,CAAC,CAAe;QACrB,OAAO,MAAM,CAAC,CAAC,EAAE,SAAS,CAAC,CAAC,OAAO,EAAE,CAAC;IACxC,CAAC;IAED,8EAA8E;IAC9E,+CAA+C;IAC/C,8EAA8E;IAE9E,oFAAoF;IAC5E,cAAc,CAAC,CAAe;QACpC,OAAO,MAAM,CAAC,CAAC,EAAE,gBAAgB,CAAC,CAAC;IACrC,CAAC;IAED,2EAA2E;IACnE,gBAAgB,CAAC,CAAS;QAChC,OAAO,IAAI,kBAAkB,CAAC,CAAC,CAAC,CAAC;IACnC,CAAC;IAED,8EAA8E;IAC9E,4CAA4C;IAC5C,8EAA8E;IAE9E;;;;;;;;;;;OAWG;IACH,KAAK,CAAC,WAAW,CACf,EAAqC,EACrC,CAAe;QAEf,IAAI,QAAwB,CAAC;QAC7B,IAAI,CAAC;YACH,0DAA0D;YAC1D,QAAQ,GAAG,MAAM,MAAM,CAAC,gCAAgC,CAA8B,CAAC;QACzF,CAAC;QAAC,OAAO,KAAK,EAAE,CAAC;YAAC,MAAM,IAAI,qBAAqB,CAAC,cAAc,EAAE,aAAa,EAAE,EAAE,KAAK,EAAE,CAAC,CAAC;QAAC,CAAC;QAE9F,2EAA2E;QAC3E,sEAAsE;QACtE,6DAA6D;QAC7D,6DAA6D;QAC7D,yCAAyC;QACzC,MAAM,OAAO,GAAG,IAAI,CAAC,cAAc,CAAC,CAAC,CAAC,CAAC;QACvC,MAAM,EAAE,KAAK,EAAE,QAAQ,EAAE,GAAG,MAAM,QAAQ,CAAC,WAAW,CAAC,EAAE,EAAE,OAAO,CAAC,CAAC;QACpE,OAAO;YACL,KAAK,EAAE,IAAI,CAAC,gBAAgB,CAAC,KAAK,CAAC;YACnC,QAAQ,EAAE,IAAI,CAAC,gBAAgB,CAAC,QAAQ,CAAC;SAC1C,CAAC;IACJ,CAAC;IAED,8EAA8E;IAC9E,4CAA4C;IAC5C,8EAA8E;IAE9E;;;;;;;;;OASG;IACH,KAAK,CAAC,WAAW,CACf,EAAqC,EACrC,CAAe,EACf,SAAwB;QAExB,IAAI,QAAwB,CAAC;QAC7B,IAAI,CAAC;YACH,0DAA0D;YAC1D,QAAQ,GAAG,MAAM,MAAM,CAAC,gCAAgC,CAA8B,CAAC;QACzF,CAAC;QAAC,OAAO,KAAK,EAAE,CAAC;YAAC,MAAM,IAAI,qBAAqB,CAAC,cAAc,EAAE,aAAa,EAAE,EAAE,KAAK,EAAE,CAAC,CAAC;QAAC,CAAC;QAE9F,+EAA+E;QAC/E,sEAAsE;QACtE,+EAA+E;QAC/E,MAAM,OAAO,GAAG,IAAI,CAAC,cAAc,CAAC,CAAC,CAAC,CAAC;QACvC,MAAM,QAAQ,GAAG,SAAS,CAAC,CAAC,CAAC,IAAI,CAAC,cAAc,CAAC,SAAS,CAAC,CAAC,CAAC,CAAC,SAAS,CAAC;QACxE,MAAM,EAAE,KAAK,EAAE,QAAQ,EAAE,GAAG,MAAM,QAAQ,CAAC,WAAW,CAAC,EAAE,EAAE,OAAO,EAAE,QAAQ,CAAC,CAAC;QAC9E,OAAO;YACL,KAAK,EAAE,IAAI,CAAC,gBAAgB,CAAC,KAAK,CAAC;YACnC,QAAQ,EAAE,IAAI,CAAC,gBAAgB,CAAC,QAAQ,CAAC;SAC1C,CAAC;IACJ,CAAC;CACF","sourcesContent":["/**\n * MathTSEngine — a TensorEngine implementation backed by\n * @danielsimonjr/mathts-tensor's rank-N Tensor. It is the only engine.\n * The MathTS packages are required dependencies; see engine-registry.ts.\n *\n * Thin adapter: it translates the TensorEngine contract onto the MathTS\n * Tensor's methods. The engine-conformance suite is what checks that\n * translation.\n *\n * @module numerical/mathts-engine\n */\nimport { Tensor } from '@danielsimonjr/mathts-tensor';\nimport type { EngineTensor, TensorEngine, EinsumSpec, ForwardGradResult, ReverseGradResult } from './tensor-engine.js';\nimport type { NestedArray } from './types.js';\nimport { NumericalBackendError } from './errors.js';\nimport { EngineCapabilityError } from './tensor-engine.js';\n\n/** EngineTensor handle wrapping a MathTS Tensor. */\nclass MathTSEngineTensor implements EngineTensor {\n  constructor(readonly inner: Tensor) {}\n  get shape(): ReadonlyArray<number> { return this.inner.shape; }\n}\n\nfunction unwrap(t: EngineTensor, op: string): Tensor {\n  if (!(t instanceof MathTSEngineTensor)) {\n    throw new NumericalBackendError(`MathTSEngine.${op}: operand is not a MathTSEngineTensor`);\n  }\n  return t.inner;\n}\n\n/**\n * Call-site shape for `@danielsimonjr/mathts-autograd`. The package is a\n * required dependency and is loaded at the call, so this interface covers\n * the dynamic-import result only.\n * @internal\n */\ninterface MathTSAutograd {\n  forwardGrad(\n    fn: (x: EngineTensor) => EngineTensor,\n    x: Tensor,\n  ): Promise<{ value: Tensor; jacobian: Tensor }>;\n  reverseGrad(\n    fn: (x: EngineTensor) => EngineTensor,\n    x: Tensor,\n    cotangent?: Tensor,\n  ): Promise<{ value: Tensor; gradient: Tensor }>;\n}\n\n/**\n * `TensorEngine` backed by `@danielsimonjr/mathts-tensor`'s rank-N Tensor.\n * The MathTS packages are required dependencies.\n *\n * @public — reachable only via the\n * `universal-physics-tensor/numerical/mathts-engine` exports subpath.\n * Intentionally NOT re-exported from the root barrel.\n */\nexport class MathTSEngine implements TensorEngine {\n  readonly name = 'MathTSEngine';\n\n  fromNested(data: NestedArray, shape: ReadonlyArray<number>): EngineTensor {\n    // UPT's NestedArray and MathTS's are the same shape; the type parameters differ.\n    return new MathTSEngineTensor(Tensor.fromNested(data as never, shape));\n  }\n  toNested(t: EngineTensor): NestedArray {\n    return unwrap(t, 'toNested').toNested();\n  }\n\n  einsum(spec: EinsumSpec, ...operands: EngineTensor[]): EngineTensor {\n    const inner = operands.map((o, i) => unwrap(o, `einsum (operand ${i})`));\n    return new MathTSEngineTensor(Tensor.einsum(spec, ...inner));\n  }\n  matMul(a: EngineTensor, b: EngineTensor): EngineTensor {\n    return new MathTSEngineTensor(unwrap(a, 'matMul').matMul(unwrap(b, 'matMul')));\n  }\n  transpose(t: EngineTensor, perm?: ReadonlyArray<number>): EngineTensor {\n    return new MathTSEngineTensor(unwrap(t, 'transpose').transpose(perm));\n  }\n  reshape(t: EngineTensor, shape: ReadonlyArray<number>): EngineTensor {\n    return new MathTSEngineTensor(unwrap(t, 'reshape').reshape(shape));\n  }\n\n  // AD dispatch. `@danielsimonjr/mathts-autograd` wraps the traced input as a\n  // DualTensor (forward mode, carries `tangent`) or a TapedTensor (reverse\n  // mode, carries `tape`) and calls these methods with it. The autograd\n  // package is loaded at the AD call, not at module load, so its classes\n  // cannot be named here; the structural casts below state the invariant\n  // the package keeps: a wrapped tensor implements add/sub/mul/scale over\n  // another wrapped tensor and returns a wrapped tensor.\n  add(a: EngineTensor, b: EngineTensor): EngineTensor {\n    if ('tangent' in a && 'tangent' in b) {\n      return (a as unknown as { add(o: unknown): EngineTensor }).add(b);\n    }\n    // AD dispatch: TapedTensor (reverse-mode)\n    if ('tape' in a && 'tape' in b) {\n      return (a as unknown as { add(o: unknown): EngineTensor }).add(b);\n    }\n    return new MathTSEngineTensor(unwrap(a, 'add').add(unwrap(b, 'add')));\n  }\n  sub(a: EngineTensor, b: EngineTensor): EngineTensor {\n    if ('tangent' in a && 'tangent' in b) {\n      return (a as unknown as { sub(o: unknown): EngineTensor }).sub(b);\n    }\n    if ('tape' in a && 'tape' in b) {\n      return (a as unknown as { sub(o: unknown): EngineTensor }).sub(b);\n    }\n    return new MathTSEngineTensor(unwrap(a, 'sub').sub(unwrap(b, 'sub')));\n  }\n  mul(a: EngineTensor, b: EngineTensor): EngineTensor {\n    if ('tangent' in a && 'tangent' in b) {\n      return (a as unknown as { mul(o: unknown): EngineTensor }).mul(b);\n    }\n    if ('tape' in a && 'tape' in b) {\n      return (a as unknown as { mul(o: unknown): EngineTensor }).mul(b);\n    }\n    return new MathTSEngineTensor(unwrap(a, 'mul').mul(unwrap(b, 'mul')));\n  }\n  scale(t: EngineTensor, k: number): EngineTensor {\n    if ('tangent' in t || 'tape' in t) {\n      return (t as unknown as { scale(k: number): EngineTensor }).scale(k);\n    }\n    return new MathTSEngineTensor(unwrap(t, 'scale').scale(k));\n  }\n\n  identity(n: number): EngineTensor {\n    return new MathTSEngineTensor(Tensor.identity(n));\n  }\n  normInf(t: EngineTensor): number {\n    return unwrap(t, 'normInf').normInf();\n  }\n\n  // ---------------------------------------------------------------------------\n  // Internal helpers for AD boundary conversions\n  // ---------------------------------------------------------------------------\n\n  /** Unwrap an EngineTensor to its underlying MathTS Tensor (for the AD boundary). */\n  private toMathTSTensor(t: EngineTensor): Tensor {\n    return unwrap(t, 'toMathTSTensor');\n  }\n\n  /** Wrap a MathTS Tensor (returned by autograd) back as an EngineTensor. */\n  private fromMathTSTensor(t: Tensor): EngineTensor {\n    return new MathTSEngineTensor(t);\n  }\n\n  // ---------------------------------------------------------------------------\n  // Forward-mode AD (Jacobian-vector product)\n  // ---------------------------------------------------------------------------\n\n  /**\n   * Forward-mode automatic differentiation via a lazy-imported\n   * `@danielsimonjr/mathts-autograd`. Throws `EngineCapabilityError` if that\n   * import fails.\n   *\n   * S1 fix: `fn` is passed UNCHANGED to `autograd.forwardGrad`. The autograd\n   * package wraps `x` as a DualTensor internally; MathTSEngine's arithmetic\n   * methods (mul/add/sub/scale) must dispatch DualTensor inputs through to\n   * mathts-autograd's dual arithmetic — there must be no fn-wrapping at the\n   * boundary, which would strip DualTensor instrumentation and silence the\n   * AD trace.\n   */\n  async forwardGrad(\n    fn: (x: EngineTensor) => EngineTensor,\n    x: EngineTensor,\n  ): Promise<ForwardGradResult> {\n    let autograd: MathTSAutograd;\n    try {\n      // Narrow the dynamic import to the local call-site shape.\n      autograd = await import('@danielsimonjr/mathts-autograd') as unknown as MathTSAutograd;\n    } catch (cause) { throw new EngineCapabilityError('MathTSEngine', 'forwardGrad', { cause }); }\n\n    // S1 fix: pass fn UNCHANGED. autograd.forwardGrad wraps x as a DualTensor;\n    // MathTSEngine's arithmetic methods (mul/add/sub/scale) MUST dispatch\n    // DualTensor inputs to mathts-autograd's dual arithmetic via\n    // `'tangent' in arg`. Wrapping fn at the boundary strips the\n    // instrumentation — the v0 sketch's bug.\n    const xMathts = this.toMathTSTensor(x);\n    const { value, jacobian } = await autograd.forwardGrad(fn, xMathts);\n    return {\n      value: this.fromMathTSTensor(value),\n      jacobian: this.fromMathTSTensor(jacobian),\n    };\n  }\n\n  // ---------------------------------------------------------------------------\n  // Reverse-mode AD (vector-Jacobian product)\n  // ---------------------------------------------------------------------------\n\n  /**\n   * Reverse-mode automatic differentiation via a lazy-imported\n   * `@danielsimonjr/mathts-autograd`. Throws `EngineCapabilityError` if that\n   * import fails.\n   *\n   * S1 fix: `fn` is passed UNCHANGED (see forwardGrad note above). The\n   * autograd package wraps `x` as a TapedTensor; MathTSEngine's op methods\n   * must dispatch TapedTensor inputs to mathts-autograd's tape arithmetic via\n   * `'tape' in arg` branching.\n   */\n  async reverseGrad(\n    fn: (x: EngineTensor) => EngineTensor,\n    x: EngineTensor,\n    cotangent?: EngineTensor,\n  ): Promise<ReverseGradResult> {\n    let autograd: MathTSAutograd;\n    try {\n      // Narrow the dynamic import to the local call-site shape.\n      autograd = await import('@danielsimonjr/mathts-autograd') as unknown as MathTSAutograd;\n    } catch (cause) { throw new EngineCapabilityError('MathTSEngine', 'reverseGrad', { cause }); }\n\n    // S1 fix: pass fn UNCHANGED (see forwardGrad note above). autograd.reverseGrad\n    // wraps x as a TapedTensor; MathTSEngine's mul/add/sub/scale dispatch\n    // TapedTensor inputs to mathts-autograd's tape arithmetic via `'tape' in arg`.\n    const xMathts = this.toMathTSTensor(x);\n    const ctMathts = cotangent ? this.toMathTSTensor(cotangent) : undefined;\n    const { value, gradient } = await autograd.reverseGrad(fn, xMathts, ctMathts);\n    return {\n      value: this.fromMathTSTensor(value),\n      gradient: this.fromMathTSTensor(gradient),\n    };\n  }\n}\n"]}
+{"version":3,"file":"mathts-engine.js","sourceRoot":"","sources":["../../src/numerical/mathts-engine.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;GAUG;AACH,OAAO,EAAE,MAAM,EAAE,MAAM,8BAA8B,CAAC;AAGtD,OAAO,EAAE,qBAAqB,EAAE,MAAM,aAAa,CAAC;AACpD,OAAO,EAAE,qBAAqB,EAAE,cAAc,EAAE,UAAU,EAAE,MAAM,oBAAoB,CAAC;AAEvF,oDAAoD;AACpD,MAAM,kBAAkB;IACD,KAAK;IAA1B,YAAqB,KAAa;qBAAb,KAAK;IAAW,CAAC;IACtC,IAAI,KAAK,KAA4B,OAAO,IAAI,CAAC,KAAK,CAAC,KAAK,CAAC,CAAC,CAAC;CAChE;AAED,SAAS,MAAM,CAAC,CAAe,EAAE,EAAU;IACzC,IAAI,CAAC,CAAC,CAAC,YAAY,kBAAkB,CAAC,EAAE,CAAC;QACvC,MAAM,IAAI,qBAAqB,CAAC,gBAAgB,EAAE,uCAAuC,CAAC,CAAC;IAC7F,CAAC;IACD,OAAO,CAAC,CAAC,KAAK,CAAC;AACjB,CAAC;AAoBD;;;;;;;GAOG;AACH,MAAM,OAAO,YAAY;IACd,IAAI,GAAG,cAAc,CAAC;IAE/B,UAAU,CAAC,IAAiB,EAAE,KAA4B;QACxD,iFAAiF;QACjF,OAAO,IAAI,kBAAkB,CAAC,MAAM,CAAC,UAAU,CAAC,IAAa,EAAE,KAAK,CAAC,CAAC,CAAC;IACzE,CAAC;IACD,QAAQ,CAAC,CAAe;QACtB,OAAO,MAAM,CAAC,CAAC,EAAE,UAAU,CAAC,CAAC,QAAQ,EAAE,CAAC;IAC1C,CAAC;IAED,MAAM,CAAC,IAAgB,EAAE,GAAG,QAAwB;QAClD,MAAM,KAAK,GAAG,QAAQ,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,MAAM,CAAC,CAAC,EAAE,mBAAmB,CAAC,GAAG,CAAC,CAAC,CAAC;QACzE,OAAO,IAAI,kBAAkB,CAAC,MAAM,CAAC,MAAM,CAAC,IAAI,EAAE,GAAG,KAAK,CAAC,CAAC,CAAC;IAC/D,CAAC;IACD,MAAM,CAAC,CAAe,EAAE,CAAe;QACrC,OAAO,IAAI,kBAAkB,CAAC,MAAM,CAAC,CAAC,EAAE,QAAQ,CAAC,CAAC,MAAM,CAAC,MAAM,CAAC,CAAC,EAAE,QAAQ,CAAC,CAAC,CAAC,CAAC;IACjF,CAAC;IACD,SAAS,CAAC,CAAe,EAAE,IAA4B;QACrD,OAAO,IAAI,kBAAkB,CAAC,MAAM,CAAC,CAAC,EAAE,WAAW,CAAC,CAAC,SAAS,CAAC,IAAI,CAAC,CAAC,CAAC;IACxE,CAAC;IACD,OAAO,CAAC,CAAe,EAAE,KAA4B;QACnD,OAAO,IAAI,kBAAkB,CAAC,MAAM,CAAC,CAAC,EAAE,SAAS,CAAC,CAAC,OAAO,CAAC,KAAK,CAAC,CAAC,CAAC;IACrE,CAAC;IAED,4EAA4E;IAC5E,yEAAyE;IACzE,0EAA0E;IAC1E,4EAA4E;IAC5E,0EAA0E;IAC1E,wEAAwE;IACxE,wDAAwD;IACxD,GAAG,CAAC,CAAe,EAAE,CAAe;QAClC,IAAI,cAAc,CAAC,CAAC,CAAC,IAAI,cAAc,CAAC,CAAC,CAAC,IAAI,UAAU,CAAC,CAAC,CAAC,KAAK,UAAU,CAAC,CAAC,CAAC;YAAE,OAAO,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;QAC/F,OAAO,IAAI,kBAAkB,CAAC,MAAM,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,GAAG,CAAC,MAAM,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC,CAAC;IACxE,CAAC;IACD,GAAG,CAAC,CAAe,EAAE,CAAe;QAClC,IAAI,cAAc,CAAC,CAAC,CAAC,IAAI,cAAc,CAAC,CAAC,CAAC,IAAI,UAAU,CAAC,CAAC,CAAC,KAAK,UAAU,CAAC,CAAC,CAAC;YAAE,OAAO,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;QAC/F,OAAO,IAAI,kBAAkB,CAAC,MAAM,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,GAAG,CAAC,MAAM,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC,CAAC;IACxE,CAAC;IACD,GAAG,CAAC,CAAe,EAAE,CAAe;QAClC,IAAI,cAAc,CAAC,CAAC,CAAC,IAAI,cAAc,CAAC,CAAC,CAAC,IAAI,UAAU,CAAC,CAAC,CAAC,KAAK,UAAU,CAAC,CAAC,CAAC;YAAE,OAAO,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;QAC/F,OAAO,IAAI,kBAAkB,CAAC,MAAM,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,GAAG,CAAC,MAAM,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC,CAAC;IACxE,CAAC;IACD,KAAK,CAAC,CAAe,EAAE,CAAS;QAC9B,IAAI,cAAc,CAAC,CAAC,CAAC;YAAE,OAAO,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC;QACzC,OAAO,IAAI,kBAAkB,CAAC,MAAM,CAAC,CAAC,EAAE,OAAO,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC;IAC7D,CAAC;IAED,QAAQ,CAAC,CAAS;QAChB,OAAO,IAAI,kBAAkB,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC;IACpD,CAAC;IACD,OAAO,CAAC,CAAe;QACrB,OAAO,MAAM,CAAC,CAAC,EAAE,SAAS,CAAC,CAAC,OAAO,EAAE,CAAC;IACxC,CAAC;IAED,8EAA8E;IAC9E,+CAA+C;IAC/C,8EAA8E;IAE9E,oFAAoF;IAC5E,cAAc,CAAC,CAAe;QACpC,OAAO,MAAM,CAAC,CAAC,EAAE,gBAAgB,CAAC,CAAC;IACrC,CAAC;IAED,2EAA2E;IACnE,gBAAgB,CAAC,CAAS;QAChC,OAAO,IAAI,kBAAkB,CAAC,CAAC,CAAC,CAAC;IACnC,CAAC;IAED,8EAA8E;IAC9E,4CAA4C;IAC5C,8EAA8E;IAE9E;;;;;;;;;;;OAWG;IACH,KAAK,CAAC,WAAW,CACf,EAAqC,EACrC,CAAe;QAEf,IAAI,QAAwB,CAAC;QAC7B,IAAI,CAAC;YACH,0DAA0D;YAC1D,QAAQ,GAAG,MAAM,MAAM,CAAC,gCAAgC,CAA8B,CAAC;QACzF,CAAC;QAAC,OAAO,KAAK,EAAE,CAAC;YAAC,MAAM,IAAI,qBAAqB,CAAC,cAAc,EAAE,aAAa,EAAE,EAAE,KAAK,EAAE,CAAC,CAAC;QAAC,CAAC;QAE9F,2EAA2E;QAC3E,sEAAsE;QACtE,6DAA6D;QAC7D,6DAA6D;QAC7D,yCAAyC;QACzC,MAAM,OAAO,GAAG,IAAI,CAAC,cAAc,CAAC,CAAC,CAAC,CAAC;QACvC,MAAM,EAAE,KAAK,EAAE,QAAQ,EAAE,GAAG,MAAM,QAAQ,CAAC,WAAW,CAAC,EAAE,EAAE,OAAO,CAAC,CAAC;QACpE,OAAO;YACL,KAAK,EAAE,IAAI,CAAC,gBAAgB,CAAC,KAAK,CAAC;YACnC,QAAQ,EAAE,IAAI,CAAC,gBAAgB,CAAC,QAAQ,CAAC;SAC1C,CAAC;IACJ,CAAC;IAED,8EAA8E;IAC9E,4CAA4C;IAC5C,8EAA8E;IAE9E;;;;;;;;;OASG;IACH,KAAK,CAAC,WAAW,CACf,EAAqC,EACrC,CAAe,EACf,SAAwB;QAExB,IAAI,QAAwB,CAAC;QAC7B,IAAI,CAAC;YACH,0DAA0D;YAC1D,QAAQ,GAAG,MAAM,MAAM,CAAC,gCAAgC,CAA8B,CAAC;QACzF,CAAC;QAAC,OAAO,KAAK,EAAE,CAAC;YAAC,MAAM,IAAI,qBAAqB,CAAC,cAAc,EAAE,aAAa,EAAE,EAAE,KAAK,EAAE,CAAC,CAAC;QAAC,CAAC;QAE9F,+EAA+E;QAC/E,sEAAsE;QACtE,+EAA+E;QAC/E,MAAM,OAAO,GAAG,IAAI,CAAC,cAAc,CAAC,CAAC,CAAC,CAAC;QACvC,MAAM,QAAQ,GAAG,SAAS,CAAC,CAAC,CAAC,IAAI,CAAC,cAAc,CAAC,SAAS,CAAC,CAAC,CAAC,CAAC,SAAS,CAAC;QACxE,MAAM,EAAE,KAAK,EAAE,QAAQ,EAAE,GAAG,MAAM,QAAQ,CAAC,WAAW,CAAC,EAAE,EAAE,OAAO,EAAE,QAAQ,CAAC,CAAC;QAC9E,OAAO;YACL,KAAK,EAAE,IAAI,CAAC,gBAAgB,CAAC,KAAK,CAAC;YACnC,QAAQ,EAAE,IAAI,CAAC,gBAAgB,CAAC,QAAQ,CAAC;SAC1C,CAAC;IACJ,CAAC;CACF","sourcesContent":["/**\n * MathTSEngine — a TensorEngine implementation backed by\n * @danielsimonjr/mathts-tensor's rank-N Tensor. It is the only engine.\n * The MathTS packages are required dependencies; see engine-registry.ts.\n *\n * Thin adapter: it translates the TensorEngine contract onto the MathTS\n * Tensor's methods. The engine-conformance suite is what checks that\n * translation.\n *\n * @module numerical/mathts-engine\n */\nimport { Tensor } from '@danielsimonjr/mathts-tensor';\nimport type { EngineTensor, TensorEngine, EinsumSpec, ForwardGradResult, ReverseGradResult } from './tensor-engine.js';\nimport type { NestedArray } from './types.js';\nimport { NumericalBackendError } from './errors.js';\nimport { EngineCapabilityError, isTracedTensor, tracedMode } from './tensor-engine.js';\n\n/** EngineTensor handle wrapping a MathTS Tensor. */\nclass MathTSEngineTensor implements EngineTensor {\n  constructor(readonly inner: Tensor) {}\n  get shape(): ReadonlyArray<number> { return this.inner.shape; }\n}\n\nfunction unwrap(t: EngineTensor, op: string): Tensor {\n  if (!(t instanceof MathTSEngineTensor)) {\n    throw new NumericalBackendError(`MathTSEngine.${op}: operand is not a MathTSEngineTensor`);\n  }\n  return t.inner;\n}\n\n/**\n * Call-site shape for `@danielsimonjr/mathts-autograd`. The package is a\n * required dependency and is loaded at the call, so this interface covers\n * the dynamic-import result only.\n * @internal\n */\ninterface MathTSAutograd {\n  forwardGrad(\n    fn: (x: EngineTensor) => EngineTensor,\n    x: Tensor,\n  ): Promise<{ value: Tensor; jacobian: Tensor }>;\n  reverseGrad(\n    fn: (x: EngineTensor) => EngineTensor,\n    x: Tensor,\n    cotangent?: Tensor,\n  ): Promise<{ value: Tensor; gradient: Tensor }>;\n}\n\n/**\n * `TensorEngine` backed by `@danielsimonjr/mathts-tensor`'s rank-N Tensor.\n * The MathTS packages are required dependencies.\n *\n * @public — reachable only via the\n * `universal-physics-tensor/numerical/mathts-engine` exports subpath.\n * Intentionally NOT re-exported from the root barrel.\n */\nexport class MathTSEngine implements TensorEngine {\n  readonly name = 'MathTSEngine';\n\n  fromNested(data: NestedArray, shape: ReadonlyArray<number>): EngineTensor {\n    // UPT's NestedArray and MathTS's are the same shape; the type parameters differ.\n    return new MathTSEngineTensor(Tensor.fromNested(data as never, shape));\n  }\n  toNested(t: EngineTensor): NestedArray {\n    return unwrap(t, 'toNested').toNested();\n  }\n\n  einsum(spec: EinsumSpec, ...operands: EngineTensor[]): EngineTensor {\n    const inner = operands.map((o, i) => unwrap(o, `einsum (operand ${i})`));\n    return new MathTSEngineTensor(Tensor.einsum(spec, ...inner));\n  }\n  matMul(a: EngineTensor, b: EngineTensor): EngineTensor {\n    return new MathTSEngineTensor(unwrap(a, 'matMul').matMul(unwrap(b, 'matMul')));\n  }\n  transpose(t: EngineTensor, perm?: ReadonlyArray<number>): EngineTensor {\n    return new MathTSEngineTensor(unwrap(t, 'transpose').transpose(perm));\n  }\n  reshape(t: EngineTensor, shape: ReadonlyArray<number>): EngineTensor {\n    return new MathTSEngineTensor(unwrap(t, 'reshape').reshape(shape));\n  }\n\n  // AD dispatch. `@danielsimonjr/mathts-autograd` wraps the traced input as a\n  // DualTensor (forward mode, carries `tangent`) or a TapedTensor (reverse\n  // mode, carries `tape`) and calls these methods with it. The operations a\n  // wrapped tensor implements are declared on `TracedEngineTensor`; the guard\n  // `isTracedTensor` is the dispatch. Two operands trace together only when\n  // both are wrapped in the same mode; a mixed or half-wrapped pair falls\n  // through to the plain path, whose `unwrap` refuses it.\n  add(a: EngineTensor, b: EngineTensor): EngineTensor {\n    if (isTracedTensor(a) && isTracedTensor(b) && tracedMode(a) === tracedMode(b)) return a.add(b);\n    return new MathTSEngineTensor(unwrap(a, 'add').add(unwrap(b, 'add')));\n  }\n  sub(a: EngineTensor, b: EngineTensor): EngineTensor {\n    if (isTracedTensor(a) && isTracedTensor(b) && tracedMode(a) === tracedMode(b)) return a.sub(b);\n    return new MathTSEngineTensor(unwrap(a, 'sub').sub(unwrap(b, 'sub')));\n  }\n  mul(a: EngineTensor, b: EngineTensor): EngineTensor {\n    if (isTracedTensor(a) && isTracedTensor(b) && tracedMode(a) === tracedMode(b)) return a.mul(b);\n    return new MathTSEngineTensor(unwrap(a, 'mul').mul(unwrap(b, 'mul')));\n  }\n  scale(t: EngineTensor, k: number): EngineTensor {\n    if (isTracedTensor(t)) return t.scale(k);\n    return new MathTSEngineTensor(unwrap(t, 'scale').scale(k));\n  }\n\n  identity(n: number): EngineTensor {\n    return new MathTSEngineTensor(Tensor.identity(n));\n  }\n  normInf(t: EngineTensor): number {\n    return unwrap(t, 'normInf').normInf();\n  }\n\n  // ---------------------------------------------------------------------------\n  // Internal helpers for AD boundary conversions\n  // ---------------------------------------------------------------------------\n\n  /** Unwrap an EngineTensor to its underlying MathTS Tensor (for the AD boundary). */\n  private toMathTSTensor(t: EngineTensor): Tensor {\n    return unwrap(t, 'toMathTSTensor');\n  }\n\n  /** Wrap a MathTS Tensor (returned by autograd) back as an EngineTensor. */\n  private fromMathTSTensor(t: Tensor): EngineTensor {\n    return new MathTSEngineTensor(t);\n  }\n\n  // ---------------------------------------------------------------------------\n  // Forward-mode AD (Jacobian-vector product)\n  // ---------------------------------------------------------------------------\n\n  /**\n   * Forward-mode automatic differentiation via a lazy-imported\n   * `@danielsimonjr/mathts-autograd`. Throws `EngineCapabilityError` if that\n   * import fails.\n   *\n   * S1 fix: `fn` is passed UNCHANGED to `autograd.forwardGrad`. The autograd\n   * package wraps `x` as a DualTensor internally; MathTSEngine's arithmetic\n   * methods (mul/add/sub/scale) must dispatch DualTensor inputs through to\n   * mathts-autograd's dual arithmetic — there must be no fn-wrapping at the\n   * boundary, which would strip DualTensor instrumentation and silence the\n   * AD trace.\n   */\n  async forwardGrad(\n    fn: (x: EngineTensor) => EngineTensor,\n    x: EngineTensor,\n  ): Promise<ForwardGradResult> {\n    let autograd: MathTSAutograd;\n    try {\n      // Narrow the dynamic import to the local call-site shape.\n      autograd = await import('@danielsimonjr/mathts-autograd') as unknown as MathTSAutograd;\n    } catch (cause) { throw new EngineCapabilityError('MathTSEngine', 'forwardGrad', { cause }); }\n\n    // S1 fix: pass fn UNCHANGED. autograd.forwardGrad wraps x as a DualTensor;\n    // MathTSEngine's arithmetic methods (mul/add/sub/scale) MUST dispatch\n    // DualTensor inputs to mathts-autograd's dual arithmetic via\n    // `'tangent' in arg`. Wrapping fn at the boundary strips the\n    // instrumentation — the v0 sketch's bug.\n    const xMathts = this.toMathTSTensor(x);\n    const { value, jacobian } = await autograd.forwardGrad(fn, xMathts);\n    return {\n      value: this.fromMathTSTensor(value),\n      jacobian: this.fromMathTSTensor(jacobian),\n    };\n  }\n\n  // ---------------------------------------------------------------------------\n  // Reverse-mode AD (vector-Jacobian product)\n  // ---------------------------------------------------------------------------\n\n  /**\n   * Reverse-mode automatic differentiation via a lazy-imported\n   * `@danielsimonjr/mathts-autograd`. Throws `EngineCapabilityError` if that\n   * import fails.\n   *\n   * S1 fix: `fn` is passed UNCHANGED (see forwardGrad note above). The\n   * autograd package wraps `x` as a TapedTensor; MathTSEngine's op methods\n   * must dispatch TapedTensor inputs to mathts-autograd's tape arithmetic via\n   * `'tape' in arg` branching.\n   */\n  async reverseGrad(\n    fn: (x: EngineTensor) => EngineTensor,\n    x: EngineTensor,\n    cotangent?: EngineTensor,\n  ): Promise<ReverseGradResult> {\n    let autograd: MathTSAutograd;\n    try {\n      // Narrow the dynamic import to the local call-site shape.\n      autograd = await import('@danielsimonjr/mathts-autograd') as unknown as MathTSAutograd;\n    } catch (cause) { throw new EngineCapabilityError('MathTSEngine', 'reverseGrad', { cause }); }\n\n    // S1 fix: pass fn UNCHANGED (see forwardGrad note above). autograd.reverseGrad\n    // wraps x as a TapedTensor; MathTSEngine's mul/add/sub/scale dispatch\n    // TapedTensor inputs to mathts-autograd's tape arithmetic via `'tape' in arg`.\n    const xMathts = this.toMathTSTensor(x);\n    const ctMathts = cotangent ? this.toMathTSTensor(cotangent) : undefined;\n    const { value, gradient } = await autograd.reverseGrad(fn, xMathts, ctMathts);\n    return {\n      value: this.fromMathTSTensor(value),\n      gradient: this.fromMathTSTensor(gradient),\n    };\n  }\n}\n"]}

package/dist/numerical/metric-inverse.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"metric-inverse.d.ts","sourceRoot":"","sources":["../../src/numerical/metric-inverse.ts"],"names":[],"mappings":"AAAA;;;;;;;;;GASG;AACH,OAAO,KAAK,EAAE,QAAQ,EAAE,SAAS,EAAE,MAAM,6BAA6B,CAAC;AACvE,OAAO,KAAK,EAAE,gBAAgB,EAAE,MAAM,qCAAqC,CAAC;AAE5E,OAAO,KAAK,EAAE,YAAY,EAAE,MAAM,oBAAoB,CAAC;AACvD,OAAO,KAAK,EAAE,eAAe,EAAE,MAAM,YAAY,CAAC;AAMlD;cACc;AACd,wBAAsB,qBAAqB,CACzC,QAAQ,EAAE,gBAAgB,EAC1B,CAAC,EAAE,gBAAgB,EACnB,MAAM,EAAE,eAAe,EACvB,SAAS,GAAE,MAA0B,EACrC,OAAO,CAAC,EAAE;IAAE,MAAM,CAAC,EAAE,YAAY,CAAA;CAAE,GAClC,OAAO,CAAC;IAAE,YAAY,EAAE,MAAM,CAAC;IAAC,OAAO,CAAC,EAAE,SAAS,CAAA;CAAE,CAAC,CAkCxD;AAED;;;2EAG2E;AAC3E,wBAAgB,iBAAiB,CAC/B,IAAI,EAAE,QAAQ,GACb;IAAE,MAAM,EAAE,gBAAgB,CAAC;IAAC,MAAM,EAAE,gBAAgB,CAAA;CAAE,GAAG,IAAI,CA6C/D"}
+{"version":3,"file":"metric-inverse.d.ts","sourceRoot":"","sources":["../../src/numerical/metric-inverse.ts"],"names":[],"mappings":"AAAA;;;;;;;;;GASG;AACH,OAAO,KAAK,EAAE,QAAQ,EAAE,SAAS,EAAE,MAAM,6BAA6B,CAAC;AACvE,OAAO,KAAK,EAAE,gBAAgB,EAAE,MAAM,qCAAqC,CAAC;AAE5E,OAAO,KAAK,EAAE,YAAY,EAAE,MAAM,oBAAoB,CAAC;AACvD,OAAO,KAAK,EAAE,eAAe,EAAE,MAAM,YAAY,CAAC;AAMlD;cACc;AACd,wBAAsB,qBAAqB,CACzC,QAAQ,EAAE,gBAAgB,EAC1B,CAAC,EAAE,gBAAgB,EACnB,MAAM,EAAE,eAAe,EACvB,SAAS,GAAE,MAA0B,EACrC,OAAO,CAAC,EAAE;IAAE,MAAM,CAAC,EAAE,YAAY,CAAA;CAAE,GAClC,OAAO,CAAC;IAAE,YAAY,EAAE,MAAM,CAAC;IAAC,OAAO,CAAC,EAAE,SAAS,CAAA;CAAE,CAAC,CAkCxD;AAED;;;2EAG2E;AAC3E,wBAAgB,iBAAiB,CAC/B,IAAI,EAAE,QAAQ,GACb;IAAE,MAAM,EAAE,gBAAgB,CAAC;IAAC,MAAM,EAAE,gBAAgB,CAAA;CAAE,GAAG,IAAI,CA0F/D"}

package/dist/numerical/metric-inverse.js

@@ -39,54 +39,91 @@ export async function evaluateMetricInverse(gInverse, g, inputs, tolerance = DEF
  *  @internal — cross-module use only; not part of the consumer surface. */
 export function scanForMetricPair(node) {
     const metrics = [];
+    // An exhaustive switch: a node kind added to `ExprNode` without an arm here
+    // fails to compile at the `never` arm, instead of becoming a silent leaf (the
+    // way `kretschmann-scalar` and `einstein-equation` were, so the warning never
+    // fired on them). A kind with no metric inside is listed as a leaf on purpose.
     const walk = (n) => {
-        if (n.kind === 'metric-tensor')
-            metrics.push(n);
-        else if (n.kind === 'op')
-            n.args.forEach(walk);
-        else if (n.kind === 'tensor-product')
-            n.args.forEach(walk);
-        else if (n.kind === 'integral') {
-            walk(n.over);
-            walk(n.integrand);
+        switch (n.kind) {
+            case 'metric-tensor':
+                metrics.push(n);
+                break;
+            case 'op':
+            case 'tensor-product':
+                n.args.forEach(walk);
+                break;
+            case 'integral':
+                walk(n.over);
+                walk(n.integrand);
+                if (n.lower !== undefined)
+                    walk(n.lower);
+                if (n.upper !== undefined)
+                    walk(n.upper);
+                break;
+            case 'derivative':
+                walk(n.of);
+                walk(n.wrt);
+                break;
+            case 'variational-derivative':
+                walk(n.functional);
+                walk(n.field);
+                walk(n.over);
+                break;
+            case 'dirac-delta':
+            case 'transcendental':
+            case 'abs':
+                walk(n.arg);
+                break;
+            case 'tensor-partial-derivative':
+                walk(n.of);
+                walk(n.wrt);
+                break;
+            case 'covariant-derivative':
+                walk(n.of);
+                walk(n.gLower);
+                walk(n.gInverse);
+                break;
+            case 'riemann-tensor':
+                walk(n.gLower);
+                walk(n.gInverse);
+                break;
+            case 'ricci-tensor':
+            case 'bianchi-residual':
+                walk(n.riemann); // the inner Riemann holds the pair
+                break;
+            case 'einstein-tensor':
+                walk(n.riemann);
+                walk(n.gLower);
+                walk(n.gInverse);
+                break;
+            case 'weyl-tensor':
+                walk(n.metric);
+                break;
+            case 'kretschmann-scalar':
+                walk(n.riemann);
+                walk(n.metric);
+                break;
+            case 'killing-vector':
+                walk(n.metric);
+                break;
+            case 'conserved-charge':
+                walk(n.killing);
+                break;
+            case 'einstein-equation':
+                walk(n.lhs);
+                walk(n.metric);
+                break;
+            case 'symbol':
+            case 'tensor-symbol':
+            case 'kronecker-delta':
+            case 'stress-energy':
+            case 'cosmological-constant':
+                break; // leaves: no metric inside
+            default: {
+                const unreachable = n;
+                throw new NumericalBackendError(`scanForMetricPair: unhandled node kind ${String(unreachable.kind)}`);
+            }
         }
-        else if (n.kind === 'derivative') {
-            walk(n.of);
-            walk(n.wrt);
-        }
-        else if (n.kind === 'tensor-partial-derivative') {
-            walk(n.of);
-            walk(n.wrt);
-        }
-        // v0.5.1 PC-3: walk the v0.4.0+v0.5.0 curvature node kinds so that
-        // InverseMetricInconsistencyWarning fires on the GR pipeline. Prior to
-        // this, scan stopped at riemann/ricci/einstein/bianchi/covariant-derivative
-        // and the warning was silently never emitted on the whole curvature stack
-        // — "a validation mechanism that doesn't run is worse than no validation
-        // at all" (Adam, v0.5.1 audit).
-        //
-        // FUTURE: when a new curvature kind is added, this scan needs an arm.
-        else if (n.kind === 'covariant-derivative') {
-            walk(n.of);
-            walk(n.gLower);
-            walk(n.gInverse);
-        }
-        else if (n.kind === 'riemann-tensor') {
-            walk(n.gLower);
-            walk(n.gInverse);
-        }
-        else if (n.kind === 'ricci-tensor') {
-            walk(n.riemann); // inner Riemann walks its own gLower/gInverse
-        }
-        else if (n.kind === 'einstein-tensor') {
-            walk(n.riemann);
-            walk(n.gLower);
-            walk(n.gInverse);
-        }
-        else if (n.kind === 'bianchi-residual') {
-            walk(n.riemann); // recurse via inner Riemann
-        }
-        // tensor-symbol / kronecker-delta / symbol: leaves, no metric inside
     };
     walk(node);
     const lower = metrics.find((m) => m.indices.every((i) => i.variance === 'lower'));

package/dist/numerical/metric-inverse.js.map

@@ -1 +1 @@
-{"version":3,"file":"metric-inverse.js","sourceRoot":"","sources":["../../src/numerical/metric-inverse.ts"],"names":[],"mappings":"AAYA,OAAO,EAAE,aAAa,EAAE,MAAM,yBAAyB,CAAC;AAGxD,OAAO,EAAE,eAAe,EAAE,MAAM,sBAAsB,CAAC;AACvD,OAAO,EAAE,qBAAqB,EAAE,MAAM,aAAa,CAAC;AAEpD,MAAM,iBAAiB,GAAG,KAAK,CAAC;AAEhC;cACc;AACd,MAAM,CAAC,KAAK,UAAU,qBAAqB,CACzC,QAA0B,EAC1B,CAAmB,EACnB,MAAuB,EACvB,SAAS,GAAW,iBAAiB,EACrC,OAAmC;IAEnC,MAAM,MAAM,GAAG,OAAO,EAAE,MAAM,IAAI,MAAM,eAAe,EAAE,CAAC;IAC1D,MAAM,CAAC,GAAG,MAAM,CAAC,SAAS,IAAI,CAAC,CAAC;IAEhC,MAAM,QAAQ,GAAG,MAAM,CAAC,OAAO,CAAC,GAAG,CAAC,QAAQ,CAAC,IAAI,CAAC,CAAC;IACnD,MAAM,KAAK,GAAG,MAAM,CAAC,OAAO,CAAC,GAAG,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC;IACzC,IAAI,QAAQ,KAAK,SAAS,EAAE,CAAC;QAC3B,MAAM,IAAI,qBAAqB,CAC7B,yDAAyD,QAAQ,CAAC,IAAI,qBAAqB,CAC5F,CAAC;IACJ,CAAC;IACD,IAAI,KAAK,KAAK,SAAS,EAAE,CAAC;QACxB,MAAM,IAAI,qBAAqB,CAC7B,yDAAyD,CAAC,CAAC,IAAI,qBAAqB,CACrF,CAAC;IACJ,CAAC;IACD,MAAM,KAAK,GAAG,MAAM,CAAC,UAAU,CAAC,QAAQ,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;IAClD,MAAM,EAAE,GAAG,MAAM,CAAC,UAAU,CAAC,KAAK,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;IAC5C,MAAM,QAAQ,GAAG,MAAM,CAAC,GAAG,CAAC,MAAM,CAAC,MAAM,CAAC,KAAK,EAAE,EAAE,CAAC,EAAE,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC;IAC1E,MAAM,YAAY,GAAG,MAAM,CAAC,OAAO,CAAC,QAAQ,CAAC,CAAC;IAE9C,IAAI,YAAY,GAAG,SAAS,EAAE,CAAC;QAC7B,MAAM,OAAO,GAAc;YACzB,QAAQ,EAAE,GAAG,QAAQ,CAAC,IAAI,IAAI,CAAC,CAAC,IAAI,EAAE;YACtC,QAAQ,EAAE,aAAa;YACvB,MAAM,EAAE,aAAa;YACrB,IAAI,EAAE,uCAAuC,QAAQ,CAAC,IAAI,IAAI,CAAC,CAAC,IAAI,YAAY;kBAC5E,GAAG,YAAY,CAAC,aAAa,CAAC,CAAC,CAAC,sBAAsB,SAAS,CAAC,aAAa,CAAC,CAAC,CAAC,KAAK;kBACrF,iDAAiD;YACrD,QAAQ,EAAE,SAAS;SACpB,CAAC;QACF,OAAO,EAAE,YAAY,EAAE,OAAO,EAAE,CAAC;IACnC,CAAC;IACD,OAAO,EAAE,YAAY,EAAE,CAAC;AAC1B,CAAC;AAED;;;2EAG2E;AAC3E,MAAM,UAAU,iBAAiB,CAC/B,IAAc;IAEd,MAAM,OAAO,GAAuB,EAAE,CAAC;IACvC,MAAM,IAAI,GAAG,CAAC,CAAW,EAAQ,EAAE;QACjC,IAAI,CAAC,CAAC,IAAI,KAAK,eAAe;YAAE,OAAO,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;aAC3C,IAAI,CAAC,CAAC,IAAI,KAAK,IAAI;YAAE,CAAC,CAAC,IAAI,CAAC,OAAO,CAAC,IAAI,CAAC,CAAC;aAC1C,IAAI,CAAC,CAAC,IAAI,KAAK,gBAAgB;YAAE,CAAC,CAAC,IAAI,CAAC,OAAO,CAAC,IAAI,CAAC,CAAC;aACtD,IAAI,CAAC,CAAC,IAAI,KAAK,UAAU,EAAE,CAAC;YAAC,IAAI,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC;YAAC,IAAI,CAAC,CAAC,CAAC,SAAS,CAAC,CAAC;QAAC,CAAC;aAC/D,IAAI,CAAC,CAAC,IAAI,KAAK,YAAY,EAAE,CAAC;YAAC,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC;YAAC,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;QAAC,CAAC;aACzD,IAAI,CAAC,CAAC,IAAI,KAAK,2BAA2B,EAAE,CAAC;YAChD,IAAI,CAAC,CAAC,CAAC,EAAc,CAAC,CAAC;YAAC,IAAI,CAAC,CAAC,CAAC,GAAe,CAAC,CAAC;QAClD,CAAC;QACD,mEAAmE;QACnE,uEAAuE;QACvE,4EAA4E;QAC5E,0EAA0E;QAC1E,yEAAyE;QACzE,gCAAgC;QAChC,EAAE;QACF,sEAAsE;aACjE,IAAI,CAAC,CAAC,IAAI,KAAK,sBAAsB,EAAE,CAAC;YAC3C,IAAI,CAAC,CAAC,CAAC,EAAc,CAAC,CAAC;YACvB,IAAI,CAAC,CAAC,CAAC,MAAM,CAAC,CAAC;YACf,IAAI,CAAC,CAAC,CAAC,QAAQ,CAAC,CAAC;QACnB,CAAC;aACI,IAAI,CAAC,CAAC,IAAI,KAAK,gBAAgB,EAAE,CAAC;YACrC,IAAI,CAAC,CAAC,CAAC,MAAM,CAAC,CAAC;YACf,IAAI,CAAC,CAAC,CAAC,QAAQ,CAAC,CAAC;QACnB,CAAC;aACI,IAAI,CAAC,CAAC,IAAI,KAAK,cAAc,EAAE,CAAC;YACnC,IAAI,CAAC,CAAC,CAAC,OAA8B,CAAC,CAAC,CAAE,8CAA8C;QACzF,CAAC;aACI,IAAI,CAAC,CAAC,IAAI,KAAK,iBAAiB,EAAE,CAAC;YACtC,IAAI,CAAC,CAAC,CAAC,OAA8B,CAAC,CAAC;YACvC,IAAI,CAAC,CAAC,CAAC,MAAM,CAAC,CAAC;YACf,IAAI,CAAC,CAAC,CAAC,QAAQ,CAAC,CAAC;QACnB,CAAC;aACI,IAAI,CAAC,CAAC,IAAI,KAAK,kBAAkB,EAAE,CAAC;YACvC,IAAI,CAAC,CAAC,CAAC,OAA8B,CAAC,CAAC,CAAE,4BAA4B;QACvE,CAAC;QACD,qEAAqE;IACvE,CAAC,CAAC;IACF,IAAI,CAAC,IAAI,CAAC,CAAC;IACX,MAAM,KAAK,GAAG,OAAO,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,OAAO,CAAC,KAAK,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,QAAQ,KAAK,OAAO,CAAC,CAAC,CAAC;IAClF,MAAM,KAAK,GAAG,OAAO,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,OAAO,CAAC,KAAK,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,QAAQ,KAAK,OAAO,CAAC,CAAC,CAAC;IAClF,OAAO,KAAK,IAAI,KAAK,CAAC,CAAC,CAAC,EAAE,MAAM,EAAE,KAAK,EAAE,MAAM,EAAE,KAAK,EAAE,CAAC,CAAC,CAAC,IAAI,CAAC;AAClE,CAAC","sourcesContent":["/**\n * InverseMetricInconsistencyWarning — numerical path. Builds g⁻¹ and g as\n * EngineTensors, computes ‖matMul(g⁻¹, g) − I‖_∞, and emits a\n * 'warning'-severity Violation when it exceeds tolerance. v0.3.5-Design.md\n * §7. Resolves the v0.3.0 §13 Q2 deferral. `scanForMetricPair` finds an\n * identifiable lower/upper metric pair in an AST so evaluateNumerical can\n * run the check automatically.\n *\n * @module numerical/metric-inverse\n */\nimport type { ExprNode, Violation } from '../dimensional/validator.js';\nimport type { MetricTensorNode } from '../dimensional/metric-validators.js';\nimport { DIMENSIONLESS } from '../dimensional/types.js';\nimport type { TensorEngine } from './tensor-engine.js';\nimport type { NumericalInputs } from './types.js';\nimport { getActiveEngine } from './engine-registry.js';\nimport { NumericalBackendError } from './errors.js';\n\nconst DEFAULT_TOLERANCE = 1e-10;\n\n/** Compute ‖g⁻¹g − I‖_∞ and, if it exceeds tolerance, return a warning.\n *  @public */\nexport async function evaluateMetricInverse(\n  gInverse: MetricTensorNode,\n  g: MetricTensorNode,\n  inputs: NumericalInputs,\n  tolerance: number = DEFAULT_TOLERANCE,\n  options?: { engine?: TensorEngine },\n): Promise<{ residualNorm: number; warning?: Violation }> {\n  const engine = options?.engine ?? await getActiveEngine();\n  const N = inputs.dimension ?? 4;\n\n  const gInvData = inputs.tensors.get(gInverse.name);\n  const gData = inputs.tensors.get(g.name);\n  if (gInvData === undefined) {\n    throw new NumericalBackendError(\n      `evaluateMetricInverse: missing tensor components for \"${gInverse.name}\" in inputs.tensors`,\n    );\n  }\n  if (gData === undefined) {\n    throw new NumericalBackendError(\n      `evaluateMetricInverse: missing tensor components for \"${g.name}\" in inputs.tensors`,\n    );\n  }\n  const gInvT = engine.fromNested(gInvData, [N, N]);\n  const gT = engine.fromNested(gData, [N, N]);\n  const residual = engine.sub(engine.matMul(gInvT, gT), engine.identity(N));\n  const residualNorm = engine.normInf(residual);\n\n  if (residualNorm > tolerance) {\n    const warning: Violation = {\n      location: `${gInverse.name}·${g.name}`,\n      expected: DIMENSIONLESS,\n      actual: DIMENSIONLESS,\n      note: `InverseMetricInconsistencyWarning: ‖${gInverse.name}·${g.name} − I‖_∞ = `\n        + `${residualNorm.toExponential(3)} exceeds tolerance ${tolerance.toExponential(1)} — `\n        + `the supplied metrics are not numerical inverses`,\n      severity: 'warning',\n    };\n    return { residualNorm, warning };\n  }\n  return { residualNorm };\n}\n\n/** Walk an ExprNode tree; return the first all-lower / all-upper\n *  metric-tensor pair found, or null. Cheap — evaluateNumerical already\n *  walks the tree to lower it.\n *  @internal — cross-module use only; not part of the consumer surface. */\nexport function scanForMetricPair(\n  node: ExprNode,\n): { gLower: MetricTensorNode; gUpper: MetricTensorNode } | null {\n  const metrics: MetricTensorNode[] = [];\n  const walk = (n: ExprNode): void => {\n    if (n.kind === 'metric-tensor') metrics.push(n);\n    else if (n.kind === 'op') n.args.forEach(walk);\n    else if (n.kind === 'tensor-product') n.args.forEach(walk);\n    else if (n.kind === 'integral') { walk(n.over); walk(n.integrand); }\n    else if (n.kind === 'derivative') { walk(n.of); walk(n.wrt); }\n    else if (n.kind === 'tensor-partial-derivative') {\n      walk(n.of as ExprNode); walk(n.wrt as ExprNode);\n    }\n    // v0.5.1 PC-3: walk the v0.4.0+v0.5.0 curvature node kinds so that\n    // InverseMetricInconsistencyWarning fires on the GR pipeline. Prior to\n    // this, scan stopped at riemann/ricci/einstein/bianchi/covariant-derivative\n    // and the warning was silently never emitted on the whole curvature stack\n    // — \"a validation mechanism that doesn't run is worse than no validation\n    // at all\" (Adam, v0.5.1 audit).\n    //\n    // FUTURE: when a new curvature kind is added, this scan needs an arm.\n    else if (n.kind === 'covariant-derivative') {\n      walk(n.of as ExprNode);\n      walk(n.gLower);\n      walk(n.gInverse);\n    }\n    else if (n.kind === 'riemann-tensor') {\n      walk(n.gLower);\n      walk(n.gInverse);\n    }\n    else if (n.kind === 'ricci-tensor') {\n      walk(n.riemann as unknown as ExprNode);  // inner Riemann walks its own gLower/gInverse\n    }\n    else if (n.kind === 'einstein-tensor') {\n      walk(n.riemann as unknown as ExprNode);\n      walk(n.gLower);\n      walk(n.gInverse);\n    }\n    else if (n.kind === 'bianchi-residual') {\n      walk(n.riemann as unknown as ExprNode);  // recurse via inner Riemann\n    }\n    // tensor-symbol / kronecker-delta / symbol: leaves, no metric inside\n  };\n  walk(node);\n  const lower = metrics.find((m) => m.indices.every((i) => i.variance === 'lower'));\n  const upper = metrics.find((m) => m.indices.every((i) => i.variance === 'upper'));\n  return lower && upper ? { gLower: lower, gUpper: upper } : null;\n}\n"]}
+{"version":3,"file":"metric-inverse.js","sourceRoot":"","sources":["../../src/numerical/metric-inverse.ts"],"names":[],"mappings":"AAYA,OAAO,EAAE,aAAa,EAAE,MAAM,yBAAyB,CAAC;AAGxD,OAAO,EAAE,eAAe,EAAE,MAAM,sBAAsB,CAAC;AACvD,OAAO,EAAE,qBAAqB,EAAE,MAAM,aAAa,CAAC;AAEpD,MAAM,iBAAiB,GAAG,KAAK,CAAC;AAEhC;cACc;AACd,MAAM,CAAC,KAAK,UAAU,qBAAqB,CACzC,QAA0B,EAC1B,CAAmB,EACnB,MAAuB,EACvB,SAAS,GAAW,iBAAiB,EACrC,OAAmC;IAEnC,MAAM,MAAM,GAAG,OAAO,EAAE,MAAM,IAAI,MAAM,eAAe,EAAE,CAAC;IAC1D,MAAM,CAAC,GAAG,MAAM,CAAC,SAAS,IAAI,CAAC,CAAC;IAEhC,MAAM,QAAQ,GAAG,MAAM,CAAC,OAAO,CAAC,GAAG,CAAC,QAAQ,CAAC,IAAI,CAAC,CAAC;IACnD,MAAM,KAAK,GAAG,MAAM,CAAC,OAAO,CAAC,GAAG,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC;IACzC,IAAI,QAAQ,KAAK,SAAS,EAAE,CAAC;QAC3B,MAAM,IAAI,qBAAqB,CAC7B,yDAAyD,QAAQ,CAAC,IAAI,qBAAqB,CAC5F,CAAC;IACJ,CAAC;IACD,IAAI,KAAK,KAAK,SAAS,EAAE,CAAC;QACxB,MAAM,IAAI,qBAAqB,CAC7B,yDAAyD,CAAC,CAAC,IAAI,qBAAqB,CACrF,CAAC;IACJ,CAAC;IACD,MAAM,KAAK,GAAG,MAAM,CAAC,UAAU,CAAC,QAAQ,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;IAClD,MAAM,EAAE,GAAG,MAAM,CAAC,UAAU,CAAC,KAAK,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;IAC5C,MAAM,QAAQ,GAAG,MAAM,CAAC,GAAG,CAAC,MAAM,CAAC,MAAM,CAAC,KAAK,EAAE,EAAE,CAAC,EAAE,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC;IAC1E,MAAM,YAAY,GAAG,MAAM,CAAC,OAAO,CAAC,QAAQ,CAAC,CAAC;IAE9C,IAAI,YAAY,GAAG,SAAS,EAAE,CAAC;QAC7B,MAAM,OAAO,GAAc;YACzB,QAAQ,EAAE,GAAG,QAAQ,CAAC,IAAI,IAAI,CAAC,CAAC,IAAI,EAAE;YACtC,QAAQ,EAAE,aAAa;YACvB,MAAM,EAAE,aAAa;YACrB,IAAI,EAAE,uCAAuC,QAAQ,CAAC,IAAI,IAAI,CAAC,CAAC,IAAI,YAAY;kBAC5E,GAAG,YAAY,CAAC,aAAa,CAAC,CAAC,CAAC,sBAAsB,SAAS,CAAC,aAAa,CAAC,CAAC,CAAC,KAAK;kBACrF,iDAAiD;YACrD,QAAQ,EAAE,SAAS;SACpB,CAAC;QACF,OAAO,EAAE,YAAY,EAAE,OAAO,EAAE,CAAC;IACnC,CAAC;IACD,OAAO,EAAE,YAAY,EAAE,CAAC;AAC1B,CAAC;AAED;;;2EAG2E;AAC3E,MAAM,UAAU,iBAAiB,CAC/B,IAAc;IAEd,MAAM,OAAO,GAAuB,EAAE,CAAC;IACvC,4EAA4E;IAC5E,8EAA8E;IAC9E,8EAA8E;IAC9E,+EAA+E;IAC/E,MAAM,IAAI,GAAG,CAAC,CAAW,EAAQ,EAAE;QACjC,QAAQ,CAAC,CAAC,IAAI,EAAE,CAAC;YACf,KAAK,eAAe;gBAClB,OAAO,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;gBAChB,MAAM;YACR,KAAK,IAAI,CAAC;YACV,KAAK,gBAAgB;gBACnB,CAAC,CAAC,IAAI,CAAC,OAAO,CAAC,IAAI,CAAC,CAAC;gBACrB,MAAM;YACR,KAAK,UAAU;gBACb,IAAI,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC;gBACb,IAAI,CAAC,CAAC,CAAC,SAAS,CAAC,CAAC;gBAClB,IAAI,CAAC,CAAC,KAAK,KAAK,SAAS;oBAAE,IAAI,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC;gBACzC,IAAI,CAAC,CAAC,KAAK,KAAK,SAAS;oBAAE,IAAI,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC;gBACzC,MAAM;YACR,KAAK,YAAY;gBACf,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC;gBACX,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;gBACZ,MAAM;YACR,KAAK,wBAAwB;gBAC3B,IAAI,CAAC,CAAC,CAAC,UAAU,CAAC,CAAC;gBACnB,IAAI,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC;gBACd,IAAI,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC;gBACb,MAAM;YACR,KAAK,aAAa,CAAC;YACnB,KAAK,gBAAgB,CAAC;YACtB,KAAK,KAAK;gBACR,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;gBACZ,MAAM;YACR,KAAK,2BAA2B;gBAC9B,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC;gBACX,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;gBACZ,MAAM;YACR,KAAK,sBAAsB;gBACzB,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC;gBACX,IAAI,CAAC,CAAC,CAAC,MAAM,CAAC,CAAC;gBACf,IAAI,CAAC,CAAC,CAAC,QAAQ,CAAC,CAAC;gBACjB,MAAM;YACR,KAAK,gBAAgB;gBACnB,IAAI,CAAC,CAAC,CAAC,MAAM,CAAC,CAAC;gBACf,IAAI,CAAC,CAAC,CAAC,QAAQ,CAAC,CAAC;gBACjB,MAAM;YACR,KAAK,cAAc,CAAC;YACpB,KAAK,kBAAkB;gBACrB,IAAI,CAAC,CAAC,CAAC,OAAO,CAAC,CAAC,CAAC,mCAAmC;gBACpD,MAAM;YACR,KAAK,iBAAiB;gBACpB,IAAI,CAAC,CAAC,CAAC,OAAO,CAAC,CAAC;gBAChB,IAAI,CAAC,CAAC,CAAC,MAAM,CAAC,CAAC;gBACf,IAAI,CAAC,CAAC,CAAC,QAAQ,CAAC,CAAC;gBACjB,MAAM;YACR,KAAK,aAAa;gBAChB,IAAI,CAAC,CAAC,CAAC,MAAM,CAAC,CAAC;gBACf,MAAM;YACR,KAAK,oBAAoB;gBACvB,IAAI,CAAC,CAAC,CAAC,OAAO,CAAC,CAAC;gBAChB,IAAI,CAAC,CAAC,CAAC,MAAM,CAAC,CAAC;gBACf,MAAM;YACR,KAAK,gBAAgB;gBACnB,IAAI,CAAC,CAAC,CAAC,MAAM,CAAC,CAAC;gBACf,MAAM;YACR,KAAK,kBAAkB;gBACrB,IAAI,CAAC,CAAC,CAAC,OAAO,CAAC,CAAC;gBAChB,MAAM;YACR,KAAK,mBAAmB;gBACtB,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;gBACZ,IAAI,CAAC,CAAC,CAAC,MAAM,CAAC,CAAC;gBACf,MAAM;YACR,KAAK,QAAQ,CAAC;YACd,KAAK,eAAe,CAAC;YACrB,KAAK,iBAAiB,CAAC;YACvB,KAAK,eAAe,CAAC;YACrB,KAAK,uBAAuB;gBAC1B,MAAM,CAAC,2BAA2B;YACpC,SAAS,CAAC;gBACR,MAAM,WAAW,GAAU,CAAC,CAAC;gBAC7B,MAAM,IAAI,qBAAqB,CAAC,0CAA0C,MAAM,CAAE,WAAiC,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC;YAC/H,CAAC;QACH,CAAC;IACH,CAAC,CAAC;IACF,IAAI,CAAC,IAAI,CAAC,CAAC;IACX,MAAM,KAAK,GAAG,OAAO,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,OAAO,CAAC,KAAK,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,QAAQ,KAAK,OAAO,CAAC,CAAC,CAAC;IAClF,MAAM,KAAK,GAAG,OAAO,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,OAAO,CAAC,KAAK,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,QAAQ,KAAK,OAAO,CAAC,CAAC,CAAC;IAClF,OAAO,KAAK,IAAI,KAAK,CAAC,CAAC,CAAC,EAAE,MAAM,EAAE,KAAK,EAAE,MAAM,EAAE,KAAK,EAAE,CAAC,CAAC,CAAC,IAAI,CAAC;AAClE,CAAC","sourcesContent":["/**\n * InverseMetricInconsistencyWarning — numerical path. Builds g⁻¹ and g as\n * EngineTensors, computes ‖matMul(g⁻¹, g) − I‖_∞, and emits a\n * 'warning'-severity Violation when it exceeds tolerance. v0.3.5-Design.md\n * §7. Resolves the v0.3.0 §13 Q2 deferral. `scanForMetricPair` finds an\n * identifiable lower/upper metric pair in an AST so evaluateNumerical can\n * run the check automatically.\n *\n * @module numerical/metric-inverse\n */\nimport type { ExprNode, Violation } from '../dimensional/validator.js';\nimport type { MetricTensorNode } from '../dimensional/metric-validators.js';\nimport { DIMENSIONLESS } from '../dimensional/types.js';\nimport type { TensorEngine } from './tensor-engine.js';\nimport type { NumericalInputs } from './types.js';\nimport { getActiveEngine } from './engine-registry.js';\nimport { NumericalBackendError } from './errors.js';\n\nconst DEFAULT_TOLERANCE = 1e-10;\n\n/** Compute ‖g⁻¹g − I‖_∞ and, if it exceeds tolerance, return a warning.\n *  @public */\nexport async function evaluateMetricInverse(\n  gInverse: MetricTensorNode,\n  g: MetricTensorNode,\n  inputs: NumericalInputs,\n  tolerance: number = DEFAULT_TOLERANCE,\n  options?: { engine?: TensorEngine },\n): Promise<{ residualNorm: number; warning?: Violation }> {\n  const engine = options?.engine ?? await getActiveEngine();\n  const N = inputs.dimension ?? 4;\n\n  const gInvData = inputs.tensors.get(gInverse.name);\n  const gData = inputs.tensors.get(g.name);\n  if (gInvData === undefined) {\n    throw new NumericalBackendError(\n      `evaluateMetricInverse: missing tensor components for \"${gInverse.name}\" in inputs.tensors`,\n    );\n  }\n  if (gData === undefined) {\n    throw new NumericalBackendError(\n      `evaluateMetricInverse: missing tensor components for \"${g.name}\" in inputs.tensors`,\n    );\n  }\n  const gInvT = engine.fromNested(gInvData, [N, N]);\n  const gT = engine.fromNested(gData, [N, N]);\n  const residual = engine.sub(engine.matMul(gInvT, gT), engine.identity(N));\n  const residualNorm = engine.normInf(residual);\n\n  if (residualNorm > tolerance) {\n    const warning: Violation = {\n      location: `${gInverse.name}·${g.name}`,\n      expected: DIMENSIONLESS,\n      actual: DIMENSIONLESS,\n      note: `InverseMetricInconsistencyWarning: ‖${gInverse.name}·${g.name} − I‖_∞ = `\n        + `${residualNorm.toExponential(3)} exceeds tolerance ${tolerance.toExponential(1)} — `\n        + `the supplied metrics are not numerical inverses`,\n      severity: 'warning',\n    };\n    return { residualNorm, warning };\n  }\n  return { residualNorm };\n}\n\n/** Walk an ExprNode tree; return the first all-lower / all-upper\n *  metric-tensor pair found, or null. Cheap — evaluateNumerical already\n *  walks the tree to lower it.\n *  @internal — cross-module use only; not part of the consumer surface. */\nexport function scanForMetricPair(\n  node: ExprNode,\n): { gLower: MetricTensorNode; gUpper: MetricTensorNode } | null {\n  const metrics: MetricTensorNode[] = [];\n  // An exhaustive switch: a node kind added to `ExprNode` without an arm here\n  // fails to compile at the `never` arm, instead of becoming a silent leaf (the\n  // way `kretschmann-scalar` and `einstein-equation` were, so the warning never\n  // fired on them). A kind with no metric inside is listed as a leaf on purpose.\n  const walk = (n: ExprNode): void => {\n    switch (n.kind) {\n      case 'metric-tensor':\n        metrics.push(n);\n        break;\n      case 'op':\n      case 'tensor-product':\n        n.args.forEach(walk);\n        break;\n      case 'integral':\n        walk(n.over);\n        walk(n.integrand);\n        if (n.lower !== undefined) walk(n.lower);\n        if (n.upper !== undefined) walk(n.upper);\n        break;\n      case 'derivative':\n        walk(n.of);\n        walk(n.wrt);\n        break;\n      case 'variational-derivative':\n        walk(n.functional);\n        walk(n.field);\n        walk(n.over);\n        break;\n      case 'dirac-delta':\n      case 'transcendental':\n      case 'abs':\n        walk(n.arg);\n        break;\n      case 'tensor-partial-derivative':\n        walk(n.of);\n        walk(n.wrt);\n        break;\n      case 'covariant-derivative':\n        walk(n.of);\n        walk(n.gLower);\n        walk(n.gInverse);\n        break;\n      case 'riemann-tensor':\n        walk(n.gLower);\n        walk(n.gInverse);\n        break;\n      case 'ricci-tensor':\n      case 'bianchi-residual':\n        walk(n.riemann); // the inner Riemann holds the pair\n        break;\n      case 'einstein-tensor':\n        walk(n.riemann);\n        walk(n.gLower);\n        walk(n.gInverse);\n        break;\n      case 'weyl-tensor':\n        walk(n.metric);\n        break;\n      case 'kretschmann-scalar':\n        walk(n.riemann);\n        walk(n.metric);\n        break;\n      case 'killing-vector':\n        walk(n.metric);\n        break;\n      case 'conserved-charge':\n        walk(n.killing);\n        break;\n      case 'einstein-equation':\n        walk(n.lhs);\n        walk(n.metric);\n        break;\n      case 'symbol':\n      case 'tensor-symbol':\n      case 'kronecker-delta':\n      case 'stress-energy':\n      case 'cosmological-constant':\n        break; // leaves: no metric inside\n      default: {\n        const unreachable: never = n;\n        throw new NumericalBackendError(`scanForMetricPair: unhandled node kind ${String((unreachable as { kind: unknown }).kind)}`);\n      }\n    }\n  };\n  walk(node);\n  const lower = metrics.find((m) => m.indices.every((i) => i.variance === 'lower'));\n  const upper = metrics.find((m) => m.indices.every((i) => i.variance === 'upper'));\n  return lower && upper ? { gLower: lower, gUpper: upper } : null;\n}\n"]}

package/dist/numerical/null-ray-integrator.d.ts

@@ -5,7 +5,8 @@ type ODESystem = (lambda: number, y: ReadonlyArray<number>) => number[];
 /**
  * Integrate `system` from affine parameter `lambda0` to `lambda1` in
  * `steps` fixed RK4 steps, starting from state `y0`. Returns the final
- * state vector. Classical 4th-order Runge-Kutta — global error O(h⁴).
+ * state vector. Classical 4th-order Runge-Kutta — global error O(h⁴). It takes
+ * exactly `steps` steps of (lambda1 − lambda0)/steps (`rk4Trajectory`).
  *
  * @internal — cross-module/test use only; not part of the consumer surface.
  */

package/dist/numerical/null-ray-integrator.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"null-ray-integrator.d.ts","sourceRoot":"","sources":["../../src/numerical/null-ray-integrator.ts"],"names":[],"mappings":"AAcA;;+EAE+E;AAC/E,KAAK,SAAS,GAAG,CAAC,MAAM,EAAE,MAAM,EAAE,CAAC,EAAE,aAAa,CAAC,MAAM,CAAC,KAAK,MAAM,EAAE,CAAC;AAUxE;;;;;;GAMG;AACH,wBAAgB,YAAY,CAC1B,MAAM,EAAE,SAAS,EACjB,EAAE,EAAE,aAAa,CAAC,MAAM,CAAC,EACzB,OAAO,EAAE,MAAM,EACf,OAAO,EAAE,MAAM,EACf,KAAK,EAAE,MAAM,GACZ,MAAM,EAAE,CAiBV"}
+{"version":3,"file":"null-ray-integrator.d.ts","sourceRoot":"","sources":["../../src/numerical/null-ray-integrator.ts"],"names":[],"mappings":"AAcA;;+EAE+E;AAC/E,KAAK,SAAS,GAAG,CAAC,MAAM,EAAE,MAAM,EAAE,CAAC,EAAE,aAAa,CAAC,MAAM,CAAC,KAAK,MAAM,EAAE,CAAC;AAUxE;;;;;;;GAOG;AACH,wBAAgB,YAAY,CAC1B,MAAM,EAAE,SAAS,EACjB,EAAE,EAAE,aAAa,CAAC,MAAM,CAAC,EACzB,OAAO,EAAE,MAAM,EACf,OAAO,EAAE,MAAM,EACf,KAAK,EAAE,MAAM,GACZ,MAAM,EAAE,CAiBV"}