universal-physics-tensor 0.45.1 → 0.46.0

package/dist/atlas/diffusion/numerics.d.ts

@@ -0,0 +1,181 @@
+/**
+ * Atlas Phase 4, S4.4 — the numerics behind the diffusion family's witnesses.
+ *
+ * Each function evaluates one claim at a RESOLUTION, so the numeric witness
+ * runner can evaluate it twice and record how the error shrank. None of them
+ * asserts an order of convergence; the runner reports the ratio and a reviewer
+ * reads it (Phase 4 design note §2.3).
+ *
+ * @module atlas/diffusion/numerics
+ */
+/**
+ * The Gaussian solution of `∂c/∂t = D ∂²c/∂x²` from a point source of unit
+ * mass, evaluated at `x`: `(4πDt)^{-1/2} exp(−x²/(4Dt))`.
+ *
+ * @internal
+ */
+export declare function diffusionKernel(x: number, t: number, D: number): number;
+/**
+ * WD1 — the probability DENSITY at the origin of a symmetric lattice random
+ * walk after `steps` steps, under the closure `D = Δx²/(2Δt)` at fixed `D`
+ * and `t`.
+ *
+ * With `Δt = t/steps` and `Δx = √(2DΔt)`, the walker sits at the origin with
+ * probability `C(steps, steps/2)/2^steps`. Only every other site is reachable
+ * after an even number of steps, so the probability is spread over a cell of
+ * width `2Δx`, and the density is that probability divided by `2Δx`.
+ * As `steps → ∞` it converges to {@link diffusionKernel}`(0, t, D)`.
+ *
+ * @param steps - the resolution; rounded up to the next even integer.
+ * @internal
+ */
+export declare function randomWalkCentralDensity(steps: number, D: number, t: number): number;
+/**
+ * WD1b — the random walk's density at lattice site `site` after `steps` steps
+ * (same closure and grid as {@link randomWalkCentralDensity}). Exactly zero
+ * outside the light cone `|site| > steps` and on sites of the wrong parity.
+ *
+ * @internal
+ */
+export declare function randomWalkDensity(site: number, steps: number, D: number, t: number): number;
+/** Parameters of the WD2 heat-conduction fixture. @internal */
+export interface HeatFixture {
+    /** Thermal conductivity κ. */
+    readonly kappa: number;
+    /** Density ρ. */
+    readonly rho: number;
+    /** Specific heat capacity c_p. */
+    readonly cp: number;
+    /** Initial Gaussian width parameter: `T(x, 0) = exp(−x²/(4 s0))`. */
+    readonly s0: number;
+    /** Final time. */
+    readonly tEnd: number;
+    /**
+     * Half-width of the domain. The edges are held at ZERO, not at the exact
+     * solution: the exact solution is the Fick side of the claim, and feeding it
+     * to the heat solver would leak the dictionary under test into the side that
+     * is supposed to be independent of it. The true edge value is below
+     * `exp(−halfWidth²/(4(s0 + D tEnd)))`, far under the tolerance.
+     */
+    readonly halfWidth: number;
+}
+/**
+ * The exact solution of the DIFFUSION equation with coefficient `D` from the
+ * WD2 initial condition: `√(s0/(s0 + D t)) · exp(−x²/(4(s0 + D t)))`.
+ *
+ * @internal
+ */
+export declare function gaussianSpread(x: number, t: number, D: number, s0: number): number;
+/**
+ * WD2 — the temperature at `x = 0, t = tEnd` from an explicit FTCS solution of
+ * the HEAT equation `ρ c_p ∂T/∂t = κ ∂²T/∂x²` on `cells` intervals.
+ *
+ * The heat equation is solved in its OWN variables: the step uses `κ`, `ρ` and
+ * `c_p`, never a pre-combined `D`. Comparing the result against
+ * {@link gaussianSpread} with `D = κ/(ρ c_p)` is what tests the dictionary.
+ * The mesh ratio `κΔt/(ρ c_p Δx²)` is held at or below 1/4, inside the FTCS
+ * stability limit of 1/2, so refinement halves `Δx` and quarters `Δt`.
+ *
+ * @param cells - the resolution: the number of spatial intervals (made even so
+ * `x = 0` is a node).
+ * @internal
+ */
+export declare function heatFtcsCentre(cells: number, f: HeatFixture): number;
+/** Parameters of the WD3 Wick-rotation fixture. @internal */
+export interface WickFixture {
+    readonly hbar: number;
+    readonly m: number;
+    /** Initial width of the α = 0 Gaussian. */
+    readonly s: number;
+    /** The point `(x, τ)` the residual is evaluated at. */
+    readonly x: number;
+    readonly tau: number;
+    /** Base finite-difference step; the resolution divides it. */
+    readonly h0: number;
+}
+/**
+ * WD3 — the finite-difference residual `|∂_τ φ − (ħ/2m) ∂²_x φ|` of the
+ * Wick-rotated free kernel `φ(x, τ)` (`witnesses/quantum-support.ts`, W5) at
+ * one point, with step `h = h0/resolution` in both variables.
+ *
+ * The residual measures whether the Wick rotation `t = −iτ` of the free
+ * Schrödinger equation really lands on the DIFFUSION equation with
+ * `D = ħ/(2m)`. Its exact value is zero; the finite-difference estimate
+ * approaches zero as `h` shrinks, and the runner records how fast.
+ *
+ * @internal
+ */
+export declare function wickHeatResidual(resolution: number, f: WickFixture, kernel: (x: number, tau: number, s: number, hbar: number, m: number) => number): number;
+/**
+ * WD3b — `∫ φ(x, τ)² dx` of the Wick-rotated free kernel, by the trapezoid
+ * rule over `[−halfWidth, halfWidth]`. The Schrödinger norm `∫|ψ|² dx` is
+ * conserved; this one DECAYS with τ, which is the unitarity the continuation
+ * does not preserve.
+ *
+ * @internal
+ */
+export declare function wickKernelSquaredNorm(tau: number, f: WickFixture, kernel: (x: number, tau: number, s: number, hbar: number, m: number) => number, halfWidth?: number, cells?: number): number;
+/** Parameters of the Langevin (Ornstein–Uhlenbeck) fixture. @internal */
+export interface LangevinFixture {
+    /** Particle mass m. */
+    readonly m: number;
+    /** Friction coefficient γ (force per velocity). */
+    readonly gamma: number;
+    /** Thermal energy k_B T. */
+    readonly kT: number;
+}
+/**
+ * WD4 — `⟨x²⟩(t) / (2Dt)` for the Langevin model `m v̇ = −γv + ξ`,
+ * `⟨ξ(t)ξ(t′)⟩ = 2γ k_B T δ(t − t′)`, started from x = 0 with a thermalized
+ * velocity, at `t = resolution · τ_p` where `τ_p = m/γ`.
+ *
+ * The second moments are integrated from the LANGEVIN model's own moment
+ * equations, `d⟨x²⟩/dt = 2⟨xv⟩`, `d⟨xv⟩/dt = ⟨v²⟩ − (γ/m)⟨xv⟩`, with
+ * `⟨v²⟩ = k_B T/m` held at equipartition, by RK4 on a fine fixed grid. The
+ * diffusion coefficient `D = k_B T/γ` (Einstein) enters only the DENOMINATOR:
+ * the ratio tends to 1 in the diffusive regime `t ≫ τ_p`.
+ *
+ * @internal
+ */
+export declare function langevinMsdRatio(resolution: number, f: LangevinFixture): number;
+/**
+ * The slow decay rate of a Fourier mode of the TELEGRAPH equation
+ * `τ u_tt + u_t = D u_xx`, divided by the diffusion rate `D q²`:
+ * `(1 − √(1 − 4ε)) / (2ε)` with `ε = τ D q²`. Tends to 1 as τ → 0; complex
+ * (oscillatory) for `ε > 1/4`, where it returns NaN.
+ *
+ * @internal
+ */
+export declare function telegraphSlowRateRatio(tau: number, D: number, q: number): number;
+/**
+ * The oscillation frequency of a telegraph mode in its underdamped range,
+ * divided by the undamped wave frequency `c q` with `c² = D/τ`:
+ * `√(1 − 1/(4ε))`, `ε = τ D q²`. Tends to 1 as ε grows.
+ *
+ * @internal
+ */
+export declare function telegraphWaveFrequencyRatio(tau: number, D: number, q: number): number;
+/** Parameters of the WD8 steady-state fixture. @internal */
+export interface SteadyStateFixture {
+    /** Thermal diffusivity α (= κ/(ρ c_p)). */
+    readonly alpha: number;
+    /** Rod length ℓ. */
+    readonly ell: number;
+    /** End temperatures. */
+    readonly tLeft: number;
+    readonly tRight: number;
+    /** Base time; the resolution multiplies it. */
+    readonly t0: number;
+    readonly cells: number;
+}
+/**
+ * WD8 — the largest deviation from the LINEAR steady profile (the Laplace
+ * equation's solution) of an FTCS heat solution with fixed end temperatures,
+ * started from a profile with a sin(πx/ℓ) bump, at `t = resolution · t0`.
+ * Tends to zero as t → ∞: the steady state is the heat equation's attractor.
+ * `reference` replaces the linear profile the deviation is measured against.
+ *
+ * @internal
+ */
+export declare function heatSteadyDeviation(resolution: number, f: SteadyStateFixture, reference?: (i: number, cells: number) => number): number;
+//# sourceMappingURL=numerics.d.ts.map

package/dist/atlas/diffusion/numerics.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"numerics.d.ts","sourceRoot":"","sources":["../../../src/atlas/diffusion/numerics.ts"],"names":[],"mappings":"AAAA;;;;;;;;;GASG;AAEH;;;;;GAKG;AACH,wBAAgB,eAAe,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,GAAG,MAAM,CAEvE;AASD;;;;;;;;;;;;;GAaG;AACH,wBAAgB,wBAAwB,CAAC,KAAK,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,GAAG,MAAM,CAMpF;AAED;;;;;;GAMG;AACH,wBAAgB,iBAAiB,CAAC,IAAI,EAAE,MAAM,EAAE,KAAK,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,GAAG,MAAM,CAM3F;AAED,+DAA+D;AAC/D,MAAM,WAAW,WAAW;IAC1B,8BAA8B;IAC9B,QAAQ,CAAC,KAAK,EAAE,MAAM,CAAC;IACvB,iBAAiB;IACjB,QAAQ,CAAC,GAAG,EAAE,MAAM,CAAC;IACrB,kCAAkC;IAClC,QAAQ,CAAC,EAAE,EAAE,MAAM,CAAC;IACpB,qEAAqE;IACrE,QAAQ,CAAC,EAAE,EAAE,MAAM,CAAC;IACpB,kBAAkB;IAClB,QAAQ,CAAC,IAAI,EAAE,MAAM,CAAC;IACtB;;;;;;OAMG;IACH,QAAQ,CAAC,SAAS,EAAE,MAAM,CAAC;CAC5B;AAED;;;;;GAKG;AACH,wBAAgB,cAAc,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,EAAE,EAAE,MAAM,GAAG,MAAM,CAGlF;AAED;;;;;;;;;;;;;GAaG;AACH,wBAAgB,cAAc,CAAC,KAAK,EAAE,MAAM,EAAE,CAAC,EAAE,WAAW,GAAG,MAAM,CA2BpE;AAED,6DAA6D;AAC7D,MAAM,WAAW,WAAW;IAC1B,QAAQ,CAAC,IAAI,EAAE,MAAM,CAAC;IACtB,QAAQ,CAAC,CAAC,EAAE,MAAM,CAAC;IACnB,2CAA2C;IAC3C,QAAQ,CAAC,CAAC,EAAE,MAAM,CAAC;IACnB,uDAAuD;IACvD,QAAQ,CAAC,CAAC,EAAE,MAAM,CAAC;IACnB,QAAQ,CAAC,GAAG,EAAE,MAAM,CAAC;IACrB,8DAA8D;IAC9D,QAAQ,CAAC,EAAE,EAAE,MAAM,CAAC;CACrB;AAED;;;;;;;;;;;GAWG;AACH,wBAAgB,gBAAgB,CAC9B,UAAU,EAAE,MAAM,EAClB,CAAC,EAAE,WAAW,EACd,MAAM,EAAE,CAAC,CAAC,EAAE,MAAM,EAAE,GAAG,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,IAAI,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,KAAK,MAAM,GAC7E,MAAM,CAMR;AAED;;;;;;;GAOG;AACH,wBAAgB,qBAAqB,CACnC,GAAG,EAAE,MAAM,EACX,CAAC,EAAE,WAAW,EACd,MAAM,EAAE,CAAC,CAAC,EAAE,MAAM,EAAE,GAAG,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,IAAI,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,KAAK,MAAM,EAC9E,SAAS,SAAK,EACd,KAAK,SAAO,GACX,MAAM,CAQR;AAID,yEAAyE;AACzE,MAAM,WAAW,eAAe;IAC9B,uBAAuB;IACvB,QAAQ,CAAC,CAAC,EAAE,MAAM,CAAC;IACnB,mDAAmD;IACnD,QAAQ,CAAC,KAAK,EAAE,MAAM,CAAC;IACvB,4BAA4B;IAC5B,QAAQ,CAAC,EAAE,EAAE,MAAM,CAAC;CACrB;AAED;;;;;;;;;;;;GAYG;AACH,wBAAgB,gBAAgB,CAAC,UAAU,EAAE,MAAM,EAAE,CAAC,EAAE,eAAe,GAAG,MAAM,CAmB/E;AAED;;;;;;;GAOG;AACH,wBAAgB,sBAAsB,CAAC,GAAG,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,GAAG,MAAM,CAMhF;AAED;;;;;;GAMG;AACH,wBAAgB,2BAA2B,CAAC,GAAG,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,GAAG,MAAM,CAIrF;AAED,4DAA4D;AAC5D,MAAM,WAAW,kBAAkB;IACjC,2CAA2C;IAC3C,QAAQ,CAAC,KAAK,EAAE,MAAM,CAAC;IACvB,oBAAoB;IACpB,QAAQ,CAAC,GAAG,EAAE,MAAM,CAAC;IACrB,wBAAwB;IACxB,QAAQ,CAAC,KAAK,EAAE,MAAM,CAAC;IACvB,QAAQ,CAAC,MAAM,EAAE,MAAM,CAAC;IACxB,+CAA+C;IAC/C,QAAQ,CAAC,EAAE,EAAE,MAAM,CAAC;IACpB,QAAQ,CAAC,KAAK,EAAE,MAAM,CAAC;CACxB;AAED;;;;;;;;GAQG;AACH,wBAAgB,mBAAmB,CACjC,UAAU,EAAE,MAAM,EAClB,CAAC,EAAE,kBAAkB,EACrB,SAAS,CAAC,EAAE,CAAC,CAAC,EAAE,MAAM,EAAE,KAAK,EAAE,MAAM,KAAK,MAAM,GAC/C,MAAM,CAyBR"}

package/dist/atlas/diffusion/numerics.js

@@ -0,0 +1,251 @@
+/**
+ * Atlas Phase 4, S4.4 — the numerics behind the diffusion family's witnesses.
+ *
+ * Each function evaluates one claim at a RESOLUTION, so the numeric witness
+ * runner can evaluate it twice and record how the error shrank. None of them
+ * asserts an order of convergence; the runner reports the ratio and a reviewer
+ * reads it (Phase 4 design note §2.3).
+ *
+ * @module atlas/diffusion/numerics
+ */
+/**
+ * The Gaussian solution of `∂c/∂t = D ∂²c/∂x²` from a point source of unit
+ * mass, evaluated at `x`: `(4πDt)^{-1/2} exp(−x²/(4Dt))`.
+ *
+ * @internal
+ */
+export function diffusionKernel(x, t, D) {
+    return Math.exp(-(x * x) / (4 * D * t)) / Math.sqrt(4 * Math.PI * D * t);
+}
+/** `ln C(n, k)` by summing logs — exact enough for the n ≤ 10⁵ used here. */
+function logBinomial(n, k) {
+    let s = 0;
+    for (let i = 1; i <= k; i++)
+        s += Math.log(n - k + i) - Math.log(i);
+    return s;
+}
+/**
+ * WD1 — the probability DENSITY at the origin of a symmetric lattice random
+ * walk after `steps` steps, under the closure `D = Δx²/(2Δt)` at fixed `D`
+ * and `t`.
+ *
+ * With `Δt = t/steps` and `Δx = √(2DΔt)`, the walker sits at the origin with
+ * probability `C(steps, steps/2)/2^steps`. Only every other site is reachable
+ * after an even number of steps, so the probability is spread over a cell of
+ * width `2Δx`, and the density is that probability divided by `2Δx`.
+ * As `steps → ∞` it converges to {@link diffusionKernel}`(0, t, D)`.
+ *
+ * @param steps - the resolution; rounded up to the next even integer.
+ * @internal
+ */
+export function randomWalkCentralDensity(steps, D, t) {
+    const n = 2 * Math.ceil(steps / 2);
+    const dt = t / n;
+    const dx = Math.sqrt(2 * D * dt);
+    const logP = logBinomial(n, n / 2) - n * Math.LN2;
+    return Math.exp(logP) / (2 * dx);
+}
+/**
+ * WD1b — the random walk's density at lattice site `site` after `steps` steps
+ * (same closure and grid as {@link randomWalkCentralDensity}). Exactly zero
+ * outside the light cone `|site| > steps` and on sites of the wrong parity.
+ *
+ * @internal
+ */
+export function randomWalkDensity(site, steps, D, t) {
+    const n = 2 * Math.ceil(steps / 2);
+    if (Math.abs(site) > n || (site + n) % 2 !== 0)
+        return 0;
+    const dx = Math.sqrt((2 * D * t) / n);
+    const k = (n + site) / 2;
+    return Math.exp(logBinomial(n, k) - n * Math.LN2) / (2 * dx);
+}
+/**
+ * The exact solution of the DIFFUSION equation with coefficient `D` from the
+ * WD2 initial condition: `√(s0/(s0 + D t)) · exp(−x²/(4(s0 + D t)))`.
+ *
+ * @internal
+ */
+export function gaussianSpread(x, t, D, s0) {
+    const s = s0 + D * t;
+    return Math.sqrt(s0 / s) * Math.exp(-(x * x) / (4 * s));
+}
+/**
+ * WD2 — the temperature at `x = 0, t = tEnd` from an explicit FTCS solution of
+ * the HEAT equation `ρ c_p ∂T/∂t = κ ∂²T/∂x²` on `cells` intervals.
+ *
+ * The heat equation is solved in its OWN variables: the step uses `κ`, `ρ` and
+ * `c_p`, never a pre-combined `D`. Comparing the result against
+ * {@link gaussianSpread} with `D = κ/(ρ c_p)` is what tests the dictionary.
+ * The mesh ratio `κΔt/(ρ c_p Δx²)` is held at or below 1/4, inside the FTCS
+ * stability limit of 1/2, so refinement halves `Δx` and quarters `Δt`.
+ *
+ * @param cells - the resolution: the number of spatial intervals (made even so
+ * `x = 0` is a node).
+ * @internal
+ */
+export function heatFtcsCentre(cells, f) {
+    const nx = 2 * Math.ceil(cells / 2);
+    const dx = (2 * f.halfWidth) / nx;
+    const capacity = f.rho * f.cp;
+    const dtMax = (0.25 * capacity * dx * dx) / f.kappa;
+    const steps = Math.ceil(f.tEnd / dtMax);
+    const dt = f.tEnd / steps;
+    const r = (f.kappa * dt) / (capacity * dx * dx);
+    let u = new Float64Array(nx + 1);
+    let next = new Float64Array(nx + 1);
+    // The initial condition exp(−x²/(4 s0)) involves no diffusivity at all.
+    for (let i = 0; i <= nx; i++) {
+        const x = -f.halfWidth + i * dx;
+        u[i] = Math.exp(-(x * x) / (4 * f.s0));
+    }
+    u[0] = 0;
+    u[nx] = 0;
+    for (let n = 1; n <= steps; n++) {
+        for (let i = 1; i < nx; i++) {
+            next[i] = u[i] + r * (u[i + 1] - 2 * u[i] + u[i - 1]);
+        }
+        next[0] = 0;
+        next[nx] = 0;
+        [u, next] = [next, u];
+    }
+    return u[nx / 2];
+}
+/**
+ * WD3 — the finite-difference residual `|∂_τ φ − (ħ/2m) ∂²_x φ|` of the
+ * Wick-rotated free kernel `φ(x, τ)` (`witnesses/quantum-support.ts`, W5) at
+ * one point, with step `h = h0/resolution` in both variables.
+ *
+ * The residual measures whether the Wick rotation `t = −iτ` of the free
+ * Schrödinger equation really lands on the DIFFUSION equation with
+ * `D = ħ/(2m)`. Its exact value is zero; the finite-difference estimate
+ * approaches zero as `h` shrinks, and the runner records how fast.
+ *
+ * @internal
+ */
+export function wickHeatResidual(resolution, f, kernel) {
+    const h = f.h0 / resolution;
+    const phi = (x, tau) => kernel(x, tau, f.s, f.hbar, f.m);
+    const dTau = (phi(f.x, f.tau + h) - phi(f.x, f.tau - h)) / (2 * h);
+    const dXX = (phi(f.x + h, f.tau) - 2 * phi(f.x, f.tau) + phi(f.x - h, f.tau)) / (h * h);
+    return Math.abs(dTau - (f.hbar / (2 * f.m)) * dXX);
+}
+/**
+ * WD3b — `∫ φ(x, τ)² dx` of the Wick-rotated free kernel, by the trapezoid
+ * rule over `[−halfWidth, halfWidth]`. The Schrödinger norm `∫|ψ|² dx` is
+ * conserved; this one DECAYS with τ, which is the unitarity the continuation
+ * does not preserve.
+ *
+ * @internal
+ */
+export function wickKernelSquaredNorm(tau, f, kernel, halfWidth = 20, cells = 4000) {
+    const dx = (2 * halfWidth) / cells;
+    let sum = 0;
+    for (let i = 0; i <= cells; i++) {
+        const v = kernel(-halfWidth + i * dx, tau, f.s, f.hbar, f.m);
+        sum += (i === 0 || i === cells ? 0.5 : 1) * v * v;
+    }
+    return sum * dx;
+}
+/**
+ * WD4 — `⟨x²⟩(t) / (2Dt)` for the Langevin model `m v̇ = −γv + ξ`,
+ * `⟨ξ(t)ξ(t′)⟩ = 2γ k_B T δ(t − t′)`, started from x = 0 with a thermalized
+ * velocity, at `t = resolution · τ_p` where `τ_p = m/γ`.
+ *
+ * The second moments are integrated from the LANGEVIN model's own moment
+ * equations, `d⟨x²⟩/dt = 2⟨xv⟩`, `d⟨xv⟩/dt = ⟨v²⟩ − (γ/m)⟨xv⟩`, with
+ * `⟨v²⟩ = k_B T/m` held at equipartition, by RK4 on a fine fixed grid. The
+ * diffusion coefficient `D = k_B T/γ` (Einstein) enters only the DENOMINATOR:
+ * the ratio tends to 1 in the diffusive regime `t ≫ τ_p`.
+ *
+ * @internal
+ */
+export function langevinMsdRatio(resolution, f) {
+    const tauP = f.m / f.gamma;
+    const tEnd = resolution * tauP;
+    const v2 = f.kT / f.m;
+    const steps = Math.max(2000, Math.ceil(200 * resolution));
+    const h = tEnd / steps;
+    let x2 = 0;
+    let xv = 0;
+    const rhs = (xvNow) => [2 * xvNow, v2 - xvNow / tauP];
+    for (let n = 0; n < steps; n++) {
+        const [a1, b1] = rhs(xv);
+        const [a2, b2] = rhs(xv + 0.5 * h * b1);
+        const [a3, b3] = rhs(xv + 0.5 * h * b2);
+        const [a4, b4] = rhs(xv + h * b3);
+        x2 += (h / 6) * (a1 + 2 * a2 + 2 * a3 + a4);
+        xv += (h / 6) * (b1 + 2 * b2 + 2 * b3 + b4);
+    }
+    const D = f.kT / f.gamma;
+    return x2 / (2 * D * tEnd);
+}
+/**
+ * The slow decay rate of a Fourier mode of the TELEGRAPH equation
+ * `τ u_tt + u_t = D u_xx`, divided by the diffusion rate `D q²`:
+ * `(1 − √(1 − 4ε)) / (2ε)` with `ε = τ D q²`. Tends to 1 as τ → 0; complex
+ * (oscillatory) for `ε > 1/4`, where it returns NaN.
+ *
+ * @internal
+ */
+export function telegraphSlowRateRatio(tau, D, q) {
+    const eps = tau * D * q * q;
+    if (eps === 0)
+        return 1;
+    const disc = 1 - 4 * eps;
+    if (disc < 0)
+        return Number.NaN;
+    return (1 - Math.sqrt(disc)) / (2 * eps);
+}
+/**
+ * The oscillation frequency of a telegraph mode in its underdamped range,
+ * divided by the undamped wave frequency `c q` with `c² = D/τ`:
+ * `√(1 − 1/(4ε))`, `ε = τ D q²`. Tends to 1 as ε grows.
+ *
+ * @internal
+ */
+export function telegraphWaveFrequencyRatio(tau, D, q) {
+    const eps = tau * D * q * q;
+    if (eps <= 0.25)
+        return Number.NaN;
+    return Math.sqrt(1 - 1 / (4 * eps));
+}
+/**
+ * WD8 — the largest deviation from the LINEAR steady profile (the Laplace
+ * equation's solution) of an FTCS heat solution with fixed end temperatures,
+ * started from a profile with a sin(πx/ℓ) bump, at `t = resolution · t0`.
+ * Tends to zero as t → ∞: the steady state is the heat equation's attractor.
+ * `reference` replaces the linear profile the deviation is measured against.
+ *
+ * @internal
+ */
+export function heatSteadyDeviation(resolution, f, reference) {
+    const nx = f.cells;
+    const dx = f.ell / nx;
+    const dtMax = (0.25 * dx * dx) / f.alpha;
+    const tEnd = resolution * f.t0;
+    const steps = Math.ceil(tEnd / dtMax);
+    const dt = tEnd / steps;
+    const r = (f.alpha * dt) / (dx * dx);
+    const linear = (i) => f.tLeft + ((f.tRight - f.tLeft) * i) / nx;
+    let u = new Float64Array(nx + 1);
+    let next = new Float64Array(nx + 1);
+    for (let i = 0; i <= nx; i++)
+        u[i] = linear(i) + Math.sin((Math.PI * i) / nx);
+    next[0] = f.tLeft;
+    next[nx] = f.tRight;
+    for (let n = 0; n < steps; n++) {
+        for (let i = 1; i < nx; i++)
+            next[i] = u[i] + r * (u[i + 1] - 2 * u[i] + u[i - 1]);
+        [u, next] = [next, u];
+    }
+    // The deviation is measured against `reference` when given — the negative
+    // control uses a WRONG steady state here — and against the linear Laplace
+    // profile otherwise. The solver itself never sees the reference.
+    const ref = reference ?? ((i) => linear(i));
+    let worst = 0;
+    for (let i = 0; i <= nx; i++)
+        worst = Math.max(worst, Math.abs(u[i] - ref(i, nx)));
+    return worst;
+}
+//# sourceMappingURL=numerics.js.map

package/dist/atlas/diffusion/numerics.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"numerics.js","sourceRoot":"","sources":["../../../src/atlas/diffusion/numerics.ts"],"names":[],"mappings":"AAAA;;;;;;;;;GASG;AAEH;;;;;GAKG;AACH,MAAM,UAAU,eAAe,CAAC,CAAS,EAAE,CAAS,EAAE,CAAS;IAC7D,OAAO,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC;AAC3E,CAAC;AAED,6EAA6E;AAC7E,SAAS,WAAW,CAAC,CAAS,EAAE,CAAS;IACvC,IAAI,CAAC,GAAG,CAAC,CAAC;IACV,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,CAAC,EAAE,CAAC,EAAE;QAAE,CAAC,IAAI,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;IACpE,OAAO,CAAC,CAAC;AACX,CAAC;AAED;;;;;;;;;;;;;GAaG;AACH,MAAM,UAAU,wBAAwB,CAAC,KAAa,EAAE,CAAS,EAAE,CAAS;IAC1E,MAAM,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,KAAK,GAAG,CAAC,CAAC,CAAC;IACnC,MAAM,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC;IACjB,MAAM,EAAE,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,GAAG,EAAE,CAAC,CAAC;IACjC,MAAM,IAAI,GAAG,WAAW,CAAC,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC;IAClD,OAAO,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC;AACnC,CAAC;AAED;;;;;;GAMG;AACH,MAAM,UAAU,iBAAiB,CAAC,IAAY,EAAE,KAAa,EAAE,CAAS,EAAE,CAAS;IACjF,MAAM,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,KAAK,GAAG,CAAC,CAAC,CAAC;IACnC,IAAI,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,IAAI,GAAG,CAAC,CAAC,GAAG,CAAC,KAAK,CAAC;QAAE,OAAO,CAAC,CAAC;IACzD,MAAM,EAAE,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;IACtC,MAAM,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC;IACzB,OAAO,IAAI,CAAC,GAAG,CAAC,WAAW,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC;AAC/D,CAAC;AAwBD;;;;;GAKG;AACH,MAAM,UAAU,cAAc,CAAC,CAAS,EAAE,CAAS,EAAE,CAAS,EAAE,EAAU;IACxE,MAAM,CAAC,GAAG,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC;IACrB,OAAO,IAAI,CAAC,IAAI,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;AAC1D,CAAC;AAED;;;;;;;;;;;;;GAaG;AACH,MAAM,UAAU,cAAc,CAAC,KAAa,EAAE,CAAc;IAC1D,MAAM,EAAE,GAAG,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,KAAK,GAAG,CAAC,CAAC,CAAC;IACpC,MAAM,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,SAAS,CAAC,GAAG,EAAE,CAAC;IAClC,MAAM,QAAQ,GAAG,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,EAAE,CAAC;IAC9B,MAAM,KAAK,GAAG,CAAC,IAAI,GAAG,QAAQ,GAAG,EAAE,GAAG,EAAE,CAAC,GAAG,CAAC,CAAC,KAAK,CAAC;IACpD,MAAM,KAAK,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,IAAI,GAAG,KAAK,CAAC,CAAC;IACxC,MAAM,EAAE,GAAG,CAAC,CAAC,IAAI,GAAG,KAAK,CAAC;IAC1B,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,KAAK,GAAG,EAAE,CAAC,GAAG,CAAC,QAAQ,GAAG,EAAE,GAAG,EAAE,CAAC,CAAC;IAEhD,IAAI,CAAC,GAAG,IAAI,YAAY,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC;IACjC,IAAI,IAAI,GAAG,IAAI,YAAY,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC;IACpC,wEAAwE;IACxE,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,EAAE,EAAE,CAAC,EAAE,EAAE,CAAC;QAC7B,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,SAAS,GAAG,CAAC,GAAG,EAAE,CAAC;QAChC,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;IACzC,CAAC;IACD,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;IACT,CAAC,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC;IACV,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,KAAK,EAAE,CAAC,EAAE,EAAE,CAAC;QAChC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,EAAE,EAAE,CAAC,EAAE,EAAE,CAAC;YAC5B,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAE,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAE,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAE,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAE,CAAC,CAAC;QAC5D,CAAC;QACD,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;QACZ,IAAI,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC;QACb,CAAC,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,IAAI,EAAE,CAAC,CAAC,CAAC;IACxB,CAAC;IACD,OAAO,CAAC,CAAC,EAAE,GAAG,CAAC,CAAE,CAAC;AACpB,CAAC;AAeD;;;;;;;;;;;GAWG;AACH,MAAM,UAAU,gBAAgB,CAC9B,UAAkB,EAClB,CAAc,EACd,MAA8E;IAE9E,MAAM,CAAC,GAAG,CAAC,CAAC,EAAE,GAAG,UAAU,CAAC;IAC5B,MAAM,GAAG,GAAG,CAAC,CAAS,EAAE,GAAW,EAAU,EAAE,CAAC,MAAM,CAAC,CAAC,EAAE,GAAG,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,IAAI,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;IACjF,MAAM,IAAI,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;IACnE,MAAM,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;IACxF,OAAO,IAAI,CAAC,GAAG,CAAC,IAAI,GAAG,CAAC,CAAC,CAAC,IAAI,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC;AACrD,CAAC;AAED;;;;;;;GAOG;AACH,MAAM,UAAU,qBAAqB,CACnC,GAAW,EACX,CAAc,EACd,MAA8E,EAC9E,SAAS,GAAG,EAAE,EACd,KAAK,GAAG,IAAI;IAEZ,MAAM,EAAE,GAAG,CAAC,CAAC,GAAG,SAAS,CAAC,GAAG,KAAK,CAAC;IACnC,IAAI,GAAG,GAAG,CAAC,CAAC;IACZ,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,KAAK,EAAE,CAAC,EAAE,EAAE,CAAC;QAChC,MAAM,CAAC,GAAG,MAAM,CAAC,CAAC,SAAS,GAAG,CAAC,GAAG,EAAE,EAAE,GAAG,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,IAAI,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;QAC7D,GAAG,IAAI,CAAC,CAAC,KAAK,CAAC,IAAI,CAAC,KAAK,KAAK,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;IACpD,CAAC;IACD,OAAO,GAAG,GAAG,EAAE,CAAC;AAClB,CAAC;AAcD;;;;;;;;;;;;GAYG;AACH,MAAM,UAAU,gBAAgB,CAAC,UAAkB,EAAE,CAAkB;IACrE,MAAM,IAAI,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,KAAK,CAAC;IAC3B,MAAM,IAAI,GAAG,UAAU,GAAG,IAAI,CAAC;IAC/B,MAAM,EAAE,GAAG,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC;IACtB,MAAM,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,EAAE,IAAI,CAAC,IAAI,CAAC,GAAG,GAAG,UAAU,CAAC,CAAC,CAAC;IAC1D,MAAM,CAAC,GAAG,IAAI,GAAG,KAAK,CAAC;IACvB,IAAI,EAAE,GAAG,CAAC,CAAC;IACX,IAAI,EAAE,GAAG,CAAC,CAAC;IACX,MAAM,GAAG,GAAG,CAAC,KAAa,EAAoB,EAAE,CAAC,CAAC,CAAC,GAAG,KAAK,EAAE,EAAE,GAAG,KAAK,GAAG,IAAI,CAAC,CAAC;IAChF,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,KAAK,EAAE,CAAC,EAAE,EAAE,CAAC;QAC/B,MAAM,CAAC,EAAE,EAAE,EAAE,CAAC,GAAG,GAAG,CAAC,EAAE,CAAC,CAAC;QACzB,MAAM,CAAC,EAAE,EAAE,EAAE,CAAC,GAAG,GAAG,CAAC,EAAE,GAAG,GAAG,GAAG,CAAC,GAAG,EAAE,CAAC,CAAC;QACxC,MAAM,CAAC,EAAE,EAAE,EAAE,CAAC,GAAG,GAAG,CAAC,EAAE,GAAG,GAAG,GAAG,CAAC,GAAG,EAAE,CAAC,CAAC;QACxC,MAAM,CAAC,EAAE,EAAE,EAAE,CAAC,GAAG,GAAG,CAAC,EAAE,GAAG,CAAC,GAAG,EAAE,CAAC,CAAC;QAClC,EAAE,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,EAAE,GAAG,CAAC,GAAG,EAAE,GAAG,CAAC,GAAG,EAAE,GAAG,EAAE,CAAC,CAAC;QAC5C,EAAE,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,EAAE,GAAG,CAAC,GAAG,EAAE,GAAG,CAAC,GAAG,EAAE,GAAG,EAAE,CAAC,CAAC;IAC9C,CAAC;IACD,MAAM,CAAC,GAAG,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,KAAK,CAAC;IACzB,OAAO,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,CAAC;AAC7B,CAAC;AAED;;;;;;;GAOG;AACH,MAAM,UAAU,sBAAsB,CAAC,GAAW,EAAE,CAAS,EAAE,CAAS;IACtE,MAAM,GAAG,GAAG,GAAG,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;IAC5B,IAAI,GAAG,KAAK,CAAC;QAAE,OAAO,CAAC,CAAC;IACxB,MAAM,IAAI,GAAG,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC;IACzB,IAAI,IAAI,GAAG,CAAC;QAAE,OAAO,MAAM,CAAC,GAAG,CAAC;IAChC,OAAO,CAAC,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC;AAC3C,CAAC;AAED;;;;;;GAMG;AACH,MAAM,UAAU,2BAA2B,CAAC,GAAW,EAAE,CAAS,EAAE,CAAS;IAC3E,MAAM,GAAG,GAAG,GAAG,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;IAC5B,IAAI,GAAG,IAAI,IAAI;QAAE,OAAO,MAAM,CAAC,GAAG,CAAC;IACnC,OAAO,IAAI,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC;AACtC,CAAC;AAgBD;;;;;;;;GAQG;AACH,MAAM,UAAU,mBAAmB,CACjC,UAAkB,EAClB,CAAqB,EACrB,SAAgD;IAEhD,MAAM,EAAE,GAAG,CAAC,CAAC,KAAK,CAAC;IACnB,MAAM,EAAE,GAAG,CAAC,CAAC,GAAG,GAAG,EAAE,CAAC;IACtB,MAAM,KAAK,GAAG,CAAC,IAAI,GAAG,EAAE,GAAG,EAAE,CAAC,GAAG,CAAC,CAAC,KAAK,CAAC;IACzC,MAAM,IAAI,GAAG,UAAU,GAAG,CAAC,CAAC,EAAE,CAAC;IAC/B,MAAM,KAAK,GAAG,IAAI,CAAC,IAAI,CAAC,IAAI,GAAG,KAAK,CAAC,CAAC;IACtC,MAAM,EAAE,GAAG,IAAI,GAAG,KAAK,CAAC;IACxB,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,KAAK,GAAG,EAAE,CAAC,GAAG,CAAC,EAAE,GAAG,EAAE,CAAC,CAAC;IACrC,MAAM,MAAM,GAAG,CAAC,CAAS,EAAU,EAAE,CAAC,CAAC,CAAC,KAAK,GAAG,CAAC,CAAC,CAAC,CAAC,MAAM,GAAG,CAAC,CAAC,KAAK,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC;IAChF,IAAI,CAAC,GAAG,IAAI,YAAY,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC;IACjC,IAAI,IAAI,GAAG,IAAI,YAAY,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC;IACpC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,EAAE,EAAE,CAAC,EAAE;QAAE,CAAC,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,IAAI,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC;IAC9E,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,KAAK,CAAC;IAClB,IAAI,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,MAAM,CAAC;IACpB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,KAAK,EAAE,CAAC,EAAE,EAAE,CAAC;QAC/B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,EAAE,EAAE,CAAC,EAAE;YAAE,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAE,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAE,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAE,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAE,CAAC,CAAC;QACvF,CAAC,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,IAAI,EAAE,CAAC,CAAC,CAAC;IACxB,CAAC;IACD,0EAA0E;IAC1E,0EAA0E;IAC1E,iEAAiE;IACjE,MAAM,GAAG,GAAG,SAAS,IAAI,CAAC,CAAC,CAAS,EAAE,EAAE,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC;IACpD,IAAI,KAAK,GAAG,CAAC,CAAC;IACd,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,EAAE,EAAE,CAAC,EAAE;QAAE,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,EAAE,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAE,GAAG,GAAG,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC;IACpF,OAAO,KAAK,CAAC;AACf,CAAC"}

package/dist/atlas/error-algebra.d.ts

@@ -0,0 +1,52 @@
+/**
+ * The Phase 0 error algebra: composition of Lipschitz-plus-offset bounds.
+ *
+ * With `∘` read as OUTER-AFTER-INNER, a bound `(K, δ)` states
+ * `|f(x) − f̃(x̃)| ≤ K|x − x̃| + δ` in the bound's own norm. Composing an outer
+ * map onto an inner one therefore gives
+ *
+ *     (K_o, δ_o) ∘ (K_i, δ_i) = (K_o·K_i, K_o·δ_i + δ_o)
+ *
+ * which is associative and NOT commutative. `IDENTITY_BOUND` is the two-sided
+ * identity, so `(bounds, ∘, IDENTITY_BOUND)` is a monoid.
+ *
+ * An exact equivalence contributes `IDENTITY_BOUND` **in the norms these
+ * Phase 0 bridges state** — not in every norm (design note §9 YELLOW (a)).
+ *
+ * @module atlas/error-algebra
+ */
+/** The `K`/`delta` pair an `ApproximationBound` composes through. @public */
+export interface BoundPair {
+    readonly K: number;
+    readonly delta: number;
+}
+/** The two-sided identity of {@link composeBounds}: an exact, error-free map. @public */
+export declare const IDENTITY_BOUND: BoundPair;
+/**
+ * Compose two bounds, `outer` after `inner`.
+ *
+ * @public
+ */
+export declare function composeBounds(outer: BoundPair, inner: BoundPair): BoundPair;
+/** Result of {@link composeBoundPath}. @public */
+export interface ComposedPath {
+    readonly bound: BoundPair;
+    /**
+     * True when the path ended on a step with no Lipschitz constant: `bound`
+     * covers the prefix, and nothing is claimed past the final step.
+     */
+    readonly terminal: boolean;
+}
+/**
+ * Fold a path of bounds in TRAVERSAL order — `[b1, b2]` means b1 then b2, so
+ * the composite is `b2 ∘ b1`.
+ *
+ * A `null` entry is a step whose Lipschitz constant is unknown. It is tolerated
+ * only as the LAST entry, where it terminates the claim; anywhere else the
+ * composite would be unbounded, and {@link MissingLipschitzError} is thrown
+ * rather than a number invented for it.
+ *
+ * @public
+ */
+export declare function composeBoundPath(bounds: readonly (BoundPair | null)[]): ComposedPath;
+//# sourceMappingURL=error-algebra.d.ts.map

package/dist/atlas/error-algebra.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"error-algebra.d.ts","sourceRoot":"","sources":["../../src/atlas/error-algebra.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;GAgBG;AAIH,6EAA6E;AAC7E,MAAM,WAAW,SAAS;IACxB,QAAQ,CAAC,CAAC,EAAE,MAAM,CAAC;IACnB,QAAQ,CAAC,KAAK,EAAE,MAAM,CAAC;CACxB;AAED,yFAAyF;AACzF,eAAO,MAAM,cAAc,EAAE,SAA8B,CAAC;AAE5D;;;;GAIG;AACH,wBAAgB,aAAa,CAAC,KAAK,EAAE,SAAS,EAAE,KAAK,EAAE,SAAS,GAAG,SAAS,CAE3E;AAED,kDAAkD;AAClD,MAAM,WAAW,YAAY;IAC3B,QAAQ,CAAC,KAAK,EAAE,SAAS,CAAC;IAC1B;;;OAGG;IACH,QAAQ,CAAC,QAAQ,EAAE,OAAO,CAAC;CAC5B;AAED;;;;;;;;;;GAUG;AACH,wBAAgB,gBAAgB,CAAC,MAAM,EAAE,SAAS,CAAC,SAAS,GAAG,IAAI,CAAC,EAAE,GAAG,YAAY,CAgBpF"}

package/dist/atlas/error-algebra.js

@@ -0,0 +1,55 @@
+/**
+ * The Phase 0 error algebra: composition of Lipschitz-plus-offset bounds.
+ *
+ * With `∘` read as OUTER-AFTER-INNER, a bound `(K, δ)` states
+ * `|f(x) − f̃(x̃)| ≤ K|x − x̃| + δ` in the bound's own norm. Composing an outer
+ * map onto an inner one therefore gives
+ *
+ *     (K_o, δ_o) ∘ (K_i, δ_i) = (K_o·K_i, K_o·δ_i + δ_o)
+ *
+ * which is associative and NOT commutative. `IDENTITY_BOUND` is the two-sided
+ * identity, so `(bounds, ∘, IDENTITY_BOUND)` is a monoid.
+ *
+ * An exact equivalence contributes `IDENTITY_BOUND` **in the norms these
+ * Phase 0 bridges state** — not in every norm (design note §9 YELLOW (a)).
+ *
+ * @module atlas/error-algebra
+ */
+import { MissingLipschitzError } from './types.js';
+/** The two-sided identity of {@link composeBounds}: an exact, error-free map. @public */
+export const IDENTITY_BOUND = { K: 1, delta: 0 };
+/**
+ * Compose two bounds, `outer` after `inner`.
+ *
+ * @public
+ */
+export function composeBounds(outer, inner) {
+    return { K: outer.K * inner.K, delta: outer.K * inner.delta + outer.delta };
+}
+/**
+ * Fold a path of bounds in TRAVERSAL order — `[b1, b2]` means b1 then b2, so
+ * the composite is `b2 ∘ b1`.
+ *
+ * A `null` entry is a step whose Lipschitz constant is unknown. It is tolerated
+ * only as the LAST entry, where it terminates the claim; anywhere else the
+ * composite would be unbounded, and {@link MissingLipschitzError} is thrown
+ * rather than a number invented for it.
+ *
+ * @public
+ */
+export function composeBoundPath(bounds) {
+    let acc = IDENTITY_BOUND;
+    for (let i = 0; i < bounds.length; i++) {
+        const step = bounds[i];
+        if (step === null) {
+            if (i !== bounds.length - 1) {
+                throw new MissingLipschitzError(`bound path has no Lipschitz constant at step ${i} of ${bounds.length}; ` +
+                    'an unknown constant is tolerated only as the final step');
+            }
+            return { bound: acc, terminal: true };
+        }
+        acc = composeBounds(step, acc);
+    }
+    return { bound: acc, terminal: false };
+}
+//# sourceMappingURL=error-algebra.js.map

package/dist/atlas/error-algebra.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"error-algebra.js","sourceRoot":"","sources":["../../src/atlas/error-algebra.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;GAgBG;AAEH,OAAO,EAAE,qBAAqB,EAAE,MAAM,YAAY,CAAC;AAQnD,yFAAyF;AACzF,MAAM,CAAC,MAAM,cAAc,GAAc,EAAE,CAAC,EAAE,CAAC,EAAE,KAAK,EAAE,CAAC,EAAE,CAAC;AAE5D;;;;GAIG;AACH,MAAM,UAAU,aAAa,CAAC,KAAgB,EAAE,KAAgB;IAC9D,OAAO,EAAE,CAAC,EAAE,KAAK,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,EAAE,KAAK,EAAE,KAAK,CAAC,CAAC,GAAG,KAAK,CAAC,KAAK,GAAG,KAAK,CAAC,KAAK,EAAE,CAAC;AAC9E,CAAC;AAYD;;;;;;;;;;GAUG;AACH,MAAM,UAAU,gBAAgB,CAAC,MAAqC;IACpE,IAAI,GAAG,GAAc,cAAc,CAAC;IACpC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,MAAM,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;QACvC,MAAM,IAAI,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC;QACvB,IAAI,IAAI,KAAK,IAAI,EAAE,CAAC;YAClB,IAAI,CAAC,KAAK,MAAM,CAAC,MAAM,GAAG,CAAC,EAAE,CAAC;gBAC5B,MAAM,IAAI,qBAAqB,CAC7B,gDAAgD,CAAC,OAAO,MAAM,CAAC,MAAM,IAAI;oBACvE,yDAAyD,CAC5D,CAAC;YACJ,CAAC;YACD,OAAO,EAAE,KAAK,EAAE,GAAG,EAAE,QAAQ,EAAE,IAAI,EAAE,CAAC;QACxC,CAAC;QACD,GAAG,GAAG,aAAa,CAAC,IAAI,EAAE,GAAG,CAAC,CAAC;IACjC,CAAC;IACD,OAAO,EAAE,KAAK,EAAE,GAAG,EAAE,QAAQ,EAAE,KAAK,EAAE,CAAC;AACzC,CAAC"}

package/dist/atlas/export.d.ts

@@ -0,0 +1,48 @@
+/**
+ * Atlas Phase 6, S6.4 — the versioned export: one combined JSON artifact for
+ * every family, and a JSON-LD projection with stable ids, QUDT quantity kinds
+ * where they genuinely resolve, and PROV-O-shaped provenance.
+ *
+ * Design note: `docs/planning/Atlas-Phase-6-Design.md` §2.
+ *
+ * Pure: no I/O. The QUDT resolution table (`data/atlas/qudt-resolution.json`) is
+ * passed IN by the emitter, so this module cannot drift from the table it was
+ * given and never reads the tree.
+ *
+ * ## Why the QUDT table is keyed by (model, parameter), never by name
+ *
+ * `kappa` is a SPRING CONSTANT in `model-chain` and a THERMAL CONDUCTIVITY in
+ * `model-heat`. A name-keyed table would assign one of them the other's
+ * quantity kind — a notation collision in the atlas's own data.
+ *
+ * ## Schema version
+ *
+ * The combined artifact keeps `schemaVersion: '0'`. Phase 4 added one OPTIONAL
+ * field (`formalRef`); an additive optional field is not a breaking change, and
+ * the plan bumps to `v1` only on one.
+ *
+ * @module atlas/export
+ * @internal
+ */
+import type { AtlasFamily } from './oscillators/index.js';
+/** The QUDT resolution table: `"<model-id>/<parameter>"` → IRI, or null when none resolves. @internal */
+export interface QudtResolution {
+    /** When and how the IRIs were checked. */
+    readonly probedAt: string;
+    readonly method: string;
+    readonly entries: Readonly<Record<string, string | null>>;
+}
+/** The stable id namespace. A URN, so no id depends on a host that could move. @internal */
+export declare const ATLAS_ID_PREFIX = "urn:upt:atlas:";
+/** The combined artifact: every family's record under one version stamp. @internal */
+export declare function toCombinedAtlasJson(families: readonly AtlasFamily[], packageVersion: string): unknown;
+/**
+ * The JSON-LD projection.
+ *
+ * @throws Error when a model parameter has NO entry in the resolution table — an
+ * unlisted parameter is a gap in the table, and exporting it silently blank
+ * would be indistinguishable from a deliberate "no QUDT kind resolves".
+ * @internal
+ */
+export declare function toAtlasJsonLd(families: readonly AtlasFamily[], packageVersion: string, qudt: QudtResolution): unknown;
+//# sourceMappingURL=export.d.ts.map

package/dist/atlas/export.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"export.d.ts","sourceRoot":"","sources":["../../src/atlas/export.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;;;;;;;GAyBG;AAEH,OAAO,KAAK,EAAE,WAAW,EAAE,MAAM,wBAAwB,CAAC;AAG1D,yGAAyG;AACzG,MAAM,WAAW,cAAc;IAC7B,0CAA0C;IAC1C,QAAQ,CAAC,QAAQ,EAAE,MAAM,CAAC;IAC1B,QAAQ,CAAC,MAAM,EAAE,MAAM,CAAC;IACxB,QAAQ,CAAC,OAAO,EAAE,QAAQ,CAAC,MAAM,CAAC,MAAM,EAAE,MAAM,GAAG,IAAI,CAAC,CAAC,CAAC;CAC3D;AAED,4FAA4F;AAC5F,eAAO,MAAM,eAAe,mBAAmB,CAAC;AAKhD,sFAAsF;AACtF,wBAAgB,mBAAmB,CAAC,QAAQ,EAAE,SAAS,WAAW,EAAE,EAAE,cAAc,EAAE,MAAM,GAAG,OAAO,CAYrG;AAED;;;;;;;GAOG;AACH,wBAAgB,aAAa,CAC3B,QAAQ,EAAE,SAAS,WAAW,EAAE,EAChC,cAAc,EAAE,MAAM,EACtB,IAAI,EAAE,cAAc,GACnB,OAAO,CAsDT"}