universal-physics-tensor 0.45.1 → 0.46.0

package/dist/atlas/waves/bridges.js

@@ -0,0 +1,265 @@
+/**
+ * Atlas Phase 4, S4.5 — the four bridges of the wave family.
+ *
+ * Design note: `docs/planning/Atlas-Phase-4-Design.md` §4.
+ *
+ * | Bridge | Relation | Premises | Witnesses |
+ * |---|---|---|---|
+ * | string → 1-D wave | `restriction` | 1 | WS1 (convergence) |
+ * | 1-D wave → d'Alembert | `derivation` | 1 | WS2 (convergence) |
+ * | Euler + adiabatic EOS → sound | `derivation` | **2 (hyperedge)** | WS3 (convergence), WS3b (Newton) |
+ * | Klein–Gordon → 1-D wave | `approximation` | 1 | WS4 (convergence), WS4b (long waves) |
+ *
+ * The chain → wave coarse-graining the plan also lists already exists in the
+ * oscillator family (`ab-chain-wave`) and is not duplicated.
+ *
+ * @module atlas/waves/bridges
+ */
+import { FORCE, FREQUENCY, VELOCITY } from '../../dimensional/types.js';
+import { DENSITY } from '../diffusion/dimensions.js';
+import { makeApproximation } from '../oscillators/bridges-limits.js';
+import { deriveRegimeGroups } from '../regime.js';
+import { LINEAR_DENSITY, PRESSURE, WAVENUMBER, WAVES_FAMILY_NAME } from './models.js';
+import { kleinGordonPhaseError } from './numerics.js';
+const TEST = 'tests/atlas/waves.test.ts';
+/** The perturbation ceiling of the sound derivation: `p₁/p₀ ≤ 1/100`. @internal */
+export const SOUND_MAX_PERTURBATION = 0.01;
+/** The dispersion-free domain edge: `ω₀/(c k) ≤ 0.1`. @internal */
+export const KG_MAX_DISPERSION_RATIO = 0.1;
+/** Bridge: string → 1-D wave equation, a restriction with c² = F/μ. @internal */
+export const BRIDGE_STRING_WAVE = {
+    id: 'ab-string-wave',
+    relation: 'restriction',
+    premises: ['model-string'],
+    conclusion: 'model-wave-1d',
+    transformation: 'y ↦ u, with c² = F/μ',
+    preserves: ['the solution operator for small slopes', 'normal-mode frequencies nπc/ℓ'],
+    doesNotPreserve: [
+        'the separate values of F and μ — only F/μ survives',
+        'the transverse geometry: the wave equation does not know u is a displacement',
+    ],
+    sideConditions: ['small slopes |y_x| ≪ 1', 'uniform tension and density', 'perfectly flexible'],
+    regime: {
+        family: WAVES_FAMILY_NAME,
+        inequalities: [],
+        groupDefinitions: deriveRegimeGroups(WAVES_FAMILY_NAME, [
+            { name: 'F', dim: FORCE },
+            { name: 'mu', dim: LINEAR_DENSITY },
+            { name: 'c', dim: VELOCITY },
+        ], []),
+    },
+    counterexamples: [],
+    evidence: new Set(['proposed', 'numerically-supported']),
+    witnesses: [
+        {
+            id: 'WS1',
+            kind: 'numeric',
+            test: TEST,
+            tolerance: 'leapfrog string solution at the midpoint within 2e-4 of sin(πx)cos(πct), c = √(F/μ), at 80 cells',
+        },
+    ],
+    citations: [
+        'Morse & Ingard, Theoretical Acoustics - the flexible string and its wave equation with c² = T/ρ',
+    ],
+    reviewStatus: 'proposed',
+};
+/** Bridge: 1-D wave equation → d'Alembert's general solution. @internal */
+export const BRIDGE_WAVE_DALEMBERT = {
+    id: 'ab-wave-dalembert',
+    relation: 'derivation',
+    premises: ['model-wave-1d'],
+    conclusion: 'model-dalembert',
+    transformation: 'characteristics ξ = x − ct, η = x + ct turn u_tt = c²u_xx into u_ξη = 0',
+    preserves: ['every C² solution on the whole line', 'the speed c'],
+    doesNotPreserve: [
+        'boundary conditions: on a bounded domain f and g are fixed only up to reflections',
+    ],
+    sideConditions: ['infinite line, or boundaries handled by the method of images', 'u ∈ C²'],
+    regime: {
+        family: WAVES_FAMILY_NAME,
+        inequalities: [],
+        groupDefinitions: deriveRegimeGroups(WAVES_FAMILY_NAME, [{ name: 'c', dim: VELOCITY }], []),
+    },
+    counterexamples: [],
+    evidence: new Set(['proposed', 'numerically-supported']),
+    witnesses: [
+        {
+            id: 'WS2',
+            kind: 'numeric',
+            test: TEST,
+            tolerance: 'finite-difference residual |u_tt − c²u_xx| below 2e-3 at h = 0.0125',
+        },
+    ],
+    citations: [
+        "d'Alembert, Recherches sur la courbe que forme une corde tendue mise en vibration (1747)",
+        'Strauss, Partial Differential Equations: An Introduction - §2.1, the wave equation on the line',
+    ],
+    reviewStatus: 'proposed',
+};
+/**
+ * Bridge: linearized Euler + adiabatic EOS → the acoustic wave equation.
+ *
+ * A HYPEREDGE: neither premise alone yields sound. The Euler pair has no
+ * closure from ρ′ to p′, and the equation of state has no dynamics.
+ *
+ * @internal
+ */
+export const BRIDGE_SOUND_SPEED = {
+    id: 'ab-sound-speed',
+    relation: 'derivation',
+    premises: ['model-euler-linear', 'model-adiabatic-eos'],
+    conclusion: 'model-sound',
+    transformation: 'close the Euler pair with p′ = (dp/dρ)|ρ₀ ρ′ = (γp₀/ρ₀) ρ′, then eliminate v: p′_tt = (γp₀/ρ₀) p′_xx',
+    preserves: ['linearity', 'the sound speed c_s² = γp₀/ρ₀'],
+    doesNotPreserve: [
+        'nonlinear steepening and shocks at finite amplitude',
+        'viscous and thermal attenuation',
+    ],
+    sideConditions: [
+        'small perturbations: |p′| ≪ p₀',
+        'adiabatic compression (no heat exchange within a wavelength)',
+        'fluid at rest, no ambient flow',
+    ],
+    regime: {
+        family: WAVES_FAMILY_NAME,
+        inequalities: [
+            {
+                group: 'p0 · p1^-1',
+                op: '>=',
+                bound: 1 / SOUND_MAX_PERTURBATION,
+                alias: '|p′| ≪ p₀ (machine form p₀/p₁ ≥ 100)',
+            },
+        ],
+        groupDefinitions: deriveRegimeGroups(WAVES_FAMILY_NAME, [
+            { name: 'rho0', dim: DENSITY },
+            { name: 'p0', dim: PRESSURE },
+            { name: 'p1', dim: PRESSURE },
+            { name: 'c', dim: VELOCITY },
+        ], ['gamma']),
+    },
+    counterexamples: [
+        {
+            description: 'The adiabatic premise is necessary. Closing the same Euler pair ISOTHERMALLY (Newton) ' +
+                'gives √(p₀/ρ₀) = 290.1 m/s for air at 20 °C, 15% below the measured ≈343 m/s; the ' +
+                'adiabatic closure gives 343.2 m/s.',
+            witness: 'WS3b',
+        },
+    ],
+    evidence: new Set(['proposed', 'numerically-supported']),
+    witnesses: [
+        {
+            id: 'WS3',
+            kind: 'numeric',
+            test: TEST,
+            tolerance: 'staggered Euler + adiabatic-closure solution within 2e-4 of the sound model’s sin(2πx)cos(2πc_s t) at 128 cells',
+        },
+        {
+            id: 'WS3b',
+            kind: 'numeric',
+            test: TEST,
+            tolerance: 'air, 101325 Pa, 1.204 kg/m³, γ = 1.4: adiabatic 343.25 m/s, isothermal 290.10 m/s',
+        },
+    ],
+    citations: [
+        'Laplace, Ann. Chim. Phys. 3 (1816) 238 - the adiabatic correction to Newton’s speed of sound',
+        'Landau & Lifshitz, Fluid Mechanics §64 - sound waves from the linearized equations of motion',
+    ],
+    reviewStatus: 'proposed',
+};
+/**
+ * Bridge: Klein–Gordon → 1-D wave, the dispersion-free limit ω₀/(ck) → 0.
+ *
+ * `delta` is the EXACT relative phase-velocity error at the domain edge, not a
+ * series term: {@link kleinGordonPhaseError} is increasing in ω₀/(ck), so its
+ * edge value is the supremum. The pendulum record learned the hard way that a
+ * truncated series can sit BELOW the error it claims to bound.
+ *
+ * @internal
+ */
+export const BRIDGE_KLEIN_GORDON_WAVE = {
+    id: 'ab-klein-gordon-wave',
+    relation: 'approximation',
+    premises: ['model-klein-gordon'],
+    conclusion: 'model-wave-1d',
+    transformation: 'drop the mass term ω₀²u; ω(k) = √(c²k² + ω₀²) → ck',
+    preserves: ['the wave speed c for short wavelengths', 'linearity'],
+    doesNotPreserve: [
+        'dispersion: wave packets spread under Klein–Gordon and not under the wave equation',
+        'the gap: Klein–Gordon has no mode below ω₀',
+    ],
+    sideConditions: ['short wavelengths: ω₀/(ck) ≤ 0.1', 'the bound is a PHASE-VELOCITY error, not uniform in time'],
+    bound: makeApproximation({
+        K: 1,
+        delta: kleinGordonPhaseError(KG_MAX_DISPERSION_RATIO, 1, 1),
+        deltaAt: (params) => {
+            const { omega0, c, k } = params;
+            return kleinGordonPhaseError(omega0 ?? Number.NaN, c ?? Number.NaN, k ?? Number.NaN);
+        },
+        norm: 'relative phase-velocity error of a Fourier mode',
+        domain: 'ω₀/(c k) ≤ 0.1',
+        horizon: 't ≪ π/(2 c k δ): the phase drift reaches π/2; machine form t < π/(2 c k δ(ω₀, c, k))',
+        horizonHolds: (t, params) => {
+            const { omega0, c, k } = params;
+            if (omega0 === undefined || c === undefined || k === undefined)
+                return false;
+            const delta = kleinGordonPhaseError(omega0, c, k);
+            if (!Number.isFinite(delta) || delta <= 0)
+                return delta === 0;
+            return t < Math.PI / (2 * c * k * delta);
+        },
+        parameterRange: 'ω₀/(c k) ≤ 0.1',
+        limitCharacter: 'regular',
+        // A phase-velocity error. The record says it is not uniform in time.
+        uniformity: ['phase velocity of one Fourier mode, for ω₀/(ck) ≤ 0.1'],
+    }),
+    regime: {
+        family: WAVES_FAMILY_NAME,
+        inequalities: [
+            {
+                group: 'c · omega0^-1 · k',
+                op: '>=',
+                bound: 1 / KG_MAX_DISPERSION_RATIO,
+                alias: 'ω₀/(c k) ≤ 0.1',
+            },
+        ],
+        groupDefinitions: deriveRegimeGroups(WAVES_FAMILY_NAME, [
+            { name: 'c', dim: VELOCITY },
+            { name: 'omega0', dim: FREQUENCY },
+            { name: 'k', dim: WAVENUMBER },
+        ], []),
+    },
+    counterexamples: [
+        {
+            description: 'At ω₀/(c k) = 1 the phase velocity is √2 c: a 41% error. The dispersion-free limit is a ' +
+                'short-wavelength statement and fails for the long waves the domain excludes.',
+            witness: 'WS4b',
+        },
+    ],
+    evidence: new Set(['proposed', 'numerically-supported']),
+    witnesses: [
+        {
+            id: 'WS4',
+            kind: 'numeric',
+            test: TEST,
+            tolerance: 'phase velocity within 2e-3 of c at k = 20 (ω₀ = c = 1); error falls ≈4× per doubling of k',
+        },
+        {
+            id: 'WS4b',
+            kind: 'numeric',
+            test: TEST,
+            tolerance: 'relative phase error √2 − 1 at ω₀/(ck) = 1',
+        },
+    ],
+    citations: [
+        'Whitham, Linear and Nonlinear Waves - §11.1, dispersive waves and the Klein-Gordon equation',
+    ],
+    reviewStatus: 'proposed',
+};
+/** The wave family's bridges, in design-note order. @internal */
+export const WAVE_BRIDGES = [
+    BRIDGE_STRING_WAVE,
+    BRIDGE_WAVE_DALEMBERT,
+    BRIDGE_SOUND_SPEED,
+    BRIDGE_KLEIN_GORDON_WAVE,
+];
+//# sourceMappingURL=bridges.js.map

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

@@ -0,0 +1 @@
+{"version":3,"file":"bridges.js","sourceRoot":"","sources":["../../../src/atlas/waves/bridges.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;GAgBG;AAEH,OAAO,EAAE,KAAK,EAAE,SAAS,EAAE,QAAQ,EAAE,MAAM,4BAA4B,CAAC;AACxE,OAAO,EAAE,OAAO,EAAE,MAAM,4BAA4B,CAAC;AACrD,OAAO,EAAE,iBAAiB,EAAE,MAAM,kCAAkC,CAAC;AACrE,OAAO,EAAE,kBAAkB,EAAE,MAAM,cAAc,CAAC;AAElD,OAAO,EAAE,cAAc,EAAE,QAAQ,EAAE,UAAU,EAAE,iBAAiB,EAAE,MAAM,aAAa,CAAC;AACtF,OAAO,EAAE,qBAAqB,EAAE,MAAM,eAAe,CAAC;AAEtD,MAAM,IAAI,GAAG,2BAA2B,CAAC;AAEzC,mFAAmF;AACnF,MAAM,CAAC,MAAM,sBAAsB,GAAG,IAAI,CAAC;AAE3C,mEAAmE;AACnE,MAAM,CAAC,MAAM,uBAAuB,GAAG,GAAG,CAAC;AAE3C,iFAAiF;AACjF,MAAM,CAAC,MAAM,kBAAkB,GAAgB;IAC7C,EAAE,EAAE,gBAAgB;IACpB,QAAQ,EAAE,aAAa;IACvB,QAAQ,EAAE,CAAC,cAAc,CAAC;IAC1B,UAAU,EAAE,eAAe;IAC3B,cAAc,EAAE,sBAAsB;IACtC,SAAS,EAAE,CAAC,wCAAwC,EAAE,+BAA+B,CAAC;IACtF,eAAe,EAAE;QACf,oDAAoD;QACpD,8EAA8E;KAC/E;IACD,cAAc,EAAE,CAAC,wBAAwB,EAAE,6BAA6B,EAAE,oBAAoB,CAAC;IAC/F,MAAM,EAAE;QACN,MAAM,EAAE,iBAAiB;QACzB,YAAY,EAAE,EAAE;QAChB,gBAAgB,EAAE,kBAAkB,CAClC,iBAAiB,EACjB;YACE,EAAE,IAAI,EAAE,GAAG,EAAE,GAAG,EAAE,KAAK,EAAE;YACzB,EAAE,IAAI,EAAE,IAAI,EAAE,GAAG,EAAE,cAAc,EAAE;YACnC,EAAE,IAAI,EAAE,GAAG,EAAE,GAAG,EAAE,QAAQ,EAAE;SAC7B,EACD,EAAE,CACH;KACF;IACD,eAAe,EAAE,EAAE;IACnB,QAAQ,EAAE,IAAI,GAAG,CAAC,CAAC,UAAU,EAAE,uBAAuB,CAAC,CAAC;IACxD,SAAS,EAAE;QACT;YACE,EAAE,EAAE,KAAK;YACT,IAAI,EAAE,SAAS;YACf,IAAI,EAAE,IAAI;YACV,SAAS,EACP,kGAAkG;SACrG;KACF;IACD,SAAS,EAAE;QACT,iGAAiG;KAClG;IACD,YAAY,EAAE,UAAU;CACzB,CAAC;AAEF,2EAA2E;AAC3E,MAAM,CAAC,MAAM,qBAAqB,GAAgB;IAChD,EAAE,EAAE,mBAAmB;IACvB,QAAQ,EAAE,YAAY;IACtB,QAAQ,EAAE,CAAC,eAAe,CAAC;IAC3B,UAAU,EAAE,iBAAiB;IAC7B,cAAc,EAAE,yEAAyE;IACzF,SAAS,EAAE,CAAC,qCAAqC,EAAE,aAAa,CAAC;IACjE,eAAe,EAAE;QACf,mFAAmF;KACpF;IACD,cAAc,EAAE,CAAC,8DAA8D,EAAE,QAAQ,CAAC;IAC1F,MAAM,EAAE;QACN,MAAM,EAAE,iBAAiB;QACzB,YAAY,EAAE,EAAE;QAChB,gBAAgB,EAAE,kBAAkB,CAAC,iBAAiB,EAAE,CAAC,EAAE,IAAI,EAAE,GAAG,EAAE,GAAG,EAAE,QAAQ,EAAE,CAAC,EAAE,EAAE,CAAC;KAC5F;IACD,eAAe,EAAE,EAAE;IACnB,QAAQ,EAAE,IAAI,GAAG,CAAC,CAAC,UAAU,EAAE,uBAAuB,CAAC,CAAC;IACxD,SAAS,EAAE;QACT;YACE,EAAE,EAAE,KAAK;YACT,IAAI,EAAE,SAAS;YACf,IAAI,EAAE,IAAI;YACV,SAAS,EAAE,qEAAqE;SACjF;KACF;IACD,SAAS,EAAE;QACT,0FAA0F;QAC1F,gGAAgG;KACjG;IACD,YAAY,EAAE,UAAU;CACzB,CAAC;AAEF;;;;;;;GAOG;AACH,MAAM,CAAC,MAAM,kBAAkB,GAAgB;IAC7C,EAAE,EAAE,gBAAgB;IACpB,QAAQ,EAAE,YAAY;IACtB,QAAQ,EAAE,CAAC,oBAAoB,EAAE,qBAAqB,CAAC;IACvD,UAAU,EAAE,aAAa;IACzB,cAAc,EACZ,sGAAsG;IACxG,SAAS,EAAE,CAAC,WAAW,EAAE,+BAA+B,CAAC;IACzD,eAAe,EAAE;QACf,qDAAqD;QACrD,iCAAiC;KAClC;IACD,cAAc,EAAE;QACd,gCAAgC;QAChC,8DAA8D;QAC9D,gCAAgC;KACjC;IACD,MAAM,EAAE;QACN,MAAM,EAAE,iBAAiB;QACzB,YAAY,EAAE;YACZ;gBACE,KAAK,EAAE,YAAY;gBACnB,EAAE,EAAE,IAAI;gBACR,KAAK,EAAE,CAAC,GAAG,sBAAsB;gBACjC,KAAK,EAAE,sCAAsC;aAC9C;SACF;QACD,gBAAgB,EAAE,kBAAkB,CAClC,iBAAiB,EACjB;YACE,EAAE,IAAI,EAAE,MAAM,EAAE,GAAG,EAAE,OAAO,EAAE;YAC9B,EAAE,IAAI,EAAE,IAAI,EAAE,GAAG,EAAE,QAAQ,EAAE;YAC7B,EAAE,IAAI,EAAE,IAAI,EAAE,GAAG,EAAE,QAAQ,EAAE;YAC7B,EAAE,IAAI,EAAE,GAAG,EAAE,GAAG,EAAE,QAAQ,EAAE;SAC7B,EACD,CAAC,OAAO,CAAC,CACV;KACF;IACD,eAAe,EAAE;QACf;YACE,WAAW,EACT,wFAAwF;gBACxF,oFAAoF;gBACpF,oCAAoC;YACtC,OAAO,EAAE,MAAM;SAChB;KACF;IACD,QAAQ,EAAE,IAAI,GAAG,CAAC,CAAC,UAAU,EAAE,uBAAuB,CAAC,CAAC;IACxD,SAAS,EAAE;QACT;YACE,EAAE,EAAE,KAAK;YACT,IAAI,EAAE,SAAS;YACf,IAAI,EAAE,IAAI;YACV,SAAS,EACP,iHAAiH;SACpH;QACD;YACE,EAAE,EAAE,MAAM;YACV,IAAI,EAAE,SAAS;YACf,IAAI,EAAE,IAAI;YACV,SAAS,EAAE,mFAAmF;SAC/F;KACF;IACD,SAAS,EAAE;QACT,8FAA8F;QAC9F,8FAA8F;KAC/F;IACD,YAAY,EAAE,UAAU;CACzB,CAAC;AAEF;;;;;;;;;GASG;AACH,MAAM,CAAC,MAAM,wBAAwB,GAAgB;IACnD,EAAE,EAAE,sBAAsB;IAC1B,QAAQ,EAAE,eAAe;IACzB,QAAQ,EAAE,CAAC,oBAAoB,CAAC;IAChC,UAAU,EAAE,eAAe;IAC3B,cAAc,EAAE,oDAAoD;IACpE,SAAS,EAAE,CAAC,wCAAwC,EAAE,WAAW,CAAC;IAClE,eAAe,EAAE;QACf,oFAAoF;QACpF,4CAA4C;KAC7C;IACD,cAAc,EAAE,CAAC,kCAAkC,EAAE,0DAA0D,CAAC;IAChH,KAAK,EAAE,iBAAiB,CAAC;QACvB,CAAC,EAAE,CAAC;QACJ,KAAK,EAAE,qBAAqB,CAAC,uBAAuB,EAAE,CAAC,EAAE,CAAC,CAAC;QAC3D,OAAO,EAAE,CAAC,MAAM,EAAE,EAAE;YAClB,MAAM,EAAE,MAAM,EAAE,CAAC,EAAE,CAAC,EAAE,GAAG,MAAM,CAAC;YAChC,OAAO,qBAAqB,CAAC,MAAM,IAAI,MAAM,CAAC,GAAG,EAAE,CAAC,IAAI,MAAM,CAAC,GAAG,EAAE,CAAC,IAAI,MAAM,CAAC,GAAG,CAAC,CAAC;QACvF,CAAC;QACD,IAAI,EAAE,iDAAiD;QACvD,MAAM,EAAE,gBAAgB;QACxB,OAAO,EACL,sFAAsF;QACxF,YAAY,EAAE,CAAC,CAAC,EAAE,MAAM,EAAE,EAAE;YAC1B,MAAM,EAAE,MAAM,EAAE,CAAC,EAAE,CAAC,EAAE,GAAG,MAAM,CAAC;YAChC,IAAI,MAAM,KAAK,SAAS,IAAI,CAAC,KAAK,SAAS,IAAI,CAAC,KAAK,SAAS;gBAAE,OAAO,KAAK,CAAC;YAC7E,MAAM,KAAK,GAAG,qBAAqB,CAAC,MAAM,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;YAClD,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,KAAK,CAAC,IAAI,KAAK,IAAI,CAAC;gBAAE,OAAO,KAAK,KAAK,CAAC,CAAC;YAC9D,OAAO,CAAC,GAAG,IAAI,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,KAAK,CAAC,CAAC;QAC3C,CAAC;QACD,cAAc,EAAE,gBAAgB;QAChC,cAAc,EAAE,SAAS;QACzB,qEAAqE;QACrE,UAAU,EAAE,CAAC,uDAAuD,CAAC;KACtE,CAAC;IACF,MAAM,EAAE;QACN,MAAM,EAAE,iBAAiB;QACzB,YAAY,EAAE;YACZ;gBACE,KAAK,EAAE,mBAAmB;gBAC1B,EAAE,EAAE,IAAI;gBACR,KAAK,EAAE,CAAC,GAAG,uBAAuB;gBAClC,KAAK,EAAE,gBAAgB;aACxB;SACF;QACD,gBAAgB,EAAE,kBAAkB,CAClC,iBAAiB,EACjB;YACE,EAAE,IAAI,EAAE,GAAG,EAAE,GAAG,EAAE,QAAQ,EAAE;YAC5B,EAAE,IAAI,EAAE,QAAQ,EAAE,GAAG,EAAE,SAAS,EAAE;YAClC,EAAE,IAAI,EAAE,GAAG,EAAE,GAAG,EAAE,UAAU,EAAE;SAC/B,EACD,EAAE,CACH;KACF;IACD,eAAe,EAAE;QACf;YACE,WAAW,EACT,0FAA0F;gBAC1F,8EAA8E;YAChF,OAAO,EAAE,MAAM;SAChB;KACF;IACD,QAAQ,EAAE,IAAI,GAAG,CAAC,CAAC,UAAU,EAAE,uBAAuB,CAAC,CAAC;IACxD,SAAS,EAAE;QACT;YACE,EAAE,EAAE,KAAK;YACT,IAAI,EAAE,SAAS;YACf,IAAI,EAAE,IAAI;YACV,SAAS,EAAE,2FAA2F;SACvG;QACD;YACE,EAAE,EAAE,MAAM;YACV,IAAI,EAAE,SAAS;YACf,IAAI,EAAE,IAAI;YACV,SAAS,EAAE,4CAA4C;SACxD;KACF;IACD,SAAS,EAAE;QACT,6FAA6F;KAC9F;IACD,YAAY,EAAE,UAAU;CACzB,CAAC;AAEF,iEAAiE;AACjE,MAAM,CAAC,MAAM,YAAY,GAA2B;IAClD,kBAAkB;IAClB,qBAAqB;IACrB,kBAAkB;IAClB,wBAAwB;CACzB,CAAC"}

package/dist/atlas/waves/index.d.ts

@@ -0,0 +1,16 @@
+/**
+ * Wave-family assembly (Atlas Phase 4, S4.5).
+ *
+ * Assembly only: this module defines no physics and no numbers of its own.
+ *
+ * @module atlas/waves
+ */
+import type { AtlasFamily } from '../oscillators/index.js';
+/**
+ * The wave family: seven models and seven bridges (four from S4.5, three from
+ * the Sprint 4 closure), no rejections. Several bridges end in other families.
+ *
+ * @internal
+ */
+export declare const WAVES_FAMILY: AtlasFamily;
+//# sourceMappingURL=index.d.ts.map

package/dist/atlas/waves/index.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"index.d.ts","sourceRoot":"","sources":["../../../src/atlas/waves/index.ts"],"names":[],"mappings":"AAAA;;;;;;GAMG;AAEH,OAAO,KAAK,EAAE,WAAW,EAAE,MAAM,yBAAyB,CAAC;AAK3D;;;;;GAKG;AACH,eAAO,MAAM,YAAY,EAAE,WAK1B,CAAC"}

package/dist/atlas/waves/index.js

@@ -0,0 +1,23 @@
+/**
+ * Wave-family assembly (Atlas Phase 4, S4.5).
+ *
+ * Assembly only: this module defines no physics and no numbers of its own.
+ *
+ * @module atlas/waves
+ */
+import { WAVE_BRIDGES } from './bridges.js';
+import { WAVE_CLOSURE_BRIDGES } from './bridges-closure.js';
+import { WAVE_MODELS, WAVES_FAMILY_NAME } from './models.js';
+/**
+ * The wave family: seven models and seven bridges (four from S4.5, three from
+ * the Sprint 4 closure), no rejections. Several bridges end in other families.
+ *
+ * @internal
+ */
+export const WAVES_FAMILY = {
+    family: WAVES_FAMILY_NAME,
+    models: WAVE_MODELS,
+    bridges: [...WAVE_BRIDGES, ...WAVE_CLOSURE_BRIDGES],
+    rejections: [],
+};
+//# sourceMappingURL=index.js.map

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

@@ -0,0 +1 @@
+{"version":3,"file":"index.js","sourceRoot":"","sources":["../../../src/atlas/waves/index.ts"],"names":[],"mappings":"AAAA;;;;;;GAMG;AAGH,OAAO,EAAE,YAAY,EAAE,MAAM,cAAc,CAAC;AAC5C,OAAO,EAAE,oBAAoB,EAAE,MAAM,sBAAsB,CAAC;AAC5D,OAAO,EAAE,WAAW,EAAE,iBAAiB,EAAE,MAAM,aAAa,CAAC;AAE7D;;;;;GAKG;AACH,MAAM,CAAC,MAAM,YAAY,GAAgB;IACvC,MAAM,EAAE,iBAAiB;IACzB,MAAM,EAAE,WAAW;IACnB,OAAO,EAAE,CAAC,GAAG,YAAY,EAAE,GAAG,oBAAoB,CAAC;IACnD,UAAU,EAAE,EAAE;CACf,CAAC"}

package/dist/atlas/waves/models.d.ts

@@ -0,0 +1,29 @@
+/**
+ * Atlas Phase 4, S4.5 — the models of the wave family.
+ *
+ * Design note: `docs/planning/Atlas-Phase-4-Design.md` §4.
+ *
+ * The 1-D wave equation itself is NOT redefined here. It is `model-wave-1d` in
+ * the oscillator family, where Phase 0 put it as the conclusion of the chain →
+ * wave coarse-graining. The wave family's bridges point at that same record,
+ * because a second definition of one model would be a second source of truth.
+ * `tests/atlas/families.test.ts` pins that model ids are unique across the
+ * whole atlas and that every bridge endpoint resolves in SOME registered family.
+ *
+ * @module atlas/waves/models
+ */
+import type { Dimension } from '../../dimensional/types.js';
+import type { AtlasModel } from '../model.js';
+/** The family name every wave record carries. @internal */
+export declare const WAVES_FAMILY_NAME = "waves";
+/** kg/m, `L^-1 M` — `src/canonical/entries/dimensional-classics.ts` (`LIN_DENSITY`). @internal */
+export declare const LINEAR_DENSITY: Dimension;
+/** Pa, `L^-1 M T^-2` — `src/canonical/entries/fluids-waves.ts` (`PRESSURE`). @internal */
+export declare const PRESSURE: Dimension;
+/** Flexural rigidity E·I, N·m², `L^3 M T^-2`. @internal */
+export declare const FLEXURAL_RIGIDITY: Dimension;
+/** rad/m, `L^-1`. @internal */
+export declare const WAVENUMBER: Dimension;
+/** The wave-family models, in design-note order. @internal */
+export declare const WAVE_MODELS: readonly AtlasModel[];
+//# sourceMappingURL=models.d.ts.map

package/dist/atlas/waves/models.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"models.d.ts","sourceRoot":"","sources":["../../../src/atlas/waves/models.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;GAaG;AAIH,OAAO,KAAK,EAAE,SAAS,EAAE,MAAM,4BAA4B,CAAC;AAG5D,OAAO,KAAK,EAAE,UAAU,EAAE,MAAM,aAAa,CAAC;AAG9C,2DAA2D;AAC3D,eAAO,MAAM,iBAAiB,UAAU,CAAC;AAEzC,kGAAkG;AAClG,eAAO,MAAM,cAAc,EAAE,SAAsB,CAAC;AAEpD,0FAA0F;AAC1F,eAAO,MAAM,QAAQ,EAAE,SAA0B,CAAC;AAElD,2DAA2D;AAC3D,eAAO,MAAM,iBAAiB,EAAE,SAAyB,CAAC;AAE1D,+BAA+B;AAC/B,eAAO,MAAM,UAAU,EAAE,SAAmB,CAAC;AA+B7C,8DAA8D;AAC9D,eAAO,MAAM,WAAW,EAAE,SAAS,UAAU,EA8E5C,CAAC"}

package/dist/atlas/waves/models.js

@@ -0,0 +1,127 @@
+/**
+ * Atlas Phase 4, S4.5 — the models of the wave family.
+ *
+ * Design note: `docs/planning/Atlas-Phase-4-Design.md` §4.
+ *
+ * The 1-D wave equation itself is NOT redefined here. It is `model-wave-1d` in
+ * the oscillator family, where Phase 0 put it as the conclusion of the chain →
+ * wave coarse-graining. The wave family's bridges point at that same record,
+ * because a second definition of one model would be a second source of truth.
+ * `tests/atlas/families.test.ts` pins that model ids are unique across the
+ * whole atlas and that every bridge endpoint resolves in SOME registered family.
+ *
+ * @module atlas/waves/models
+ */
+import { dim } from '../../dimensional/ast-builders.js';
+import { FORCE, FREQUENCY, VELOCITY } from '../../dimensional/types.js';
+import { DENSITY } from '../diffusion/dimensions.js';
+import { deriveRegimeGroups } from '../regime.js';
+/** The family name every wave record carries. @internal */
+export const WAVES_FAMILY_NAME = 'waves';
+/** kg/m, `L^-1 M` — `src/canonical/entries/dimensional-classics.ts` (`LIN_DENSITY`). @internal */
+export const LINEAR_DENSITY = dim(-1, 1);
+/** Pa, `L^-1 M T^-2` — `src/canonical/entries/fluids-waves.ts` (`PRESSURE`). @internal */
+export const PRESSURE = dim(-1, 1, -2);
+/** Flexural rigidity E·I, N·m², `L^3 M T^-2`. @internal */
+export const FLEXURAL_RIGIDITY = dim(3, 1, -2);
+/** rad/m, `L^-1`. @internal */
+export const WAVENUMBER = dim(-1);
+function model(spec) {
+    const dimensionless = spec.dimensionlessInputs ?? [];
+    return {
+        id: spec.id,
+        family: WAVES_FAMILY_NAME,
+        stateSpace: spec.stateSpace,
+        dynamics: spec.dynamics,
+        observables: spec.observables,
+        parameters: spec.parameters,
+        dimensionlessInputs: dimensionless,
+        canonicalRefs: spec.canonicalRefs,
+        regime: {
+            family: WAVES_FAMILY_NAME,
+            inequalities: [],
+            groupDefinitions: deriveRegimeGroups(WAVES_FAMILY_NAME, spec.parameters, dimensionless),
+        },
+    };
+}
+/** The wave-family models, in design-note order. @internal */
+export const WAVE_MODELS = [
+    model({
+        id: 'model-string',
+        stateSpace: 'y(x, t), transverse displacement of a string on [0, ℓ]',
+        dynamics: 'μ y_tt = F y_xx (small slopes, |y_x| ≪ 1)',
+        observables: ['y', 'y_t', 'normal-mode frequencies'],
+        parameters: [
+            { name: 'F', dim: FORCE },
+            { name: 'mu', dim: LINEAR_DENSITY },
+        ],
+        canonicalRefs: ['CE-string-wave-speed'],
+    }),
+    model({
+        id: 'model-dalembert',
+        stateSpace: 'u(x, t) = f(x − ct) + g(x + ct), f and g twice differentiable',
+        dynamics: 'two profiles translating rigidly at ±c',
+        observables: ['u', 'the two travelling profiles'],
+        parameters: [{ name: 'c', dim: VELOCITY }],
+        canonicalRefs: [],
+    }),
+    model({
+        id: 'model-euler-linear',
+        stateSpace: '(ρ′, v)(x, t), small perturbations of a fluid at rest',
+        dynamics: '∂ρ′/∂t = −ρ₀ ∂v/∂x,  ρ₀ ∂v/∂t = −∂p′/∂x',
+        observables: ['ρ′', 'v', 'p′'],
+        parameters: [
+            { name: 'rho0', dim: DENSITY },
+            { name: 'p0', dim: PRESSURE },
+        ],
+        canonicalRefs: [],
+    }),
+    model({
+        id: 'model-adiabatic-eos',
+        stateSpace: 'p(ρ), a barotropic closure',
+        dynamics: 'p = p₀ (ρ/ρ₀)^γ',
+        observables: ['p', 'dp/dρ'],
+        parameters: [
+            { name: 'rho0', dim: DENSITY },
+            { name: 'p0', dim: PRESSURE },
+        ],
+        dimensionlessInputs: ['gamma'],
+        canonicalRefs: [],
+    }),
+    model({
+        id: 'model-sound',
+        stateSpace: 'p′(x, t), the acoustic pressure',
+        dynamics: 'p′_tt = c_s² p′_xx with c_s² = γ p₀/ρ₀',
+        observables: ['p′', 'sound speed c_s'],
+        parameters: [
+            { name: 'rho0', dim: DENSITY },
+            { name: 'p0', dim: PRESSURE },
+        ],
+        dimensionlessInputs: ['gamma'],
+        canonicalRefs: ['CE-sound-speed'],
+    }),
+    model({
+        id: 'model-klein-gordon',
+        stateSpace: 'u(x, t), a field on ℝ × ℝ',
+        dynamics: 'u_tt = c² u_xx − ω₀² u (dispersive: ω² = c²k² + ω₀²)',
+        observables: ['u', 'dispersion ω(k)', 'phase velocity'],
+        parameters: [
+            { name: 'c', dim: VELOCITY },
+            { name: 'omega0', dim: FREQUENCY },
+        ],
+        canonicalRefs: [],
+    }),
+    model({
+        id: 'model-stiff-string',
+        stateSpace: 'y(x, t), transverse displacement of a string with bending stiffness',
+        dynamics: 'μ y_tt = F y_xx − EI y_xxxx (dispersive: ω² = (F/μ)k² + (EI/μ)k⁴)',
+        observables: ['y', 'partial frequencies (inharmonicity)'],
+        parameters: [
+            { name: 'F', dim: FORCE },
+            { name: 'mu', dim: LINEAR_DENSITY },
+            { name: 'EI', dim: FLEXURAL_RIGIDITY },
+        ],
+        canonicalRefs: [],
+    }),
+];
+//# sourceMappingURL=models.js.map

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

@@ -0,0 +1 @@
+{"version":3,"file":"models.js","sourceRoot":"","sources":["../../../src/atlas/waves/models.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;GAaG;AAEH,OAAO,EAAE,GAAG,EAAE,MAAM,mCAAmC,CAAC;AAGxD,OAAO,EAAE,KAAK,EAAE,SAAS,EAAE,QAAQ,EAAE,MAAM,4BAA4B,CAAC;AACxE,OAAO,EAAE,OAAO,EAAE,MAAM,4BAA4B,CAAC;AAErD,OAAO,EAAE,kBAAkB,EAAE,MAAM,cAAc,CAAC;AAElD,2DAA2D;AAC3D,MAAM,CAAC,MAAM,iBAAiB,GAAG,OAAO,CAAC;AAEzC,kGAAkG;AAClG,MAAM,CAAC,MAAM,cAAc,GAAc,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;AAEpD,0FAA0F;AAC1F,MAAM,CAAC,MAAM,QAAQ,GAAc,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;AAElD,2DAA2D;AAC3D,MAAM,CAAC,MAAM,iBAAiB,GAAc,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;AAE1D,+BAA+B;AAC/B,MAAM,CAAC,MAAM,UAAU,GAAc,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;AAY7C,SAAS,KAAK,CAAC,IAAe;IAC5B,MAAM,aAAa,GAAG,IAAI,CAAC,mBAAmB,IAAI,EAAE,CAAC;IACrD,OAAO;QACL,EAAE,EAAE,IAAI,CAAC,EAAE;QACX,MAAM,EAAE,iBAAiB;QACzB,UAAU,EAAE,IAAI,CAAC,UAAU;QAC3B,QAAQ,EAAE,IAAI,CAAC,QAAQ;QACvB,WAAW,EAAE,IAAI,CAAC,WAAW;QAC7B,UAAU,EAAE,IAAI,CAAC,UAAU;QAC3B,mBAAmB,EAAE,aAAa;QAClC,aAAa,EAAE,IAAI,CAAC,aAAa;QACjC,MAAM,EAAE;YACN,MAAM,EAAE,iBAAiB;YACzB,YAAY,EAAE,EAAE;YAChB,gBAAgB,EAAE,kBAAkB,CAAC,iBAAiB,EAAE,IAAI,CAAC,UAAU,EAAE,aAAa,CAAC;SACxF;KACF,CAAC;AACJ,CAAC;AAED,8DAA8D;AAC9D,MAAM,CAAC,MAAM,WAAW,GAA0B;IAChD,KAAK,CAAC;QACJ,EAAE,EAAE,cAAc;QAClB,UAAU,EAAE,wDAAwD;QACpE,QAAQ,EAAE,2CAA2C;QACrD,WAAW,EAAE,CAAC,GAAG,EAAE,KAAK,EAAE,yBAAyB,CAAC;QACpD,UAAU,EAAE;YACV,EAAE,IAAI,EAAE,GAAG,EAAE,GAAG,EAAE,KAAK,EAAE;YACzB,EAAE,IAAI,EAAE,IAAI,EAAE,GAAG,EAAE,cAAc,EAAE;SACpC;QACD,aAAa,EAAE,CAAC,sBAAsB,CAAC;KACxC,CAAC;IACF,KAAK,CAAC;QACJ,EAAE,EAAE,iBAAiB;QACrB,UAAU,EAAE,+DAA+D;QAC3E,QAAQ,EAAE,wCAAwC;QAClD,WAAW,EAAE,CAAC,GAAG,EAAE,6BAA6B,CAAC;QACjD,UAAU,EAAE,CAAC,EAAE,IAAI,EAAE,GAAG,EAAE,GAAG,EAAE,QAAQ,EAAE,CAAC;QAC1C,aAAa,EAAE,EAAE;KAClB,CAAC;IACF,KAAK,CAAC;QACJ,EAAE,EAAE,oBAAoB;QACxB,UAAU,EAAE,uDAAuD;QACnE,QAAQ,EAAE,yCAAyC;QACnD,WAAW,EAAE,CAAC,IAAI,EAAE,GAAG,EAAE,IAAI,CAAC;QAC9B,UAAU,EAAE;YACV,EAAE,IAAI,EAAE,MAAM,EAAE,GAAG,EAAE,OAAO,EAAE;YAC9B,EAAE,IAAI,EAAE,IAAI,EAAE,GAAG,EAAE,QAAQ,EAAE;SAC9B;QACD,aAAa,EAAE,EAAE;KAClB,CAAC;IACF,KAAK,CAAC;QACJ,EAAE,EAAE,qBAAqB;QACzB,UAAU,EAAE,4BAA4B;QACxC,QAAQ,EAAE,iBAAiB;QAC3B,WAAW,EAAE,CAAC,GAAG,EAAE,OAAO,CAAC;QAC3B,UAAU,EAAE;YACV,EAAE,IAAI,EAAE,MAAM,EAAE,GAAG,EAAE,OAAO,EAAE;YAC9B,EAAE,IAAI,EAAE,IAAI,EAAE,GAAG,EAAE,QAAQ,EAAE;SAC9B;QACD,mBAAmB,EAAE,CAAC,OAAO,CAAC;QAC9B,aAAa,EAAE,EAAE;KAClB,CAAC;IACF,KAAK,CAAC;QACJ,EAAE,EAAE,aAAa;QACjB,UAAU,EAAE,iCAAiC;QAC7C,QAAQ,EAAE,wCAAwC;QAClD,WAAW,EAAE,CAAC,IAAI,EAAE,iBAAiB,CAAC;QACtC,UAAU,EAAE;YACV,EAAE,IAAI,EAAE,MAAM,EAAE,GAAG,EAAE,OAAO,EAAE;YAC9B,EAAE,IAAI,EAAE,IAAI,EAAE,GAAG,EAAE,QAAQ,EAAE;SAC9B;QACD,mBAAmB,EAAE,CAAC,OAAO,CAAC;QAC9B,aAAa,EAAE,CAAC,gBAAgB,CAAC;KAClC,CAAC;IACF,KAAK,CAAC;QACJ,EAAE,EAAE,oBAAoB;QACxB,UAAU,EAAE,2BAA2B;QACvC,QAAQ,EAAE,sDAAsD;QAChE,WAAW,EAAE,CAAC,GAAG,EAAE,iBAAiB,EAAE,gBAAgB,CAAC;QACvD,UAAU,EAAE;YACV,EAAE,IAAI,EAAE,GAAG,EAAE,GAAG,EAAE,QAAQ,EAAE;YAC5B,EAAE,IAAI,EAAE,QAAQ,EAAE,GAAG,EAAE,SAAS,EAAE;SACnC;QACD,aAAa,EAAE,EAAE;KAClB,CAAC;IACF,KAAK,CAAC;QACJ,EAAE,EAAE,oBAAoB;QACxB,UAAU,EAAE,qEAAqE;QACjF,QAAQ,EAAE,mEAAmE;QAC7E,WAAW,EAAE,CAAC,GAAG,EAAE,qCAAqC,CAAC;QACzD,UAAU,EAAE;YACV,EAAE,IAAI,EAAE,GAAG,EAAE,GAAG,EAAE,KAAK,EAAE;YACzB,EAAE,IAAI,EAAE,IAAI,EAAE,GAAG,EAAE,cAAc,EAAE;YACnC,EAAE,IAAI,EAAE,IAAI,EAAE,GAAG,EAAE,iBAAiB,EAAE;SACvC;QACD,aAAa,EAAE,EAAE;KAClB,CAAC;CACH,CAAC"}

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

@@ -0,0 +1,149 @@
+/**
+ * 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
+ */
+/** Parameters of the WS1 string fixture. @internal */
+export interface StringFixture {
+    /** Tension F. */
+    readonly tension: number;
+    /** Linear mass density μ. */
+    readonly mu: number;
+    /** Final time. */
+    readonly tEnd: number;
+}
+/**
+ * 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 declare function stringLeapfrogMidpoint(cells: number, f: StringFixture): number;
+/** Parameters of the WS2 d'Alembert fixture. @internal */
+export interface DalembertFixture {
+    readonly c: number;
+    /** The point `(x, t)` the residual is evaluated at. */
+    readonly x: number;
+    readonly t: number;
+    /** Base finite-difference step; the resolution divides it. */
+    readonly h0: number;
+}
+/**
+ * The d'Alembert form `u = f(x − ct) + g(x + ct)` with the WS2 profiles.
+ *
+ * @internal
+ */
+export declare function dalembert(x: number, t: number, c: number): number;
+/**
+ * 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 declare function dalembertResidual(resolution: number, f: DalembertFixture): number;
+/** Parameters of the WS3 acoustic fixture. @internal */
+export interface AcousticFixture {
+    /** Ambient pressure p₀. */
+    readonly p0: number;
+    /** Ambient density ρ₀. */
+    readonly rho0: number;
+    /** Adiabatic index γ. */
+    readonly gamma: number;
+    readonly tEnd: number;
+}
+/**
+ * The adiabatic equation of state `p(ρ) = p₀ (ρ/ρ₀)^γ` — the SECOND premise of
+ * the sound-speed hyperedge — and its derivative at `ρ`.
+ *
+ * @internal
+ */
+export declare function adiabaticSlope(rho: number, f: AcousticFixture): number;
+/**
+ * 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 declare function acousticLeapfrogQuarter(cells: number, f: AcousticFixture): number;
+/**
+ * Speed of sound from the adiabatic (Laplace) and isothermal (Newton) closures.
+ *
+ * @internal
+ */
+export declare function soundSpeeds(p0: number, rho0: number, gamma: number): {
+    readonly adiabatic: number;
+    readonly isothermal: number;
+};
+/**
+ * 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 declare function kleinGordonPhaseError(omega0: number, c: number, k: number): number;
+/**
+ * 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 declare function kleinGordonPhaseVelocity(resolution: number, omega0: number, c: number, k0: number): number;
+/**
+ * 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 declare function kgNonrelativisticError(x: number): number;
+/**
+ * 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 declare function kgUniformModeValue(steps: number, c: number, omega0: number, tEnd: number, cells?: number): number;
+/**
+ * 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 declare function stiffStringPhaseError(F: number, EI: number, k: number): number;
+/**
+ * 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 declare function stiffStringPhaseVelocity(resolution: number, F: number, mu: number, EI: number, k0: number): number;
+//# sourceMappingURL=numerics.d.ts.map

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

@@ -0,0 +1 @@
+{"version":3,"file":"numerics.d.ts","sourceRoot":"","sources":["../../../src/atlas/waves/numerics.ts"],"names":[],"mappings":"AAAA;;;;;;;;GAQG;AAEH,sDAAsD;AACtD,MAAM,WAAW,aAAa;IAC5B,iBAAiB;IACjB,QAAQ,CAAC,OAAO,EAAE,MAAM,CAAC;IACzB,6BAA6B;IAC7B,QAAQ,CAAC,EAAE,EAAE,MAAM,CAAC;IACpB,kBAAkB;IAClB,QAAQ,CAAC,IAAI,EAAE,MAAM,CAAC;CACvB;AAED;;;;;;;;;;;GAWG;AACH,wBAAgB,sBAAsB,CAAC,KAAK,EAAE,MAAM,EAAE,CAAC,EAAE,aAAa,GAAG,MAAM,CAwB9E;AAED,0DAA0D;AAC1D,MAAM,WAAW,gBAAgB;IAC/B,QAAQ,CAAC,CAAC,EAAE,MAAM,CAAC;IACnB,uDAAuD;IACvD,QAAQ,CAAC,CAAC,EAAE,MAAM,CAAC;IACnB,QAAQ,CAAC,CAAC,EAAE,MAAM,CAAC;IACnB,8DAA8D;IAC9D,QAAQ,CAAC,EAAE,EAAE,MAAM,CAAC;CACrB;AAMD;;;;GAIG;AACH,wBAAgB,SAAS,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,GAAG,MAAM,CAEjE;AAED;;;;;;GAMG;AACH,wBAAgB,iBAAiB,CAAC,UAAU,EAAE,MAAM,EAAE,CAAC,EAAE,gBAAgB,GAAG,MAAM,CAMjF;AAED,wDAAwD;AACxD,MAAM,WAAW,eAAe;IAC9B,2BAA2B;IAC3B,QAAQ,CAAC,EAAE,EAAE,MAAM,CAAC;IACpB,0BAA0B;IAC1B,QAAQ,CAAC,IAAI,EAAE,MAAM,CAAC;IACtB,yBAAyB;IACzB,QAAQ,CAAC,KAAK,EAAE,MAAM,CAAC;IACvB,QAAQ,CAAC,IAAI,EAAE,MAAM,CAAC;CACvB;AAED;;;;;GAKG;AACH,wBAAgB,cAAc,CAAC,GAAG,EAAE,MAAM,EAAE,CAAC,EAAE,eAAe,GAAG,MAAM,CAEtE;AAED;;;;;;;;;;;;;GAaG;AACH,wBAAgB,uBAAuB,CAAC,KAAK,EAAE,MAAM,EAAE,CAAC,EAAE,eAAe,GAAG,MAAM,CA0BjF;AAED;;;;GAIG;AACH,wBAAgB,WAAW,CAAC,EAAE,EAAE,MAAM,EAAE,IAAI,EAAE,MAAM,EAAE,KAAK,EAAE,MAAM,GAAG;IACpE,QAAQ,CAAC,SAAS,EAAE,MAAM,CAAC;IAC3B,QAAQ,CAAC,UAAU,EAAE,MAAM,CAAC;CAC7B,CAEA;AAED;;;;;;;;;GASG;AACH,wBAAgB,qBAAqB,CAAC,MAAM,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,GAAG,MAAM,CAGlF;AAED;;;;;;GAMG;AACH,wBAAgB,wBAAwB,CACtC,UAAU,EAAE,MAAM,EAClB,MAAM,EAAE,MAAM,EACd,CAAC,EAAE,MAAM,EACT,EAAE,EAAE,MAAM,GACT,MAAM,CAGR;AAID;;;;;;;GAOG;AACH,wBAAgB,sBAAsB,CAAC,CAAC,EAAE,MAAM,GAAG,MAAM,CAGxD;AAED;;;;;;;GAOG;AACH,wBAAgB,kBAAkB,CAChC,KAAK,EAAE,MAAM,EACb,CAAC,EAAE,MAAM,EACT,MAAM,EAAE,MAAM,EACd,IAAI,EAAE,MAAM,EACZ,KAAK,SAAK,GACT,MAAM,CAkBR;AAED;;;;;;GAMG;AACH,wBAAgB,qBAAqB,CAAC,CAAC,EAAE,MAAM,EAAE,EAAE,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,GAAG,MAAM,CAE9E;AAED;;;;;;GAMG;AACH,wBAAgB,wBAAwB,CACtC,UAAU,EAAE,MAAM,EAClB,CAAC,EAAE,MAAM,EACT,EAAE,EAAE,MAAM,EACV,EAAE,EAAE,MAAM,EACV,EAAE,EAAE,MAAM,GACT,MAAM,CAGR"}