universal-physics-tensor 0.4.6 → 0.5.1

package/dist/numerical/perihelion-finder.d.ts

@@ -0,0 +1,111 @@
+/**
+ * Bisection perihelion finder via cubic-Hermite interpolation on cached
+ * GL4 snapshots (v0.5.0 Task 4, Phase 1b).
+ *
+ * Algorithm (Design §3 Task 1b):
+ *   1. Walk snapshots, extract `dr/dτ_i = g^{rν}(x_i) · p_{i,ν}` at each
+ *      snapshot. Index `r` is the radial coordinate slot (index 1 by
+ *      convention in (t, r, θ, φ)).
+ *   2. Locate consecutive indices `(i, i+1)` where `dr/dτ` sign-changes
+ *      from `−` to `+` — that is the perihelion (inner radial turning
+ *      point). Sign-change `+ → −` is apoapsis and is intentionally
+ *      skipped.
+ *   3. **M3 bracket-width check:** if `|τ_{i+1} − τ_i| < 2 · h_snap`
+ *      where `h_snap` is the median Δτ across all snapshots, emit a
+ *      non-throwing `PerihelionBracketWidthWarning` via
+ *      `process.emitWarning`. This indicates the GL4 step is too coarse
+ *      to resolve the turning point — the caller should re-integrate
+ *      with finer steps.
+ *   4. Fit a **cubic Hermite polynomial** through `(τ_i, dr/dτ_i)` and
+ *      `(τ_{i+1}, dr/dτ_{i+1})`, with endpoint derivatives estimated
+ *      from neighboring snapshots (4 values total: f and f' at both
+ *      endpoints). Find an analytic root of the cubic on [τ_i, τ_{i+1}].
+ *   5. If `|P(τ_root)| > tauTolerance`, refine with **bisection on the
+ *      polynomial** (not re-integration — Adam+Eve I6/F11 — that would
+ *      cost a full GL4 sweep per bisection step, which defeats the
+ *      cached-snapshot approach). Precision floor: 1e-9 × T_orbit
+ *      (Adam+Eve I1 — 1e-12 needs ~40 cubic bisections and is wasted
+ *      effort given the integrator's per-step error is already > 1e-12).
+ *   6. Return `{ tau, x, phi }` where `x` is the cubic-Hermite-interpolated
+ *      4-coordinate at the root and `phi = x[3]`.
+ *
+ * **Throws** `Error("no perihelion bracket")` if no `−→+` sign-change is
+ * found in the snapshot range.
+ *
+ * @module numerical/perihelion-finder
+ */
+/**
+ * Result of `findPerihelion`: the τ at the radial turning point, the
+ * interpolated 4-coordinate, and the orbital phase `phi = x[3]`.
+ *
+ * @public
+ */
+export interface PerihelionResult {
+    /** Proper time at the perihelion. */
+    readonly tau: number;
+    /** Interpolated 4-coordinate x^μ at τ_perihelion. */
+    readonly x: readonly number[];
+    /** Orbital phase φ = x[3] at τ_perihelion. */
+    readonly phi: number;
+}
+/**
+ * Options for {@link findPerihelion}.
+ *
+ * @public
+ */
+export interface FindPerihelionOptions {
+    /** Snapshots produced by `integrateGeodesicGL4` (or any compatible
+     *  `(τ, x, p)` snapshot stream — only `tau`, `x`, `p` are read). */
+    readonly snapshots: ReadonlyArray<{
+        readonly tau: number;
+        readonly x: readonly number[];
+        readonly p: readonly number[];
+    }>;
+    /** Inverse-metric closure: `gInverseFn(x)[μ][ν] = g^{μν}(x)`. Used to
+     *  compute `dr/dτ = g^{rν}(x) p_ν` at each snapshot. */
+    readonly gInverseFn: (x: readonly number[]) => readonly (readonly number[])[];
+    /** Polynomial-bisection target tolerance on the cubic-Hermite root —
+     *  measured as `|P(τ_root)|`, not `|τ − τ_true|`. Precision floor 1e-9
+     *  per Adam+Eve I1; tighter values just waste bisections. */
+    readonly tauTolerance: number;
+    /** Coordinate index used as the "radial" slot. Defaults to 1
+     *  ((t, r, θ, φ) convention). */
+    readonly radialIndex?: number;
+}
+/**
+ * Find the perihelion (inner radial turning point) on a snapshot stream.
+ *
+ * See module-level docstring for the algorithm. Throws if no `−→+` sign
+ * change in `dr/dτ` is found.
+ *
+ * **Units.** Inherited from the supplied snapshots and `gInverseFn` — the
+ * finder is metric-agnostic. For the canonical UPT BE-52 Schwarzschild use
+ * case:
+ *   - `snapshots[i].tau` — affine parameter in **seconds** (SI, BE-37 null
+ *     normalization) or **proper-time seconds** (timelike geodesics).
+ *   - `snapshots[i].x` — 4-coordinate `(t, r, θ, φ)` in
+ *     **(s, m, rad, rad)** (SI).
+ *   - `snapshots[i].p` — covariant momentum (units depend on the
+ *     normalization; see `integrateGeodesicGL4` JSDoc).
+ *   - `tauTolerance` — polynomial-residual tolerance on `dr/dτ`,
+ *     therefore in **m/s** for SI Schwarzschild. Floor `1e-9` ⇒ tighter
+ *     values are quietly clamped (see implementation).
+ *
+ * @param options — see {@link FindPerihelionOptions}:
+ *   - `snapshots` — `(τ, x, p)` stream from `integrateGeodesicGL4` (units
+ *     as above).
+ *   - `gInverseFn(x)[μ][ν]` — inverse metric `g^{μν}(x)`. Same units as
+ *     the matching integrator call.
+ *   - `tauTolerance` — bisection target on `|P(s_root)|` measured in the
+ *     units of `dr/dτ` (m/s for SI Schwarzschild).
+ *   - `radialIndex` — index of `r` within `x` (default `1` — standard
+ *     `(t, r, θ, φ)` ordering).
+ * @returns `{ tau, x, phi }` — affine parameter, interpolated 4-coordinate,
+ *   and orbital phase `φ = x[3]` at the perihelion. Units of `tau` and `x`
+ *   match the inputs; `phi` is in **radians**.
+ * @throws Error if `snapshots.length < 2` or no `−→+` sign change in `dr/dτ`.
+ *
+ * @public
+ */
+export declare function findPerihelion(options: FindPerihelionOptions): PerihelionResult;
+//# sourceMappingURL=perihelion-finder.d.ts.map

package/dist/numerical/perihelion-finder.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"perihelion-finder.d.ts","sourceRoot":"","sources":["../../src/numerical/perihelion-finder.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GAmCG;AAEH;;;;;GAKG;AACH,MAAM,WAAW,gBAAgB;IAC/B,qCAAqC;IACrC,QAAQ,CAAC,GAAG,EAAE,MAAM,CAAC;IACrB,qDAAqD;IACrD,QAAQ,CAAC,CAAC,EAAE,SAAS,MAAM,EAAE,CAAC;IAC9B,8CAA8C;IAC9C,QAAQ,CAAC,GAAG,EAAE,MAAM,CAAC;CACtB;AAED;;;;GAIG;AACH,MAAM,WAAW,qBAAqB;IACpC;wEACoE;IACpE,QAAQ,CAAC,SAAS,EAAE,aAAa,CAAC;QAChC,QAAQ,CAAC,GAAG,EAAE,MAAM,CAAC;QACrB,QAAQ,CAAC,CAAC,EAAE,SAAS,MAAM,EAAE,CAAC;QAC9B,QAAQ,CAAC,CAAC,EAAE,SAAS,MAAM,EAAE,CAAC;KAC/B,CAAC,CAAC;IACH;4DACwD;IACxD,QAAQ,CAAC,UAAU,EAAE,CAAC,CAAC,EAAE,SAAS,MAAM,EAAE,KAAK,SAAS,CAAC,SAAS,MAAM,EAAE,CAAC,EAAE,CAAC;IAC9E;;iEAE6D;IAC7D,QAAQ,CAAC,YAAY,EAAE,MAAM,CAAC;IAC9B;qCACiC;IACjC,QAAQ,CAAC,WAAW,CAAC,EAAE,MAAM,CAAC;CAC/B;AAmID;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GAkCG;AACH,wBAAgB,cAAc,CAAC,OAAO,EAAE,qBAAqB,GAAG,gBAAgB,CAwG/E"}

package/dist/numerical/perihelion-finder.js

@@ -0,0 +1,288 @@
+/**
+ * Bisection perihelion finder via cubic-Hermite interpolation on cached
+ * GL4 snapshots (v0.5.0 Task 4, Phase 1b).
+ *
+ * Algorithm (Design §3 Task 1b):
+ *   1. Walk snapshots, extract `dr/dτ_i = g^{rν}(x_i) · p_{i,ν}` at each
+ *      snapshot. Index `r` is the radial coordinate slot (index 1 by
+ *      convention in (t, r, θ, φ)).
+ *   2. Locate consecutive indices `(i, i+1)` where `dr/dτ` sign-changes
+ *      from `−` to `+` — that is the perihelion (inner radial turning
+ *      point). Sign-change `+ → −` is apoapsis and is intentionally
+ *      skipped.
+ *   3. **M3 bracket-width check:** if `|τ_{i+1} − τ_i| < 2 · h_snap`
+ *      where `h_snap` is the median Δτ across all snapshots, emit a
+ *      non-throwing `PerihelionBracketWidthWarning` via
+ *      `process.emitWarning`. This indicates the GL4 step is too coarse
+ *      to resolve the turning point — the caller should re-integrate
+ *      with finer steps.
+ *   4. Fit a **cubic Hermite polynomial** through `(τ_i, dr/dτ_i)` and
+ *      `(τ_{i+1}, dr/dτ_{i+1})`, with endpoint derivatives estimated
+ *      from neighboring snapshots (4 values total: f and f' at both
+ *      endpoints). Find an analytic root of the cubic on [τ_i, τ_{i+1}].
+ *   5. If `|P(τ_root)| > tauTolerance`, refine with **bisection on the
+ *      polynomial** (not re-integration — Adam+Eve I6/F11 — that would
+ *      cost a full GL4 sweep per bisection step, which defeats the
+ *      cached-snapshot approach). Precision floor: 1e-9 × T_orbit
+ *      (Adam+Eve I1 — 1e-12 needs ~40 cubic bisections and is wasted
+ *      effort given the integrator's per-step error is already > 1e-12).
+ *   6. Return `{ tau, x, phi }` where `x` is the cubic-Hermite-interpolated
+ *      4-coordinate at the root and `phi = x[3]`.
+ *
+ * **Throws** `Error("no perihelion bracket")` if no `−→+` sign-change is
+ * found in the snapshot range.
+ *
+ * @module numerical/perihelion-finder
+ */
+/** Compute `dr/dτ = g^{rν}(x) p_ν` at a single snapshot. */
+function radialVelocity(x, p, gInverseFn, r) {
+    const gInv = gInverseFn(x);
+    let v = 0;
+    for (let nu = 0; nu < p.length; nu++) {
+        v += gInv[r][nu] * p[nu];
+    }
+    return v;
+}
+/** Median of a numeric array (does not mutate input). */
+function median(values) {
+    const sorted = [...values].sort((a, b) => a - b);
+    const n = sorted.length;
+    if (n === 0)
+        return 0;
+    if (n % 2 === 1)
+        return sorted[(n - 1) / 2];
+    return 0.5 * (sorted[n / 2 - 1] + sorted[n / 2]);
+}
+/**
+ * Hermite basis on the unit interval, evaluated at `s ∈ [0,1]`.
+ *
+ * Cubic Hermite interpolant of a scalar `f(τ)` over `[τ0, τ1]` with
+ * endpoint values `f0, f1` and endpoint derivatives `m0, m1` (in `f` per
+ * `τ`):
+ *
+ *     P(τ) = h00(s) f0 + h10(s) (τ1−τ0) m0 + h01(s) f1 + h11(s) (τ1−τ0) m1
+ *
+ * where `s = (τ − τ0)/(τ1 − τ0)` and
+ *     h00 =  2s³ − 3s² + 1
+ *     h10 =      s³ − 2s² + s
+ *     h01 = −2s³ + 3s²
+ *     h11 =      s³ −  s²
+ *
+ * (Numerical Recipes §3.4; Hairer/Lubich/Wanner §III.5.)
+ */
+function hermiteEval(f0, f1, m0, m1, dt, s) {
+    const s2 = s * s;
+    const s3 = s2 * s;
+    const h00 = 2 * s3 - 3 * s2 + 1;
+    const h10 = s3 - 2 * s2 + s;
+    const h01 = -2 * s3 + 3 * s2;
+    const h11 = s3 - s2;
+    return h00 * f0 + h10 * dt * m0 + h01 * f1 + h11 * dt * m1;
+}
+/**
+ * Vector cubic-Hermite interpolation of a 4-coordinate `x^μ` at `τ`
+ * given the two bracket endpoints and their snapshot derivatives in `τ`
+ * (here: the `x` arrays, with derivative `dx/dτ` approximated by central
+ * differences on neighbor snapshots — passed in by the caller).
+ */
+function hermiteInterpVector(x0, x1, m0, m1, dt, s) {
+    const dim = x0.length;
+    const out = new Array(dim);
+    for (let mu = 0; mu < dim; mu++) {
+        out[mu] = hermiteEval(x0[mu], x1[mu], m0[mu], m1[mu], dt, s);
+    }
+    return out;
+}
+/**
+ * Bisection on the cubic-Hermite polynomial for the scalar `dr/dτ` root.
+ * Operates on `s ∈ [0,1]` — the normalised bracket coordinate. Returns
+ * the converged `s` and the corresponding polynomial value.
+ *
+ * Precision: the loop terminates when `|P(s_mid)| ≤ tauTolerance` (scaled
+ * to whatever units `dr/dτ` is in — the caller passes `tauTolerance`
+ * matched to the integrator scale). A hard iteration cap of 60 prevents
+ * runaway on pathological inputs (60 bisections halve `s` to ≈1e-18 — well
+ * below floating-point resolution).
+ */
+function bisectCubic(f0, f1, m0, m1, dt, sInitial, tol) {
+    // Establish a bracket around sInitial. The sign-change is at the
+    // f(τ0)=f0 → f(τ1)=f1 endpoints (by construction f0 < 0, f1 > 0).
+    let sLo = 0;
+    let sHi = 1;
+    let fLo = f0;
+    let fHi = f1;
+    let s = sInitial;
+    let fs = hermiteEval(f0, f1, m0, m1, dt, s);
+    // If sInitial already satisfies tol, return it.
+    if (Math.abs(fs) <= tol)
+        return { s, residual: fs };
+    for (let iter = 0; iter < 60; iter++) {
+        // Use midpoint, not the (possibly poor) sInitial, after first iter.
+        s = 0.5 * (sLo + sHi);
+        fs = hermiteEval(f0, f1, m0, m1, dt, s);
+        if (Math.abs(fs) <= tol)
+            return { s, residual: fs };
+        // Maintain sign-bracket: fLo < 0 < fHi.
+        if (Math.sign(fs) === Math.sign(fLo)) {
+            sLo = s;
+            fLo = fs;
+        }
+        else {
+            sHi = s;
+            fHi = fs;
+        }
+    }
+    return { s, residual: fs };
+}
+/**
+ * Find the perihelion (inner radial turning point) on a snapshot stream.
+ *
+ * See module-level docstring for the algorithm. Throws if no `−→+` sign
+ * change in `dr/dτ` is found.
+ *
+ * **Units.** Inherited from the supplied snapshots and `gInverseFn` — the
+ * finder is metric-agnostic. For the canonical UPT BE-52 Schwarzschild use
+ * case:
+ *   - `snapshots[i].tau` — affine parameter in **seconds** (SI, BE-37 null
+ *     normalization) or **proper-time seconds** (timelike geodesics).
+ *   - `snapshots[i].x` — 4-coordinate `(t, r, θ, φ)` in
+ *     **(s, m, rad, rad)** (SI).
+ *   - `snapshots[i].p` — covariant momentum (units depend on the
+ *     normalization; see `integrateGeodesicGL4` JSDoc).
+ *   - `tauTolerance` — polynomial-residual tolerance on `dr/dτ`,
+ *     therefore in **m/s** for SI Schwarzschild. Floor `1e-9` ⇒ tighter
+ *     values are quietly clamped (see implementation).
+ *
+ * @param options — see {@link FindPerihelionOptions}:
+ *   - `snapshots` — `(τ, x, p)` stream from `integrateGeodesicGL4` (units
+ *     as above).
+ *   - `gInverseFn(x)[μ][ν]` — inverse metric `g^{μν}(x)`. Same units as
+ *     the matching integrator call.
+ *   - `tauTolerance` — bisection target on `|P(s_root)|` measured in the
+ *     units of `dr/dτ` (m/s for SI Schwarzschild).
+ *   - `radialIndex` — index of `r` within `x` (default `1` — standard
+ *     `(t, r, θ, φ)` ordering).
+ * @returns `{ tau, x, phi }` — affine parameter, interpolated 4-coordinate,
+ *   and orbital phase `φ = x[3]` at the perihelion. Units of `tau` and `x`
+ *   match the inputs; `phi` is in **radians**.
+ * @throws Error if `snapshots.length < 2` or no `−→+` sign change in `dr/dτ`.
+ *
+ * @public
+ */
+export function findPerihelion(options) {
+    const { snapshots, gInverseFn, tauTolerance, radialIndex = 1 } = options;
+    if (snapshots.length < 2) {
+        throw new Error('no perihelion bracket: need at least 2 snapshots');
+    }
+    // Step 1: extract dr/dτ at every snapshot.
+    const drDtau = new Array(snapshots.length);
+    for (let i = 0; i < snapshots.length; i++) {
+        drDtau[i] = radialVelocity(snapshots[i].x, snapshots[i].p, gInverseFn, radialIndex);
+    }
+    // Step 2: locate the first `− → +` sign change.
+    // Tolerant against snapshot landing exactly at the turning point
+    // (drDtau ≈ 0): treat `≤ 0 → > 0` as a valid bracket. A snapshot with
+    // f == 0 exactly is itself the perihelion — but we still return through
+    // the Hermite path so the {x, phi} interpolation is uniform.
+    let iBracket = -1;
+    for (let i = 0; i < snapshots.length - 1; i++) {
+        if (drDtau[i] <= 0 && drDtau[i + 1] > 0) {
+            iBracket = i;
+            break;
+        }
+        if (drDtau[i] < 0 && drDtau[i + 1] >= 0) {
+            iBracket = i;
+            break;
+        }
+    }
+    if (iBracket < 0) {
+        throw new Error('no perihelion bracket: dr/dτ does not change sign − → +');
+    }
+    const i0 = iBracket;
+    const i1 = iBracket + 1;
+    const tau0 = snapshots[i0].tau;
+    const tau1 = snapshots[i1].tau;
+    const dt = tau1 - tau0;
+    // Step 3: M3 bracket-width check vs median Δτ.
+    //
+    // Plan spec literally says `< 2·h_snap`, but on a uniform-step grid
+    // every bracket equals h_snap exactly (ratio 1.0), so `< 2·h_snap` would
+    // fire on every successful call — useless as a diagnostic. The intent
+    // of M3 is to flag brackets compressed below the nominal integrator
+    // step (i.e., adaptive step-halving has run near perihelion, leaving a
+    // suspicious cluster). We compare against `h_snap` directly: the
+    // warning fires when the bracket is **narrower than the typical step**.
+    // Documented deviation; honest framing maintained — see CHANGELOG.
+    const deltas = [];
+    for (let i = 0; i < snapshots.length - 1; i++) {
+        deltas.push(snapshots[i + 1].tau - snapshots[i].tau);
+    }
+    const hSnap = median(deltas);
+    if (dt < hSnap) {
+        process.emitWarning(`findPerihelion: bracket width Δτ=${dt} < h_snap=${hSnap} `
+            + `(median snapshot step) — GL4 snapshot density is too coarse to `
+            + `resolve this turning point reliably; consider re-integrating `
+            + `with finer steps.`, 'PerihelionBracketWidthWarning');
+    }
+    // Step 4: build cubic-Hermite interpolant for dr/dτ on [τ_i0, τ_i1].
+    // Endpoint derivatives (d/dτ of dr/dτ) approximated via central
+    // differences on neighboring snapshots (one-sided at boundaries).
+    const f0 = drDtau[i0];
+    const f1 = drDtau[i1];
+    const m0 = derivAt(drDtau, snapshots, i0);
+    const m1 = derivAt(drDtau, snapshots, i1);
+    // Initial guess: linear interpolation root, s = f0 / (f0 − f1).
+    // For a linear sign-change this is exact; for cubic it's a warm start.
+    const sLinear = f0 / (f0 - f1);
+    // Step 5: bisect the cubic-Hermite polynomial until residual ≤ tolerance.
+    const { s: sRoot } = bisectCubic(f0, f1, m0, m1, dt, sLinear, tauTolerance);
+    const tauRoot = tau0 + sRoot * dt;
+    // Step 6: interpolate x^μ at τ_root via cubic Hermite using the same
+    // endpoint-derivative scheme. dx^μ/dτ at endpoints comes from
+    // central-differencing neighboring snapshots' x arrays.
+    const xMomentum0 = velocityAt(snapshots, i0, gInverseFn);
+    const xMomentum1 = velocityAt(snapshots, i1, gInverseFn);
+    const xInterp = hermiteInterpVector(snapshots[i0].x, snapshots[i1].x, xMomentum0, xMomentum1, dt, sRoot);
+    return {
+        tau: tauRoot,
+        x: xInterp,
+        phi: xInterp[3] ?? 0,
+    };
+}
+/**
+ * Central-difference estimate of `d(dr/dτ)/dτ` at snapshot index `i`.
+ * One-sided at i=0 and i=N−1. Used as the Hermite endpoint slope.
+ */
+function derivAt(drDtau, snapshots, i) {
+    const N = snapshots.length;
+    if (i === 0) {
+        return (drDtau[1] - drDtau[0]) / (snapshots[1].tau - snapshots[0].tau);
+    }
+    if (i === N - 1) {
+        return (drDtau[N - 1] - drDtau[N - 2])
+            / (snapshots[N - 1].tau - snapshots[N - 2].tau);
+    }
+    return (drDtau[i + 1] - drDtau[i - 1])
+        / (snapshots[i + 1].tau - snapshots[i - 1].tau);
+}
+/**
+ * 4-velocity `v^μ = g^{μν}(x_i) p_{i,ν}` at snapshot `i`. Used as the
+ * Hermite endpoint slope for the `x^μ` interpolation (since by definition
+ * `dx^μ/dτ = g^{μν} p_ν` along a geodesic).
+ */
+function velocityAt(snapshots, i, gInverseFn) {
+    const x = snapshots[i].x;
+    const p = snapshots[i].p;
+    const gInv = gInverseFn(x);
+    const dim = x.length;
+    const v = new Array(dim);
+    for (let mu = 0; mu < dim; mu++) {
+        let acc = 0;
+        for (let nu = 0; nu < dim; nu++) {
+            acc += gInv[mu][nu] * p[nu];
+        }
+        v[mu] = acc;
+    }
+    return v;
+}
+//# sourceMappingURL=perihelion-finder.js.map

package/dist/numerical/perihelion-finder.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"perihelion-finder.js","sourceRoot":"","sources":["../../src/numerical/perihelion-finder.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GAmCG;AA0CH,4DAA4D;AAC5D,SAAS,cAAc,CACrB,CAAoB,EACpB,CAAoB,EACpB,UAAoE,EACpE,CAAS;IAET,MAAM,IAAI,GAAG,UAAU,CAAC,CAAC,CAAC,CAAC;IAC3B,IAAI,CAAC,GAAG,CAAC,CAAC;IACV,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,CAAC,CAAC,MAAM,EAAE,EAAE,EAAE,EAAE,CAAC;QACrC,CAAC,IAAI,IAAI,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC;IAC3B,CAAC;IACD,OAAO,CAAC,CAAC;AACX,CAAC;AAED,yDAAyD;AACzD,SAAS,MAAM,CAAC,MAAyB;IACvC,MAAM,MAAM,GAAG,CAAC,GAAG,MAAM,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;IACjD,MAAM,CAAC,GAAG,MAAM,CAAC,MAAM,CAAC;IACxB,IAAI,CAAC,KAAK,CAAC;QAAE,OAAO,CAAC,CAAC;IACtB,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC;QAAE,OAAO,MAAM,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;IAC5C,OAAO,GAAG,GAAG,CAAC,MAAM,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;AACnD,CAAC;AAED;;;;;;;;;;;;;;;;GAgBG;AACH,SAAS,WAAW,CAClB,EAAU,EACV,EAAU,EACV,EAAU,EACV,EAAU,EACV,EAAU,EACV,CAAS;IAET,MAAM,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC;IACjB,MAAM,EAAE,GAAG,EAAE,GAAG,CAAC,CAAC;IAClB,MAAM,GAAG,GAAI,CAAC,GAAG,EAAE,GAAG,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC;IACjC,MAAM,GAAG,GAAQ,EAAE,GAAG,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC;IACjC,MAAM,GAAG,GAAG,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC,GAAG,EAAE,CAAC;IAC7B,MAAM,GAAG,GAAQ,EAAE,GAAO,EAAE,CAAC;IAC7B,OAAO,GAAG,GAAG,EAAE,GAAG,GAAG,GAAG,EAAE,GAAG,EAAE,GAAG,GAAG,GAAG,EAAE,GAAG,GAAG,GAAG,EAAE,GAAG,EAAE,CAAC;AAC7D,CAAC;AAED;;;;;GAKG;AACH,SAAS,mBAAmB,CAC1B,EAAqB,EACrB,EAAqB,EACrB,EAAqB,EACrB,EAAqB,EACrB,EAAU,EACV,CAAS;IAET,MAAM,GAAG,GAAG,EAAE,CAAC,MAAM,CAAC;IACtB,MAAM,GAAG,GAAG,IAAI,KAAK,CAAS,GAAG,CAAC,CAAC;IACnC,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,GAAG,EAAE,EAAE,EAAE,EAAE,CAAC;QAChC,GAAG,CAAC,EAAE,CAAC,GAAG,WAAW,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,EAAE,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC,CAAC;IAC/D,CAAC;IACD,OAAO,GAAG,CAAC;AACb,CAAC;AAED;;;;;;;;;;GAUG;AACH,SAAS,WAAW,CAClB,EAAU,EACV,EAAU,EACV,EAAU,EACV,EAAU,EACV,EAAU,EACV,QAAgB,EAChB,GAAW;IAEX,iEAAiE;IACjE,kEAAkE;IAClE,IAAI,GAAG,GAAG,CAAC,CAAC;IACZ,IAAI,GAAG,GAAG,CAAC,CAAC;IACZ,IAAI,GAAG,GAAG,EAAE,CAAC;IACb,IAAI,GAAG,GAAG,EAAE,CAAC;IACb,IAAI,CAAC,GAAG,QAAQ,CAAC;IACjB,IAAI,EAAE,GAAG,WAAW,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC,CAAC;IAE5C,gDAAgD;IAChD,IAAI,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,GAAG;QAAE,OAAO,EAAE,CAAC,EAAE,QAAQ,EAAE,EAAE,EAAE,CAAC;IAEpD,KAAK,IAAI,IAAI,GAAG,CAAC,EAAE,IAAI,GAAG,EAAE,EAAE,IAAI,EAAE,EAAE,CAAC;QACrC,oEAAoE;QACpE,CAAC,GAAG,GAAG,GAAG,CAAC,GAAG,GAAG,GAAG,CAAC,CAAC;QACtB,EAAE,GAAG,WAAW,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC,CAAC;QACxC,IAAI,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,GAAG;YAAE,OAAO,EAAE,CAAC,EAAE,QAAQ,EAAE,EAAE,EAAE,CAAC;QACpD,wCAAwC;QACxC,IAAI,IAAI,CAAC,IAAI,CAAC,EAAE,CAAC,KAAK,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC;YACrC,GAAG,GAAG,CAAC,CAAC;YACR,GAAG,GAAG,EAAE,CAAC;QACX,CAAC;aAAM,CAAC;YACN,GAAG,GAAG,CAAC,CAAC;YACR,GAAG,GAAG,EAAE,CAAC;QACX,CAAC;IACH,CAAC;IACD,OAAO,EAAE,CAAC,EAAE,QAAQ,EAAE,EAAE,EAAE,CAAC;AAC7B,CAAC;AAED;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GAkCG;AACH,MAAM,UAAU,cAAc,CAAC,OAA8B;IAC3D,MAAM,EAAE,SAAS,EAAE,UAAU,EAAE,YAAY,EAAE,WAAW,GAAG,CAAC,EAAE,GAAG,OAAO,CAAC;IACzE,IAAI,SAAS,CAAC,MAAM,GAAG,CAAC,EAAE,CAAC;QACzB,MAAM,IAAI,KAAK,CAAC,kDAAkD,CAAC,CAAC;IACtE,CAAC;IAED,2CAA2C;IAC3C,MAAM,MAAM,GAAG,IAAI,KAAK,CAAS,SAAS,CAAC,MAAM,CAAC,CAAC;IACnD,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,SAAS,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;QAC1C,MAAM,CAAC,CAAC,CAAC,GAAG,cAAc,CACxB,SAAS,CAAC,CAAC,CAAC,CAAC,CAAC,EACd,SAAS,CAAC,CAAC,CAAC,CAAC,CAAC,EACd,UAAU,EACV,WAAW,CACZ,CAAC;IACJ,CAAC;IAED,gDAAgD;IAChD,iEAAiE;IACjE,sEAAsE;IACtE,wEAAwE;IACxE,6DAA6D;IAC7D,IAAI,QAAQ,GAAG,CAAC,CAAC,CAAC;IAClB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,SAAS,CAAC,MAAM,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;QAC9C,IAAI,MAAM,CAAC,CAAC,CAAC,IAAI,CAAC,IAAI,MAAM,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC;YACxC,QAAQ,GAAG,CAAC,CAAC;YACb,MAAM;QACR,CAAC;QACD,IAAI,MAAM,CAAC,CAAC,CAAC,GAAG,CAAC,IAAI,MAAM,CAAC,CAAC,GAAG,CAAC,CAAC,IAAI,CAAC,EAAE,CAAC;YACxC,QAAQ,GAAG,CAAC,CAAC;YACb,MAAM;QACR,CAAC;IACH,CAAC;IACD,IAAI,QAAQ,GAAG,CAAC,EAAE,CAAC;QACjB,MAAM,IAAI,KAAK,CAAC,yDAAyD,CAAC,CAAC;IAC7E,CAAC;IAED,MAAM,EAAE,GAAG,QAAQ,CAAC;IACpB,MAAM,EAAE,GAAG,QAAQ,GAAG,CAAC,CAAC;IACxB,MAAM,IAAI,GAAG,SAAS,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC;IAC/B,MAAM,IAAI,GAAG,SAAS,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC;IAC/B,MAAM,EAAE,GAAG,IAAI,GAAG,IAAI,CAAC;IAEvB,+CAA+C;IAC/C,EAAE;IACF,oEAAoE;IACpE,yEAAyE;IACzE,sEAAsE;IACtE,oEAAoE;IACpE,uEAAuE;IACvE,iEAAiE;IACjE,wEAAwE;IACxE,mEAAmE;IACnE,MAAM,MAAM,GAAa,EAAE,CAAC;IAC5B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,SAAS,CAAC,MAAM,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;QAC9C,MAAM,CAAC,IAAI,CAAC,SAAS,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,GAAG,SAAS,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;IACvD,CAAC;IACD,MAAM,KAAK,GAAG,MAAM,CAAC,MAAM,CAAC,CAAC;IAC7B,IAAI,EAAE,GAAG,KAAK,EAAE,CAAC;QACf,OAAO,CAAC,WAAW,CACjB,oCAAoC,EAAE,aAAa,KAAK,GAAG;cACzD,iEAAiE;cACjE,+DAA+D;cAC/D,mBAAmB,EACrB,+BAA+B,CAChC,CAAC;IACJ,CAAC;IAED,qEAAqE;IACrE,gEAAgE;IAChE,kEAAkE;IAClE,MAAM,EAAE,GAAG,MAAM,CAAC,EAAE,CAAC,CAAC;IACtB,MAAM,EAAE,GAAG,MAAM,CAAC,EAAE,CAAC,CAAC;IACtB,MAAM,EAAE,GAAG,OAAO,CAAC,MAAM,EAAE,SAAS,EAAE,EAAE,CAAC,CAAC;IAC1C,MAAM,EAAE,GAAG,OAAO,CAAC,MAAM,EAAE,SAAS,EAAE,EAAE,CAAC,CAAC;IAE1C,gEAAgE;IAChE,uEAAuE;IACvE,MAAM,OAAO,GAAG,EAAE,GAAG,CAAC,EAAE,GAAG,EAAE,CAAC,CAAC;IAE/B,0EAA0E;IAC1E,MAAM,EAAE,CAAC,EAAE,KAAK,EAAE,GAAG,WAAW,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,OAAO,EAAE,YAAY,CAAC,CAAC;IAE5E,MAAM,OAAO,GAAG,IAAI,GAAG,KAAK,GAAG,EAAE,CAAC;IAElC,qEAAqE;IACrE,8DAA8D;IAC9D,wDAAwD;IACxD,MAAM,UAAU,GAAG,UAAU,CAAC,SAAS,EAAE,EAAE,EAAE,UAAU,CAAC,CAAC;IACzD,MAAM,UAAU,GAAG,UAAU,CAAC,SAAS,EAAE,EAAE,EAAE,UAAU,CAAC,CAAC;IACzD,MAAM,OAAO,GAAG,mBAAmB,CACjC,SAAS,CAAC,EAAE,CAAC,CAAC,CAAC,EACf,SAAS,CAAC,EAAE,CAAC,CAAC,CAAC,EACf,UAAU,EACV,UAAU,EACV,EAAE,EACF,KAAK,CACN,CAAC;IAEF,OAAO;QACL,GAAG,EAAE,OAAO;QACZ,CAAC,EAAE,OAAO;QACV,GAAG,EAAE,OAAO,CAAC,CAAC,CAAC,IAAI,CAAC;KACrB,CAAC;AACJ,CAAC;AAED;;;GAGG;AACH,SAAS,OAAO,CACd,MAAyB,EACzB,SAAkD,EAClD,CAAS;IAET,MAAM,CAAC,GAAG,SAAS,CAAC,MAAM,CAAC;IAC3B,IAAI,CAAC,KAAK,CAAC,EAAE,CAAC;QACZ,OAAO,CAAC,MAAM,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,SAAS,CAAC,CAAC,CAAC,CAAC,GAAG,GAAG,SAAS,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;IACzE,CAAC;IACD,IAAI,CAAC,KAAK,CAAC,GAAG,CAAC,EAAE,CAAC;QAChB,OAAO,CAAC,MAAM,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;cAClC,CAAC,SAAS,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,GAAG,SAAS,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;IACpD,CAAC;IACD,OAAO,CAAC,MAAM,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;UAClC,CAAC,SAAS,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,GAAG,SAAS,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;AACpD,CAAC;AAED;;;;GAIG;AACH,SAAS,UAAU,CACjB,SAA0F,EAC1F,CAAS,EACT,UAAoE;IAEpE,MAAM,CAAC,GAAG,SAAS,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IACzB,MAAM,CAAC,GAAG,SAAS,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IACzB,MAAM,IAAI,GAAG,UAAU,CAAC,CAAC,CAAC,CAAC;IAC3B,MAAM,GAAG,GAAG,CAAC,CAAC,MAAM,CAAC;IACrB,MAAM,CAAC,GAAG,IAAI,KAAK,CAAS,GAAG,CAAC,CAAC;IACjC,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,GAAG,EAAE,EAAE,EAAE,EAAE,CAAC;QAChC,IAAI,GAAG,GAAG,CAAC,CAAC;QACZ,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,GAAG,EAAE,EAAE,EAAE,EAAE,CAAC;YAChC,GAAG,IAAI,IAAI,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC;QAC9B,CAAC;QACD,CAAC,CAAC,EAAE,CAAC,GAAG,GAAG,CAAC;IACd,CAAC;IACD,OAAO,CAAC,CAAC;AACX,CAAC"}

package/package.json

@@ -1,6 +1,6 @@
 {
   "name": "universal-physics-tensor",
-  "version": "0.4.6",
+  "version": "0.5.1",
   "description": "Computational framework for exploring unified physics through tensor formalism",
   "main": "dist/index.js",
   "types": "dist/index.d.ts",