universal-physics-tensor 9.0.0 → 11.0.0

package/dist/atlas/oscillators/bridges-limits.d.ts

@@ -35,11 +35,11 @@ import type { ApproximationBound, AtlasBridge, RelationContract } from '../types
  * @internal
  */
 export declare function makeApproximation(bound: ApproximationBound): ApproximationBound;
-/** Arithmetic–geometric mean, the standard evaluation route for `K`. @internal */
-export declare function agm(a0: number, b0: number): number;
 /**
  * The EXACT relative period error of the pendulum at amplitude `theta0`:
- * `(2/π) K(sin(θ0/2)) − 1 = 1/AGM(1, cos(θ0/2)) − 1`.
+ * `(2/π) K(sin(θ0/2)) − 1`, with `K` MathTS `ellipticK`, which takes the
+ * PARAMETER `m = k² = sin²(θ0/2)` and evaluates it by the AGM
+ * (`K = π / (2·AGM(1, √(1−m)))`, so this is `1/AGM(1, cos(θ0/2)) − 1`).
  *
  * This is the machine form of `ab-pendulum-linear`'s `delta`. It is SHARP —
  * the bound equals the error it bounds — and it is strictly increasing in
@@ -102,19 +102,6 @@ export declare const AB_PENDULUM_LINEAR: AtlasBridge;
 export declare const AB_DAMPED_MASSLESS: AtlasBridge;
 /** The two approximation bridges of this module, in design-note order. @internal */
 export declare const LIMIT_BRIDGES: readonly AtlasBridge[];
-/**
- * The two approximation bridges re-registered through `RelationContract`.
- *
- * The union REQUIRES `bound` on an `approximation`, and `ApproximationBound`
- * requires its own machine `horizonHolds`, so these two lines are what turn
- * Sprint 0's doc-comment requirement into a compile-and-throw one. Both bounds
- * are the records' own, measured ones — nothing is restated here.
- *
- * @internal
- */
-export declare const CONTRACT_PENDULUM_LINEAR: RelationContract;
-/** @internal */
-export declare const CONTRACT_DAMPED_MASSLESS: RelationContract;
 /** The two limit contracts, in design-note order. @internal */
 export declare const LIMIT_CONTRACTS: readonly RelationContract[];
 //# sourceMappingURL=bridges-limits.d.ts.map

package/dist/atlas/oscillators/bridges-limits.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"bridges-limits.d.ts","sourceRoot":"","sources":["../../../src/atlas/oscillators/bridges-limits.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;;;GAqBG;AAKH,OAAO,KAAK,EAAE,kBAAkB,EAAE,WAAW,EAAE,gBAAgB,EAAU,MAAM,aAAa,CAAC;AAS7F;;;;;;;;;;;;GAYG;AACH,wBAAgB,iBAAiB,CAAC,KAAK,EAAE,kBAAkB,GAAG,kBAAkB,CAO/E;AAED,kFAAkF;AAClF,wBAAgB,GAAG,CAAC,EAAE,EAAE,MAAM,EAAE,EAAE,EAAE,MAAM,GAAG,MAAM,CAOlD;AAED;;;;;;;;;;;;;;;;;GAiBG;AACH,wBAAgB,qBAAqB,CAAC,MAAM,EAAE,QAAQ,CAAC,MAAM,CAAC,MAAM,EAAE,MAAM,CAAC,CAAC,GAAG,MAAM,CAItF;AAED;;;;;;;;;;;;GAYG;AACH,wBAAgB,mBAAmB,CAAC,MAAM,EAAE,QAAQ,CAAC,MAAM,CAAC,MAAM,EAAE,MAAM,CAAC,CAAC,GAAG,MAAM,CAKpF;AAuCD;;;;;;;;;;;;;;;;GAgBG;AACH,eAAO,MAAM,kBAAkB,EAAE,WAiFhC,CAAC;AAEF;;;;;;;;;;GAUG;AACH,eAAO,MAAM,kBAAkB,EAAE,WAoFhC,CAAC;AAEF,oFAAoF;AACpF,eAAO,MAAM,aAAa,EAAE,SAAS,WAAW,EAA6C,CAAC;AAE9F;;;;;;;;;GASG;AACH,eAAO,MAAM,wBAAwB,EAAE,gBACC,CAAC;AAEzC,gBAAgB;AAChB,eAAO,MAAM,wBAAwB,EAAE,gBACC,CAAC;AAEzC,+DAA+D;AAC/D,eAAO,MAAM,eAAe,EAAE,SAAS,gBAAgB,EAGtD,CAAC"}
+{"version":3,"file":"bridges-limits.d.ts","sourceRoot":"","sources":["../../../src/atlas/oscillators/bridges-limits.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;;;GAqBG;AAOH,OAAO,KAAK,EAAE,kBAAkB,EAAE,WAAW,EAAE,gBAAgB,EAAU,MAAM,aAAa,CAAC;AAS7F;;;;;;;;;;;;GAYG;AACH,wBAAgB,iBAAiB,CAAC,KAAK,EAAE,kBAAkB,GAAG,kBAAkB,CAO/E;AAED;;;;;;;;;;;;;;;;;;;GAmBG;AACH,wBAAgB,qBAAqB,CAAC,MAAM,EAAE,QAAQ,CAAC,MAAM,CAAC,MAAM,EAAE,MAAM,CAAC,CAAC,GAAG,MAAM,CAMtF;AAED;;;;;;;;;;;;GAYG;AACH,wBAAgB,mBAAmB,CAAC,MAAM,EAAE,QAAQ,CAAC,MAAM,CAAC,MAAM,EAAE,MAAM,CAAC,CAAC,GAAG,MAAM,CAKpF;AAiCD;;;;;;;;;;;;;;;;GAgBG;AACH,eAAO,MAAM,kBAAkB,EAAE,WAgFhC,CAAC;AAEF;;;;;;;;;;GAUG;AACH,eAAO,MAAM,kBAAkB,EAAE,WAmFhC,CAAC;AAEF,oFAAoF;AACpF,eAAO,MAAM,aAAa,EAAE,SAAS,WAAW,EAA6C,CAAC;AAmB9F,+DAA+D;AAC/D,eAAO,MAAM,eAAe,EAAE,SAAS,gBAAgB,EAGtD,CAAC"}

package/dist/atlas/oscillators/bridges-limits.js

@@ -20,8 +20,10 @@
  *
  * @module atlas/oscillators/bridges-limits
  */
+import { ellipticK } from '@danielsimonjr/mathts-functions';
+import { numberFrom } from '../../core/mathts-read.js';
 import { physjsFormalRef } from '../physjs-ref.js';
-import { deriveRegimeGroups } from '../regime.js';
+import { makeRegime } from '../regime.js';
 import { MissingHorizonError } from '../types.js';
 import { relationContractOf } from './bridges-exact.js';
 import { DAMPING, SPRING_CONSTANT } from './dimensions.js';
@@ -48,18 +50,11 @@ export function makeApproximation(bound) {
     }
     return bound;
 }
-/** Arithmetic–geometric mean, the standard evaluation route for `K`. @internal */
-export function agm(a0, b0) {
-    let a = a0;
-    let b = b0;
-    for (let i = 0; i < 60 && a !== b; i++) {
-        [a, b] = [(a + b) / 2, Math.sqrt(a * b)];
-    }
-    return a;
-}
 /**
  * The EXACT relative period error of the pendulum at amplitude `theta0`:
- * `(2/π) K(sin(θ0/2)) − 1 = 1/AGM(1, cos(θ0/2)) − 1`.
+ * `(2/π) K(sin(θ0/2)) − 1`, with `K` MathTS `ellipticK`, which takes the
+ * PARAMETER `m = k² = sin²(θ0/2)` and evaluates it by the AGM
+ * (`K = π / (2·AGM(1, √(1−m)))`, so this is `1/AGM(1, cos(θ0/2)) − 1`).
  *
  * This is the machine form of `ab-pendulum-linear`'s `delta`. It is SHARP —
  * the bound equals the error it bounds — and it is strictly increasing in
@@ -79,7 +74,9 @@ export function pendulumPeriodErrorAt(params) {
     const { theta0 } = params;
     if (!Number.isFinite(theta0) || Math.abs(theta0) >= Math.PI)
         return Infinity;
-    return 1 / agm(1, Math.cos(theta0 / 2)) - 1;
+    const m = Math.sin(theta0 / 2) ** 2;
+    const k = numberFrom('ellipticK', ellipticK(m));
+    return (2 / Math.PI) * k - 1;
 }
 /**
  * The position-offset bound of `ab-damped-massless` at a point:
@@ -103,14 +100,10 @@ export function dampedOffsetBoundAt(params) {
     return (2 * (1 + Math.abs(v0)) * m) / b;
 }
 /** Pendulum regime: `θ0 ≤ 0.5 rad`, stated on the dimensionless input. @internal */
-const PENDULUM_REGIME = {
-    family: FAMILY,
-    inequalities: [{ group: 'theta0', op: '<=', bound: 0.5, alias: 'θ0 ≤ 0.5 rad' }],
-    groupDefinitions: deriveRegimeGroups(FAMILY, [
-        { name: 'g', dim: ACCELERATION },
-        { name: 'ell', dim: LENGTH },
-    ], ['theta0']),
-};
+const PENDULUM_REGIME = makeRegime(FAMILY, [
+    { name: 'g', dim: ACCELERATION },
+    { name: 'ell', dim: LENGTH },
+], ['theta0'], [{ group: 'theta0', op: '<=', bound: 0.5, alias: 'θ0 ≤ 0.5 rad' }]);
 /**
  * Damped-spring regime: overdamped, `m k / b² < 1/4`.
  *
@@ -120,15 +113,11 @@ const PENDULUM_REGIME = {
  *
  * @internal
  */
-const DAMPED_REGIME = {
-    family: FAMILY,
-    inequalities: [{ group: 'm · b^-2 · k', op: '<', bound: 0.25, alias: 'ζ > 1' }],
-    groupDefinitions: deriveRegimeGroups(FAMILY, [
-        { name: 'm', dim: MASS },
-        { name: 'b', dim: DAMPING },
-        { name: 'k', dim: SPRING_CONSTANT },
-    ], []),
-};
+const DAMPED_REGIME = makeRegime(FAMILY, [
+    { name: 'm', dim: MASS },
+    { name: 'b', dim: DAMPING },
+    { name: 'k', dim: SPRING_CONSTANT },
+], [], [{ group: 'm · b^-2 · k', op: '<', bound: 0.25, alias: 'ζ > 1' }]);
 /**
  * `ab-pendulum-linear` — small-angle pendulum → harmonic spring.
  *
@@ -158,8 +147,8 @@ export const AB_PENDULUM_LINEAR = {
     bound: makeApproximation({
         K: 1,
         // Sup over `θ0 ≤ 0.5`: the error is strictly increasing in |θ0|, so the
-        // edge value IS the supremum. 0.0158525311014... (AGM), independently
-        // reproduced by quadrature of the complete elliptic integral.
+        // edge value IS the supremum. 0.0158525311014... (MathTS ellipticK, by the
+        // AGM), independently reproduced by quadrature of the complete elliptic integral.
         delta: pendulumPeriodErrorAt({ theta0: 0.5 }),
         deltaAt: pendulumPeriodErrorAt,
         deltaAtBasis: 'closed-form',
@@ -194,7 +183,6 @@ export const AB_PENDULUM_LINEAR = {
             witness: 'W7b',
         },
     ],
-    evidence: new Set(['numerically-supported']),
     witnesses: [
         {
             id: 'W7',
@@ -225,7 +213,7 @@ export const AB_PENDULUM_LINEAR = {
     // certifies the TRANSFORMATION (mass mℓ², spring constant mgℓ, ω0² = g/ℓ),
     // not `bound.delta`. Physlib's period results concern `periodFormula`, which
     // its own TODO has not yet tied to the motion. The sanity lemmas instantiate
-    // that imported statement. Keyed by `formal/physjs/manifest.json`.
+    // that imported statement. Keyed by `formal/lean/manifest.json`.
     formalRef: physjsFormalRef('ab-pendulum-linear'),
 };
 /**
@@ -302,7 +290,6 @@ export const AB_DAMPED_MASSLESS = {
             witness: 'W8b',
         },
     ],
-    evidence: new Set(['numerically-supported']),
     witnesses: [
         {
             id: 'W8',
@@ -335,9 +322,9 @@ export const LIMIT_BRIDGES = [AB_PENDULUM_LINEAR, AB_DAMPED_MASSLESS];
  *
  * @internal
  */
-export const CONTRACT_PENDULUM_LINEAR = relationContractOf(AB_PENDULUM_LINEAR);
+const CONTRACT_PENDULUM_LINEAR = relationContractOf(AB_PENDULUM_LINEAR);
 /** @internal */
-export const CONTRACT_DAMPED_MASSLESS = relationContractOf(AB_DAMPED_MASSLESS);
+const CONTRACT_DAMPED_MASSLESS = relationContractOf(AB_DAMPED_MASSLESS);
 /** The two limit contracts, in design-note order. @internal */
 export const LIMIT_CONTRACTS = [
     CONTRACT_PENDULUM_LINEAR,

package/dist/atlas/oscillators/bridges-limits.js.map

@@ -1 +1 @@
-{"version":3,"file":"bridges-limits.js","sourceRoot":"","sources":["../../../src/atlas/oscillators/bridges-limits.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;;;GAqBG;AAEH,OAAO,EAAE,eAAe,EAAE,MAAM,kBAAkB,CAAC;AACnD,OAAO,EAAE,kBAAkB,EAAE,MAAM,cAAc,CAAC;AAClD,OAAO,EAAE,mBAAmB,EAAE,MAAM,aAAa,CAAC;AAElD,OAAO,EAAE,kBAAkB,EAAE,MAAM,oBAAoB,CAAC;AACxD,OAAO,EAAE,OAAO,EAAE,eAAe,EAAE,MAAM,iBAAiB,CAAC;AAC3D,OAAO,EAAE,oBAAoB,EAAE,MAAM,YAAY,CAAC;AAClD,OAAO,EAAE,YAAY,EAAE,MAAM,EAAE,IAAI,EAAE,MAAM,4BAA4B,CAAC;AAExE,MAAM,MAAM,GAAG,aAAa,CAAC;AAC7B,MAAM,SAAS,GAAG,wCAAwC,CAAC;AAE3D;;;;;;;;;;;;GAYG;AACH,MAAM,UAAU,iBAAiB,CAAC,KAAyB;IACzD,IAAI,KAAK,CAAC,OAAO,CAAC,IAAI,EAAE,KAAK,EAAE,EAAE,CAAC;QAChC,MAAM,IAAI,mBAAmB,CAC3B,+EAA+E,CAChF,CAAC;IACJ,CAAC;IACD,OAAO,KAAK,CAAC;AACf,CAAC;AAED,kFAAkF;AAClF,MAAM,UAAU,GAAG,CAAC,EAAU,EAAE,EAAU;IACxC,IAAI,CAAC,GAAG,EAAE,CAAC;IACX,IAAI,CAAC,GAAG,EAAE,CAAC;IACX,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,EAAE,IAAI,CAAC,KAAK,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;QACvC,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,EAAE,IAAI,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;IAC3C,CAAC;IACD,OAAO,CAAC,CAAC;AACX,CAAC;AAED;;;;;;;;;;;;;;;;;GAiBG;AACH,MAAM,UAAU,qBAAqB,CAAC,MAAwC;IAC5E,MAAM,EAAE,MAAM,EAAE,GAAG,MAAM,CAAC;IAC1B,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,MAAM,CAAC,IAAI,IAAI,CAAC,GAAG,CAAC,MAAM,CAAC,IAAI,IAAI,CAAC,EAAE;QAAE,OAAO,QAAQ,CAAC;IAC7E,OAAO,CAAC,GAAG,GAAG,CAAC,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;AAC9C,CAAC;AAED;;;;;;;;;;;;GAYG;AACH,MAAM,UAAU,mBAAmB,CAAC,MAAwC;IAC1E,MAAM,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,GAAG,MAAM,CAAC;IACnC,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC,KAAK,CAAC,MAAM,CAAC,QAAQ,CAAC;QAAE,OAAO,QAAQ,CAAC;IAC/D,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,KAAK,CAAC,IAAI,EAAE,KAAK,CAAC;QAAE,OAAO,QAAQ,CAAC;IACpD,OAAO,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC;AAC1C,CAAC;AAED,oFAAoF;AACpF,MAAM,eAAe,GAAW;IAC9B,MAAM,EAAE,MAAM;IACd,YAAY,EAAE,CAAC,EAAE,KAAK,EAAE,QAAQ,EAAE,EAAE,EAAE,IAAI,EAAE,KAAK,EAAE,GAAG,EAAE,KAAK,EAAE,cAAc,EAAE,CAAC;IAChF,gBAAgB,EAAE,kBAAkB,CAClC,MAAM,EACN;QACE,EAAE,IAAI,EAAE,GAAG,EAAE,GAAG,EAAE,YAAY,EAAE;QAChC,EAAE,IAAI,EAAE,KAAK,EAAE,GAAG,EAAE,MAAM,EAAE;KAC7B,EACD,CAAC,QAAQ,CAAC,CACX;CACF,CAAC;AAEF;;;;;;;;GAQG;AACH,MAAM,aAAa,GAAW;IAC5B,MAAM,EAAE,MAAM;IACd,YAAY,EAAE,CAAC,EAAE,KAAK,EAAE,cAAc,EAAE,EAAE,EAAE,GAAG,EAAE,KAAK,EAAE,IAAI,EAAE,KAAK,EAAE,OAAO,EAAE,CAAC;IAC/E,gBAAgB,EAAE,kBAAkB,CAClC,MAAM,EACN;QACE,EAAE,IAAI,EAAE,GAAG,EAAE,GAAG,EAAE,IAAI,EAAE;QACxB,EAAE,IAAI,EAAE,GAAG,EAAE,GAAG,EAAE,OAAO,EAAE;QAC3B,EAAE,IAAI,EAAE,GAAG,EAAE,GAAG,EAAE,eAAe,EAAE;KACpC,EACD,EAAE,CACH;CACF,CAAC;AAEF;;;;;;;;;;;;;;;;GAgBG;AACH,MAAM,CAAC,MAAM,kBAAkB,GAAgB;IAC7C,EAAE,EAAE,oBAAoB;IACxB,QAAQ,EAAE,eAAe;IACzB,QAAQ,EAAE,CAAC,gBAAgB,CAAC;IAC5B,UAAU,EAAE,cAAc;IAC1B,cAAc,EAAE,uCAAuC;IACvD,SAAS,EAAE,CAAC,iCAAiC,EAAE,qBAAqB,EAAE,wBAAwB,CAAC;IAC/F,eAAe,EAAE,CAAC,wCAAwC,EAAE,gCAAgC,CAAC;IAC7F,cAAc,EAAE,CAAC,cAAc,EAAE,wDAAwD,CAAC;IAC1F,KAAK,EAAE,iBAAiB,CAAC;QACvB,CAAC,EAAE,CAAC;QACJ,wEAAwE;QACxE,sEAAsE;QACtE,8DAA8D;QAC9D,KAAK,EAAE,qBAAqB,CAAC,EAAE,MAAM,EAAE,GAAG,EAAE,CAAC;QAC7C,OAAO,EAAE,qBAAqB;QAC9B,YAAY,EAAE,aAAa;QAC3B,IAAI,EAAE,oBAAoB;QAC1B,MAAM,EAAE,cAAc;QACtB,OAAO,EAAE,8DAA8D;QACvE,sEAAsE;QACtE,yEAAyE;QACzE,yEAAyE;QACzE,wEAAwE;QACxE,uEAAuE;QACvE,uEAAuE;QACvE,kEAAkE;QAClE,+DAA+D;QAC/D,YAAY,EAAE,CAAC,CAAC,EAAE,MAAM,EAAE,EAAE;YAC1B,MAAM,EAAE,EAAE,EAAE,MAAM,EAAE,GAAG,MAAM,CAAC;YAC9B,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,EAAE,CAAC,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,MAAM,CAAC,IAAI,MAAM,KAAK,CAAC;gBAAE,OAAO,KAAK,CAAC;YACnF,OAAO,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC,GAAG,MAAM,IAAI,CAAC,CAAC;QACpC,CAAC;QACD,cAAc,EAAE,cAAc;QAC9B,cAAc,EAAE,SAAS;QACzB,wEAAwE;QACxE,cAAc;QACd,UAAU,EAAE,CAAC,yCAAyC,CAAC;KACxD,CAAC;IACF,MAAM,EAAE,eAAe;IACvB,eAAe,EAAE;QACf;YACE,WAAW,EACT,oFAAoF;gBACpF,yDAAyD;YAC3D,OAAO,EAAE,KAAK;SACf;KACF;IACD,QAAQ,EAAE,IAAI,GAAG,CAAC,CAAC,uBAAuB,CAAC,CAAC;IAC5C,SAAS,EAAE;QACT;YACE,EAAE,EAAE,IAAI;YACR,IAAI,EAAE,SAAS;YACf,IAAI,EAAE,SAAS;YACf,SAAS,EAAE,4EAA4E;SACxF;QACD;YACE,EAAE,EAAE,KAAK;YACT,IAAI,EAAE,SAAS;YACf,IAAI,EAAE,SAAS;YACf,SAAS,EAAE,6CAA6C;SACzD;QACD;YACE,EAAE,EAAE,KAAK;YACT,IAAI,EAAE,SAAS;YACf,IAAI,EAAE,SAAS;YACf,SAAS,EAAE,sEAAsE;SAClF;KACF;IACD,SAAS,EAAE;QACT,oFAAoF;QACpF,iDAAiD;KAClD;IACD,YAAY,EAAE,UAAU;IACxB,4EAA4E;IAC5E,yEAAyE;IACzE,2EAA2E;IAC3E,6EAA6E;IAC7E,6EAA6E;IAC7E,mEAAmE;IACnE,SAAS,EAAE,eAAe,CAAC,oBAAoB,CAAC;CACjD,CAAC;AAEF;;;;;;;;;;GAUG;AACH,MAAM,CAAC,MAAM,kBAAkB,GAAgB;IAC7C,EAAE,EAAE,oBAAoB;IACxB,QAAQ,EAAE,eAAe;IACzB,QAAQ,EAAE,CAAC,qBAAqB,CAAC;IACjC,UAAU,EAAE,mBAAmB;IAC/B,cAAc,EAAE,+BAA+B;IAC/C,SAAS,EAAE,CAAC,sCAAsC,EAAE,wCAAwC,CAAC;IAC7F,eAAe,EAAE;QACf,qCAAqC;QACrC,qEAAqE;QACrE,wBAAwB;KACzB;IACD,cAAc,EAAE;QACd,4BAA4B;QAC5B,4DAA4D;KAC7D;IACD,KAAK,EAAE,iBAAiB,CAAC;QACvB,CAAC,EAAE,CAAC;QACJ,yEAAyE;QACzE,oEAAoE;QACpE,0EAA0E;QAC1E,0EAA0E;QAC1E,uEAAuE;QACvE,0EAA0E;QAC1E,mEAAmE;QACnE,kBAAkB;QAClB,wEAAwE;QACxE,4EAA4E;QAC5E,yEAAyE;QACzE,kEAAkE;QAClE,2EAA2E;QAC3E,6EAA6E;QAC7E,2DAA2D;QAC3D,yDAAyD;QACzD,yDAAyD;QACzD,KAAK,EAAE,mBAAmB,CAAC,EAAE,CAAC,EAAE,IAAI,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,EAAE,EAAE,CAAC,EAAE,EAAE,EAAE,CAAC,EAAE,CAAC;QACjE,OAAO,EAAE,mBAAmB;QAC5B,sFAAsF;QACtF,YAAY,EAAE,uBAAuB;QACrC,IAAI,EAAE,mCAAmC;QACzC,MAAM,EAAE,+DAA+D;QACvE,OAAO,EAAE,wCAAwC;QACjD,YAAY,EAAE,CAAC,CAAC,EAAE,MAAM,EAAE,EAAE;YAC1B,MAAM,EAAE,CAAC,EAAE,CAAC,EAAE,GAAG,MAAM,CAAC;YACxB,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,IAAI,CAAC,KAAK,CAAC;gBAAE,OAAO,KAAK,CAAC;YACxE,OAAO,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC;QAC1B,CAAC;QACD,cAAc,EAAE,0DAA0D;QAC1E,cAAc,EAAE,UAAU;QAC1B,UAAU,EAAE;YACV,+DAA+D;YAC/D,0BAA0B;SAC3B;KACF,CAAC;IACF,MAAM,EAAE,aAAa;IACrB,eAAe,EAAE;QACf;YACE,WAAW,EACT,8EAA8E;gBAC9E,wEAAwE;YAC1E,OAAO,EAAE,KAAK;SACf;KACF;IACD,QAAQ,EAAE,IAAI,GAAG,CAAC,CAAC,uBAAuB,CAAC,CAAC;IAC5C,SAAS,EAAE;QACT;YACE,EAAE,EAAE,IAAI;YACR,IAAI,EAAE,SAAS;YACf,IAAI,EAAE,SAAS;YACf,SAAS,EAAE,wEAAwE;SACpF;QACD;YACE,EAAE,EAAE,KAAK;YACT,IAAI,EAAE,SAAS;YACf,IAAI,EAAE,SAAS;YACf,SAAS,EACP,2FAA2F;SAC9F;KACF;IACD,SAAS,EAAE;QACT,4FAA4F;QAC5F,mFAAmF;KACpF;IACD,YAAY,EAAE,UAAU;CACzB,CAAC;AAEF,oFAAoF;AACpF,MAAM,CAAC,MAAM,aAAa,GAA2B,CAAC,kBAAkB,EAAE,kBAAkB,CAAC,CAAC;AAE9F;;;;;;;;;GASG;AACH,MAAM,CAAC,MAAM,wBAAwB,GACnC,kBAAkB,CAAC,kBAAkB,CAAC,CAAC;AAEzC,gBAAgB;AAChB,MAAM,CAAC,MAAM,wBAAwB,GACnC,kBAAkB,CAAC,kBAAkB,CAAC,CAAC;AAEzC,+DAA+D;AAC/D,MAAM,CAAC,MAAM,eAAe,GAAgC;IAC1D,wBAAwB;IACxB,wBAAwB;CACzB,CAAC"}
+{"version":3,"file":"bridges-limits.js","sourceRoot":"","sources":["../../../src/atlas/oscillators/bridges-limits.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;;;GAqBG;AAEH,OAAO,EAAE,SAAS,EAAE,MAAM,iCAAiC,CAAC;AAC5D,OAAO,EAAE,UAAU,EAAE,MAAM,2BAA2B,CAAC;AACvD,OAAO,EAAE,eAAe,EAAE,MAAM,kBAAkB,CAAC;AACnD,OAAO,EAAE,UAAU,EAAE,MAAM,cAAc,CAAC;AAC1C,OAAO,EAAE,mBAAmB,EAAE,MAAM,aAAa,CAAC;AAElD,OAAO,EAAE,kBAAkB,EAAE,MAAM,oBAAoB,CAAC;AACxD,OAAO,EAAE,OAAO,EAAE,eAAe,EAAE,MAAM,iBAAiB,CAAC;AAC3D,OAAO,EAAE,oBAAoB,EAAE,MAAM,YAAY,CAAC;AAClD,OAAO,EAAE,YAAY,EAAE,MAAM,EAAE,IAAI,EAAE,MAAM,4BAA4B,CAAC;AAExE,MAAM,MAAM,GAAG,aAAa,CAAC;AAC7B,MAAM,SAAS,GAAG,wCAAwC,CAAC;AAE3D;;;;;;;;;;;;GAYG;AACH,MAAM,UAAU,iBAAiB,CAAC,KAAyB;IACzD,IAAI,KAAK,CAAC,OAAO,CAAC,IAAI,EAAE,KAAK,EAAE,EAAE,CAAC;QAChC,MAAM,IAAI,mBAAmB,CAC3B,+EAA+E,CAChF,CAAC;IACJ,CAAC;IACD,OAAO,KAAK,CAAC;AACf,CAAC;AAED;;;;;;;;;;;;;;;;;;;GAmBG;AACH,MAAM,UAAU,qBAAqB,CAAC,MAAwC;IAC5E,MAAM,EAAE,MAAM,EAAE,GAAG,MAAM,CAAC;IAC1B,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,MAAM,CAAC,IAAI,IAAI,CAAC,GAAG,CAAC,MAAM,CAAC,IAAI,IAAI,CAAC,EAAE;QAAE,OAAO,QAAQ,CAAC;IAC7E,MAAM,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,MAAM,GAAG,CAAC,CAAC,IAAI,CAAC,CAAC;IACpC,MAAM,CAAC,GAAG,UAAU,CAAC,WAAW,EAAE,SAAS,CAAC,CAAC,CAAC,CAAC,CAAC;IAChD,OAAO,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;AAC/B,CAAC;AAED;;;;;;;;;;;;GAYG;AACH,MAAM,UAAU,mBAAmB,CAAC,MAAwC;IAC1E,MAAM,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,GAAG,MAAM,CAAC;IACnC,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC,KAAK,CAAC,MAAM,CAAC,QAAQ,CAAC;QAAE,OAAO,QAAQ,CAAC;IAC/D,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,KAAK,CAAC,IAAI,EAAE,KAAK,CAAC;QAAE,OAAO,QAAQ,CAAC;IACpD,OAAO,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC;AAC1C,CAAC;AAED,oFAAoF;AACpF,MAAM,eAAe,GAAW,UAAU,CACxC,MAAM,EACN;IACE,EAAE,IAAI,EAAE,GAAG,EAAE,GAAG,EAAE,YAAY,EAAE;IAChC,EAAE,IAAI,EAAE,KAAK,EAAE,GAAG,EAAE,MAAM,EAAE;CAC7B,EACD,CAAC,QAAQ,CAAC,EACV,CAAC,EAAE,KAAK,EAAE,QAAQ,EAAE,EAAE,EAAE,IAAI,EAAE,KAAK,EAAE,GAAG,EAAE,KAAK,EAAE,cAAc,EAAE,CAAC,CACnE,CAAC;AAEF;;;;;;;;GAQG;AACH,MAAM,aAAa,GAAW,UAAU,CACtC,MAAM,EACN;IACE,EAAE,IAAI,EAAE,GAAG,EAAE,GAAG,EAAE,IAAI,EAAE;IACxB,EAAE,IAAI,EAAE,GAAG,EAAE,GAAG,EAAE,OAAO,EAAE;IAC3B,EAAE,IAAI,EAAE,GAAG,EAAE,GAAG,EAAE,eAAe,EAAE;CACpC,EACD,EAAE,EACF,CAAC,EAAE,KAAK,EAAE,cAAc,EAAE,EAAE,EAAE,GAAG,EAAE,KAAK,EAAE,IAAI,EAAE,KAAK,EAAE,OAAO,EAAE,CAAC,CAClE,CAAC;AAEF;;;;;;;;;;;;;;;;GAgBG;AACH,MAAM,CAAC,MAAM,kBAAkB,GAAgB;IAC7C,EAAE,EAAE,oBAAoB;IACxB,QAAQ,EAAE,eAAe;IACzB,QAAQ,EAAE,CAAC,gBAAgB,CAAC;IAC5B,UAAU,EAAE,cAAc;IAC1B,cAAc,EAAE,uCAAuC;IACvD,SAAS,EAAE,CAAC,iCAAiC,EAAE,qBAAqB,EAAE,wBAAwB,CAAC;IAC/F,eAAe,EAAE,CAAC,wCAAwC,EAAE,gCAAgC,CAAC;IAC7F,cAAc,EAAE,CAAC,cAAc,EAAE,wDAAwD,CAAC;IAC1F,KAAK,EAAE,iBAAiB,CAAC;QACvB,CAAC,EAAE,CAAC;QACJ,wEAAwE;QACxE,2EAA2E;QAC3E,kFAAkF;QAClF,KAAK,EAAE,qBAAqB,CAAC,EAAE,MAAM,EAAE,GAAG,EAAE,CAAC;QAC7C,OAAO,EAAE,qBAAqB;QAC9B,YAAY,EAAE,aAAa;QAC3B,IAAI,EAAE,oBAAoB;QAC1B,MAAM,EAAE,cAAc;QACtB,OAAO,EAAE,8DAA8D;QACvE,sEAAsE;QACtE,yEAAyE;QACzE,yEAAyE;QACzE,wEAAwE;QACxE,uEAAuE;QACvE,uEAAuE;QACvE,kEAAkE;QAClE,+DAA+D;QAC/D,YAAY,EAAE,CAAC,CAAC,EAAE,MAAM,EAAE,EAAE;YAC1B,MAAM,EAAE,EAAE,EAAE,MAAM,EAAE,GAAG,MAAM,CAAC;YAC9B,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,EAAE,CAAC,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,MAAM,CAAC,IAAI,MAAM,KAAK,CAAC;gBAAE,OAAO,KAAK,CAAC;YACnF,OAAO,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC,GAAG,MAAM,IAAI,CAAC,CAAC;QACpC,CAAC;QACD,cAAc,EAAE,cAAc;QAC9B,cAAc,EAAE,SAAS;QACzB,wEAAwE;QACxE,cAAc;QACd,UAAU,EAAE,CAAC,yCAAyC,CAAC;KACxD,CAAC;IACF,MAAM,EAAE,eAAe;IACvB,eAAe,EAAE;QACf;YACE,WAAW,EACT,oFAAoF;gBACpF,yDAAyD;YAC3D,OAAO,EAAE,KAAK;SACf;KACF;IACD,SAAS,EAAE;QACT;YACE,EAAE,EAAE,IAAI;YACR,IAAI,EAAE,SAAS;YACf,IAAI,EAAE,SAAS;YACf,SAAS,EAAE,4EAA4E;SACxF;QACD;YACE,EAAE,EAAE,KAAK;YACT,IAAI,EAAE,SAAS;YACf,IAAI,EAAE,SAAS;YACf,SAAS,EAAE,6CAA6C;SACzD;QACD;YACE,EAAE,EAAE,KAAK;YACT,IAAI,EAAE,SAAS;YACf,IAAI,EAAE,SAAS;YACf,SAAS,EAAE,sEAAsE;SAClF;KACF;IACD,SAAS,EAAE;QACT,oFAAoF;QACpF,iDAAiD;KAClD;IACD,YAAY,EAAE,UAAU;IACxB,4EAA4E;IAC5E,yEAAyE;IACzE,2EAA2E;IAC3E,6EAA6E;IAC7E,6EAA6E;IAC7E,iEAAiE;IACjE,SAAS,EAAE,eAAe,CAAC,oBAAoB,CAAC;CACjD,CAAC;AAEF;;;;;;;;;;GAUG;AACH,MAAM,CAAC,MAAM,kBAAkB,GAAgB;IAC7C,EAAE,EAAE,oBAAoB;IACxB,QAAQ,EAAE,eAAe;IACzB,QAAQ,EAAE,CAAC,qBAAqB,CAAC;IACjC,UAAU,EAAE,mBAAmB;IAC/B,cAAc,EAAE,+BAA+B;IAC/C,SAAS,EAAE,CAAC,sCAAsC,EAAE,wCAAwC,CAAC;IAC7F,eAAe,EAAE;QACf,qCAAqC;QACrC,qEAAqE;QACrE,wBAAwB;KACzB;IACD,cAAc,EAAE;QACd,4BAA4B;QAC5B,4DAA4D;KAC7D;IACD,KAAK,EAAE,iBAAiB,CAAC;QACvB,CAAC,EAAE,CAAC;QACJ,yEAAyE;QACzE,oEAAoE;QACpE,0EAA0E;QAC1E,0EAA0E;QAC1E,uEAAuE;QACvE,0EAA0E;QAC1E,mEAAmE;QACnE,kBAAkB;QAClB,wEAAwE;QACxE,4EAA4E;QAC5E,yEAAyE;QACzE,kEAAkE;QAClE,2EAA2E;QAC3E,6EAA6E;QAC7E,2DAA2D;QAC3D,yDAAyD;QACzD,yDAAyD;QACzD,KAAK,EAAE,mBAAmB,CAAC,EAAE,CAAC,EAAE,IAAI,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,EAAE,EAAE,CAAC,EAAE,EAAE,EAAE,CAAC,EAAE,CAAC;QACjE,OAAO,EAAE,mBAAmB;QAC5B,sFAAsF;QACtF,YAAY,EAAE,uBAAuB;QACrC,IAAI,EAAE,mCAAmC;QACzC,MAAM,EAAE,+DAA+D;QACvE,OAAO,EAAE,wCAAwC;QACjD,YAAY,EAAE,CAAC,CAAC,EAAE,MAAM,EAAE,EAAE;YAC1B,MAAM,EAAE,CAAC,EAAE,CAAC,EAAE,GAAG,MAAM,CAAC;YACxB,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,IAAI,CAAC,KAAK,CAAC;gBAAE,OAAO,KAAK,CAAC;YACxE,OAAO,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC;QAC1B,CAAC;QACD,cAAc,EAAE,0DAA0D;QAC1E,cAAc,EAAE,UAAU;QAC1B,UAAU,EAAE;YACV,+DAA+D;YAC/D,0BAA0B;SAC3B;KACF,CAAC;IACF,MAAM,EAAE,aAAa;IACrB,eAAe,EAAE;QACf;YACE,WAAW,EACT,8EAA8E;gBAC9E,wEAAwE;YAC1E,OAAO,EAAE,KAAK;SACf;KACF;IACD,SAAS,EAAE;QACT;YACE,EAAE,EAAE,IAAI;YACR,IAAI,EAAE,SAAS;YACf,IAAI,EAAE,SAAS;YACf,SAAS,EAAE,wEAAwE;SACpF;QACD;YACE,EAAE,EAAE,KAAK;YACT,IAAI,EAAE,SAAS;YACf,IAAI,EAAE,SAAS;YACf,SAAS,EACP,2FAA2F;SAC9F;KACF;IACD,SAAS,EAAE;QACT,4FAA4F;QAC5F,mFAAmF;KACpF;IACD,YAAY,EAAE,UAAU;CACzB,CAAC;AAEF,oFAAoF;AACpF,MAAM,CAAC,MAAM,aAAa,GAA2B,CAAC,kBAAkB,EAAE,kBAAkB,CAAC,CAAC;AAE9F;;;;;;;;;GASG;AACH,MAAM,wBAAwB,GAC5B,kBAAkB,CAAC,kBAAkB,CAAC,CAAC;AAEzC,gBAAgB;AAChB,MAAM,wBAAwB,GAC5B,kBAAkB,CAAC,kBAAkB,CAAC,CAAC;AAEzC,+DAA+D;AAC/D,MAAM,CAAC,MAAM,eAAe,GAAgC;IAC1D,wBAAwB;IACxB,wBAAwB;CACzB,CAAC","sourcesContent":["/**\n * The two Phase 0 APPROXIMATION bridges — design note §3, bridges 3 and 4.\n *\n * - `ab-pendulum-linear` — a REGULAR limit. `θ'' + ω0² sin θ = 0` reduces to\n *   `θ'' + ω0² θ = 0` as `θ0 → 0`. The order of the system does not drop, the\n *   error is `O(θ0²)` in the relative period, and the approximation is\n *   NON-UNIFORM in time: a fixed relative period error accumulates phase\n *   without bound, so the bound carries the horizon `t ≪ 16 T0/θ0²`.\n * - `ab-damped-massless` — a SINGULAR limit. `m x'' + b x' + k x = 0` reduces\n *   to `b x' + k x = 0` as `m → 0`. The order DROPS from two to one, one\n *   initial condition is lost, and the reduced model cannot satisfy `x'(0)`.\n *   The bound therefore holds only OUTSIDE the boundary layer, `t ≥ 5 m/b`.\n *\n * On the position offset of the singular bridge: `x_full − x_red` is an\n * `O((1 + |v0|) m)` deviation that PERSISTS THROUGHOUT THE OUTER REGION and\n * decays only on the slow time scale `b/k`. It is the matched-asymptotics\n * correction, not a boundary-layer transient, and waiting out the fast mode\n * does not remove it. (Design note §9 YELLOW (5) — \"non-vanishing\" was the\n * imprecise wording this replaces: the offset does go to zero as `t → ∞`.)\n *\n * @module atlas/oscillators/bridges-limits\n */\n\nimport { ellipticK } from '@danielsimonjr/mathts-functions';\nimport { numberFrom } from '../../core/mathts-read.js';\nimport { physjsFormalRef } from '../physjs-ref.js';\nimport { makeRegime } from '../regime.js';\nimport { MissingHorizonError } from '../types.js';\nimport type { ApproximationBound, AtlasBridge, RelationContract, Regime } from '../types.js';\nimport { relationContractOf } from './bridges-exact.js';\nimport { DAMPING, SPRING_CONSTANT } from './dimensions.js';\nimport { RELATIVE_PERIOD_NORM } from './norms.js';\nimport { ACCELERATION, LENGTH, MASS } from '../../dimensional/types.js';\n\nconst FAMILY = 'oscillators';\nconst TEST_FILE = 'tests/atlas/oscillators-limits.test.ts';\n\n/**\n * Build an `ApproximationBound`, rejecting one with no prose horizon.\n *\n * Every `approximation` bridge in this module is constructed through here, so\n * the mandatory-horizon rule of `ApproximationBound` has something that FAILS\n * rather than a comment asking for compliance. `horizonHolds` is mandatory by\n * the type; the prose `horizon` is the field a caller can satisfy with `''`,\n * so it is the one checked at run time. `uniformity: null` is accepted: not\n * yet analysed is a real state, and the refusal lives in `boundPath`.\n *\n * @throws MissingHorizonError if `horizon` is empty or whitespace.\n * @internal\n */\nexport function makeApproximation(bound: ApproximationBound): ApproximationBound {\n  if (bound.horizon.trim() === '') {\n    throw new MissingHorizonError(\n      'an approximation bound requires a non-empty prose horizon beside horizonHolds',\n    );\n  }\n  return bound;\n}\n\n/**\n * The EXACT relative period error of the pendulum at amplitude `theta0`:\n * `(2/π) K(sin(θ0/2)) − 1`, with `K` MathTS `ellipticK`, which takes the\n * PARAMETER `m = k² = sin²(θ0/2)` and evaluates it by the AGM\n * (`K = π / (2·AGM(1, √(1−m)))`, so this is `1/AGM(1, cos(θ0/2)) − 1`).\n *\n * This is the machine form of `ab-pendulum-linear`'s `delta`. It is SHARP —\n * the bound equals the error it bounds — and it is strictly increasing in\n * `|θ0|`, so its value at the edge of a declared amplitude range is the\n * supremum over that range. That last property is what lets the record state\n * one scalar for the whole domain honestly.\n *\n * The series `θ0²/16` is the LEADING TERM of this quantity, not an upper bound\n * on it: at the domain edge `θ0 = 0.5` the series gives `0.0156250` while the\n * exact error is `0.0158525311`, so the series UNDERSTATES by 1.456%. A record\n * that declares the series value is violated at its own boundary.\n *\n * @returns `Infinity` outside `|θ0| < π`, where no finite period error exists.\n * @internal\n */\nexport function pendulumPeriodErrorAt(params: Readonly<Record<string, number>>): number {\n  const { theta0 } = params;\n  if (!Number.isFinite(theta0) || Math.abs(theta0) >= Math.PI) return Infinity;\n  const m = Math.sin(theta0 / 2) ** 2;\n  const k = numberFrom('ellipticK', ellipticK(m));\n  return (2 / Math.PI) * k - 1;\n}\n\n/**\n * The position-offset bound of `ab-damped-massless` at a point:\n * `2 (1 + |v0|) m / b`, the formula this record's docstring has always\n * carried. It is a bound ONLY at the witness normalisation `b = k = x0 = 1`:\n * the leading offset is `(m/b)(v0 + k x0/b)`, and the formula has no `k` or\n * `x0` in it. At `m = 1e-3, b = 1, k = 100` (inside `m k/b² < 1/4`) the true\n * error is more than 20 times the formula.\n *\n * @returns `Infinity` unless `m`, `b`, `k`, `x0` and `v0` are all finite and\n *   `b = k = x0 = 1` — a missing parameter would otherwise silently return a\n *   bound for a normalisation the caller did not ask about.\n * @internal\n */\nexport function dampedOffsetBoundAt(params: Readonly<Record<string, number>>): number {\n  const { m, b, k, x0, v0 } = params;\n  if (![m, b, k, x0, v0].every(Number.isFinite)) return Infinity;\n  if (b !== 1 || k !== 1 || x0 !== 1) return Infinity;\n  return (2 * (1 + Math.abs(v0)) * m) / b;\n}\n\n/** Pendulum regime: `θ0 ≤ 0.5 rad`, stated on the dimensionless input. @internal */\nconst PENDULUM_REGIME: Regime = makeRegime(\n  FAMILY,\n  [\n    { name: 'g', dim: ACCELERATION },\n    { name: 'ell', dim: LENGTH },\n  ],\n  ['theta0'],\n  [{ group: 'theta0', op: '<=', bound: 0.5, alias: 'θ0 ≤ 0.5 rad' }],\n);\n\n/**\n * Damped-spring regime: overdamped, `m k / b² < 1/4`.\n *\n * The inequality is written on the π-group, not on `ζ = ½ (m k / b²)^(−½)`,\n * which is a derived quantity and not a group (`RegimeInequality` doc).\n * `m k / b² < 1/4` is exactly `ζ > 1`, recorded as the display alias.\n *\n * @internal\n */\nconst DAMPED_REGIME: Regime = makeRegime(\n  FAMILY,\n  [\n    { name: 'm', dim: MASS },\n    { name: 'b', dim: DAMPING },\n    { name: 'k', dim: SPRING_CONSTANT },\n  ],\n  [],\n  [{ group: 'm · b^-2 · k', op: '<', bound: 0.25, alias: 'ζ > 1' }],\n);\n\n/**\n * `ab-pendulum-linear` — small-angle pendulum → harmonic spring.\n *\n * `θ0²/16` is the leading term of the exact relative period error\n * `(2/π) K(sin(θ0/2)) − 1`. At `θ0 = 0.2` the residual is `5.744e-6`, which\n * the next series term `11 θ0⁴/3072 = 5.729e-6` accounts for to `1.5e-8`.\n *\n * `delta` is the EXACT error at the edge of the declared range, not that\n * series term. The record previously declared `0.5²/16 = 0.0156250`, which is\n * the series evaluated at `θ0 = 0.5` — and the exact error there is\n * `0.0158525311`, so the declared bound was VIOLATED AT ITS OWN BOUNDARY by\n * 1.456%. A truncated series is an approximation of the error, not a bound on\n * it, and the sign of its omitted tail decides which. Here the tail is\n * positive, so the series sits below what it was asked to cover.\n *\n * @internal\n */\nexport const AB_PENDULUM_LINEAR: AtlasBridge = {\n  id: 'ab-pendulum-linear',\n  relation: 'approximation',\n  premises: ['model-pendulum'],\n  conclusion: 'model-spring',\n  transformation: 'sin θ → θ, with x ↔ ℓ θ and ω0² = g/ℓ',\n  preserves: ['harmonic frequency ω0 to O(θ0²)', 'energy conservation', 'time-reversal symmetry'],\n  doesNotPreserve: ['the amplitude dependence of the period', 'phase over times t ≳ 16 T0/θ0²'],\n  sideConditions: ['θ0 ≤ 0.5 rad', 'the bound is a PERIOD error and is not uniform in time'],\n  bound: makeApproximation({\n    K: 1,\n    // Sup over `θ0 ≤ 0.5`: the error is strictly increasing in |θ0|, so the\n    // edge value IS the supremum. 0.0158525311014... (MathTS ellipticK, by the\n    // AGM), independently reproduced by quadrature of the complete elliptic integral.\n    delta: pendulumPeriodErrorAt({ theta0: 0.5 }),\n    deltaAt: pendulumPeriodErrorAt,\n    deltaAtBasis: 'closed-form',\n    norm: RELATIVE_PERIOD_NORM,\n    domain: 'θ0 ≤ 0.5 rad',\n    horizon: 't ≪ 16 T0/θ0²; machine form t < 4 T0/θ0², the π/2-drift time',\n    // ⚠ The machine form is a QUARTER of the prose scale, and it must be.\n    // `16 T0/θ0²` is the time at which the accumulated phase drift reaches a\n    // FULL 2π (measured: 400 T0 at θ0 = 0.2, drift 6.2976 rad) — by then the\n    // approximation has lapped, so `t < 16 T0/θ0²` is not a horizon at all.\n    // The `≪` is what carries the whole claim, and `<` does not render it.\n    // `4 T0/θ0²` is the π/2-drift time, which is exactly what W7b measures\n    // independently (99.77 cycles at θ0 = 0.2). It is also what makes\n    // `horizonHolds(200 T0)` false and `horizonHolds(10 T0)` true.\n    horizonHolds: (t, params) => {\n      const { T0, theta0 } = params;\n      if (!Number.isFinite(T0) || !Number.isFinite(theta0) || theta0 === 0) return false;\n      return t < (4 * T0) / theta0 ** 2;\n    },\n    parameterRange: 'θ0 ≤ 0.5 rad',\n    limitCharacter: 'regular',\n    // A period error. The record says it is not uniform in time, so time is\n    // not listed.\n    uniformity: ['one period, for θ0 in the stated domain'],\n  }),\n  regime: PENDULUM_REGIME,\n  counterexamples: [\n    {\n      description:\n        'the approximation is non-uniform in time: at θ0 = 0.2 the phase drift reaches π/2 ' +\n        'after ~100 cycles, however small the per-cycle error is',\n      witness: 'W7b',\n    },\n  ],\n  witnesses: [\n    {\n      id: 'W7',\n      kind: 'numeric',\n      test: TEST_FILE,\n      tolerance: 'T/T0 − 1 ∈ [0.002505, 0.002507]; residual ∈ [5.70e-6, 5.76e-6] at θ0 = 0.2',\n    },\n    {\n      id: 'W7b',\n      kind: 'numeric',\n      test: TEST_FILE,\n      tolerance: 'cycles to π/2 drift ∈ [99, 101] at θ0 = 0.2',\n    },\n    {\n      id: 'W7c',\n      kind: 'numeric',\n      test: TEST_FILE,\n      tolerance: 'RK4 zero-crossing lag within 0.05° of the elliptic prediction 89.98°',\n    },\n  ],\n  citations: [\n    'Landau & Lifshitz, Mechanics §11 (pendulum period as a complete elliptic integral)',\n    'Abramowitz & Stegun §17.6 (AGM evaluation of K)',\n  ],\n  reviewStatus: 'proposed',\n  // Milestone 1. The reference is the PhysJS theorem, which imports Physlib's\n  // `ClassicalMechanics.SimplePendulum.linearizedEquationOfMotion_iff`. It\n  // certifies the TRANSFORMATION (mass mℓ², spring constant mgℓ, ω0² = g/ℓ),\n  // not `bound.delta`. Physlib's period results concern `periodFormula`, which\n  // its own TODO has not yet tied to the motion. The sanity lemmas instantiate\n  // that imported statement. Keyed by `formal/lean/manifest.json`.\n  formalRef: physjsFormalRef('ab-pendulum-linear'),\n};\n\n/**\n * `ab-damped-massless` — the SINGULAR `m → 0` limit of the damped spring.\n *\n * `deltaAt = 2(1 + |v0|) m/b` bounds the position offset outside the boundary\n * layer, and `delta` is its supremum over the declared range rather than its\n * value at one fixture. The lost initial condition is visible in the VELOCITY: inside the\n * layer the reduced model's `x'` is wrong by O(1/m·…), and it is only after\n * `t ≈ 5 m/b` that the fast mode has decayed.\n *\n * @internal\n */\nexport const AB_DAMPED_MASSLESS: AtlasBridge = {\n  id: 'ab-damped-massless',\n  relation: 'approximation',\n  premises: ['model-damped-spring'],\n  conclusion: 'model-first-order',\n  transformation: 'm → 0, dropping the m x″ term',\n  preserves: ['the slow relaxation rate k/b to O(m)', 'the sign and monotonicity of the decay'],\n  doesNotPreserve: [\n    'the order of the system (two → one)',\n    \"the initial condition x'(0), which the reduced model cannot satisfy\",\n    'the fast mode r ≈ −b/m',\n  ],\n  sideConditions: [\n    'overdamped: m k / b² < 1/4',\n    'the bound holds only outside the boundary layer, t ≥ 5 m/b',\n  ],\n  bound: makeApproximation({\n    K: 1,\n    // Sup over the DECLARED range, not the fixture. `2(1+|v0|)m/b` increases\n    // in both m and |v0|, and the declared range is `m k/b² < 1/4` with\n    // `|v0| ≤ 5`; at the witness normalisation `b = k = 1` that is `m < 1/4`,\n    // so the supremum is `2·(1+5)·(1/4) = 3`, approached at the open edge and\n    // not attained. The record previously declared `2·(1+5)·1e-3 = 0.012`,\n    // the same formula frozen at the ONE fixture mass `m = 1e-3` — true error\n    // there is 5.96e-3, but at m = 0.24 it is 5.19e-1, forty times the\n    // declared bound.\n    // KNOWN LOOSENESS (0.47.0 persona re-test, L9 residual; documented, not\n    // tightened): for x0 = 1 the measured sup of |x − x_reduced| over the whole\n    // declared range (m < 1/4, |v0| ≤ 5, t ≥ 5m/b) is 0.524, at m ≈ 0.21 and\n    // v0 = +5, and it does not grow toward critical damping (0.514 at\n    // m = 0.249999). So δ = 3 holds but is 5.7× loose, and it exceeds the unit\n    // signal it bounds. A tighter δ needs a bound with its own basis, not a grid\n    // maximum; until then `upt path` withholds the point bound\n    // (deltaAtBasis below). Independent check: the persona's\n    // checks/recheck_047.py and a 115 × 2 grid over (m, ±5).\n    delta: dampedOffsetBoundAt({ m: 0.25, b: 1, k: 1, x0: 1, v0: 5 }),\n    deltaAt: dampedOffsetBoundAt,\n    // Witness W8b supports it numerically for m ∈ {1e-1, 1e-2, 1e-3}; no proof covers it.\n    deltaAtBasis: 'numerically-supported',\n    norm: 'sup |x − x_reduced| for t ≥ 5 m/b',\n    domain: 't ≥ 5 m/b, overdamped, at the witness normalisation b = k = 1',\n    horizon: 't ≥ 5 m/b (outside the boundary layer)',\n    horizonHolds: (t, params) => {\n      const { m, b } = params;\n      if (!Number.isFinite(m) || !Number.isFinite(b) || b === 0) return false;\n      return t >= (5 * m) / b;\n    },\n    parameterRange: 'm k / b² < 1/4, |v0| ≤ 5; with b = k = 1 this is m < 1/4',\n    limitCharacter: 'singular',\n    uniformity: [\n      't ≥ 5 m/b, overdamped, at the witness normalisation b = k = 1',\n      'm k / b² < 1/4, |v0| ≤ 5',\n    ],\n  }),\n  regime: DAMPED_REGIME,\n  counterexamples: [\n    {\n      description:\n        \"inside the boundary layer the reduced model's velocity is wrong by O(1): at \" +\n        \"t = 0.5 m/b with v0 = 5 the velocity error is 3.6, against x'_red = −1\",\n      witness: 'W8b',\n    },\n  ],\n  witnesses: [\n    {\n      id: 'W8',\n      kind: 'numeric',\n      test: TEST_FILE,\n      tolerance: '|r_slow + k/b| < 2m and |r_fast·m + b| < 2m for m ∈ {1e-1, 1e-2, 1e-3}',\n    },\n    {\n      id: 'W8b',\n      kind: 'numeric',\n      test: TEST_FILE,\n      tolerance:\n        \"|x_full − x_red| < 2(1+v0)·m for t ≥ 5 m/b; |x'| error > 1 at 0.5 m/b and < 0.05 at 5 m/b\",\n    },\n  ],\n  citations: [\n    \"Kevorkian & Cole, Multiple Scale and Singular Perturbation Methods §2 (Tikhonov's theorem)\",\n    'Verhulst, Methods and Applications of Singular Perturbations §8 (boundary layers)',\n  ],\n  reviewStatus: 'proposed',\n};\n\n/** The two approximation bridges of this module, in design-note order. @internal */\nexport const LIMIT_BRIDGES: readonly AtlasBridge[] = [AB_PENDULUM_LINEAR, AB_DAMPED_MASSLESS];\n\n/**\n * The two approximation bridges re-registered through `RelationContract`.\n *\n * The union REQUIRES `bound` on an `approximation`, and `ApproximationBound`\n * requires its own machine `horizonHolds`, so these two lines are what turn\n * Sprint 0's doc-comment requirement into a compile-and-throw one. Both bounds\n * are the records' own, measured ones — nothing is restated here.\n *\n * @internal\n */\nconst CONTRACT_PENDULUM_LINEAR: RelationContract =\n  relationContractOf(AB_PENDULUM_LINEAR);\n\n/** @internal */\nconst CONTRACT_DAMPED_MASSLESS: RelationContract =\n  relationContractOf(AB_DAMPED_MASSLESS);\n\n/** The two limit contracts, in design-note order. @internal */\nexport const LIMIT_CONTRACTS: readonly RelationContract[] = [\n  CONTRACT_PENDULUM_LINEAR,\n  CONTRACT_DAMPED_MASSLESS,\n];\n"]}

package/dist/atlas/oscillators/dimensions.d.ts

@@ -1,32 +1,20 @@
 /**
- * Dimension constants the oscillator pilot needs that the canonical entry
- * files declare MODULE-LOCAL and do not export.
+ * Dimension constants the oscillator pilot needs.
  *
- * They are redefined here with the same `dim(...)` arguments, quoted from
- * source in design note §10 Q2. `tests/atlas/models.test.ts` asserts the
- * values, so a future divergence in an entry file fails a test.
+ * Capacitance, inductance, resistance and the spring constant are the named
+ * shapes of `src/dimensional/types.ts`, re-exported here so the family's
+ * modules keep one import site; `tests/atlas/models.test.ts` asserts they are
+ * those objects, not copies. `DAMPING` and `CUBIC_STIFFNESS` are this family's
+ * own and are defined here.
  *
  * `FREQUENCY`, `MASS`, `LENGTH`, `VELOCITY` and `ACCELERATION` are importable
- * from `src/dimensional/types.ts` and are NOT redefined here.
+ * from `src/dimensional/types.ts` and are NOT re-exported here.
  *
  * @module atlas/oscillators/dimensions
  */
 import type { Dimension } from '../../dimensional/types.js';
-/** farad, `L^-2 M^-1 T^4 I^2` — `src/canonical/entries/electromagnetism.ts:42`. @internal */
-export declare const CAPACITANCE: Dimension;
-/** henry, `L^2 M T^-2 I^-2` — `src/canonical/entries/electromagnetism.ts:43`. @internal */
-export declare const INDUCTANCE: Dimension;
-/** ohm, `L^2 M T^-3 I^-2` — `src/canonical/entries/electromagnetism.ts:40`. @internal */
-export declare const RESISTANCE: Dimension;
-/**
- * `M T^-2` — defined TWICE, identically, at
- * `src/canonical/entries/mechanics.ts:35` and
- * `src/canonical/entries/fluids-waves.ts:48`. One definition here; the test
- * asserts agreement with both.
- *
- * @internal
- */
-export declare const SPRING_CONSTANT: Dimension;
+/** Farad, henry, ohm and newton per metre: the named shapes of `dimensional/types.ts`. @internal */
+export { CAPACITANCE, INDUCTANCE, RESISTANCE, SPRING_CONSTANT } from '../../dimensional/types.js';
 /** `M T^-1` — the `b` of `m x'' + b x' + k x = 0`. @internal */
 export declare const DAMPING: Dimension;
 /** `M L^-2 T^-2` — the `β` of `m x'' + k x + β x³ = 0`. @internal */

package/dist/atlas/oscillators/dimensions.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"dimensions.d.ts","sourceRoot":"","sources":["../../../src/atlas/oscillators/dimensions.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;GAYG;AAGH,OAAO,KAAK,EAAE,SAAS,EAAE,MAAM,4BAA4B,CAAC;AAE5D,6FAA6F;AAC7F,eAAO,MAAM,WAAW,EAAE,SAA6B,CAAC;AAExD,2FAA2F;AAC3F,eAAO,MAAM,UAAU,EAAE,SAA6B,CAAC;AAEvD,yFAAyF;AACzF,eAAO,MAAM,UAAU,EAAE,SAA6B,CAAC;AAEvD;;;;;;;GAOG;AACH,eAAO,MAAM,eAAe,EAAE,SAAyB,CAAC;AAExD,gEAAgE;AAChE,eAAO,MAAM,OAAO,EAAE,SAAyB,CAAC;AAEhD,qEAAqE;AACrE,eAAO,MAAM,eAAe,EAAE,SAA0B,CAAC"}
+{"version":3,"file":"dimensions.d.ts","sourceRoot":"","sources":["../../../src/atlas/oscillators/dimensions.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;GAaG;AAGH,OAAO,KAAK,EAAE,SAAS,EAAE,MAAM,4BAA4B,CAAC;AAE5D,oGAAoG;AACpG,OAAO,EAAE,WAAW,EAAE,UAAU,EAAE,UAAU,EAAE,eAAe,EAAE,MAAM,4BAA4B,CAAC;AAElG,gEAAgE;AAChE,eAAO,MAAM,OAAO,EAAE,SAAyB,CAAC;AAEhD,qEAAqE;AACrE,eAAO,MAAM,eAAe,EAAE,SAA0B,CAAC"}

package/dist/atlas/oscillators/dimensions.js

@@ -1,32 +1,20 @@
 /**
- * Dimension constants the oscillator pilot needs that the canonical entry
- * files declare MODULE-LOCAL and do not export.
+ * Dimension constants the oscillator pilot needs.
  *
- * They are redefined here with the same `dim(...)` arguments, quoted from
- * source in design note §10 Q2. `tests/atlas/models.test.ts` asserts the
- * values, so a future divergence in an entry file fails a test.
+ * Capacitance, inductance, resistance and the spring constant are the named
+ * shapes of `src/dimensional/types.ts`, re-exported here so the family's
+ * modules keep one import site; `tests/atlas/models.test.ts` asserts they are
+ * those objects, not copies. `DAMPING` and `CUBIC_STIFFNESS` are this family's
+ * own and are defined here.
  *
  * `FREQUENCY`, `MASS`, `LENGTH`, `VELOCITY` and `ACCELERATION` are importable
- * from `src/dimensional/types.ts` and are NOT redefined here.
+ * from `src/dimensional/types.ts` and are NOT re-exported here.
  *
  * @module atlas/oscillators/dimensions
  */
 import { dim } from '../../dimensional/ast-builders.js';
-/** farad, `L^-2 M^-1 T^4 I^2` — `src/canonical/entries/electromagnetism.ts:42`. @internal */
-export const CAPACITANCE = dim(-2, -1, 4, 2);
-/** henry, `L^2 M T^-2 I^-2` — `src/canonical/entries/electromagnetism.ts:43`. @internal */
-export const INDUCTANCE = dim(2, 1, -2, -2);
-/** ohm, `L^2 M T^-3 I^-2` — `src/canonical/entries/electromagnetism.ts:40`. @internal */
-export const RESISTANCE = dim(2, 1, -3, -2);
-/**
- * `M T^-2` — defined TWICE, identically, at
- * `src/canonical/entries/mechanics.ts:35` and
- * `src/canonical/entries/fluids-waves.ts:48`. One definition here; the test
- * asserts agreement with both.
- *
- * @internal
- */
-export const SPRING_CONSTANT = dim(0, 1, -2);
+/** Farad, henry, ohm and newton per metre: the named shapes of `dimensional/types.ts`. @internal */
+export { CAPACITANCE, INDUCTANCE, RESISTANCE, SPRING_CONSTANT } from '../../dimensional/types.js';
 /** `M T^-1` — the `b` of `m x'' + b x' + k x = 0`. @internal */
 export const DAMPING = dim(0, 1, -1);
 /** `M L^-2 T^-2` — the `β` of `m x'' + k x + β x³ = 0`. @internal */

package/dist/atlas/oscillators/dimensions.js.map

@@ -1 +1 @@
-{"version":3,"file":"dimensions.js","sourceRoot":"","sources":["../../../src/atlas/oscillators/dimensions.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;GAYG;AAEH,OAAO,EAAE,GAAG,EAAE,MAAM,mCAAmC,CAAC;AAGxD,6FAA6F;AAC7F,MAAM,CAAC,MAAM,WAAW,GAAc,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;AAExD,2FAA2F;AAC3F,MAAM,CAAC,MAAM,UAAU,GAAc,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;AAEvD,yFAAyF;AACzF,MAAM,CAAC,MAAM,UAAU,GAAc,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;AAEvD;;;;;;;GAOG;AACH,MAAM,CAAC,MAAM,eAAe,GAAc,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;AAExD,gEAAgE;AAChE,MAAM,CAAC,MAAM,OAAO,GAAc,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;AAEhD,qEAAqE;AACrE,MAAM,CAAC,MAAM,eAAe,GAAc,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC"}
+{"version":3,"file":"dimensions.js","sourceRoot":"","sources":["../../../src/atlas/oscillators/dimensions.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;GAaG;AAEH,OAAO,EAAE,GAAG,EAAE,MAAM,mCAAmC,CAAC;AAGxD,oGAAoG;AACpG,OAAO,EAAE,WAAW,EAAE,UAAU,EAAE,UAAU,EAAE,eAAe,EAAE,MAAM,4BAA4B,CAAC;AAElG,gEAAgE;AAChE,MAAM,CAAC,MAAM,OAAO,GAAc,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;AAEhD,qEAAqE;AACrE,MAAM,CAAC,MAAM,eAAe,GAAc,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC","sourcesContent":["/**\n * Dimension constants the oscillator pilot needs.\n *\n * Capacitance, inductance, resistance and the spring constant are the named\n * shapes of `src/dimensional/types.ts`, re-exported here so the family's\n * modules keep one import site; `tests/atlas/models.test.ts` asserts they are\n * those objects, not copies. `DAMPING` and `CUBIC_STIFFNESS` are this family's\n * own and are defined here.\n *\n * `FREQUENCY`, `MASS`, `LENGTH`, `VELOCITY` and `ACCELERATION` are importable\n * from `src/dimensional/types.ts` and are NOT re-exported here.\n *\n * @module atlas/oscillators/dimensions\n */\n\nimport { dim } from '../../dimensional/ast-builders.js';\nimport type { Dimension } from '../../dimensional/types.js';\n\n/** Farad, henry, ohm and newton per metre: the named shapes of `dimensional/types.ts`. @internal */\nexport { CAPACITANCE, INDUCTANCE, RESISTANCE, SPRING_CONSTANT } from '../../dimensional/types.js';\n\n/** `M T^-1` — the `b` of `m x'' + b x' + k x = 0`. @internal */\nexport const DAMPING: Dimension = dim(0, 1, -1);\n\n/** `M L^-2 T^-2` — the `β` of `m x'' + k x + β x³ = 0`. @internal */\nexport const CUBIC_STIFFNESS: Dimension = dim(-2, 1, -2);\n"]}

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

@@ -12,22 +12,7 @@
  *
  * @module atlas/oscillators
  */
-import type { AtlasBridge, AtlasRejection } from '../types.js';
-import type { AtlasModel } from '../model.js';
-/**
- * One atlas family: its models, the bridges between them, and the claimed
- * bridges that were refuted. A rejection is evidence, not an omission, so it
- * travels with the family rather than being dropped.
- *
- * @internal
- */
-export interface AtlasFamily {
-    /** `'oscillators'`. */
-    readonly family: string;
-    readonly models: readonly AtlasModel[];
-    readonly bridges: readonly AtlasBridge[];
-    readonly rejections: readonly AtlasRejection[];
-}
+import type { AtlasFamily } from '../family.js';
 /**
  * The oscillator pilot family, in design-note order: the two exact
  * equivalences, the two approximations, then the coarse-graining.

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

@@ -1 +1 @@
-{"version":3,"file":"index.d.ts","sourceRoot":"","sources":["../../../src/atlas/oscillators/index.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;GAaG;AAEH,OAAO,KAAK,EAAE,WAAW,EAAE,cAAc,EAAE,MAAM,aAAa,CAAC;AAC/D,OAAO,KAAK,EAAE,UAAU,EAAE,MAAM,aAAa,CAAC;AAO9C;;;;;;GAMG;AACH,MAAM,WAAW,WAAW;IAC1B,uBAAuB;IACvB,QAAQ,CAAC,MAAM,EAAE,MAAM,CAAC;IACxB,QAAQ,CAAC,MAAM,EAAE,SAAS,UAAU,EAAE,CAAC;IACvC,QAAQ,CAAC,OAAO,EAAE,SAAS,WAAW,EAAE,CAAC;IACzC,QAAQ,CAAC,UAAU,EAAE,SAAS,cAAc,EAAE,CAAC;CAChD;AAED;;;;;GAKG;AACH,eAAO,MAAM,iBAAiB,EAAE,WAU/B,CAAC"}
+{"version":3,"file":"index.d.ts","sourceRoot":"","sources":["../../../src/atlas/oscillators/index.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;GAaG;AAEH,OAAO,KAAK,EAAE,WAAW,EAAE,MAAM,cAAc,CAAC;AAOhD;;;;;GAKG;AACH,eAAO,MAAM,iBAAiB,EAAE,WAU/B,CAAC"}

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

@@ -1 +1 @@
-{"version":3,"file":"index.js","sourceRoot":"","sources":["../../../src/atlas/oscillators/index.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;GAaG;AAIH,OAAO,EAAE,YAAY,EAAE,MAAM,aAAa,CAAC;AAC3C,OAAO,EAAE,gBAAgB,EAAE,iBAAiB,EAAE,MAAM,oBAAoB,CAAC;AACzE,OAAO,EAAE,aAAa,EAAE,MAAM,qBAAqB,CAAC;AACpD,OAAO,EAAE,iBAAiB,EAAE,MAAM,qBAAqB,CAAC;AACxD,OAAO,EAAE,gBAAgB,EAAE,MAAM,iBAAiB,CAAC;AAiBnD;;;;;GAKG;AACH,MAAM,CAAC,MAAM,iBAAiB,GAAgB;IAC5C,MAAM,EAAE,aAAa;IACrB,MAAM,EAAE,YAAY;IACpB,OAAO,EAAE;QACP,gBAAgB;QAChB,iBAAiB;QACjB,GAAG,aAAa;QAChB,iBAAiB;KAClB;IACD,UAAU,EAAE,gBAAgB;CAC7B,CAAC"}
+{"version":3,"file":"index.js","sourceRoot":"","sources":["../../../src/atlas/oscillators/index.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;GAaG;AAGH,OAAO,EAAE,YAAY,EAAE,MAAM,aAAa,CAAC;AAC3C,OAAO,EAAE,gBAAgB,EAAE,iBAAiB,EAAE,MAAM,oBAAoB,CAAC;AACzE,OAAO,EAAE,aAAa,EAAE,MAAM,qBAAqB,CAAC;AACpD,OAAO,EAAE,iBAAiB,EAAE,MAAM,qBAAqB,CAAC;AACxD,OAAO,EAAE,gBAAgB,EAAE,MAAM,iBAAiB,CAAC;AAEnD;;;;;GAKG;AACH,MAAM,CAAC,MAAM,iBAAiB,GAAgB;IAC5C,MAAM,EAAE,aAAa;IACrB,MAAM,EAAE,YAAY;IACpB,OAAO,EAAE;QACP,gBAAgB;QAChB,iBAAiB;QACjB,GAAG,aAAa;QAChB,iBAAiB;KAClB;IACD,UAAU,EAAE,gBAAgB;CAC7B,CAAC","sourcesContent":["/**\n * Oscillator-family assembly (Atlas Phase 0, S0.6).\n *\n * Collects the W1/W2 records — the nine models, the five bridges and the one\n * rejection — into the single value the JSON projection serializes. Assembly\n * only: this module defines no physics and no numbers of its own, so a change\n * to a bridge is a change in ONE place.\n *\n * The design note calls this file `family.ts` (§1 layout); the S0.6 brief\n * places it at `oscillators/index.ts`. The brief is followed, because the\n * package's `exports` map and every other barrel in this tree use `index.ts`.\n *\n * @module atlas/oscillators\n */\n\nimport type { AtlasFamily } from '../family.js';\nimport { ATLAS_MODELS } from './models.js';\nimport { BRIDGE_SPRING_LC, BRIDGE_DAMPED_RLC } from './bridges-exact.js';\nimport { LIMIT_BRIDGES } from './bridges-limits.js';\nimport { BRIDGE_CHAIN_WAVE } from './bridges-coarse.js';\nimport { ATLAS_REJECTIONS } from './rejections.js';\n\n/**\n * The oscillator pilot family, in design-note order: the two exact\n * equivalences, the two approximations, then the coarse-graining.\n *\n * @internal\n */\nexport const OSCILLATOR_FAMILY: AtlasFamily = {\n  family: 'oscillators',\n  models: ATLAS_MODELS,\n  bridges: [\n    BRIDGE_SPRING_LC,\n    BRIDGE_DAMPED_RLC,\n    ...LIMIT_BRIDGES,\n    BRIDGE_CHAIN_WAVE,\n  ],\n  rejections: ATLAS_REJECTIONS,\n};\n"]}

package/dist/atlas/oscillators/limit-witnesses.d.ts

@@ -3,8 +3,8 @@
  * bridges that had none: W7 (`ab-pendulum-linear`), W8b (`ab-damped-massless`)
  * and W9 (`ab-chain-wave`).
  *
- * Each one INTEGRATES the premise model's own equation with MathTS
- * `solveODESystem` and `dt`, and reads the
+ * Each one INTEGRATES the premise model's own equation with `rk4Trajectory`
+ * (MathTS `solveODESystem` with `dt`, every state kept) and reads the
  * claimed quantity off the motion. None evaluates the closed form its record
  * states, so a test comparing the measurement with that closed form compares
  * two methods, not one method with itself.

package/dist/atlas/oscillators/limit-witnesses.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"limit-witnesses.d.ts","sourceRoot":"","sources":["../../../src/atlas/oscillators/limit-witnesses.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;GAaG;AAKH;;;;;GAKG;AACH,eAAO,MAAM,UAAU;aAAK,MAAM,EAAE,GAAG;aAAE,cAAc,EAAE,KAAM;CAAW,CAAC;AAE3E;;;;;;GAMG;AACH,wBAAgB,0BAA0B,CAAC,MAAM,EAAE,MAAM,EAAE,cAAc,GAAE,MAAkC,GAAG,MAAM,CAKrH;AAED;;;;;;;GAOG;AACH,eAAO,MAAM,WAAW;aAAK,CAAC,EAAE,IAAI;aAAE,CAAC,EAAE,CAAC;aAAE,CAAC,EAAE,CAAC;aAAE,EAAE,EAAE,CAAC;aAAE,EAAE,EAAE,CAAC;aAAE,IAAI,EAAE,EAAE;aAAE,YAAY,EAAE,EAAE;CAAW,CAAC;AAEtG;;;;;;;GAOG;AACH,wBAAgB,qBAAqB,CACnC,CAAC,EAAE,MAAM,EACT,OAAO,GAAE;IAAE,QAAQ,CAAC,CAAC,EAAE,MAAM,CAAC;IAAC,QAAQ,CAAC,CAAC,EAAE,MAAM,CAAC;IAAC,QAAQ,CAAC,EAAE,EAAE,MAAM,CAAC;IAAC,QAAQ,CAAC,EAAE,EAAE,MAAM,CAAC;IAAC,QAAQ,CAAC,IAAI,EAAE,MAAM,CAAC;IAAC,QAAQ,CAAC,YAAY,EAAE,MAAM,CAAA;CAAgB,GAChK,MAAM,CAmBR;AAED;;;;;;GAMG;AACH,eAAO,MAAM,UAAU;aAAK,MAAM,EAAE,EAAE;aAAE,KAAK,EAAE,CAAC;aAAE,CAAC,EAAE,CAAC;aAAE,CAAC,EAAE,CAAC;aAAE,gBAAgB,EAAE,IAAI;CAAW,CAAC;AAEhG;;;;;;;GAOG;AACH,wBAAgB,2BAA2B,CACzC,MAAM,EAAE,MAAM,EACd,OAAO,GAAE;IAAE,QAAQ,CAAC,KAAK,EAAE,MAAM,CAAC;IAAC,QAAQ,CAAC,CAAC,EAAE,MAAM,CAAC;IAAC,QAAQ,CAAC,CAAC,EAAE,MAAM,CAAC;IAAC,QAAQ,CAAC,gBAAgB,EAAE,MAAM,CAAA;CAAe,GAC1H,MAAM,CAmDR"}
+{"version":3,"file":"limit-witnesses.d.ts","sourceRoot":"","sources":["../../../src/atlas/oscillators/limit-witnesses.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;GAaG;AAKH;;;;;GAKG;AACH,eAAO,MAAM,UAAU;aAAK,MAAM,EAAE,GAAG;aAAE,cAAc,EAAE,KAAM;CAAW,CAAC;AAE3E;;;;;;GAMG;AACH,wBAAgB,0BAA0B,CAAC,MAAM,EAAE,MAAM,EAAE,cAAc,GAAE,MAAkC,GAAG,MAAM,CAKrH;AAED;;;;;;;GAOG;AACH,eAAO,MAAM,WAAW;aAAK,CAAC,EAAE,IAAI;aAAE,CAAC,EAAE,CAAC;aAAE,CAAC,EAAE,CAAC;aAAE,EAAE,EAAE,CAAC;aAAE,EAAE,EAAE,CAAC;aAAE,IAAI,EAAE,EAAE;aAAE,YAAY,EAAE,EAAE;CAAW,CAAC;AAEtG;;;;;;;GAOG;AACH,wBAAgB,qBAAqB,CACnC,CAAC,EAAE,MAAM,EACT,OAAO,GAAE;IAAE,QAAQ,CAAC,CAAC,EAAE,MAAM,CAAC;IAAC,QAAQ,CAAC,CAAC,EAAE,MAAM,CAAC;IAAC,QAAQ,CAAC,EAAE,EAAE,MAAM,CAAC;IAAC,QAAQ,CAAC,EAAE,EAAE,MAAM,CAAC;IAAC,QAAQ,CAAC,IAAI,EAAE,MAAM,CAAC;IAAC,QAAQ,CAAC,YAAY,EAAE,MAAM,CAAA;CAAgB,GAChK,MAAM,CAcR;AAED;;;;;;GAMG;AACH,eAAO,MAAM,UAAU;aAAK,MAAM,EAAE,EAAE;aAAE,KAAK,EAAE,CAAC;aAAE,CAAC,EAAE,CAAC;aAAE,CAAC,EAAE,CAAC;aAAE,gBAAgB,EAAE,IAAI;CAAW,CAAC;AAEhG;;;;;;;GAOG;AACH,wBAAgB,2BAA2B,CACzC,MAAM,EAAE,MAAM,EACd,OAAO,GAAE;IAAE,QAAQ,CAAC,KAAK,EAAE,MAAM,CAAC;IAAC,QAAQ,CAAC,CAAC,EAAE,MAAM,CAAC;IAAC,QAAQ,CAAC,CAAC,EAAE,MAAM,CAAC;IAAC,QAAQ,CAAC,gBAAgB,EAAE,MAAM,CAAA;CAAe,GAC1H,MAAM,CAsDR"}

package/dist/atlas/oscillators/limit-witnesses.js

@@ -3,8 +3,8 @@
  * bridges that had none: W7 (`ab-pendulum-linear`), W8b (`ab-damped-massless`)
  * and W9 (`ab-chain-wave`).
  *
- * Each one INTEGRATES the premise model's own equation with MathTS
- * `solveODESystem` and `dt`, and reads the
+ * Each one INTEGRATES the premise model's own equation with `rk4Trajectory`
+ * (MathTS `solveODESystem` with `dt`, every state kept) and reads the
  * claimed quantity off the motion. None evaluates the closed form its record
  * states, so a test comparing the measurement with that closed form compares
  * two methods, not one method with itself.
@@ -12,7 +12,7 @@
  * @module atlas/oscillators/limit-witnesses
  * @internal
  */
-import { solveODESystem } from '@danielsimonjr/mathts-functions';
+import { rk4Trajectory } from '../../numerical/rk4-trajectory.js';
 import { quarterPeriod } from './norm-transport-witness.js';
 /**
  * W7: the pendulum θ'' = −ω0² sin θ at ω0 = 1, released from rest; the
@@ -59,10 +59,10 @@ export function measureMasslessOffset(m, fixture = W8B_FIXTURE) {
     const h = m / stepsPerMass;
     const layerSteps = Math.ceil((5 * stepsPerMass) / b);
     const steps = Math.ceil(tEnd / h);
-    const sol = solveODESystem((_t, y) => [y[1], -(b * y[1] + k * y[0]) / m], [x0, v0], [0, steps * h], { dt: h });
+    const states = rk4Trajectory((_t, y) => [y[1], -(b * y[1] + k * y[0]) / m], [x0, v0], 0, steps * h, steps);
     let sup = 0;
     for (let i = layerSteps; i <= steps; i++) {
-        const x = sol.y[i]?.[0];
+        const x = states[i]?.[0];
         if (x === undefined)
             return Number.NaN;
         sup = Math.max(sup, Math.abs(x - x0 * Math.exp((-k / b) * i * h)));
@@ -98,21 +98,24 @@ export function measureChainDispersionError(masses, fixture = W9_FIXTURE) {
         return out;
     };
     const h = 1 / stepsPerUnitTime;
-    let state = [
+    const start = [
         ...Float64Array.from({ length: N }, (_, n) => Math.cos((2 * Math.PI * n) / N)),
         ...new Array(N).fill(0),
     ];
     // Twice the lattice quarter period π/(4√(κ/m) sin(π/N)): the crossing is inside.
     const maxSteps = Math.ceil(((2 * Math.PI) / (4 * Math.sqrt(w) * Math.sin(Math.PI / N))) * stepsPerUnitTime);
+    // One continuous run; the system ignores t, so the states are those of a
+    // step-by-step loop, and the first crossing of u_0 is read off them.
+    const states = rk4Trajectory((_t, y) => {
+        const pos = Float64Array.from(y.slice(0, N));
+        const vel = y.slice(N);
+        const a = accel(pos);
+        return [...vel, ...a];
+    }, start, 0, maxSteps * h, maxSteps);
     for (let i = 0; i < maxSteps; i++) {
-        const sol = solveODESystem((_t, y) => {
-            const pos = Float64Array.from(y.slice(0, N));
-            const vel = y.slice(N);
-            const a = accel(pos);
-            return [...vel, ...a];
-        }, state, [0, h], { dt: h });
-        const next = sol.y[sol.y.length - 1];
-        if (next === undefined)
+        const state = states[i];
+        const next = states[i + 1];
+        if (state === undefined || next === undefined)
             return Number.NaN;
         const u0 = state[0];
         const nu0 = next[0];
@@ -130,7 +133,6 @@ export function measureChainDispersionError(masses, fixture = W9_FIXTURE) {
             const q = (2 * Math.PI) / (N * a);
             return 1 - omega / (c * q);
         }
-        state = next.slice();
     }
     return Number.NaN;
 }

package/dist/atlas/oscillators/limit-witnesses.js.map

@@ -1 +1 @@
-{"version":3,"file":"limit-witnesses.js","sourceRoot":"","sources":["../../../src/atlas/oscillators/limit-witnesses.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;GAaG;AAEH,OAAO,EAAE,cAAc,EAAE,MAAM,iCAAiC,CAAC;AACjE,OAAO,EAAE,aAAa,EAAE,MAAM,6BAA6B,CAAC;AAE5D;;;;;GAKG;AACH,MAAM,CAAC,MAAM,UAAU,GAAG,EAAE,MAAM,EAAE,GAAG,EAAE,cAAc,EAAE,MAAM,EAAW,CAAC;AAE3E;;;;;;GAMG;AACH,MAAM,UAAU,0BAA0B,CAAC,MAAc,EAAE,cAAc,GAAW,UAAU,CAAC,cAAc;IAC3G,IAAI,CAAC,CAAC,MAAM,GAAG,CAAC,IAAI,MAAM,GAAG,IAAI,CAAC,EAAE,CAAC;QAAE,OAAO,MAAM,CAAC,GAAG,CAAC;IACzD,MAAM,EAAE,GAAG,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC;IACvB,MAAM,OAAO,GAAG,aAAa,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,MAAM,EAAE,EAAE,GAAG,cAAc,EAAE,CAAC,GAAG,cAAc,CAAC,CAAC;IACpG,OAAO,CAAC,CAAC,GAAG,OAAO,CAAC,GAAG,EAAE,CAAC;AAC5B,CAAC;AAED;;;;;;;GAOG;AACH,MAAM,CAAC,MAAM,WAAW,GAAG,EAAE,CAAC,EAAE,IAAI,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,EAAE,EAAE,CAAC,EAAE,EAAE,EAAE,CAAC,EAAE,IAAI,EAAE,EAAE,EAAE,YAAY,EAAE,EAAE,EAAW,CAAC;AAEtG;;;;;;;GAOG;AACH,MAAM,UAAU,qBAAqB,CACnC,CAAS,EACT,OAAO,GAA+I,WAAW;IAEjK,MAAM,EAAE,CAAC,EAAE,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,IAAI,EAAE,YAAY,EAAE,GAAG,OAAO,CAAC;IACrD,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC;QAAE,OAAO,MAAM,CAAC,GAAG,CAAC;IAC5E,MAAM,CAAC,GAAG,CAAC,GAAG,YAAY,CAAC;IAC3B,MAAM,UAAU,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,YAAY,CAAC,GAAG,CAAC,CAAC,CAAC;IACrD,MAAM,KAAK,GAAG,IAAI,CAAC,IAAI,CAAC,IAAI,GAAG,CAAC,CAAC,CAAC;IAClC,MAAM,GAAG,GAAG,cAAc,CACxB,CAAC,EAAE,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAE,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAE,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAE,CAAC,GAAG,CAAC,CAAC,EAChD,CAAC,EAAE,EAAE,EAAE,CAAC,EACR,CAAC,CAAC,EAAE,KAAK,GAAG,CAAC,CAAC,EACd,EAAE,EAAE,EAAE,CAAC,EAAE,CACV,CAAC;IACF,IAAI,GAAG,GAAG,CAAC,CAAC;IACZ,KAAK,IAAI,CAAC,GAAG,UAAU,EAAE,CAAC,IAAI,KAAK,EAAE,CAAC,EAAE,EAAE,CAAC;QACzC,MAAM,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;QACxB,IAAI,CAAC,KAAK,SAAS;YAAE,OAAO,MAAM,CAAC,GAAG,CAAC;QACvC,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,EAAE,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;IACrE,CAAC;IACD,OAAO,GAAG,CAAC;AACb,CAAC;AAED;;;;;;GAMG;AACH,MAAM,CAAC,MAAM,UAAU,GAAG,EAAE,MAAM,EAAE,EAAE,EAAE,KAAK,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,gBAAgB,EAAE,IAAI,EAAW,CAAC;AAEhG;;;;;;;GAOG;AACH,MAAM,UAAU,2BAA2B,CACzC,MAAc,EACd,OAAO,GAA0G,UAAU;IAE3H,MAAM,CAAC,GAAG,MAAM,CAAC;IACjB,IAAI,CAAC,MAAM,CAAC,SAAS,CAAC,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC;QAAE,OAAO,MAAM,CAAC,GAAG,CAAC;IACrD,MAAM,EAAE,KAAK,EAAE,CAAC,EAAE,CAAC,EAAE,gBAAgB,EAAE,GAAG,OAAO,CAAC;IAClD,MAAM,CAAC,GAAG,KAAK,GAAG,CAAC,CAAC;IACpB,MAAM,KAAK,GAAG,CAAC,CAAe,EAAgB,EAAE;QAC9C,MAAM,GAAG,GAAG,IAAI,YAAY,CAAC,CAAC,CAAC,CAAC;QAChC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE;YAAE,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAE,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAE,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAE,CAAC,CAAC;QAC7F,OAAO,GAAG,CAAC;IACb,CAAC,CAAC;IACF,MAAM,CAAC,GAAG,CAAC,GAAG,gBAAgB,CAAC;IAC/B,IAAI,KAAK,GAAG;QACV,GAAG,YAAY,CAAC,IAAI,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;QAC9E,GAAG,IAAI,KAAK,CAAS,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC;KAChC,CAAC;IACF,iFAAiF;IACjF,MAAM,QAAQ,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,gBAAgB,CAAC,CAAC;IAC5G,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,QAAQ,EAAE,CAAC,EAAE,EAAE,CAAC;QAClC,MAAM,GAAG,GAAG,cAAc,CACxB,CAAC,EAAE,EAAE,CAAC,EAAE,EAAE;YACR,MAAM,GAAG,GAAG,YAAY,CAAC,IAAI,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;YAC7C,MAAM,GAAG,GAAG,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC;YACvB,MAAM,CAAC,GAAG,KAAK,CAAC,GAAG,CAAC,CAAC;YACrB,OAAO,CAAC,GAAG,GAAG,EAAE,GAAG,CAAC,CAAC,CAAC;QACxB,CAAC,EACD,KAAK,EACL,CAAC,CAAC,EAAE,CAAC,CAAC,EACN,EAAE,EAAE,EAAE,CAAC,EAAE,CACV,CAAC;QACF,MAAM,IAAI,GAAG,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC;QACrC,IAAI,IAAI,KAAK,SAAS;YAAE,OAAO,MAAM,CAAC,GAAG,CAAC;QAC1C,MAAM,EAAE,GAAG,KAAK,CAAC,CAAC,CAAE,CAAC;QACrB,MAAM,GAAG,GAAG,IAAI,CAAC,CAAC,CAAE,CAAC;QACrB,MAAM,EAAE,GAAG,KAAK,CAAC,CAAC,CAAE,CAAC;QACrB,MAAM,GAAG,GAAG,IAAI,CAAC,CAAC,CAAE,CAAC;QACrB,IAAI,EAAE,GAAG,CAAC,IAAI,GAAG,IAAI,CAAC,EAAE,CAAC;YACvB,MAAM,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,GAAG,CAAC,EAAE,EAAE,EAAE,GAAG,CAAC,EAAE,GAAG,EAAE,GAAG,GAAG,CAAC,CAAC,CAAC;YACpD,MAAM,CAAC,GAAG,CAAC,CAAS,EAAU,EAAE,CAC9B,CAAC,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC,GAAG,EAAE,CAAC;YACjI,MAAM,EAAE,GAAG,CAAC,CAAS,EAAU,EAAE,CAC/B,CAAC,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC;YACrH,IAAI,CAAC,GAAG,EAAE,GAAG,CAAC,EAAE,GAAG,EAAE,CAAC,CAAC;YACvB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE;gBAAE,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC;YAC9C,MAAM,KAAK,GAAG,IAAI,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;YAC1C,MAAM,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;YAC3B,MAAM,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;YAClC,OAAO,CAAC,GAAG,KAAK,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;QAC7B,CAAC;QACD,KAAK,GAAG,IAAI,CAAC,KAAK,EAAE,CAAC;IACvB,CAAC;IACD,OAAO,MAAM,CAAC,GAAG,CAAC;AACpB,CAAC"}
+{"version":3,"file":"limit-witnesses.js","sourceRoot":"","sources":["../../../src/atlas/oscillators/limit-witnesses.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;GAaG;AAEH,OAAO,EAAE,aAAa,EAAE,MAAM,mCAAmC,CAAC;AAClE,OAAO,EAAE,aAAa,EAAE,MAAM,6BAA6B,CAAC;AAE5D;;;;;GAKG;AACH,MAAM,CAAC,MAAM,UAAU,GAAG,EAAE,MAAM,EAAE,GAAG,EAAE,cAAc,EAAE,MAAM,EAAW,CAAC;AAE3E;;;;;;GAMG;AACH,MAAM,UAAU,0BAA0B,CAAC,MAAc,EAAE,cAAc,GAAW,UAAU,CAAC,cAAc;IAC3G,IAAI,CAAC,CAAC,MAAM,GAAG,CAAC,IAAI,MAAM,GAAG,IAAI,CAAC,EAAE,CAAC;QAAE,OAAO,MAAM,CAAC,GAAG,CAAC;IACzD,MAAM,EAAE,GAAG,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC;IACvB,MAAM,OAAO,GAAG,aAAa,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,MAAM,EAAE,EAAE,GAAG,cAAc,EAAE,CAAC,GAAG,cAAc,CAAC,CAAC;IACpG,OAAO,CAAC,CAAC,GAAG,OAAO,CAAC,GAAG,EAAE,CAAC;AAC5B,CAAC;AAED;;;;;;;GAOG;AACH,MAAM,CAAC,MAAM,WAAW,GAAG,EAAE,CAAC,EAAE,IAAI,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,EAAE,EAAE,CAAC,EAAE,EAAE,EAAE,CAAC,EAAE,IAAI,EAAE,EAAE,EAAE,YAAY,EAAE,EAAE,EAAW,CAAC;AAEtG;;;;;;;GAOG;AACH,MAAM,UAAU,qBAAqB,CACnC,CAAS,EACT,OAAO,GAA+I,WAAW;IAEjK,MAAM,EAAE,CAAC,EAAE,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,IAAI,EAAE,YAAY,EAAE,GAAG,OAAO,CAAC;IACrD,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC;QAAE,OAAO,MAAM,CAAC,GAAG,CAAC;IAC5E,MAAM,CAAC,GAAG,CAAC,GAAG,YAAY,CAAC;IAC3B,MAAM,UAAU,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,YAAY,CAAC,GAAG,CAAC,CAAC,CAAC;IACrD,MAAM,KAAK,GAAG,IAAI,CAAC,IAAI,CAAC,IAAI,GAAG,CAAC,CAAC,CAAC;IAClC,MAAM,MAAM,GAAG,aAAa,CAAC,CAAC,EAAE,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAE,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAE,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAE,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,EAAE,CAAC,EAAE,CAAC,EAAE,KAAK,GAAG,CAAC,EAAE,KAAK,CAAC,CAAC;IAC9G,IAAI,GAAG,GAAG,CAAC,CAAC;IACZ,KAAK,IAAI,CAAC,GAAG,UAAU,EAAE,CAAC,IAAI,KAAK,EAAE,CAAC,EAAE,EAAE,CAAC;QACzC,MAAM,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;QACzB,IAAI,CAAC,KAAK,SAAS;YAAE,OAAO,MAAM,CAAC,GAAG,CAAC;QACvC,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,EAAE,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;IACrE,CAAC;IACD,OAAO,GAAG,CAAC;AACb,CAAC;AAED;;;;;;GAMG;AACH,MAAM,CAAC,MAAM,UAAU,GAAG,EAAE,MAAM,EAAE,EAAE,EAAE,KAAK,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,gBAAgB,EAAE,IAAI,EAAW,CAAC;AAEhG;;;;;;;GAOG;AACH,MAAM,UAAU,2BAA2B,CACzC,MAAc,EACd,OAAO,GAA0G,UAAU;IAE3H,MAAM,CAAC,GAAG,MAAM,CAAC;IACjB,IAAI,CAAC,MAAM,CAAC,SAAS,CAAC,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC;QAAE,OAAO,MAAM,CAAC,GAAG,CAAC;IACrD,MAAM,EAAE,KAAK,EAAE,CAAC,EAAE,CAAC,EAAE,gBAAgB,EAAE,GAAG,OAAO,CAAC;IAClD,MAAM,CAAC,GAAG,KAAK,GAAG,CAAC,CAAC;IACpB,MAAM,KAAK,GAAG,CAAC,CAAe,EAAgB,EAAE;QAC9C,MAAM,GAAG,GAAG,IAAI,YAAY,CAAC,CAAC,CAAC,CAAC;QAChC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE;YAAE,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAE,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAE,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAE,CAAC,CAAC;QAC7F,OAAO,GAAG,CAAC;IACb,CAAC,CAAC;IACF,MAAM,CAAC,GAAG,CAAC,GAAG,gBAAgB,CAAC;IAC/B,MAAM,KAAK,GAAG;QACZ,GAAG,YAAY,CAAC,IAAI,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;QAC9E,GAAG,IAAI,KAAK,CAAS,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC;KAChC,CAAC;IACF,iFAAiF;IACjF,MAAM,QAAQ,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,gBAAgB,CAAC,CAAC;IAC5G,yEAAyE;IACzE,qEAAqE;IACrE,MAAM,MAAM,GAAG,aAAa,CAC1B,CAAC,EAAE,EAAE,CAAC,EAAE,EAAE;QACR,MAAM,GAAG,GAAG,YAAY,CAAC,IAAI,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;QAC7C,MAAM,GAAG,GAAG,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC;QACvB,MAAM,CAAC,GAAG,KAAK,CAAC,GAAG,CAAC,CAAC;QACrB,OAAO,CAAC,GAAG,GAAG,EAAE,GAAG,CAAC,CAAC,CAAC;IACxB,CAAC,EACD,KAAK,EACL,CAAC,EACD,QAAQ,GAAG,CAAC,EACZ,QAAQ,CACT,CAAC;IACF,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,QAAQ,EAAE,CAAC,EAAE,EAAE,CAAC;QAClC,MAAM,KAAK,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC;QACxB,MAAM,IAAI,GAAG,MAAM,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;QAC3B,IAAI,KAAK,KAAK,SAAS,IAAI,IAAI,KAAK,SAAS;YAAE,OAAO,MAAM,CAAC,GAAG,CAAC;QACjE,MAAM,EAAE,GAAG,KAAK,CAAC,CAAC,CAAE,CAAC;QACrB,MAAM,GAAG,GAAG,IAAI,CAAC,CAAC,CAAE,CAAC;QACrB,MAAM,EAAE,GAAG,KAAK,CAAC,CAAC,CAAE,CAAC;QACrB,MAAM,GAAG,GAAG,IAAI,CAAC,CAAC,CAAE,CAAC;QACrB,IAAI,EAAE,GAAG,CAAC,IAAI,GAAG,IAAI,CAAC,EAAE,CAAC;YACvB,MAAM,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,GAAG,CAAC,EAAE,EAAE,EAAE,GAAG,CAAC,EAAE,GAAG,EAAE,GAAG,GAAG,CAAC,CAAC,CAAC;YACpD,MAAM,CAAC,GAAG,CAAC,CAAS,EAAU,EAAE,CAC9B,CAAC,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC,GAAG,EAAE,CAAC;YACjI,MAAM,EAAE,GAAG,CAAC,CAAS,EAAU,EAAE,CAC/B,CAAC,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC;YACrH,IAAI,CAAC,GAAG,EAAE,GAAG,CAAC,EAAE,GAAG,EAAE,CAAC,CAAC;YACvB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE;gBAAE,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC;YAC9C,MAAM,KAAK,GAAG,IAAI,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;YAC1C,MAAM,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;YAC3B,MAAM,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;YAClC,OAAO,CAAC,GAAG,KAAK,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;QAC7B,CAAC;IACH,CAAC;IACD,OAAO,MAAM,CAAC,GAAG,CAAC;AACpB,CAAC","sourcesContent":["/**\n * The measurements behind the in-process witnesses of the three oscillator\n * bridges that had none: W7 (`ab-pendulum-linear`), W8b (`ab-damped-massless`)\n * and W9 (`ab-chain-wave`).\n *\n * Each one INTEGRATES the premise model's own equation with `rk4Trajectory`\n * (MathTS `solveODESystem` with `dt`, every state kept) and reads the\n * claimed quantity off the motion. None evaluates the closed form its record\n * states, so a test comparing the measurement with that closed form compares\n * two methods, not one method with itself.\n *\n * @module atlas/oscillators/limit-witnesses\n * @internal\n */\n\nimport { rk4Trajectory } from '../../numerical/rk4-trajectory.js';\nimport { quarterPeriod } from './norm-transport-witness.js';\n\n/**\n * W7: the pendulum θ'' = −ω0² sin θ at ω0 = 1, released from rest; the\n * amplitude is `theta0 / resolution`, so refinement moves into the θ0 → 0 limit.\n *\n * @internal\n */\nexport const W7_FIXTURE = { theta0: 0.4, stepsPerPeriod: 16_384 } as const;\n\n/**\n * T/T0 of the pendulum released from rest at `theta0`, measured by RK4 with\n * `stepsPerPeriod` steps per linear period T0 = 2π (ω0 = 1).\n *\n * @returns NaN unless 0 < theta0 < π.\n * @internal\n */\nexport function measurePendulumPeriodRatio(theta0: number, stepsPerPeriod: number = W7_FIXTURE.stepsPerPeriod): number {\n  if (!(theta0 > 0 && theta0 < Math.PI)) return Number.NaN;\n  const T0 = 2 * Math.PI;\n  const quarter = quarterPeriod((x) => -Math.sin(x), theta0, T0 / stepsPerPeriod, 2 * stepsPerPeriod);\n  return (4 * quarter) / T0;\n}\n\n/**\n * W8b: the damped spring m x'' + b x' + k x = 0 at the record's witness\n * normalisation b = k = x0 = 1, with v0 = 5 (the edge of |v0| ≤ 5). The mass is\n * `m / resolution`. The RK4 step is m/20, so the fast rate b/m is resolved at\n * every mass, and 5m/b is a whole number of steps.\n *\n * @internal\n */\nexport const W8B_FIXTURE = { m: 0.02, b: 1, k: 1, x0: 1, v0: 5, tEnd: 10, stepsPerMass: 20 } as const;\n\n/**\n * sup |x(t) − x0 e^(−kt/b)| over the RK4 grid on t ∈ [5m/b, tEnd], x from the\n * full second-order model: the norm of `ab-damped-massless`'s bound. The\n * window stops at `tEnd`; both motions there are below x0·e^(−k·tEnd/b).\n *\n * @returns NaN unless m, b and k are positive and finite.\n * @internal\n */\nexport function measureMasslessOffset(\n  m: number,\n  fixture: { readonly b: number; readonly k: number; readonly x0: number; readonly v0: number; readonly tEnd: number; readonly stepsPerMass: number } = W8B_FIXTURE,\n): number {\n  const { b, k, x0, v0, tEnd, stepsPerMass } = fixture;\n  if (![m, b, k].every((v) => Number.isFinite(v) && v > 0)) return Number.NaN;\n  const h = m / stepsPerMass;\n  const layerSteps = Math.ceil((5 * stepsPerMass) / b);\n  const steps = Math.ceil(tEnd / h);\n  const states = rk4Trajectory((_t, y) => [y[1]!, -(b * y[1]! + k * y[0]!) / m], [x0, v0], 0, steps * h, steps);\n  let sup = 0;\n  for (let i = layerSteps; i <= steps; i++) {\n    const x = states[i]?.[0];\n    if (x === undefined) return Number.NaN;\n    sup = Math.max(sup, Math.abs(x - x0 * Math.exp((-k / b) * i * h)));\n  }\n  return sup;\n}\n\n/**\n * W9: a ring of `masses · resolution` masses, m u_n'' = κ(u_{n+1} − 2u_n + u_{n−1}),\n * started at rest in its longest mode u_n = cos(2π n/N). Refinement doubles N, so\n * qa = 2π/N halves.\n *\n * @internal\n */\nexport const W9_FIXTURE = { masses: 16, kappa: 1, m: 1, a: 1, stepsPerUnitTime: 1000 } as const;\n\n/**\n * 1 − ω/(c q) on a ring of `masses` masses in its longest mode, with ω measured\n * from the first zero of u_0(t) = cos ωt (RK4, cubic Hermite and Newton on the\n * bracketing step) and c = a √(κ/m), q = 2π/(N a).\n *\n * @returns NaN for fewer than 3 masses or no crossing.\n * @internal\n */\nexport function measureChainDispersionError(\n  masses: number,\n  fixture: { readonly kappa: number; readonly m: number; readonly a: number; readonly stepsPerUnitTime: number } = W9_FIXTURE,\n): number {\n  const N = masses;\n  if (!Number.isInteger(N) || N < 3) return Number.NaN;\n  const { kappa, m, a, stepsPerUnitTime } = fixture;\n  const w = kappa / m;\n  const accel = (u: Float64Array): Float64Array => {\n    const out = new Float64Array(N);\n    for (let n = 0; n < N; n++) out[n] = w * (u[(n + 1) % N]! - 2 * u[n]! + u[(n - 1 + N) % N]!);\n    return out;\n  };\n  const h = 1 / stepsPerUnitTime;\n  const start = [\n    ...Float64Array.from({ length: N }, (_, n) => Math.cos((2 * Math.PI * n) / N)),\n    ...new Array<number>(N).fill(0),\n  ];\n  // Twice the lattice quarter period π/(4√(κ/m) sin(π/N)): the crossing is inside.\n  const maxSteps = Math.ceil(((2 * Math.PI) / (4 * Math.sqrt(w) * Math.sin(Math.PI / N))) * stepsPerUnitTime);\n  // One continuous run; the system ignores t, so the states are those of a\n  // step-by-step loop, and the first crossing of u_0 is read off them.\n  const states = rk4Trajectory(\n    (_t, y) => {\n      const pos = Float64Array.from(y.slice(0, N));\n      const vel = y.slice(N);\n      const a = accel(pos);\n      return [...vel, ...a];\n    },\n    start,\n    0,\n    maxSteps * h,\n    maxSteps,\n  );\n  for (let i = 0; i < maxSteps; i++) {\n    const state = states[i];\n    const next = states[i + 1];\n    if (state === undefined || next === undefined) return Number.NaN;\n    const u0 = state[0]!;\n    const nu0 = next[0]!;\n    const v0 = state[N]!;\n    const nv0 = next[N]!;\n    if (u0 > 0 && nu0 <= 0) {\n      const [p0, m0, p1, m1] = [u0, v0 * h, nu0, nv0 * h];\n      const p = (s: number): number =>\n        (2 * s ** 3 - 3 * s ** 2 + 1) * p0 + (s ** 3 - 2 * s ** 2 + s) * m0 + (-2 * s ** 3 + 3 * s ** 2) * p1 + (s ** 3 - s ** 2) * m1;\n      const dp = (s: number): number =>\n        (6 * s ** 2 - 6 * s) * p0 + (3 * s ** 2 - 4 * s + 1) * m0 + (-6 * s ** 2 + 6 * s) * p1 + (3 * s ** 2 - 2 * s) * m1;\n      let s = p0 / (p0 - p1);\n      for (let j = 0; j < 8; j++) s -= p(s) / dp(s);\n      const omega = Math.PI / (2 * (i + s) * h);\n      const c = a * Math.sqrt(w);\n      const q = (2 * Math.PI) / (N * a);\n      return 1 - omega / (c * q);\n    }\n  }\n  return Number.NaN;\n}\n"]}

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

@@ -1 +1 @@
-{"version":3,"file":"models.d.ts","sourceRoot":"","sources":["../../../src/atlas/oscillators/models.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;GAmBG;AASH,OAAO,KAAK,EAAE,UAAU,EAAE,MAAM,aAAa,CAAC;AA6B9C,6EAA6E;AAC7E,eAAO,MAAM,YAAY,EAAE,SAAS,UAAU,EAmH7C,CAAC;AAIF,mEAAmE;AACnE,wBAAgB,aAAa,CAAC,EAAE,EAAE,MAAM,GAAG,UAAU,CAIpD"}
+{"version":3,"file":"models.d.ts","sourceRoot":"","sources":["../../../src/atlas/oscillators/models.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;GAmBG;AASH,OAAO,KAAK,EAAE,UAAU,EAAE,MAAM,aAAa,CAAC;AAyB9C,6EAA6E;AAC7E,eAAO,MAAM,YAAY,EAAE,SAAS,UAAU,EAmH7C,CAAC;AAIF,mEAAmE;AACnE,wBAAgB,aAAa,CAAC,EAAE,EAAE,MAAM,GAAG,UAAU,CAIpD"}

package/dist/atlas/oscillators/models.js

@@ -19,18 +19,14 @@
  * @module atlas/oscillators/models
  */
 import { ACCELERATION, LENGTH, MASS, VELOCITY, } from '../../dimensional/types.js';
-import { deriveRegimeGroups } from '../regime.js';
+import { makeRegime } from '../regime.js';
 import { CAPACITANCE, CUBIC_STIFFNESS, DAMPING, INDUCTANCE, RESISTANCE, SPRING_CONSTANT, } from './dimensions.js';
 const FAMILY = 'oscillators';
 function model(spec) {
     return {
         ...spec,
         family: FAMILY,
-        regime: {
-            family: FAMILY,
-            inequalities: [],
-            groupDefinitions: deriveRegimeGroups(FAMILY, spec.parameters, spec.dimensionlessInputs),
-        },
+        regime: makeRegime(FAMILY, spec.parameters, spec.dimensionlessInputs, []),
     };
 }
 /** The nine Phase 0 oscillator models, in design-note §3 order. @internal */

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

@@ -1 +1 @@
-{"version":3,"file":"models.js","sourceRoot":"","sources":["../../../src/atlas/oscillators/models.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;GAmBG;AAEH,OAAO,EACL,YAAY,EACZ,MAAM,EACN,IAAI,EACJ,QAAQ,GACT,MAAM,4BAA4B,CAAC;AACpC,OAAO,EAAE,kBAAkB,EAAE,MAAM,cAAc,CAAC;AAGlD,OAAO,EACL,WAAW,EACX,eAAe,EACf,OAAO,EACP,UAAU,EACV,UAAU,EACV,eAAe,GAChB,MAAM,iBAAiB,CAAC;AAEzB,MAAM,MAAM,GAAG,aAAa,CAAC;AAE7B,SAAS,KAAK,CACZ,IAEC;IAED,OAAO;QACL,GAAG,IAAI;QACP,MAAM,EAAE,MAAM;QACd,MAAM,EAAE;YACN,MAAM,EAAE,MAAM;YACd,YAAY,EAAE,EAAE;YAChB,gBAAgB,EAAE,kBAAkB,CAAC,MAAM,EAAE,IAAI,CAAC,UAAU,EAAE,IAAI,CAAC,mBAAmB,CAAC;SACxF;KACF,CAAC;AACJ,CAAC;AAED,6EAA6E;AAC7E,MAAM,CAAC,MAAM,YAAY,GAA0B;IACjD,KAAK,CAAC;QACJ,EAAE,EAAE,cAAc;QAClB,UAAU,EAAE,cAAc;QAC1B,QAAQ,EAAE,gBAAgB;QAC1B,WAAW,EAAE,CAAC,GAAG,EAAE,IAAI,EAAE,QAAQ,CAAC;QAClC,UAAU,EAAE;YACV,EAAE,IAAI,EAAE,GAAG,EAAE,GAAG,EAAE,IAAI,EAAE;YACxB,EAAE,IAAI,EAAE,GAAG,EAAE,GAAG,EAAE,eAAe,EAAE;SACpC;QACD,mBAAmB,EAAE,EAAE;QACvB,aAAa,EAAE;YACb,8BAA8B;YAC9B,4BAA4B;YAC5B,sBAAsB;SACvB;KACF,CAAC;IACF,KAAK,CAAC;QACJ,EAAE,EAAE,UAAU;QACd,UAAU,EAAE,cAAc;QAC1B,QAAQ,EAAE,gBAAgB;QAC1B,WAAW,EAAE,CAAC,GAAG,EAAE,IAAI,EAAE,QAAQ,CAAC;QAClC,UAAU,EAAE;YACV,EAAE,IAAI,EAAE,GAAG,EAAE,GAAG,EAAE,UAAU,EAAE;YAC9B,EAAE,IAAI,EAAE,GAAG,EAAE,GAAG,EAAE,WAAW,EAAE;SAChC;QACD,mBAAmB,EAAE,EAAE;QACvB,aAAa,EAAE,CAAC,iBAAiB,EAAE,qBAAqB,CAAC;KAC1D,CAAC;IACF,KAAK,CAAC;QACJ,EAAE,EAAE,qBAAqB;QACzB,UAAU,EAAE,cAAc;QAC1B,QAAQ,EAAE,uBAAuB;QACjC,WAAW,EAAE,CAAC,GAAG,EAAE,IAAI,CAAC;QACxB,UAAU,EAAE;YACV,EAAE,IAAI,EAAE,GAAG,EAAE,GAAG,EAAE,IAAI,EAAE;YACxB,EAAE,IAAI,EAAE,GAAG,EAAE,GAAG,EAAE,OAAO,EAAE;YAC3B,EAAE,IAAI,EAAE,GAAG,EAAE,GAAG,EAAE,eAAe,EAAE;SACpC;QACD,mBAAmB,EAAE,EAAE;QACvB,aAAa,EAAE,EAAE;KAClB,CAAC;IACF,KAAK,CAAC;QACJ,EAAE,EAAE,WAAW;QACf,UAAU,EAAE,cAAc;QAC1B,QAAQ,EAAE,uBAAuB;QACjC,WAAW,EAAE,CAAC,GAAG,EAAE,IAAI,CAAC;QACxB,UAAU,EAAE;YACV,EAAE,IAAI,EAAE,GAAG,EAAE,GAAG,EAAE,UAAU,EAAE;YAC9B,EAAE,IAAI,EAAE,GAAG,EAAE,GAAG,EAAE,UAAU,EAAE;YAC9B,EAAE,IAAI,EAAE,GAAG,EAAE,GAAG,EAAE,WAAW,EAAE;SAChC;QACD,mBAAmB,EAAE,EAAE;QACvB,aAAa,EAAE,EAAE;KAClB,CAAC;IACF,KAAK,CAAC;QACJ,EAAE,EAAE,gBAAgB;QACpB,UAAU,EAAE,kBAAkB;QAC9B,QAAQ,EAAE,sBAAsB;QAChC,WAAW,EAAE,CAAC,OAAO,EAAE,QAAQ,EAAE,QAAQ,CAAC;QAC1C,UAAU,EAAE;YACV,EAAE,IAAI,EAAE,GAAG,EAAE,GAAG,EAAE,YAAY,EAAE;YAChC,EAAE,IAAI,EAAE,KAAK,EAAE,GAAG,EAAE,MAAM,EAAE;SAC7B;QACD,mBAAmB,EAAE,CAAC,QAAQ,CAAC;QAC/B,aAAa,EAAE,CAAC,oBAAoB,CAAC;KACtC,CAAC;IACF,KAAK,CAAC;QACJ,EAAE,EAAE,aAAa;QACjB,UAAU,EAAE,yBAAyB;QACrC,QAAQ,EAAE,sCAAsC;QAChD,WAAW,EAAE,CAAC,KAAK,EAAE,MAAM,EAAE,YAAY,CAAC;QAC1C,UAAU,EAAE;YACV,EAAE,IAAI,EAAE,GAAG,EAAE,GAAG,EAAE,IAAI,EAAE;YACxB,EAAE,IAAI,EAAE,OAAO,EAAE,GAAG,EAAE,eAAe,EAAE;YACvC,EAAE,IAAI,EAAE,GAAG,EAAE,GAAG,EAAE,MAAM,EAAE;SAC3B;QACD,mBAAmB,EAAE,CAAC,IAAI,CAAC;QAC3B,aAAa,EAAE,EAAE;KAClB,CAAC;IACF,KAAK,CAAC;QACJ,EAAE,EAAE,eAAe;QACnB,UAAU,EAAE,2BAA2B;QACvC,QAAQ,EAAE,gBAAgB;QAC1B,WAAW,EAAE,CAAC,GAAG,EAAE,KAAK,EAAE,YAAY,CAAC;QACvC,UAAU,EAAE,CAAC,EAAE,IAAI,EAAE,GAAG,EAAE,GAAG,EAAE,QAAQ,EAAE,CAAC;QAC1C,mBAAmB,EAAE,EAAE;QACvB,aAAa,EAAE,CAAC,eAAe,CAAC;KACjC,CAAC;IACF,KAAK,CAAC;QACJ,EAAE,EAAE,oBAAoB;QACxB,UAAU,EAAE,cAAc;QAC1B,QAAQ,EAAE,uBAAuB;QACjC,WAAW,EAAE,CAAC,GAAG,EAAE,IAAI,EAAE,4BAA4B,CAAC;QACtD,UAAU,EAAE;YACV,EAAE,IAAI,EAAE,GAAG,EAAE,GAAG,EAAE,IAAI,EAAE;YACxB,EAAE,IAAI,EAAE,GAAG,EAAE,GAAG,EAAE,eAAe,EAAE;YACnC,EAAE,IAAI,EAAE,MAAM,EAAE,GAAG,EAAE,eAAe,EAAE;YACtC,EAAE,IAAI,EAAE,IAAI,EAAE,GAAG,EAAE,MAAM,EAAE;SAC5B;QACD,mBAAmB,EAAE,EAAE;QACvB,aAAa,EAAE,EAAE;KAClB,CAAC;IACF,KAAK,CAAC;QACJ,EAAE,EAAE,mBAAmB;QACvB,UAAU,EAAE,OAAO;QACnB,QAAQ,EAAE,gBAAgB;QAC1B,WAAW,EAAE,CAAC,GAAG,CAAC;QAClB,UAAU,EAAE;YACV,EAAE,IAAI,EAAE,GAAG,EAAE,GAAG,EAAE,OAAO,EAAE;YAC3B,EAAE,IAAI,EAAE,GAAG,EAAE,GAAG,EAAE,eAAe,EAAE;SACpC;QACD,mBAAmB,EAAE,EAAE;QACvB,aAAa,EAAE,EAAE;KAClB,CAAC;CACH,CAAC;AAEF,MAAM,KAAK,GAAoC,IAAI,GAAG,CAAC,YAAY,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;AAE3F,mEAAmE;AACnE,MAAM,UAAU,aAAa,CAAC,EAAU;IACtC,MAAM,KAAK,GAAG,KAAK,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC;IAC5B,IAAI,KAAK,KAAK,SAAS;QAAE,MAAM,IAAI,KAAK,CAAC,2BAA2B,EAAE,EAAE,CAAC,CAAC;IAC1E,OAAO,KAAK,CAAC;AACf,CAAC"}
+{"version":3,"file":"models.js","sourceRoot":"","sources":["../../../src/atlas/oscillators/models.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;GAmBG;AAEH,OAAO,EACL,YAAY,EACZ,MAAM,EACN,IAAI,EACJ,QAAQ,GACT,MAAM,4BAA4B,CAAC;AACpC,OAAO,EAAE,UAAU,EAAE,MAAM,cAAc,CAAC;AAG1C,OAAO,EACL,WAAW,EACX,eAAe,EACf,OAAO,EACP,UAAU,EACV,UAAU,EACV,eAAe,GAChB,MAAM,iBAAiB,CAAC;AAEzB,MAAM,MAAM,GAAG,aAAa,CAAC;AAE7B,SAAS,KAAK,CACZ,IAEC;IAED,OAAO;QACL,GAAG,IAAI;QACP,MAAM,EAAE,MAAM;QACd,MAAM,EAAE,UAAU,CAAC,MAAM,EAAE,IAAI,CAAC,UAAU,EAAE,IAAI,CAAC,mBAAmB,EAAE,EAAE,CAAC;KAC1E,CAAC;AACJ,CAAC;AAED,6EAA6E;AAC7E,MAAM,CAAC,MAAM,YAAY,GAA0B;IACjD,KAAK,CAAC;QACJ,EAAE,EAAE,cAAc;QAClB,UAAU,EAAE,cAAc;QAC1B,QAAQ,EAAE,gBAAgB;QAC1B,WAAW,EAAE,CAAC,GAAG,EAAE,IAAI,EAAE,QAAQ,CAAC;QAClC,UAAU,EAAE;YACV,EAAE,IAAI,EAAE,GAAG,EAAE,GAAG,EAAE,IAAI,EAAE;YACxB,EAAE,IAAI,EAAE,GAAG,EAAE,GAAG,EAAE,eAAe,EAAE;SACpC;QACD,mBAAmB,EAAE,EAAE;QACvB,aAAa,EAAE;YACb,8BAA8B;YAC9B,4BAA4B;YAC5B,sBAAsB;SACvB;KACF,CAAC;IACF,KAAK,CAAC;QACJ,EAAE,EAAE,UAAU;QACd,UAAU,EAAE,cAAc;QAC1B,QAAQ,EAAE,gBAAgB;QAC1B,WAAW,EAAE,CAAC,GAAG,EAAE,IAAI,EAAE,QAAQ,CAAC;QAClC,UAAU,EAAE;YACV,EAAE,IAAI,EAAE,GAAG,EAAE,GAAG,EAAE,UAAU,EAAE;YAC9B,EAAE,IAAI,EAAE,GAAG,EAAE,GAAG,EAAE,WAAW,EAAE;SAChC;QACD,mBAAmB,EAAE,EAAE;QACvB,aAAa,EAAE,CAAC,iBAAiB,EAAE,qBAAqB,CAAC;KAC1D,CAAC;IACF,KAAK,CAAC;QACJ,EAAE,EAAE,qBAAqB;QACzB,UAAU,EAAE,cAAc;QAC1B,QAAQ,EAAE,uBAAuB;QACjC,WAAW,EAAE,CAAC,GAAG,EAAE,IAAI,CAAC;QACxB,UAAU,EAAE;YACV,EAAE,IAAI,EAAE,GAAG,EAAE,GAAG,EAAE,IAAI,EAAE;YACxB,EAAE,IAAI,EAAE,GAAG,EAAE,GAAG,EAAE,OAAO,EAAE;YAC3B,EAAE,IAAI,EAAE,GAAG,EAAE,GAAG,EAAE,eAAe,EAAE;SACpC;QACD,mBAAmB,EAAE,EAAE;QACvB,aAAa,EAAE,EAAE;KAClB,CAAC;IACF,KAAK,CAAC;QACJ,EAAE,EAAE,WAAW;QACf,UAAU,EAAE,cAAc;QAC1B,QAAQ,EAAE,uBAAuB;QACjC,WAAW,EAAE,CAAC,GAAG,EAAE,IAAI,CAAC;QACxB,UAAU,EAAE;YACV,EAAE,IAAI,EAAE,GAAG,EAAE,GAAG,EAAE,UAAU,EAAE;YAC9B,EAAE,IAAI,EAAE,GAAG,EAAE,GAAG,EAAE,UAAU,EAAE;YAC9B,EAAE,IAAI,EAAE,GAAG,EAAE,GAAG,EAAE,WAAW,EAAE;SAChC;QACD,mBAAmB,EAAE,EAAE;QACvB,aAAa,EAAE,EAAE;KAClB,CAAC;IACF,KAAK,CAAC;QACJ,EAAE,EAAE,gBAAgB;QACpB,UAAU,EAAE,kBAAkB;QAC9B,QAAQ,EAAE,sBAAsB;QAChC,WAAW,EAAE,CAAC,OAAO,EAAE,QAAQ,EAAE,QAAQ,CAAC;QAC1C,UAAU,EAAE;YACV,EAAE,IAAI,EAAE,GAAG,EAAE,GAAG,EAAE,YAAY,EAAE;YAChC,EAAE,IAAI,EAAE,KAAK,EAAE,GAAG,EAAE,MAAM,EAAE;SAC7B;QACD,mBAAmB,EAAE,CAAC,QAAQ,CAAC;QAC/B,aAAa,EAAE,CAAC,oBAAoB,CAAC;KACtC,CAAC;IACF,KAAK,CAAC;QACJ,EAAE,EAAE,aAAa;QACjB,UAAU,EAAE,yBAAyB;QACrC,QAAQ,EAAE,sCAAsC;QAChD,WAAW,EAAE,CAAC,KAAK,EAAE,MAAM,EAAE,YAAY,CAAC;QAC1C,UAAU,EAAE;YACV,EAAE,IAAI,EAAE,GAAG,EAAE,GAAG,EAAE,IAAI,EAAE;YACxB,EAAE,IAAI,EAAE,OAAO,EAAE,GAAG,EAAE,eAAe,EAAE;YACvC,EAAE,IAAI,EAAE,GAAG,EAAE,GAAG,EAAE,MAAM,EAAE;SAC3B;QACD,mBAAmB,EAAE,CAAC,IAAI,CAAC;QAC3B,aAAa,EAAE,EAAE;KAClB,CAAC;IACF,KAAK,CAAC;QACJ,EAAE,EAAE,eAAe;QACnB,UAAU,EAAE,2BAA2B;QACvC,QAAQ,EAAE,gBAAgB;QAC1B,WAAW,EAAE,CAAC,GAAG,EAAE,KAAK,EAAE,YAAY,CAAC;QACvC,UAAU,EAAE,CAAC,EAAE,IAAI,EAAE,GAAG,EAAE,GAAG,EAAE,QAAQ,EAAE,CAAC;QAC1C,mBAAmB,EAAE,EAAE;QACvB,aAAa,EAAE,CAAC,eAAe,CAAC;KACjC,CAAC;IACF,KAAK,CAAC;QACJ,EAAE,EAAE,oBAAoB;QACxB,UAAU,EAAE,cAAc;QAC1B,QAAQ,EAAE,uBAAuB;QACjC,WAAW,EAAE,CAAC,GAAG,EAAE,IAAI,EAAE,4BAA4B,CAAC;QACtD,UAAU,EAAE;YACV,EAAE,IAAI,EAAE,GAAG,EAAE,GAAG,EAAE,IAAI,EAAE;YACxB,EAAE,IAAI,EAAE,GAAG,EAAE,GAAG,EAAE,eAAe,EAAE;YACnC,EAAE,IAAI,EAAE,MAAM,EAAE,GAAG,EAAE,eAAe,EAAE;YACtC,EAAE,IAAI,EAAE,IAAI,EAAE,GAAG,EAAE,MAAM,EAAE;SAC5B;QACD,mBAAmB,EAAE,EAAE;QACvB,aAAa,EAAE,EAAE;KAClB,CAAC;IACF,KAAK,CAAC;QACJ,EAAE,EAAE,mBAAmB;QACvB,UAAU,EAAE,OAAO;QACnB,QAAQ,EAAE,gBAAgB;QAC1B,WAAW,EAAE,CAAC,GAAG,CAAC;QAClB,UAAU,EAAE;YACV,EAAE,IAAI,EAAE,GAAG,EAAE,GAAG,EAAE,OAAO,EAAE;YAC3B,EAAE,IAAI,EAAE,GAAG,EAAE,GAAG,EAAE,eAAe,EAAE;SACpC;QACD,mBAAmB,EAAE,EAAE;QACvB,aAAa,EAAE,EAAE;KAClB,CAAC;CACH,CAAC;AAEF,MAAM,KAAK,GAAoC,IAAI,GAAG,CAAC,YAAY,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;AAE3F,mEAAmE;AACnE,MAAM,UAAU,aAAa,CAAC,EAAU;IACtC,MAAM,KAAK,GAAG,KAAK,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC;IAC5B,IAAI,KAAK,KAAK,SAAS;QAAE,MAAM,IAAI,KAAK,CAAC,2BAA2B,EAAE,EAAE,CAAC,CAAC;IAC1E,OAAO,KAAK,CAAC;AACf,CAAC","sourcesContent":["/**\n * The nine oscillator models of design note §3.\n *\n * Each model's `regime` carries the π-groups `deriveRegimeGroups` produces\n * from its own parameters and dimensionless inputs. `inequalities` is EMPTY\n * on every model: Phase 0 states its inequalities on the BRIDGES (which carry\n * the regime of validity), and the design note fixes no per-model inequality.\n * An invented one would be a fabricated constraint, so none is written.\n *\n * The Phase 3 fields `boundaryData`, `initialData` and `symmetryGroup`\n * (`../model.ts`) are ABSENT on all nine for the same reason `inequalities`\n * is empty: no source records them for these ODE models, and an absent field\n * says \"not recorded\" where a present empty one would claim \"there are none\".\n *\n * `canonicalRefs` are restricted to ids that exist in `CANONICAL_EQUATIONS`;\n * `tests/atlas/models.test.ts` resolves every one of them, and a new one needs\n * a numeric check in `tests/atlas/canonical-links.test.ts`.\n *\n * @module atlas/oscillators/models\n */\n\nimport {\n  ACCELERATION,\n  LENGTH,\n  MASS,\n  VELOCITY,\n} from '../../dimensional/types.js';\nimport { makeRegime } from '../regime.js';\nimport type { AtlasModel } from '../model.js';\nimport type { DimensionalVariable } from '../../dimensional/buckingham.js';\nimport {\n  CAPACITANCE,\n  CUBIC_STIFFNESS,\n  DAMPING,\n  INDUCTANCE,\n  RESISTANCE,\n  SPRING_CONSTANT,\n} from './dimensions.js';\n\nconst FAMILY = 'oscillators';\n\nfunction model(\n  spec: Omit<AtlasModel, 'family' | 'regime'> & {\n    readonly parameters: readonly DimensionalVariable[];\n  },\n): AtlasModel {\n  return {\n    ...spec,\n    family: FAMILY,\n    regime: makeRegime(FAMILY, spec.parameters, spec.dimensionlessInputs, []),\n  };\n}\n\n/** The nine Phase 0 oscillator models, in design-note §3 order. @internal */\nexport const ATLAS_MODELS: readonly AtlasModel[] = [\n  model({\n    id: 'model-spring',\n    stateSpace: '(x, x′) ∈ ℝ²',\n    dynamics: 'm x″ + k x = 0',\n    observables: ['x', 'x′', 'energy'],\n    parameters: [\n      { name: 'm', dim: MASS },\n      { name: 'k', dim: SPRING_CONSTANT },\n    ],\n    dimensionlessInputs: [],\n    canonicalRefs: [\n      'CE-simple-harmonic-frequency',\n      'CE-spring-potential-energy',\n      'CE-oscillator-energy',\n    ],\n  }),\n  model({\n    id: 'model-lc',\n    stateSpace: '(q, q′) ∈ ℝ²',\n    dynamics: 'L q″ + q/C = 0',\n    observables: ['q', 'q′', 'energy'],\n    parameters: [\n      { name: 'L', dim: INDUCTANCE },\n      { name: 'C', dim: CAPACITANCE },\n    ],\n    dimensionlessInputs: [],\n    canonicalRefs: ['CE-lc-resonance', 'CE-capacitor-energy'],\n  }),\n  model({\n    id: 'model-damped-spring',\n    stateSpace: '(x, x′) ∈ ℝ²',\n    dynamics: 'm x″ + b x′ + k x = 0',\n    observables: ['x', 'x′'],\n    parameters: [\n      { name: 'm', dim: MASS },\n      { name: 'b', dim: DAMPING },\n      { name: 'k', dim: SPRING_CONSTANT },\n    ],\n    dimensionlessInputs: [],\n    canonicalRefs: [],\n  }),\n  model({\n    id: 'model-rlc',\n    stateSpace: '(q, q′) ∈ ℝ²',\n    dynamics: 'L q″ + R q′ + q/C = 0',\n    observables: ['q', 'q′'],\n    parameters: [\n      { name: 'L', dim: INDUCTANCE },\n      { name: 'R', dim: RESISTANCE },\n      { name: 'C', dim: CAPACITANCE },\n    ],\n    dimensionlessInputs: [],\n    canonicalRefs: [],\n  }),\n  model({\n    id: 'model-pendulum',\n    stateSpace: '(θ, θ′) ∈ S¹ × ℝ',\n    dynamics: 'θ″ + (g/ℓ) sin θ = 0',\n    observables: ['theta', 'theta′', 'period'],\n    parameters: [\n      { name: 'g', dim: ACCELERATION },\n      { name: 'ell', dim: LENGTH },\n    ],\n    dimensionlessInputs: ['theta0'],\n    canonicalRefs: ['CE-pendulum-period'],\n  }),\n  model({\n    id: 'model-chain',\n    stateSpace: '(u_n, u_n′) ∈ ℝ^N × ℝ^N',\n    dynamics: 'm u_n″ = κ(u_{n+1} − 2u_n + u_{n−1})',\n    observables: ['u_n', 'u_n′', 'dispersion'],\n    parameters: [\n      { name: 'm', dim: MASS },\n      { name: 'kappa', dim: SPRING_CONSTANT },\n      { name: 'a', dim: LENGTH },\n    ],\n    dimensionlessInputs: ['qa'],\n    canonicalRefs: [],\n  }),\n  model({\n    id: 'model-wave-1d',\n    stateSpace: 'u(x, t), a field on ℝ × ℝ',\n    dynamics: 'u_tt = c² u_xx',\n    observables: ['u', 'u_t', 'dispersion'],\n    parameters: [{ name: 'c', dim: VELOCITY }],\n    dimensionlessInputs: [],\n    canonicalRefs: ['CE-wave-speed'],\n  }),\n  model({\n    id: 'model-cubic-spring',\n    stateSpace: '(x, x′) ∈ ℝ²',\n    dynamics: 'm x″ + k x + β x³ = 0',\n    observables: ['x', 'x′', 'amplitude-dependent period'],\n    parameters: [\n      { name: 'm', dim: MASS },\n      { name: 'k', dim: SPRING_CONSTANT },\n      { name: 'beta', dim: CUBIC_STIFFNESS },\n      { name: 'x0', dim: LENGTH },\n    ],\n    dimensionlessInputs: [],\n    canonicalRefs: [],\n  }),\n  model({\n    id: 'model-first-order',\n    stateSpace: 'x ∈ ℝ',\n    dynamics: 'b x′ + k x = 0',\n    observables: ['x'],\n    parameters: [\n      { name: 'b', dim: DAMPING },\n      { name: 'k', dim: SPRING_CONSTANT },\n    ],\n    dimensionlessInputs: [],\n    canonicalRefs: [],\n  }),\n];\n\nconst BY_ID: ReadonlyMap<string, AtlasModel> = new Map(ATLAS_MODELS.map((m) => [m.id, m]));\n\n/** Look up a model by id, throwing on an unknown one. @internal */\nexport function getAtlasModel(id: string): AtlasModel {\n  const found = BY_ID.get(id);\n  if (found === undefined) throw new Error(`unknown atlas model id: ${id}`);\n  return found;\n}\n"]}