universal-physics-tensor 0.45.2 → 0.47.0

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"}

package/dist/atlas/export.js

@@ -0,0 +1,106 @@
+/**
+ * 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 { toAtlasJson } from './serialize.js';
+/** The stable id namespace. A URN, so no id depends on a host that could move. @internal */
+export const ATLAS_ID_PREFIX = 'urn:upt:atlas:';
+const modelIri = (id) => `${ATLAS_ID_PREFIX}model:${id}`;
+const bridgeIri = (id) => `${ATLAS_ID_PREFIX}bridge:${id}`;
+/** The combined artifact: every family's record under one version stamp. @internal */
+export function toCombinedAtlasJson(families, packageVersion) {
+    return {
+        schemaVersion: '0',
+        packageVersion,
+        families: families.map((f) => {
+            const { packageVersion: _dropped, ...record } = toAtlasJson(f, packageVersion);
+            return record;
+        }),
+    };
+}
+/**
+ * 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 function toAtlasJsonLd(families, packageVersion, qudt) {
+    const graph = [];
+    for (const family of families) {
+        for (const m of family.models) {
+            graph.push({
+                '@id': modelIri(m.id),
+                '@type': 'upt:Model',
+                'upt:family': family.family,
+                'upt:dynamics': m.dynamics,
+                'upt:parameter': m.parameters.map((p) => {
+                    const key = `${m.id}/${p.name}`;
+                    if (!(key in qudt.entries)) {
+                        throw new Error(`toAtlasJsonLd: no QUDT resolution entry for ${key}`);
+                    }
+                    const kind = qudt.entries[key];
+                    return {
+                        'upt:symbol': p.name,
+                        'upt:dimension': [p.dim.L, p.dim.M, p.dim.T, p.dim.I, p.dim.Theta, p.dim.N, p.dim.J],
+                        ...(kind === null ? {} : { 'qudt:hasQuantityKind': { '@id': kind } }),
+                    };
+                }),
+                ...(m.canonicalRefs.length === 0 ? {} : { 'upt:canonicalRef': [...m.canonicalRefs] }),
+            });
+        }
+        for (const b of family.bridges) {
+            graph.push({
+                '@id': bridgeIri(b.id),
+                '@type': ['upt:Bridge', 'prov:Entity'],
+                'upt:relation': b.relation,
+                'upt:premise': b.premises.map((p) => ({ '@id': modelIri(p) })),
+                'upt:conclusion': { '@id': modelIri(b.conclusion) },
+                // PROV-O shape: the bridge's claim is derived from its premises, and
+                // its sources are the citations it states.
+                'prov:wasDerivedFrom': b.premises.map((p) => ({ '@id': modelIri(p) })),
+                'dcterms:source': [...b.citations],
+                'upt:reviewStatus': b.reviewStatus,
+                'upt:witness': b.witnesses.map((w) => ({ 'upt:id': w.id, 'upt:kind': w.kind, 'upt:test': w.test })),
+                ...(b.formalRef === undefined
+                    ? {}
+                    : { 'upt:formalRef': { 'upt:statement': b.formalRef.statement, 'upt:fidelity': b.formalRef.fidelity } }),
+            });
+        }
+    }
+    return {
+        '@context': {
+            upt: `${ATLAS_ID_PREFIX}vocab:`,
+            prov: 'http://www.w3.org/ns/prov#',
+            qudt: 'http://qudt.org/schema/qudt/',
+            dcterms: 'http://purl.org/dc/terms/',
+        },
+        '@id': `${ATLAS_ID_PREFIX}release:${packageVersion}`,
+        'upt:qudtProbedAt': qudt.probedAt,
+        '@graph': graph,
+    };
+}
+//# sourceMappingURL=export.js.map

package/dist/atlas/export.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"export.js","sourceRoot":"","sources":["../../src/atlas/export.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;;;;;;;GAyBG;AAGH,OAAO,EAAE,WAAW,EAAE,MAAM,gBAAgB,CAAC;AAU7C,4FAA4F;AAC5F,MAAM,CAAC,MAAM,eAAe,GAAG,gBAAgB,CAAC;AAEhD,MAAM,QAAQ,GAAG,CAAC,EAAU,EAAU,EAAE,CAAC,GAAG,eAAe,SAAS,EAAE,EAAE,CAAC;AACzE,MAAM,SAAS,GAAG,CAAC,EAAU,EAAU,EAAE,CAAC,GAAG,eAAe,UAAU,EAAE,EAAE,CAAC;AAE3E,sFAAsF;AACtF,MAAM,UAAU,mBAAmB,CAAC,QAAgC,EAAE,cAAsB;IAC1F,OAAO;QACL,aAAa,EAAE,GAAG;QAClB,cAAc;QACd,QAAQ,EAAE,QAAQ,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE;YAC3B,MAAM,EAAE,cAAc,EAAE,QAAQ,EAAE,GAAG,MAAM,EAAE,GAAG,WAAW,CAAC,CAAC,EAAE,cAAc,CAG5E,CAAC;YACF,OAAO,MAAM,CAAC;QAChB,CAAC,CAAC;KACH,CAAC;AACJ,CAAC;AAED;;;;;;;GAOG;AACH,MAAM,UAAU,aAAa,CAC3B,QAAgC,EAChC,cAAsB,EACtB,IAAoB;IAEpB,MAAM,KAAK,GAAc,EAAE,CAAC;IAC5B,KAAK,MAAM,MAAM,IAAI,QAAQ,EAAE,CAAC;QAC9B,KAAK,MAAM,CAAC,IAAI,MAAM,CAAC,MAAM,EAAE,CAAC;YAC9B,KAAK,CAAC,IAAI,CAAC;gBACT,KAAK,EAAE,QAAQ,CAAC,CAAC,CAAC,EAAE,CAAC;gBACrB,OAAO,EAAE,WAAW;gBACpB,YAAY,EAAE,MAAM,CAAC,MAAM;gBAC3B,cAAc,EAAE,CAAC,CAAC,QAAQ;gBAC1B,eAAe,EAAE,CAAC,CAAC,UAAU,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE;oBACtC,MAAM,GAAG,GAAG,GAAG,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,IAAI,EAAE,CAAC;oBAChC,IAAI,CAAC,CAAC,GAAG,IAAI,IAAI,CAAC,OAAO,CAAC,EAAE,CAAC;wBAC3B,MAAM,IAAI,KAAK,CAAC,+CAA+C,GAAG,EAAE,CAAC,CAAC;oBACxE,CAAC;oBACD,MAAM,IAAI,GAAG,IAAI,CAAC,OAAO,CAAC,GAAG,CAAC,CAAC;oBAC/B,OAAO;wBACL,YAAY,EAAE,CAAC,CAAC,IAAI;wBACpB,eAAe,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,KAAK,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;wBACpF,GAAG,CAAC,IAAI,KAAK,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,sBAAsB,EAAE,EAAE,KAAK,EAAE,IAAI,EAAE,EAAE,CAAC;qBACtE,CAAC;gBACJ,CAAC,CAAC;gBACF,GAAG,CAAC,CAAC,CAAC,aAAa,CAAC,MAAM,KAAK,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,kBAAkB,EAAE,CAAC,GAAG,CAAC,CAAC,aAAa,CAAC,EAAE,CAAC;aACtF,CAAC,CAAC;QACL,CAAC;QACD,KAAK,MAAM,CAAC,IAAI,MAAM,CAAC,OAAO,EAAE,CAAC;YAC/B,KAAK,CAAC,IAAI,CAAC;gBACT,KAAK,EAAE,SAAS,CAAC,CAAC,CAAC,EAAE,CAAC;gBACtB,OAAO,EAAE,CAAC,YAAY,EAAE,aAAa,CAAC;gBACtC,cAAc,EAAE,CAAC,CAAC,QAAQ;gBAC1B,aAAa,EAAE,CAAC,CAAC,QAAQ,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,EAAE,KAAK,EAAE,QAAQ,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC;gBAC9D,gBAAgB,EAAE,EAAE,KAAK,EAAE,QAAQ,CAAC,CAAC,CAAC,UAAU,CAAC,EAAE;gBACnD,qEAAqE;gBACrE,2CAA2C;gBAC3C,qBAAqB,EAAE,CAAC,CAAC,QAAQ,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,EAAE,KAAK,EAAE,QAAQ,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC;gBACtE,gBAAgB,EAAE,CAAC,GAAG,CAAC,CAAC,SAAS,CAAC;gBAClC,kBAAkB,EAAE,CAAC,CAAC,YAAY;gBAClC,aAAa,EAAE,CAAC,CAAC,SAAS,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,EAAE,QAAQ,EAAE,CAAC,CAAC,EAAE,EAAE,UAAU,EAAE,CAAC,CAAC,IAAI,EAAE,UAAU,EAAE,CAAC,CAAC,IAAI,EAAE,CAAC,CAAC;gBACnG,GAAG,CAAC,CAAC,CAAC,SAAS,KAAK,SAAS;oBAC3B,CAAC,CAAC,EAAE;oBACJ,CAAC,CAAC,EAAE,eAAe,EAAE,EAAE,eAAe,EAAE,CAAC,CAAC,SAAS,CAAC,SAAS,EAAE,cAAc,EAAE,CAAC,CAAC,SAAS,CAAC,QAAQ,EAAE,EAAE,CAAC;aAC3G,CAAC,CAAC;QACL,CAAC;IACH,CAAC;IACD,OAAO;QACL,UAAU,EAAE;YACV,GAAG,EAAE,GAAG,eAAe,QAAQ;YAC/B,IAAI,EAAE,4BAA4B;YAClC,IAAI,EAAE,8BAA8B;YACpC,OAAO,EAAE,2BAA2B;SACrC;QACD,KAAK,EAAE,GAAG,eAAe,WAAW,cAAc,EAAE;QACpD,kBAAkB,EAAE,IAAI,CAAC,QAAQ;QACjC,QAAQ,EAAE,KAAK;KAChB,CAAC;AACJ,CAAC"}

package/dist/atlas/families.d.ts

@@ -0,0 +1,14 @@
+/**
+ * Every registered atlas family, in the order they were built.
+ *
+ * Gates that must cover the WHOLE atlas (the evidence rule, admission, the
+ * JSON artifacts) iterate this list rather than naming one family, so a new
+ * family is checked the moment it is registered here. A gate that named
+ * `OSCILLATOR_FAMILY` directly would pass forever over a family it never read.
+ *
+ * @module atlas/families
+ */
+import type { AtlasFamily } from './oscillators/index.js';
+/** All atlas families. @internal */
+export declare const ATLAS_FAMILIES: readonly AtlasFamily[];
+//# sourceMappingURL=families.d.ts.map

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

@@ -0,0 +1 @@
+{"version":3,"file":"families.d.ts","sourceRoot":"","sources":["../../src/atlas/families.ts"],"names":[],"mappings":"AAAA;;;;;;;;;GASG;AAKH,OAAO,KAAK,EAAE,WAAW,EAAE,MAAM,wBAAwB,CAAC;AAE1D,oCAAoC;AACpC,eAAO,MAAM,cAAc,EAAE,SAAS,WAAW,EAIhD,CAAC"}

package/dist/atlas/families.js

@@ -0,0 +1,20 @@
+/**
+ * Every registered atlas family, in the order they were built.
+ *
+ * Gates that must cover the WHOLE atlas (the evidence rule, admission, the
+ * JSON artifacts) iterate this list rather than naming one family, so a new
+ * family is checked the moment it is registered here. A gate that named
+ * `OSCILLATOR_FAMILY` directly would pass forever over a family it never read.
+ *
+ * @module atlas/families
+ */
+import { DIFFUSION_FAMILY } from './diffusion/index.js';
+import { OSCILLATOR_FAMILY } from './oscillators/index.js';
+import { WAVES_FAMILY } from './waves/index.js';
+/** All atlas families. @internal */
+export const ATLAS_FAMILIES = [
+    OSCILLATOR_FAMILY,
+    DIFFUSION_FAMILY,
+    WAVES_FAMILY,
+];
+//# sourceMappingURL=families.js.map

package/dist/atlas/families.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"families.js","sourceRoot":"","sources":["../../src/atlas/families.ts"],"names":[],"mappings":"AAAA;;;;;;;;;GASG;AAEH,OAAO,EAAE,gBAAgB,EAAE,MAAM,sBAAsB,CAAC;AACxD,OAAO,EAAE,iBAAiB,EAAE,MAAM,wBAAwB,CAAC;AAC3D,OAAO,EAAE,YAAY,EAAE,MAAM,kBAAkB,CAAC;AAGhD,oCAAoC;AACpC,MAAM,CAAC,MAAM,cAAc,GAA2B;IACpD,iBAAiB;IACjB,gBAAgB;IAChB,YAAY;CACb,CAAC"}