universal-physics-tensor 0.45.2 → 0.46.0

package/dist/atlas/waves/numerics.js

@@ -0,0 +1,218 @@
+/**
+ * Atlas Phase 4, S4.5 — the numerics behind the wave family's witnesses.
+ *
+ * Each function evaluates one claim at a RESOLUTION, so the numeric witness
+ * runner can record how the error shrinks under refinement. No function
+ * asserts an order; the runner reports the ratio (Phase 4 design note §2.3).
+ *
+ * @module atlas/waves/numerics
+ */
+/**
+ * WS1 — the displacement at the midpoint `x = ½, t = tEnd` of a string of unit
+ * length with fixed ends, from a leapfrog solution of `μ y_tt = F y_xx` written
+ * in the string's OWN variables, starting from `y = sin(πx)` at rest.
+ *
+ * The comparison target, `sin(πx) cos(π c t)` with `c = √(F/μ)`, is the 1-D
+ * wave equation's solution, so the pair tests the restriction's dictionary.
+ * The Courant number `cΔt/Δx` is held at 1/2.
+ *
+ * @param cells - spatial intervals (made even so the midpoint is a node).
+ * @internal
+ */
+export function stringLeapfrogMidpoint(cells, f) {
+    const nx = 2 * Math.ceil(cells / 2);
+    const dx = 1 / nx;
+    // Courant 1/2 from the string's own coefficients: Δt = ½ Δx √(μ/F).
+    const dtMax = 0.5 * dx * Math.sqrt(f.mu / f.tension);
+    const steps = Math.ceil(f.tEnd / dtMax);
+    const dt = f.tEnd / steps;
+    const r2 = (f.tension * dt * dt) / (f.mu * dx * dx);
+    let prev = new Float64Array(nx + 1);
+    let cur = new Float64Array(nx + 1);
+    let next = new Float64Array(nx + 1);
+    for (let i = 0; i <= nx; i++)
+        prev[i] = Math.sin(Math.PI * i * dx);
+    // First step from rest: y¹ = y⁰ + ½ r² δ²y⁰ (second-order start).
+    for (let i = 1; i < nx; i++) {
+        cur[i] = prev[i] + 0.5 * r2 * (prev[i + 1] - 2 * prev[i] + prev[i - 1]);
+    }
+    for (let n = 1; n < steps; n++) {
+        for (let i = 1; i < nx; i++) {
+            next[i] = 2 * cur[i] - prev[i] + r2 * (cur[i + 1] - 2 * cur[i] + cur[i - 1]);
+        }
+        [prev, cur, next] = [cur, next, prev];
+    }
+    return cur[nx / 2];
+}
+/** The two profiles of WS2: deliberately different, so neither direction is trivial. */
+const profileRight = (s) => Math.exp(-s * s);
+const profileLeft = (s) => 0.5 * Math.exp(-((s - 1) * (s - 1)) / 0.5);
+/**
+ * The d'Alembert form `u = f(x − ct) + g(x + ct)` with the WS2 profiles.
+ *
+ * @internal
+ */
+export function dalembert(x, t, c) {
+    return profileRight(x - c * t) + profileLeft(x + c * t);
+}
+/**
+ * WS2 — the finite-difference residual `|u_tt − c² u_xx|` of the d'Alembert
+ * form at one point, step `h = h0/resolution` in both variables. Exactly zero
+ * in the limit; the runner records how fast the estimate gets there.
+ *
+ * @internal
+ */
+export function dalembertResidual(resolution, f) {
+    const h = f.h0 / resolution;
+    const u = (x, t) => dalembert(x, t, f.c);
+    const utt = (u(f.x, f.t + h) - 2 * u(f.x, f.t) + u(f.x, f.t - h)) / (h * h);
+    const uxx = (u(f.x + h, f.t) - 2 * u(f.x, f.t) + u(f.x - h, f.t)) / (h * h);
+    return Math.abs(utt - f.c * f.c * uxx);
+}
+/**
+ * The adiabatic equation of state `p(ρ) = p₀ (ρ/ρ₀)^γ` — the SECOND premise of
+ * the sound-speed hyperedge — and its derivative at `ρ`.
+ *
+ * @internal
+ */
+export function adiabaticSlope(rho, f) {
+    return (f.gamma * f.p0 * Math.pow(rho / f.rho0, f.gamma - 1)) / f.rho0;
+}
+/**
+ * WS3 — the pressure perturbation at `x = ¼, t = tEnd` from a staggered
+ * leapfrog solution of the LINEARIZED EULER equations
+ * `∂ρ′/∂t = −ρ₀ ∂v/∂x`, `ρ₀ ∂v/∂t = −∂p′/∂x`, closed by the ADIABATIC EOS
+ * `p′ = (dp/dρ)|ρ₀ ρ′`, on a periodic unit domain from `p′ = sin 2πx` at rest.
+ *
+ * Both premises enter in their own form: the Euler pair steps `ρ′` and `v`,
+ * and the only link from `ρ′` to `p′` is {@link adiabaticSlope}. The
+ * comparison target `sin(2πx) cos(2π c t)` with `c² = γ p₀/ρ₀` is the SOUND
+ * model's solution, so the pair tests the hyperedge's derivation.
+ *
+ * @param cells - grid intervals (made a multiple of 4 so `x = ¼` is a node).
+ * @internal
+ */
+export function acousticLeapfrogQuarter(cells, f) {
+    const nx = 4 * Math.ceil(cells / 4);
+    const dx = 1 / nx;
+    const slope = adiabaticSlope(f.rho0, f);
+    const cMax = Math.sqrt(slope);
+    const dtMax = (0.5 * dx) / cMax;
+    const steps = Math.ceil(f.tEnd / dtMax);
+    const dt = f.tEnd / steps;
+    // ρ′ at nodes i·Δx; v at half-nodes (i + ½)Δx, periodic.
+    const rho = new Float64Array(nx);
+    const v = new Float64Array(nx);
+    for (let i = 0; i < nx; i++)
+        rho[i] = Math.sin(2 * Math.PI * i * dx) / slope;
+    const pressure = (i) => slope * rho[(i + nx) % nx];
+    // Half step from rest: v^{1/2} = −(Δt/2ρ₀) ∂p′/∂x.
+    for (let i = 0; i < nx; i++)
+        v[i] = (-(0.5 * dt) / f.rho0) * ((pressure(i + 1) - pressure(i)) / dx);
+    for (let n = 0; n < steps; n++) {
+        for (let i = 0; i < nx; i++) {
+            const dv = (v[i] - v[(i - 1 + nx) % nx]) / dx;
+            rho[i] = rho[i] - dt * f.rho0 * dv;
+        }
+        if (n < steps - 1) {
+            for (let i = 0; i < nx; i++)
+                v[i] = v[i] - (dt / f.rho0) * ((pressure(i + 1) - pressure(i)) / dx);
+        }
+    }
+    return pressure(nx / 4);
+}
+/**
+ * Speed of sound from the adiabatic (Laplace) and isothermal (Newton) closures.
+ *
+ * @internal
+ */
+export function soundSpeeds(p0, rho0, gamma) {
+    return { adiabatic: Math.sqrt((gamma * p0) / rho0), isothermal: Math.sqrt(p0 / rho0) };
+}
+/**
+ * Relative phase-velocity error of the Klein–Gordon mode against the
+ * dispersion-free wave: `√(1 + (ω₀/(c k))²) − 1`.
+ *
+ * This is the exact quantity the approximation bound of `ab-klein-gordon-wave`
+ * is stated in; it is increasing in `ω₀/(c k)`, so its value at the edge of the
+ * declared range is the supremum over the range.
+ *
+ * @internal
+ */
+export function kleinGordonPhaseError(omega0, c, k) {
+    const ratio = omega0 / (c * k);
+    return Math.sqrt(1 + ratio * ratio) - 1;
+}
+/**
+ * WS4 — the Klein–Gordon phase velocity `√(c² + ω₀²/k²)` at wavenumber
+ * `k = k₀ · resolution`. It tends to the wave speed `c` as `k` grows: the
+ * dispersion-free limit, approached from above.
+ *
+ * @internal
+ */
+export function kleinGordonPhaseVelocity(resolution, omega0, c, k0) {
+    const k = k0 * resolution;
+    return Math.sqrt(c * c + (omega0 * omega0) / (k * k));
+}
+// ── Sprint 4 closure: Klein–Gordon limits, weak damping ────────────────────
+/**
+ * The relative error of the non-relativistic kinetic frequency: the
+ * Klein–Gordon mode's `ω − ω₀ = ω₀(√(1 + x²) − 1)` against the free
+ * Schrödinger particle's `ω₀ x²/2` under `ħ/m ↦ c²/ω₀`, with `x = c k/ω₀`.
+ * Equals `|(√(1 + x²) − 1)/(x²/2) − 1|`, increasing in x.
+ *
+ * @internal
+ */
+export function kgNonrelativisticError(x) {
+    if (x === 0)
+        return 0;
+    return Math.abs((Math.sqrt(1 + x * x) - 1) / ((x * x) / 2) - 1);
+}
+/**
+ * WS6 — `u(tEnd)` of the SPATIALLY UNIFORM Klein–Gordon field on a periodic
+ * grid, stepped by leapfrog in the full PDE `u_tt = c² u_xx − ω₀² u` from
+ * `u = 1` at rest, with `steps` time steps. The uniform mode's Laplacian
+ * vanishes, so the field should follow the oscillator solution `cos(ω₀ t)`.
+ *
+ * @internal
+ */
+export function kgUniformModeValue(steps, c, omega0, tEnd, cells = 16) {
+    const dt = tEnd / steps;
+    const dx = 1 / cells;
+    let prev = new Float64Array(cells).fill(1);
+    let cur = new Float64Array(cells);
+    let next = new Float64Array(cells);
+    const lap = (a, i) => (a[(i + 1) % cells] - 2 * a[i] + a[(i - 1 + cells) % cells]) / (dx * dx);
+    for (let i = 0; i < cells; i++) {
+        cur[i] = prev[i] + 0.5 * dt * dt * (c * c * lap(prev, i) - omega0 * omega0 * prev[i]);
+    }
+    for (let n = 1; n < steps; n++) {
+        for (let i = 0; i < cells; i++) {
+            next[i] = 2 * cur[i] - prev[i] + dt * dt * (c * c * lap(cur, i) - omega0 * omega0 * cur[i]);
+        }
+        [prev, cur, next] = [cur, next, prev];
+    }
+    return cur[0];
+}
+/**
+ * Relative phase-velocity error of a stiff string against the flexible string:
+ * `√(1 + β) − 1` with `β = EI k²/F`. Increasing in β, so its value at the
+ * domain edge is the supremum over the domain.
+ *
+ * @internal
+ */
+export function stiffStringPhaseError(F, EI, k) {
+    return Math.sqrt(1 + (EI * k * k) / F) - 1;
+}
+/**
+ * WS7 — the stiff string's phase velocity `√(F/μ + (EI/μ)k²)` at wavenumber
+ * `k = k₀/resolution`: it tends to the flexible-string speed `√(F/μ)` as the
+ * wavelength grows.
+ *
+ * @internal
+ */
+export function stiffStringPhaseVelocity(resolution, F, mu, EI, k0) {
+    const k = k0 / resolution;
+    return Math.sqrt(F / mu + (EI / mu) * k * k);
+}
+//# sourceMappingURL=numerics.js.map

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

@@ -0,0 +1 @@
+{"version":3,"file":"numerics.js","sourceRoot":"","sources":["../../../src/atlas/waves/numerics.ts"],"names":[],"mappings":"AAAA;;;;;;;;GAQG;AAYH;;;;;;;;;;;GAWG;AACH,MAAM,UAAU,sBAAsB,CAAC,KAAa,EAAE,CAAgB;IACpE,MAAM,EAAE,GAAG,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,KAAK,GAAG,CAAC,CAAC,CAAC;IACpC,MAAM,EAAE,GAAG,CAAC,GAAG,EAAE,CAAC;IAClB,oEAAoE;IACpE,MAAM,KAAK,GAAG,GAAG,GAAG,EAAE,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,OAAO,CAAC,CAAC;IACrD,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,EAAE,GAAG,CAAC,CAAC,CAAC,OAAO,GAAG,EAAE,GAAG,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC,CAAC;IAEpD,IAAI,IAAI,GAAG,IAAI,YAAY,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC;IACpC,IAAI,GAAG,GAAG,IAAI,YAAY,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC;IACnC,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,IAAI,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,EAAE,GAAG,CAAC,GAAG,EAAE,CAAC,CAAC;IACnE,kEAAkE;IAClE,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,EAAE,EAAE,CAAC,EAAE,EAAE,CAAC;QAC5B,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,CAAE,GAAG,GAAG,GAAG,EAAE,GAAG,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,CAAE,GAAG,CAAC,GAAG,IAAI,CAAC,CAAC,CAAE,GAAG,IAAI,CAAC,CAAC,GAAG,CAAC,CAAE,CAAC,CAAC;IAC9E,CAAC;IACD,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,EAAE,CAAC;YAC5B,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC,CAAC,CAAE,GAAG,IAAI,CAAC,CAAC,CAAE,GAAG,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAE,GAAG,CAAC,GAAG,GAAG,CAAC,CAAC,CAAE,GAAG,GAAG,CAAC,CAAC,GAAG,CAAC,CAAE,CAAC,CAAC;QACpF,CAAC;QACD,CAAC,IAAI,EAAE,GAAG,EAAE,IAAI,CAAC,GAAG,CAAC,GAAG,EAAE,IAAI,EAAE,IAAI,CAAC,CAAC;IACxC,CAAC;IACD,OAAO,GAAG,CAAC,EAAE,GAAG,CAAC,CAAE,CAAC;AACtB,CAAC;AAYD,wFAAwF;AACxF,MAAM,YAAY,GAAG,CAAC,CAAS,EAAU,EAAE,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;AAC7D,MAAM,WAAW,GAAG,CAAC,CAAS,EAAU,EAAE,CAAC,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC;AAEtF;;;;GAIG;AACH,MAAM,UAAU,SAAS,CAAC,CAAS,EAAE,CAAS,EAAE,CAAS;IACvD,OAAO,YAAY,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,WAAW,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC;AAC1D,CAAC;AAED;;;;;;GAMG;AACH,MAAM,UAAU,iBAAiB,CAAC,UAAkB,EAAE,CAAmB;IACvE,MAAM,CAAC,GAAG,CAAC,CAAC,EAAE,GAAG,UAAU,CAAC;IAC5B,MAAM,CAAC,GAAG,CAAC,CAAS,EAAE,CAAS,EAAU,EAAE,CAAC,SAAS,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;IACjE,MAAM,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;IAC5E,MAAM,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;IAC5E,OAAO,IAAI,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC;AACzC,CAAC;AAaD;;;;;GAKG;AACH,MAAM,UAAU,cAAc,CAAC,GAAW,EAAE,CAAkB;IAC5D,OAAO,CAAC,CAAC,CAAC,KAAK,GAAG,CAAC,CAAC,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC,CAAC,IAAI,EAAE,CAAC,CAAC,KAAK,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,IAAI,CAAC;AACzE,CAAC;AAED;;;;;;;;;;;;;GAaG;AACH,MAAM,UAAU,uBAAuB,CAAC,KAAa,EAAE,CAAkB;IACvE,MAAM,EAAE,GAAG,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,KAAK,GAAG,CAAC,CAAC,CAAC;IACpC,MAAM,EAAE,GAAG,CAAC,GAAG,EAAE,CAAC;IAClB,MAAM,KAAK,GAAG,cAAc,CAAC,CAAC,CAAC,IAAI,EAAE,CAAC,CAAC,CAAC;IACxC,MAAM,IAAI,GAAG,IAAI,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC;IAC9B,MAAM,KAAK,GAAG,CAAC,GAAG,GAAG,EAAE,CAAC,GAAG,IAAI,CAAC;IAChC,MAAM,KAAK,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,IAAI,GAAG,KAAK,CAAC,CAAC;IACxC,MAAM,EAAE,GAAG,CAAC,CAAC,IAAI,GAAG,KAAK,CAAC;IAE1B,yDAAyD;IACzD,MAAM,GAAG,GAAG,IAAI,YAAY,CAAC,EAAE,CAAC,CAAC;IACjC,MAAM,CAAC,GAAG,IAAI,YAAY,CAAC,EAAE,CAAC,CAAC;IAC/B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,EAAE,EAAE,CAAC,EAAE;QAAE,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,GAAG,CAAC,GAAG,EAAE,CAAC,GAAG,KAAK,CAAC;IAC7E,MAAM,QAAQ,GAAG,CAAC,CAAS,EAAU,EAAE,CAAC,KAAK,GAAG,GAAG,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,GAAG,EAAE,CAAE,CAAC;IACpE,mDAAmD;IACnD,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,EAAE,EAAE,CAAC,EAAE;QAAE,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,GAAG,EAAE,CAAC,GAAG,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,QAAQ,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,QAAQ,CAAC,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC;IACpG,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,EAAE,CAAC;YAC5B,MAAM,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,CAAE,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,EAAE,CAAC,GAAG,EAAE,CAAE,CAAC,GAAG,EAAE,CAAC;YAChD,GAAG,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,CAAE,GAAG,EAAE,GAAG,CAAC,CAAC,IAAI,GAAG,EAAE,CAAC;QACtC,CAAC;QACD,IAAI,CAAC,GAAG,KAAK,GAAG,CAAC,EAAE,CAAC;YAClB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,EAAE,EAAE,CAAC,EAAE;gBAAE,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAE,GAAG,CAAC,EAAE,GAAG,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,QAAQ,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,QAAQ,CAAC,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC;QACrG,CAAC;IACH,CAAC;IACD,OAAO,QAAQ,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC;AAC1B,CAAC;AAED;;;;GAIG;AACH,MAAM,UAAU,WAAW,CAAC,EAAU,EAAE,IAAY,EAAE,KAAa;IAIjE,OAAO,EAAE,SAAS,EAAE,IAAI,CAAC,IAAI,CAAC,CAAC,KAAK,GAAG,EAAE,CAAC,GAAG,IAAI,CAAC,EAAE,UAAU,EAAE,IAAI,CAAC,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,EAAE,CAAC;AACzF,CAAC;AAED;;;;;;;;;GASG;AACH,MAAM,UAAU,qBAAqB,CAAC,MAAc,EAAE,CAAS,EAAE,CAAS;IACxE,MAAM,KAAK,GAAG,MAAM,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;IAC/B,OAAO,IAAI,CAAC,IAAI,CAAC,CAAC,GAAG,KAAK,GAAG,KAAK,CAAC,GAAG,CAAC,CAAC;AAC1C,CAAC;AAED;;;;;;GAMG;AACH,MAAM,UAAU,wBAAwB,CACtC,UAAkB,EAClB,MAAc,EACd,CAAS,EACT,EAAU;IAEV,MAAM,CAAC,GAAG,EAAE,GAAG,UAAU,CAAC;IAC1B,OAAO,IAAI,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,MAAM,GAAG,MAAM,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;AACxD,CAAC;AAED,8EAA8E;AAE9E;;;;;;;GAOG;AACH,MAAM,UAAU,sBAAsB,CAAC,CAAS;IAC9C,IAAI,CAAC,KAAK,CAAC;QAAE,OAAO,CAAC,CAAC;IACtB,OAAO,IAAI,CAAC,GAAG,CAAC,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;AAClE,CAAC;AAED;;;;;;;GAOG;AACH,MAAM,UAAU,kBAAkB,CAChC,KAAa,EACb,CAAS,EACT,MAAc,EACd,IAAY,EACZ,KAAK,GAAG,EAAE;IAEV,MAAM,EAAE,GAAG,IAAI,GAAG,KAAK,CAAC;IACxB,MAAM,EAAE,GAAG,CAAC,GAAG,KAAK,CAAC;IACrB,IAAI,IAAI,GAAG,IAAI,YAAY,CAAC,KAAK,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;IAC3C,IAAI,GAAG,GAAG,IAAI,YAAY,CAAC,KAAK,CAAC,CAAC;IAClC,IAAI,IAAI,GAAG,IAAI,YAAY,CAAC,KAAK,CAAC,CAAC;IACnC,MAAM,GAAG,GAAG,CAAC,CAAe,EAAE,CAAS,EAAU,EAAE,CACjD,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,KAAK,CAAE,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAE,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,KAAK,CAAC,GAAG,KAAK,CAAE,CAAC,GAAG,CAAC,EAAE,GAAG,EAAE,CAAC,CAAC;IAC9E,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,KAAK,EAAE,CAAC,EAAE,EAAE,CAAC;QAC/B,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,CAAE,GAAG,GAAG,GAAG,EAAE,GAAG,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC,IAAI,EAAE,CAAC,CAAC,GAAG,MAAM,GAAG,MAAM,GAAG,IAAI,CAAC,CAAC,CAAE,CAAC,CAAC;IAC1F,CAAC;IACD,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,KAAK,EAAE,CAAC,EAAE,EAAE,CAAC;QAC/B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,KAAK,EAAE,CAAC,EAAE,EAAE,CAAC;YAC/B,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC,CAAC,CAAE,GAAG,IAAI,CAAC,CAAC,CAAE,GAAG,EAAE,GAAG,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC,GAAG,EAAE,CAAC,CAAC,GAAG,MAAM,GAAG,MAAM,GAAG,GAAG,CAAC,CAAC,CAAE,CAAC,CAAC;QACjG,CAAC;QACD,CAAC,IAAI,EAAE,GAAG,EAAE,IAAI,CAAC,GAAG,CAAC,GAAG,EAAE,IAAI,EAAE,IAAI,CAAC,CAAC;IACxC,CAAC;IACD,OAAO,GAAG,CAAC,CAAC,CAAE,CAAC;AACjB,CAAC;AAED;;;;;;GAMG;AACH,MAAM,UAAU,qBAAqB,CAAC,CAAS,EAAE,EAAU,EAAE,CAAS;IACpE,OAAO,IAAI,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC;AAC7C,CAAC;AAED;;;;;;GAMG;AACH,MAAM,UAAU,wBAAwB,CACtC,UAAkB,EAClB,CAAS,EACT,EAAU,EACV,EAAU,EACV,EAAU;IAEV,MAAM,CAAC,GAAG,EAAE,GAAG,UAAU,CAAC;IAC1B,OAAO,IAAI,CAAC,IAAI,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC,EAAE,GAAG,EAAE,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC;AAC/C,CAAC"}

package/dist/atlas/witness-artifact.d.ts

@@ -0,0 +1,65 @@
+/**
+ * Atlas Phase 4, S4.3 — the witness-results ARTIFACT: its shape, the run that
+ * produces it, and the one reader that turns it into passing witness ids.
+ *
+ * Design note: `docs/planning/Atlas-Phase-4-Design.md` §3.
+ *
+ * **The load-bearing invariant: no test writes the artifact.** The Lead runs
+ * `scripts/emit-witness-results.mjs`, which calls {@link runWitnessRegistry}
+ * and writes `data/atlas/witness-results.json`. `tests/atlas/witness-results.test.ts`
+ * deep-equals the committed file against a fresh in-process run — the
+ * `atlas-json` pattern. A test that wrote the artifact would certify whatever
+ * the last run happened to produce, which is the opposite of a check.
+ *
+ * **The artifact is deterministic by construction.** `elapsedMs` is dropped:
+ * wall-clock time differs on every run, and a pin that deep-equals it could
+ * never pass twice. Everything kept is a function of the specs and the peer.
+ *
+ * @module atlas/witness-artifact
+ * @internal
+ */
+import type { Convergence } from './witness-numeric.js';
+import type { UnresolvedReason, WitnessStatus } from './witness-result.js';
+import type { SymbolicSimplifier } from './witness-symbolic.js';
+import type { RegisteredWitness } from './witness-specs.js';
+/** One row of the artifact. @internal */
+export interface WitnessResultRecord {
+    readonly recordId: string;
+    readonly witnessId: string;
+    readonly kind: 'symbolic' | 'numeric';
+    readonly status: WitnessStatus;
+    /** Present iff `status === 'unresolved'`. */
+    readonly reason?: UnresolvedReason;
+    readonly detail: string;
+    /** Numeric witnesses only, and only when a number was produced. */
+    readonly convergence?: Convergence;
+}
+/** The committed artifact, `data/atlas/witness-results.json`. @internal */
+export interface WitnessResultsArtifact {
+    readonly schemaVersion: '0';
+    readonly results: readonly WitnessResultRecord[];
+}
+/**
+ * Run every registered witness and assemble the artifact.
+ *
+ * @param registry - the witnesses to run; defaults to {@link WITNESS_REGISTRY}.
+ * @param simplifier - forwarded to `runSymbolicWitness`. Omit it for the real
+ * peer (resolved through `isSimplifierAvailable()`); pass `null` to run every
+ * symbolic witness as peer-absent.
+ * @returns the artifact; never throws, because neither runner does.
+ * @internal
+ */
+export declare function runWitnessRegistry(registry?: readonly RegisteredWitness[], simplifier?: SymbolicSimplifier | null): Promise<WitnessResultsArtifact>;
+/**
+ * The ids of witnesses of `recordId` that the artifact records as `'checked'`
+ * — the value to pass as `deriveEvidence`'s `passingWitnessIds` for the
+ * artifact-governed kinds.
+ *
+ * Only `'checked'` counts. An `'unresolved'` row (peer absent, timeout, an
+ * irreducible difference) is not a pass, and admitting it would hand
+ * `symbolically-checked` to a claim nobody examined.
+ *
+ * @internal
+ */
+export declare function artifactPassingWitnessIds(artifact: WitnessResultsArtifact, recordId: string): ReadonlySet<string>;
+//# sourceMappingURL=witness-artifact.d.ts.map

package/dist/atlas/witness-artifact.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"witness-artifact.d.ts","sourceRoot":"","sources":["../../src/atlas/witness-artifact.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;GAmBG;AAEH,OAAO,KAAK,EAAE,WAAW,EAAE,MAAM,sBAAsB,CAAC;AAExD,OAAO,KAAK,EAAE,gBAAgB,EAAoB,aAAa,EAAE,MAAM,qBAAqB,CAAC;AAC7F,OAAO,KAAK,EAAE,kBAAkB,EAAE,MAAM,uBAAuB,CAAC;AAEhE,OAAO,KAAK,EAAE,iBAAiB,EAAE,MAAM,oBAAoB,CAAC;AAG5D,yCAAyC;AACzC,MAAM,WAAW,mBAAmB;IAClC,QAAQ,CAAC,QAAQ,EAAE,MAAM,CAAC;IAC1B,QAAQ,CAAC,SAAS,EAAE,MAAM,CAAC;IAC3B,QAAQ,CAAC,IAAI,EAAE,UAAU,GAAG,SAAS,CAAC;IACtC,QAAQ,CAAC,MAAM,EAAE,aAAa,CAAC;IAC/B,6CAA6C;IAC7C,QAAQ,CAAC,MAAM,CAAC,EAAE,gBAAgB,CAAC;IACnC,QAAQ,CAAC,MAAM,EAAE,MAAM,CAAC;IACxB,mEAAmE;IACnE,QAAQ,CAAC,WAAW,CAAC,EAAE,WAAW,CAAC;CACpC;AAED,2EAA2E;AAC3E,MAAM,WAAW,sBAAsB;IACrC,QAAQ,CAAC,aAAa,EAAE,GAAG,CAAC;IAC5B,QAAQ,CAAC,OAAO,EAAE,SAAS,mBAAmB,EAAE,CAAC;CAClD;AAED;;;;;;;;;GASG;AACH,wBAAsB,kBAAkB,CACtC,QAAQ,GAAE,SAAS,iBAAiB,EAAqB,EACzD,UAAU,CAAC,EAAE,kBAAkB,GAAG,IAAI,GACrC,OAAO,CAAC,sBAAsB,CAAC,CA0BjC;AAED;;;;;;;;;;GAUG;AACH,wBAAgB,yBAAyB,CACvC,QAAQ,EAAE,sBAAsB,EAChC,QAAQ,EAAE,MAAM,GACf,WAAW,CAAC,MAAM,CAAC,CAMrB"}

package/dist/atlas/witness-artifact.js

@@ -0,0 +1,78 @@
+/**
+ * Atlas Phase 4, S4.3 — the witness-results ARTIFACT: its shape, the run that
+ * produces it, and the one reader that turns it into passing witness ids.
+ *
+ * Design note: `docs/planning/Atlas-Phase-4-Design.md` §3.
+ *
+ * **The load-bearing invariant: no test writes the artifact.** The Lead runs
+ * `scripts/emit-witness-results.mjs`, which calls {@link runWitnessRegistry}
+ * and writes `data/atlas/witness-results.json`. `tests/atlas/witness-results.test.ts`
+ * deep-equals the committed file against a fresh in-process run — the
+ * `atlas-json` pattern. A test that wrote the artifact would certify whatever
+ * the last run happened to produce, which is the opposite of a check.
+ *
+ * **The artifact is deterministic by construction.** `elapsedMs` is dropped:
+ * wall-clock time differs on every run, and a pin that deep-equals it could
+ * never pass twice. Everything kept is a function of the specs and the peer.
+ *
+ * @module atlas/witness-artifact
+ * @internal
+ */
+import { runNumericWitness } from './witness-numeric.js';
+import { runSymbolicWitness } from './witness-symbolic.js';
+import { WITNESS_REGISTRY } from './witness-specs.js';
+/**
+ * Run every registered witness and assemble the artifact.
+ *
+ * @param registry - the witnesses to run; defaults to {@link WITNESS_REGISTRY}.
+ * @param simplifier - forwarded to `runSymbolicWitness`. Omit it for the real
+ * peer (resolved through `isSimplifierAvailable()`); pass `null` to run every
+ * symbolic witness as peer-absent.
+ * @returns the artifact; never throws, because neither runner does.
+ * @internal
+ */
+export async function runWitnessRegistry(registry = WITNESS_REGISTRY, simplifier) {
+    const results = [];
+    for (const entry of registry) {
+        // Sequential on purpose: the order of the artifact is the order of the
+        // registry, and a symbolic run's timeout budget should not compete with
+        // its neighbours for the event loop.
+        let run;
+        let convergence;
+        if (entry.kind === 'symbolic') {
+            run = await runSymbolicWitness(entry.spec, simplifier);
+        }
+        else {
+            const numeric = runNumericWitness(entry.spec);
+            run = numeric;
+            convergence = numeric.convergence;
+        }
+        results.push({
+            recordId: entry.recordId,
+            witnessId: run.witnessId,
+            kind: run.kind,
+            status: run.status,
+            ...(run.reason === undefined ? {} : { reason: run.reason }),
+            detail: run.detail,
+            ...(convergence === undefined ? {} : { convergence }),
+        });
+    }
+    return { schemaVersion: '0', results };
+}
+/**
+ * The ids of witnesses of `recordId` that the artifact records as `'checked'`
+ * — the value to pass as `deriveEvidence`'s `passingWitnessIds` for the
+ * artifact-governed kinds.
+ *
+ * Only `'checked'` counts. An `'unresolved'` row (peer absent, timeout, an
+ * irreducible difference) is not a pass, and admitting it would hand
+ * `symbolically-checked` to a claim nobody examined.
+ *
+ * @internal
+ */
+export function artifactPassingWitnessIds(artifact, recordId) {
+    return new Set(artifact.results
+        .filter((r) => r.recordId === recordId && r.status === 'checked')
+        .map((r) => r.witnessId));
+}
+//# sourceMappingURL=witness-artifact.js.map

package/dist/atlas/witness-artifact.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"witness-artifact.js","sourceRoot":"","sources":["../../src/atlas/witness-artifact.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;GAmBG;AAGH,OAAO,EAAE,iBAAiB,EAAE,MAAM,sBAAsB,CAAC;AAGzD,OAAO,EAAE,kBAAkB,EAAE,MAAM,uBAAuB,CAAC;AAE3D,OAAO,EAAE,gBAAgB,EAAE,MAAM,oBAAoB,CAAC;AAqBtD;;;;;;;;;GASG;AACH,MAAM,CAAC,KAAK,UAAU,kBAAkB,CACtC,QAAQ,GAAiC,gBAAgB,EACzD,UAAsC;IAEtC,MAAM,OAAO,GAA0B,EAAE,CAAC;IAC1C,KAAK,MAAM,KAAK,IAAI,QAAQ,EAAE,CAAC;QAC7B,uEAAuE;QACvE,wEAAwE;QACxE,qCAAqC;QACrC,IAAI,GAAqB,CAAC;QAC1B,IAAI,WAAoC,CAAC;QACzC,IAAI,KAAK,CAAC,IAAI,KAAK,UAAU,EAAE,CAAC;YAC9B,GAAG,GAAG,MAAM,kBAAkB,CAAC,KAAK,CAAC,IAAI,EAAE,UAAU,CAAC,CAAC;QACzD,CAAC;aAAM,CAAC;YACN,MAAM,OAAO,GAAG,iBAAiB,CAAC,KAAK,CAAC,IAAI,CAAC,CAAC;YAC9C,GAAG,GAAG,OAAO,CAAC;YACd,WAAW,GAAG,OAAO,CAAC,WAAW,CAAC;QACpC,CAAC;QACD,OAAO,CAAC,IAAI,CAAC;YACX,QAAQ,EAAE,KAAK,CAAC,QAAQ;YACxB,SAAS,EAAE,GAAG,CAAC,SAAS;YACxB,IAAI,EAAE,GAAG,CAAC,IAAI;YACd,MAAM,EAAE,GAAG,CAAC,MAAM;YAClB,GAAG,CAAC,GAAG,CAAC,MAAM,KAAK,SAAS,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,MAAM,EAAE,GAAG,CAAC,MAAM,EAAE,CAAC;YAC3D,MAAM,EAAE,GAAG,CAAC,MAAM;YAClB,GAAG,CAAC,WAAW,KAAK,SAAS,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,WAAW,EAAE,CAAC;SACtD,CAAC,CAAC;IACL,CAAC;IACD,OAAO,EAAE,aAAa,EAAE,GAAG,EAAE,OAAO,EAAE,CAAC;AACzC,CAAC;AAED;;;;;;;;;;GAUG;AACH,MAAM,UAAU,yBAAyB,CACvC,QAAgC,EAChC,QAAgB;IAEhB,OAAO,IAAI,GAAG,CACZ,QAAQ,CAAC,OAAO;SACb,MAAM,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,QAAQ,KAAK,QAAQ,IAAI,CAAC,CAAC,MAAM,KAAK,SAAS,CAAC;SAChE,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,SAAS,CAAC,CAC3B,CAAC;AACJ,CAAC"}

package/dist/atlas/witness-numeric.d.ts

@@ -0,0 +1,81 @@
+/**
+ * Atlas Phase 4, S4.2 — the numeric witness runner.
+ *
+ * Design note: `docs/planning/Atlas-Phase-4-Design.md` §2.3.
+ *
+ * A numeric witness evaluates a claim at TWO resolutions and reports whether
+ * refinement moved the answer toward the target. Two resolutions are the
+ * minimum that can say anything at all: one number agreeing with a target is
+ * consistent with a scheme that is simply wrong by an amount the tolerance
+ * happens to admit.
+ *
+ * ## The ratio is REPORTED, never compared against an expected order
+ *
+ * `convergence.ratio = |coarse − target| / |fine − target|` is the factor by
+ * which refinement reduced the error. A second-order scheme halving its step
+ * should give ≈ 4; a first-order one ≈ 2. **Nothing here asserts which.** The
+ * witnesses in hand do not all record their scheme's order, and asserting an
+ * order nobody wrote down would be fabrication of exactly the kind this sprint
+ * forbids. The number goes into the results artifact; a reviewer reads it.
+ *
+ * What the runner DOES decide is the weaker, checkable thing: a ratio at or
+ * below 1 means refinement did NOT improve the answer, and that earns
+ * `'unresolved'` rather than `'checked'` — a fine value that happens to land
+ * inside the tolerance while refinement is not converging is agreement by
+ * accident, and the tolerance alone cannot tell the two apart.
+ *
+ * @module atlas/witness-numeric
+ * @internal
+ */
+import type { WitnessRunResult } from './witness-result.js';
+/** The two-resolution error record a numeric witness carries. @internal */
+export interface Convergence {
+    /** |coarse − target| in the witness's own units. */
+    readonly coarse: number;
+    /** |fine − target|. */
+    readonly fine: number;
+    /**
+     * `coarse / fine`. Greater than 1 ⇒ refinement helped. `Infinity` when the
+     * fine error is exactly zero, which is a legitimate outcome for an exact
+     * scheme and must not be reported as a failure.
+     */
+    readonly ratio: number;
+}
+/** What a numeric witness asserts, and how to evaluate it. @internal */
+export interface NumericWitnessSpec {
+    /** The `Witness.id`. */
+    readonly id: string;
+    /**
+     * Evaluate the claim at a resolution. The meaning of the argument is the
+     * witness's own (a step count, a grid size, a truncation order); the runner
+     * only requires that `fineResolution` is the more refined of the two.
+     */
+    readonly evaluate: (resolution: number) => number;
+    /** The value the claim should reproduce. */
+    readonly target: number;
+    readonly coarseResolution: number;
+    readonly fineResolution: number;
+    /**
+     * Absolute tolerance the FINE evaluation must meet. Required: a default
+     * tolerance is a claim about accuracy that the caller, not this module,
+     * is entitled to make.
+     */
+    readonly tolerance: number;
+}
+/** A result that also carries the two-resolution record, when one was produced. */
+export interface NumericWitnessRunResult extends WitnessRunResult {
+    readonly kind: 'numeric';
+    /** Absent only when evaluation threw before any number existed. */
+    readonly convergence?: Convergence;
+}
+/**
+ * Run one numeric witness at two resolutions.
+ *
+ * @param spec - the claim and its two resolutions.
+ * @returns a result carrying `convergence`; **never throws** — an evaluation
+ * that throws is `'unresolved'` with `'parse-error'`, because a scheme that
+ * blew up has not refuted the physics it was evaluating.
+ * @internal
+ */
+export declare function runNumericWitness(spec: NumericWitnessSpec): NumericWitnessRunResult;
+//# sourceMappingURL=witness-numeric.d.ts.map

package/dist/atlas/witness-numeric.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"witness-numeric.d.ts","sourceRoot":"","sources":["../../src/atlas/witness-numeric.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;;;;;;;;;;GA4BG;AAEH,OAAO,KAAK,EAAE,gBAAgB,EAAE,MAAM,qBAAqB,CAAC;AAE5D,2EAA2E;AAC3E,MAAM,WAAW,WAAW;IAC1B,oDAAoD;IACpD,QAAQ,CAAC,MAAM,EAAE,MAAM,CAAC;IACxB,uBAAuB;IACvB,QAAQ,CAAC,IAAI,EAAE,MAAM,CAAC;IACtB;;;;OAIG;IACH,QAAQ,CAAC,KAAK,EAAE,MAAM,CAAC;CACxB;AAED,wEAAwE;AACxE,MAAM,WAAW,kBAAkB;IACjC,wBAAwB;IACxB,QAAQ,CAAC,EAAE,EAAE,MAAM,CAAC;IACpB;;;;OAIG;IACH,QAAQ,CAAC,QAAQ,EAAE,CAAC,UAAU,EAAE,MAAM,KAAK,MAAM,CAAC;IAClD,4CAA4C;IAC5C,QAAQ,CAAC,MAAM,EAAE,MAAM,CAAC;IACxB,QAAQ,CAAC,gBAAgB,EAAE,MAAM,CAAC;IAClC,QAAQ,CAAC,cAAc,EAAE,MAAM,CAAC;IAChC;;;;OAIG;IACH,QAAQ,CAAC,SAAS,EAAE,MAAM,CAAC;CAC5B;AAED,mFAAmF;AACnF,MAAM,WAAW,uBAAwB,SAAQ,gBAAgB;IAC/D,QAAQ,CAAC,IAAI,EAAE,SAAS,CAAC;IACzB,mEAAmE;IACnE,QAAQ,CAAC,WAAW,CAAC,EAAE,WAAW,CAAC;CACpC;AAED;;;;;;;;GAQG;AACH,wBAAgB,iBAAiB,CAAC,IAAI,EAAE,kBAAkB,GAAG,uBAAuB,CA0EnF"}

package/dist/atlas/witness-numeric.js

@@ -0,0 +1,105 @@
+/**
+ * Atlas Phase 4, S4.2 — the numeric witness runner.
+ *
+ * Design note: `docs/planning/Atlas-Phase-4-Design.md` §2.3.
+ *
+ * A numeric witness evaluates a claim at TWO resolutions and reports whether
+ * refinement moved the answer toward the target. Two resolutions are the
+ * minimum that can say anything at all: one number agreeing with a target is
+ * consistent with a scheme that is simply wrong by an amount the tolerance
+ * happens to admit.
+ *
+ * ## The ratio is REPORTED, never compared against an expected order
+ *
+ * `convergence.ratio = |coarse − target| / |fine − target|` is the factor by
+ * which refinement reduced the error. A second-order scheme halving its step
+ * should give ≈ 4; a first-order one ≈ 2. **Nothing here asserts which.** The
+ * witnesses in hand do not all record their scheme's order, and asserting an
+ * order nobody wrote down would be fabrication of exactly the kind this sprint
+ * forbids. The number goes into the results artifact; a reviewer reads it.
+ *
+ * What the runner DOES decide is the weaker, checkable thing: a ratio at or
+ * below 1 means refinement did NOT improve the answer, and that earns
+ * `'unresolved'` rather than `'checked'` — a fine value that happens to land
+ * inside the tolerance while refinement is not converging is agreement by
+ * accident, and the tolerance alone cannot tell the two apart.
+ *
+ * @module atlas/witness-numeric
+ * @internal
+ */
+/**
+ * Run one numeric witness at two resolutions.
+ *
+ * @param spec - the claim and its two resolutions.
+ * @returns a result carrying `convergence`; **never throws** — an evaluation
+ * that throws is `'unresolved'` with `'parse-error'`, because a scheme that
+ * blew up has not refuted the physics it was evaluating.
+ * @internal
+ */
+export function runNumericWitness(spec) {
+    const started = Date.now();
+    const base = { witnessId: spec.id, kind: 'numeric' };
+    const elapsed = () => Date.now() - started;
+    let coarseValue;
+    let fineValue;
+    try {
+        coarseValue = spec.evaluate(spec.coarseResolution);
+        fineValue = spec.evaluate(spec.fineResolution);
+    }
+    catch (err) {
+        return {
+            ...base,
+            status: 'unresolved',
+            reason: 'parse-error',
+            detail: `Evaluation threw rather than producing a number: ${String(err)}.`,
+            elapsedMs: elapsed(),
+        };
+    }
+    if (!Number.isFinite(coarseValue) || !Number.isFinite(fineValue)) {
+        return {
+            ...base,
+            status: 'unresolved',
+            reason: 'parse-error',
+            detail: `Evaluation produced a non-finite value (coarse ${coarseValue}, fine ` +
+                `${fineValue}). Nothing is concluded from it.`,
+            elapsedMs: elapsed(),
+        };
+    }
+    const coarse = Math.abs(coarseValue - spec.target);
+    const fine = Math.abs(fineValue - spec.target);
+    // fine === 0 ⇒ an exact scheme. Infinity is the honest ratio there, and the
+    // comparisons below treat it as "refinement helped", which it did.
+    const ratio = fine === 0 ? (coarse === 0 ? 1 : Infinity) : coarse / fine;
+    const convergence = { coarse, fine, ratio };
+    if (fine > spec.tolerance) {
+        return {
+            ...base,
+            status: 'refuted',
+            convergence,
+            detail: `The fine evaluation misses the target by ${fine}, outside the stated ` +
+                `tolerance ${spec.tolerance}.`,
+            elapsedMs: elapsed(),
+        };
+    }
+    if (ratio <= 1) {
+        return {
+            ...base,
+            status: 'unresolved',
+            reason: 'no-convergence',
+            convergence,
+            detail: `Refinement did not reduce the error (ratio ${ratio}). The fine value is ` +
+                'inside tolerance, but agreement that does not improve with refinement is ' +
+                'not evidence that the scheme converges to the target.',
+            elapsedMs: elapsed(),
+        };
+    }
+    return {
+        ...base,
+        status: 'checked',
+        convergence,
+        detail: `Fine error ${fine} is within tolerance ${spec.tolerance}, and refinement ` +
+            `reduced the error by a factor of ${ratio}.`,
+        elapsedMs: elapsed(),
+    };
+}
+//# sourceMappingURL=witness-numeric.js.map

package/dist/atlas/witness-numeric.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"witness-numeric.js","sourceRoot":"","sources":["../../src/atlas/witness-numeric.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;;;;;;;;;;GA4BG;AA+CH;;;;;;;;GAQG;AACH,MAAM,UAAU,iBAAiB,CAAC,IAAwB;IACxD,MAAM,OAAO,GAAG,IAAI,CAAC,GAAG,EAAE,CAAC;IAC3B,MAAM,IAAI,GAAG,EAAE,SAAS,EAAE,IAAI,CAAC,EAAE,EAAE,IAAI,EAAE,SAAkB,EAAE,CAAC;IAC9D,MAAM,OAAO,GAAG,GAAW,EAAE,CAAC,IAAI,CAAC,GAAG,EAAE,GAAG,OAAO,CAAC;IAEnD,IAAI,WAAmB,CAAC;IACxB,IAAI,SAAiB,CAAC;IACtB,IAAI,CAAC;QACH,WAAW,GAAG,IAAI,CAAC,QAAQ,CAAC,IAAI,CAAC,gBAAgB,CAAC,CAAC;QACnD,SAAS,GAAG,IAAI,CAAC,QAAQ,CAAC,IAAI,CAAC,cAAc,CAAC,CAAC;IACjD,CAAC;IAAC,OAAO,GAAG,EAAE,CAAC;QACb,OAAO;YACL,GAAG,IAAI;YACP,MAAM,EAAE,YAAY;YACpB,MAAM,EAAE,aAAa;YACrB,MAAM,EAAE,oDAAoD,MAAM,CAAC,GAAG,CAAC,GAAG;YAC1E,SAAS,EAAE,OAAO,EAAE;SACrB,CAAC;IACJ,CAAC;IAED,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,WAAW,CAAC,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,SAAS,CAAC,EAAE,CAAC;QACjE,OAAO;YACL,GAAG,IAAI;YACP,MAAM,EAAE,YAAY;YACpB,MAAM,EAAE,aAAa;YACrB,MAAM,EACJ,kDAAkD,WAAW,SAAS;gBACtE,GAAG,SAAS,kCAAkC;YAChD,SAAS,EAAE,OAAO,EAAE;SACrB,CAAC;IACJ,CAAC;IAED,MAAM,MAAM,GAAG,IAAI,CAAC,GAAG,CAAC,WAAW,GAAG,IAAI,CAAC,MAAM,CAAC,CAAC;IACnD,MAAM,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,SAAS,GAAG,IAAI,CAAC,MAAM,CAAC,CAAC;IAC/C,4EAA4E;IAC5E,mEAAmE;IACnE,MAAM,KAAK,GAAG,IAAI,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,MAAM,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,MAAM,GAAG,IAAI,CAAC;IACzE,MAAM,WAAW,GAAgB,EAAE,MAAM,EAAE,IAAI,EAAE,KAAK,EAAE,CAAC;IAEzD,IAAI,IAAI,GAAG,IAAI,CAAC,SAAS,EAAE,CAAC;QAC1B,OAAO;YACL,GAAG,IAAI;YACP,MAAM,EAAE,SAAS;YACjB,WAAW;YACX,MAAM,EACJ,4CAA4C,IAAI,uBAAuB;gBACvE,aAAa,IAAI,CAAC,SAAS,GAAG;YAChC,SAAS,EAAE,OAAO,EAAE;SACrB,CAAC;IACJ,CAAC;IAED,IAAI,KAAK,IAAI,CAAC,EAAE,CAAC;QACf,OAAO;YACL,GAAG,IAAI;YACP,MAAM,EAAE,YAAY;YACpB,MAAM,EAAE,gBAAgB;YACxB,WAAW;YACX,MAAM,EACJ,8CAA8C,KAAK,uBAAuB;gBAC1E,2EAA2E;gBAC3E,uDAAuD;YACzD,SAAS,EAAE,OAAO,EAAE;SACrB,CAAC;IACJ,CAAC;IAED,OAAO;QACL,GAAG,IAAI;QACP,MAAM,EAAE,SAAS;QACjB,WAAW;QACX,MAAM,EACJ,cAAc,IAAI,wBAAwB,IAAI,CAAC,SAAS,mBAAmB;YAC3E,oCAAoC,KAAK,GAAG;QAC9C,SAAS,EAAE,OAAO,EAAE;KACrB,CAAC;AACJ,CAAC"}

package/dist/atlas/witness-result.d.ts

@@ -0,0 +1,77 @@
+/**
+ * Atlas Phase 4, S4.2 — the shape a witness run reports.
+ *
+ * Design note: `docs/planning/Atlas-Phase-4-Design.md` §2.2.
+ *
+ * **The load-bearing distinction in this file: a check that DID NOT RUN is not
+ * a check that FAILED.** Peer absent, budget exhausted, and parser throw are
+ * all `'unresolved'`, each carrying the reason it could not answer. Only a
+ * check that COMPLETES and finds a disagreement is `'refuted'`.
+ *
+ * This repo has been bitten by exactly that conflation more than once — a
+ * checker that discarded its own stderr manufactured a false negative
+ * indistinguishable from a clean pass. Folding "could not run" into "failed"
+ * is the same defect with the sign flipped: it turns a missing optional peer
+ * into a refutation of the physics.
+ *
+ * Nothing here throws. A runner that threw would force every caller into a
+ * `try`/`catch` whose `catch` arm cannot tell which of the three outcomes it
+ * is looking at.
+ *
+ * @module atlas/witness-result
+ * @internal
+ */
+/**
+ * What a witness run concluded.
+ *
+ * - `'checked'` — the check RAN and the claim held.
+ * - `'refuted'` — the check RAN and the claim did not hold.
+ * - `'unresolved'` — the check did not run, or ran and could not decide.
+ *
+ * @internal
+ */
+export type WitnessStatus = 'checked' | 'refuted' | 'unresolved';
+/**
+ * Why a witness is `'unresolved'`. Required on that status so a reader never
+ * has to guess which of four very different situations produced it.
+ *
+ * @internal
+ */
+export type UnresolvedReason = 
+/** The optional peer that performs the check is not installed. */
+'peer-absent'
+/** The wall-clock budget elapsed before an answer arrived. */
+ | 'timeout'
+/** The input could not be parsed or evaluated at all. */
+ | 'parse-error'
+/** The check ran to completion and neither confirmed nor refuted. */
+ | 'not-simplified'
+/** Refinement did not improve the answer, so nothing is concluded. */
+ | 'no-convergence';
+/** What one witness run produced. @internal */
+export interface WitnessRunResult {
+    /** The `Witness.id` this run corresponds to (`'W7b'`). */
+    readonly witnessId: string;
+    readonly kind: 'symbolic' | 'numeric';
+    readonly status: WitnessStatus;
+    /** REQUIRED when `status === 'unresolved'`, absent otherwise. */
+    readonly reason?: UnresolvedReason;
+    /** Human-readable note — never the discriminator, always the explanation. */
+    readonly detail: string;
+    /** Milliseconds the runner waited. */
+    readonly elapsedMs: number;
+}
+/**
+ * The ids of runs that CHECKED, as `deriveEvidence` wants its
+ * `passingWitnessIds` argument.
+ *
+ * A convenience with one purpose: keep the mapping from run results to
+ * "passing" in ONE place. Writing `results.filter(r => r.status !== 'refuted')`
+ * at a call site would hand `symbolically-checked` to every unresolved run —
+ * the free-tag defect `derive-evidence.ts` documents at length, reintroduced
+ * one layer up.
+ *
+ * @internal
+ */
+export declare function passingWitnessIds(results: readonly WitnessRunResult[]): ReadonlySet<string>;
+//# sourceMappingURL=witness-result.d.ts.map