universal-physics-tensor 10.0.0 → 11.0.0

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","sourcesContent":["/**\n * Dimension constants the oscillator pilot needs that the canonical entry\n * files declare MODULE-LOCAL and do not export.\n *\n * They are redefined here with the same `dim(...)` arguments, quoted from\n * source in design note §10 Q2. `tests/atlas/models.test.ts` asserts the\n * values, so a future divergence in an entry file fails a test.\n *\n * `FREQUENCY`, `MASS`, `LENGTH`, `VELOCITY` and `ACCELERATION` are importable\n * from `src/dimensional/types.ts` and are NOT redefined 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, `L^-2 M^-1 T^4 I^2` — `src/canonical/entries/electromagnetism.ts:42`. @internal */\nexport const CAPACITANCE: Dimension = dim(-2, -1, 4, 2);\n\n/** henry, `L^2 M T^-2 I^-2` — `src/canonical/entries/electromagnetism.ts:43`. @internal */\nexport const INDUCTANCE: Dimension = dim(2, 1, -2, -2);\n\n/** ohm, `L^2 M T^-3 I^-2` — `src/canonical/entries/electromagnetism.ts:40`. @internal */\nexport const RESISTANCE: Dimension = dim(2, 1, -3, -2);\n\n/**\n * `M T^-2` — defined TWICE, identically, at\n * `src/canonical/entries/mechanics.ts:35` and\n * `src/canonical/entries/fluids-waves.ts:48`. One definition here; the test\n * asserts agreement with both.\n *\n * @internal\n */\nexport const SPRING_CONSTANT: Dimension = dim(0, 1, -2);\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"]}
+{"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/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","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 MathTS\n * `solveODESystem` and `dt`, 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 { solveODESystem } from '@danielsimonjr/mathts-functions';\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 sol = solveODESystem(\n    (_t, y) => [y[1]!, -(b * y[1]! + k * y[0]!) / m],\n    [x0, v0],\n    [0, steps * h],\n    { dt: h },\n  );\n  let sup = 0;\n  for (let i = layerSteps; i <= steps; i++) {\n    const x = sol.y[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  let state = [\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  for (let i = 0; i < maxSteps; i++) {\n    const sol = solveODESystem(\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      state,\n      [0, h],\n      { dt: h },\n    );\n    const next = sol.y[sol.y.length - 1];\n    if (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    state = next.slice();\n  }\n  return Number.NaN;\n}\n"]}
+{"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","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 { deriveRegimeGroups } 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: {\n      family: FAMILY,\n      inequalities: [],\n      groupDefinitions: deriveRegimeGroups(FAMILY, spec.parameters, spec.dimensionlessInputs),\n    },\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"]}
+{"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"]}

package/dist/atlas/oscillators/pendulum-motion.d.ts

@@ -9,6 +9,8 @@
  * @module atlas/oscillators/pendulum-motion
  * @internal
  */
+/** 2π: one period of the angle variable, the T0 = 1 time unit's ω0. @internal */
+export declare const TWO_PI: number;
 /** ω0² at T0 = 1. @internal */
 export declare const W2: number;
 /** @internal */
@@ -16,8 +18,8 @@ export declare const pendulumAccel: (x: number) => number;
 /** @internal */
 export declare const linearAccel: (x: number) => number;
 /**
- * One fixed step of x″ = accel(x). MathTS `solveODESystem` with `dt` is the
- * classical RK4 step.
+ * One fixed step of x″ = accel(x): `integrateRK4` (MathTS `solveODESystem`
+ * with `dt`, the classical RK4 step) over one step of length `h`.
  *
  * @internal
  */
@@ -27,17 +29,38 @@ export declare function crossingTimes(accel: (x: number) => number, theta0: numb
 /** The phase at `t`, interpolated through the anchors (0, 0) and (c_k, π/2 + kπ). @internal */
 export declare function phaseAt(crossings: readonly number[], t: number): number;
 /**
- * θ0 in [0, π) whose exact relative period error is `eps`, by bisection on the
- * strictly increasing `pendulumPeriodErrorAt`.
+ * θ0 in [0, π) whose exact relative period error is `eps`, by MathTS
+ * `findRoot` (Brent) on the strictly increasing `pendulumPeriodErrorAt`.
+ *
+ * The residual is solved as `err/eps − 1`, of order 1, so `findRoot`'s one
+ * `tol` (it stops on the residual AND on the bracket width) cannot stop on a
+ * residual that is small only because `eps` is. The upper end is just below
+ * π, where the error is finite; at π it is `Infinity`, which Brent's
+ * interpolation cannot take.
  *
  * @returns NaN for a negative or non-finite `eps`.
  * @internal
  */
 export declare function theta0OfPeriodError(eps: number): number;
 /**
- * The nome q = exp(−π K′/K) of the modulus k = sin(θ0/2), with
- * K′/K = AGM(1, cos(θ0/2)) / AGM(1, sin(θ0/2)).
+ * The nome q = exp(−π K′/K) of the modulus k = sin(θ0/2): K = K(m) and
+ * K′ = K(1 − m) at the parameter m = k², both MathTS `ellipticK` (whose
+ * route is the AGM, so K′/K = AGM(1, cos(θ0/2)) / AGM(1, sin(θ0/2))).
  * @internal
  */
 export declare function pendulumNome(theta0: number): number;
+/**
+ * The exact pendulum released from rest at `theta0`, in units of T0 = 1
+ * (ω0 = 2π): θ(t) = 2 asin(k sn(K − ω0 t, m)), with k = sin(θ0/2), the
+ * parameter m = k², and K = K(m) the quarter period in the sn argument, read
+ * from `pendulumPeriodErrorAt` as (π/2)(1 + ε) so one `ellipticK` read
+ * serves both. At t = 0, sn K = 1 and θ = θ0; the period is T = 4K/ω0 =
+ * (1 + ε) T0. MathTS `jacobiSN` is a plain number function of (u, m). A
+ * closed form, independent of every integrator: `tests/atlas/pendulum-exact.test.ts`
+ * checks it against the cn·dn identities (θ″ = −2ω0² k sn u dn u, sin θ =
+ * 2k sn u dn u) and against RK4.
+ *
+ * @internal
+ */
+export declare function exactPendulumAngle(theta0: number, t: number): number;
 //# sourceMappingURL=pendulum-motion.d.ts.map

package/dist/atlas/oscillators/pendulum-motion.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"pendulum-motion.d.ts","sourceRoot":"","sources":["../../../src/atlas/oscillators/pendulum-motion.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;GAUG;AAKH,+BAA+B;AAC/B,eAAO,MAAM,EAAE,QAAqB,CAAC;AACrC,gBAAgB;AAChB,eAAO,MAAM,aAAa,MAAO,MAAM,KAAG,MAA2B,CAAC;AACtE,gBAAgB;AAChB,eAAO,MAAM,WAAW,MAAO,MAAM,KAAG,MAAiB,CAAC;AAE1D;;;;;GAKG;AACH,wBAAgB,OAAO,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,KAAK,EAAE,CAAC,CAAC,EAAE,MAAM,KAAK,MAAM,GAAG,CAAC,MAAM,EAAE,MAAM,CAAC,CAIvG;AA4BD,kGAAkG;AAClG,wBAAgB,aAAa,CAAC,KAAK,EAAE,CAAC,CAAC,EAAE,MAAM,KAAK,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,cAAc,EAAE,MAAM,EAAE,OAAO,EAAE,MAAM,GAAG,MAAM,EAAE,CAY7H;AAED,+FAA+F;AAC/F,wBAAgB,OAAO,CAAC,SAAS,EAAE,SAAS,MAAM,EAAE,EAAE,CAAC,EAAE,MAAM,GAAG,MAAM,CAWvE;AAED;;;;;;GAMG;AACH,wBAAgB,mBAAmB,CAAC,GAAG,EAAE,MAAM,GAAG,MAAM,CAYvD;AAED;;;;GAIG;AACH,wBAAgB,YAAY,CAAC,MAAM,EAAE,MAAM,GAAG,MAAM,CAInD"}
+{"version":3,"file":"pendulum-motion.d.ts","sourceRoot":"","sources":["../../../src/atlas/oscillators/pendulum-motion.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;GAUG;AAOH,iFAAiF;AACjF,eAAO,MAAM,MAAM,QAAc,CAAC;AAClC,+BAA+B;AAC/B,eAAO,MAAM,EAAE,QAAc,CAAC;AAC9B,gBAAgB;AAChB,eAAO,MAAM,aAAa,MAAO,MAAM,KAAG,MAA2B,CAAC;AACtE,gBAAgB;AAChB,eAAO,MAAM,WAAW,MAAO,MAAM,KAAG,MAAiB,CAAC;AAE1D;;;;;GAKG;AACH,wBAAgB,OAAO,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,KAAK,EAAE,CAAC,CAAC,EAAE,MAAM,KAAK,MAAM,GAAG,CAAC,MAAM,EAAE,MAAM,CAAC,CAGvG;AAqBD,kGAAkG;AAClG,wBAAgB,aAAa,CAAC,KAAK,EAAE,CAAC,CAAC,EAAE,MAAM,KAAK,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,cAAc,EAAE,MAAM,EAAE,OAAO,EAAE,MAAM,GAAG,MAAM,EAAE,CAY7H;AAED,+FAA+F;AAC/F,wBAAgB,OAAO,CAAC,SAAS,EAAE,SAAS,MAAM,EAAE,EAAE,CAAC,EAAE,MAAM,GAAG,MAAM,CAWvE;AAED;;;;;;;;;;;;GAYG;AACH,wBAAgB,mBAAmB,CAAC,GAAG,EAAE,MAAM,GAAG,MAAM,CASvD;AAED;;;;;GAKG;AACH,wBAAgB,YAAY,CAAC,MAAM,EAAE,MAAM,GAAG,MAAM,CAMnD;AAED;;;;;;;;;;;;GAYG;AACH,wBAAgB,kBAAkB,CAAC,MAAM,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,GAAG,MAAM,CAKpE"}

package/dist/atlas/oscillators/pendulum-motion.js

@@ -9,50 +9,46 @@
  * @module atlas/oscillators/pendulum-motion
  * @internal
  */
-import { solveODESystem } from '@danielsimonjr/mathts-functions';
-import { agm, pendulumPeriodErrorAt } from './bridges-limits.js';
+import { ellipticK, findRoot, jacobiSN } from '@danielsimonjr/mathts-functions';
+import { numberFrom } from '../../core/mathts-read.js';
+import { integrateRK4 } from '../../numerical/null-ray-integrator.js';
+import { pendulumPeriodErrorAt } from './bridges-limits.js';
+/** 2π: one period of the angle variable, the T0 = 1 time unit's ω0. @internal */
+export const TWO_PI = 2 * Math.PI;
 /** ω0² at T0 = 1. @internal */
-export const W2 = (2 * Math.PI) ** 2;
+export const W2 = TWO_PI ** 2;
 /** @internal */
 export const pendulumAccel = (x) => -W2 * Math.sin(x);
 /** @internal */
 export const linearAccel = (x) => -W2 * x;
 /**
- * One fixed step of x″ = accel(x). MathTS `solveODESystem` with `dt` is the
- * classical RK4 step.
+ * One fixed step of x″ = accel(x): `integrateRK4` (MathTS `solveODESystem`
+ * with `dt`, the classical RK4 step) over one step of length `h`.
  *
  * @internal
  */
 export function rk4Step(x, v, h, accel) {
-    const sol = solveODESystem((_t, y) => [y[1], accel(y[0])], [x, v], [0, h], { dt: h });
-    const end = sol.y[sol.y.length - 1] ?? [x, v];
+    const end = integrateRK4((_t, y) => [y[1], accel(y[0])], [x, v], 0, h, 1);
     return [end[0], end[1]];
 }
 /**
  * The root in (0, h] of x(τ), the position one partial RK4 step of length τ
- * from (x, v), by the secant method from the bracket (0, x) and (h, nx).
+ * from (x, v), by MathTS `findRoot` (Brent) on the bracket (0, x), (h, nx).
  *
  * The crossing is then located to the integrator's own truncation error,
  * rather than to the O(h³) error of a straight line between the two steps.
  * That truncation error is nearly the same for the pendulum and the linear
  * oscillator at small θ0, so it cancels in their phase difference; the
- * interpolation error does not.
+ * interpolation error does not. `findRoot` evaluates both ends, and a step
+ * of length 0 is the start itself (`solveODESystem` cannot take `dt = 0`).
+ * Its `tol` stops on the residual (an angle) and on the bracket width (a
+ * time in T0), both at 1e-12, below the integrator's error either way.
  */
 function refineCrossing(x, v, nx, h, accel) {
-    let a = 0;
-    let fa = x;
-    let b = h;
-    let fb = nx;
-    for (let k = 0; k < 50 && fb !== 0 && fb !== fa; k++) {
-        const c = b - (fb * (b - a)) / (fb - fa);
-        if (!(c >= -h && c <= 2 * h) || c === b)
-            break;
-        a = b;
-        fa = fb;
-        b = c;
-        fb = rk4Step(x, v, c, accel)[0];
-    }
-    return b;
+    if (nx === 0)
+        return h;
+    const position = (tau) => (tau === 0 ? x : tau === h ? nx : rk4Step(x, v, tau, accel)[0]);
+    return findRoot(position, 0, h, { tol: 1e-12, maxIter: 200 });
 }
 /** Zero-crossing times of a release from rest at `theta0`, over `periods` of T0 = 1. @internal */
 export function crossingTimes(accel, theta0, stepsPerPeriod, periods) {
@@ -84,8 +80,14 @@ export function phaseAt(crossings, t) {
     return Number.NaN;
 }
 /**
- * θ0 in [0, π) whose exact relative period error is `eps`, by bisection on the
- * strictly increasing `pendulumPeriodErrorAt`.
+ * θ0 in [0, π) whose exact relative period error is `eps`, by MathTS
+ * `findRoot` (Brent) on the strictly increasing `pendulumPeriodErrorAt`.
+ *
+ * The residual is solved as `err/eps − 1`, of order 1, so `findRoot`'s one
+ * `tol` (it stops on the residual AND on the bracket width) cannot stop on a
+ * residual that is small only because `eps` is. The upper end is just below
+ * π, where the error is finite; at π it is `Infinity`, which Brent's
+ * interpolation cannot take.
  *
  * @returns NaN for a negative or non-finite `eps`.
  * @internal
@@ -95,28 +97,46 @@ export function theta0OfPeriodError(eps) {
         return Number.NaN;
     if (eps === 0)
         return 0;
-    let lo = 0;
-    let hi = Math.PI;
-    for (let i = 0; i < 200 && hi - lo > 0; i++) {
-        const mid = (lo + hi) / 2;
-        if (mid === lo || mid === hi)
-            break;
-        if (pendulumPeriodErrorAt({ theta0: mid }) < eps)
-            lo = mid;
-        else
-            hi = mid;
-    }
-    return (lo + hi) / 2;
+    const hi = Math.PI * (1 - Number.EPSILON);
+    if (pendulumPeriodErrorAt({ theta0: hi }) <= eps)
+        return hi;
+    return findRoot((theta0) => pendulumPeriodErrorAt({ theta0 }) / eps - 1, 0, hi, {
+        tol: 4 * Number.EPSILON,
+        maxIter: 200,
+    });
 }
 /**
- * The nome q = exp(−π K′/K) of the modulus k = sin(θ0/2), with
- * K′/K = AGM(1, cos(θ0/2)) / AGM(1, sin(θ0/2)).
+ * The nome q = exp(−π K′/K) of the modulus k = sin(θ0/2): K = K(m) and
+ * K′ = K(1 − m) at the parameter m = k², both MathTS `ellipticK` (whose
+ * route is the AGM, so K′/K = AGM(1, cos(θ0/2)) / AGM(1, sin(θ0/2))).
  * @internal
  */
 export function pendulumNome(theta0) {
     const k = Math.sin(Math.abs(theta0) / 2);
     if (k === 0)
         return 0;
-    return Math.exp((-Math.PI * agm(1, Math.cos(theta0 / 2))) / agm(1, k));
+    const K = numberFrom('ellipticK', ellipticK(k * k));
+    const Kp = numberFrom('ellipticK', ellipticK(Math.cos(theta0 / 2) ** 2));
+    return Math.exp((-Math.PI * Kp) / K);
+}
+/**
+ * The exact pendulum released from rest at `theta0`, in units of T0 = 1
+ * (ω0 = 2π): θ(t) = 2 asin(k sn(K − ω0 t, m)), with k = sin(θ0/2), the
+ * parameter m = k², and K = K(m) the quarter period in the sn argument, read
+ * from `pendulumPeriodErrorAt` as (π/2)(1 + ε) so one `ellipticK` read
+ * serves both. At t = 0, sn K = 1 and θ = θ0; the period is T = 4K/ω0 =
+ * (1 + ε) T0. MathTS `jacobiSN` is a plain number function of (u, m). A
+ * closed form, independent of every integrator: `tests/atlas/pendulum-exact.test.ts`
+ * checks it against the cn·dn identities (θ″ = −2ω0² k sn u dn u, sin θ =
+ * 2k sn u dn u) and against RK4.
+ *
+ * @internal
+ */
+export function exactPendulumAngle(theta0, t) {
+    const k = Math.sin(theta0 / 2);
+    if (k === 0)
+        return 0;
+    const K = (Math.PI / 2) * (1 + pendulumPeriodErrorAt({ theta0 }));
+    return 2 * Math.asin(k * jacobiSN(K - TWO_PI * t, k * k));
 }
 //# sourceMappingURL=pendulum-motion.js.map

package/dist/atlas/oscillators/pendulum-motion.js.map

@@ -1 +1 @@
-{"version":3,"file":"pendulum-motion.js","sourceRoot":"","sources":["../../../src/atlas/oscillators/pendulum-motion.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;GAUG;AAEH,OAAO,EAAE,cAAc,EAAE,MAAM,iCAAiC,CAAC;AACjE,OAAO,EAAE,GAAG,EAAE,qBAAqB,EAAE,MAAM,qBAAqB,CAAC;AAEjE,+BAA+B;AAC/B,MAAM,CAAC,MAAM,EAAE,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC,IAAI,CAAC,CAAC;AACrC,gBAAgB;AAChB,MAAM,CAAC,MAAM,aAAa,GAAG,CAAC,CAAS,EAAU,EAAE,CAAC,CAAC,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;AACtE,gBAAgB;AAChB,MAAM,CAAC,MAAM,WAAW,GAAG,CAAC,CAAS,EAAU,EAAE,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC;AAE1D;;;;;GAKG;AACH,MAAM,UAAU,OAAO,CAAC,CAAS,EAAE,CAAS,EAAE,CAAS,EAAE,KAA4B;IACnF,MAAM,GAAG,GAAG,cAAc,CAAC,CAAC,EAAE,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAE,EAAE,KAAK,CAAC,CAAC,CAAC,CAAC,CAAE,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,EAAE,EAAE,EAAE,CAAC,EAAE,CAAC,CAAC;IACxF,MAAM,GAAG,GAAG,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,MAAM,GAAG,CAAC,CAAC,IAAI,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;IAC9C,OAAO,CAAC,GAAG,CAAC,CAAC,CAAE,EAAE,GAAG,CAAC,CAAC,CAAE,CAAC,CAAC;AAC5B,CAAC;AAED;;;;;;;;;GASG;AACH,SAAS,cAAc,CAAC,CAAS,EAAE,CAAS,EAAE,EAAU,EAAE,CAAS,EAAE,KAA4B;IAC/F,IAAI,CAAC,GAAG,CAAC,CAAC;IACV,IAAI,EAAE,GAAG,CAAC,CAAC;IACX,IAAI,CAAC,GAAG,CAAC,CAAC;IACV,IAAI,EAAE,GAAG,EAAE,CAAC;IACZ,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,EAAE,IAAI,EAAE,KAAK,CAAC,IAAI,EAAE,KAAK,EAAE,EAAE,CAAC,EAAE,EAAE,CAAC;QACrD,MAAM,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,EAAE,GAAG,EAAE,CAAC,CAAC;QACzC,IAAI,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,IAAI,CAAC,KAAK,CAAC;YAAE,MAAM;QAC/C,CAAC,GAAG,CAAC,CAAC;QACN,EAAE,GAAG,EAAE,CAAC;QACR,CAAC,GAAG,CAAC,CAAC;QACN,EAAE,GAAG,OAAO,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC;IAClC,CAAC;IACD,OAAO,CAAC,CAAC;AACX,CAAC;AAED,kGAAkG;AAClG,MAAM,UAAU,aAAa,CAAC,KAA4B,EAAE,MAAc,EAAE,cAAsB,EAAE,OAAe;IACjH,MAAM,CAAC,GAAG,CAAC,GAAG,cAAc,CAAC;IAC7B,MAAM,GAAG,GAAa,EAAE,CAAC;IACzB,IAAI,CAAC,GAAG,MAAM,CAAC;IACf,IAAI,CAAC,GAAG,CAAC,CAAC;IACV,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,cAAc,GAAG,OAAO,EAAE,CAAC,EAAE,EAAE,CAAC;QAClD,MAAM,CAAC,EAAE,EAAE,EAAE,CAAC,GAAG,OAAO,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,KAAK,CAAC,CAAC;QACzC,IAAI,EAAE,KAAK,CAAC,IAAI,CAAC,GAAG,CAAC,KAAK,EAAE,GAAG,CAAC;YAAE,GAAG,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,GAAG,cAAc,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,EAAE,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC;QACvF,CAAC,GAAG,EAAE,CAAC;QACP,CAAC,GAAG,EAAE,CAAC;IACT,CAAC;IACD,OAAO,GAAG,CAAC;AACb,CAAC;AAED,+FAA+F;AAC/F,MAAM,UAAU,OAAO,CAAC,SAA4B,EAAE,CAAS;IAC7D,IAAI,KAAK,GAAG,CAAC,CAAC;IACd,IAAI,SAAS,GAAG,CAAC,CAAC;IAClB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,SAAS,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;QAC1C,MAAM,CAAC,GAAG,SAAS,CAAC,CAAC,CAAE,CAAC;QACxB,MAAM,KAAK,GAAG,IAAI,CAAC,EAAE,GAAG,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC;QACxC,IAAI,CAAC,IAAI,CAAC;YAAE,OAAO,SAAS,GAAG,CAAC,CAAC,KAAK,GAAG,SAAS,CAAC,GAAG,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC;QACjF,KAAK,GAAG,CAAC,CAAC;QACV,SAAS,GAAG,KAAK,CAAC;IACpB,CAAC;IACD,OAAO,MAAM,CAAC,GAAG,CAAC;AACpB,CAAC;AAED;;;;;;GAMG;AACH,MAAM,UAAU,mBAAmB,CAAC,GAAW;IAC7C,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,GAAG,CAAC,IAAI,GAAG,GAAG,CAAC;QAAE,OAAO,MAAM,CAAC,GAAG,CAAC;IACxD,IAAI,GAAG,KAAK,CAAC;QAAE,OAAO,CAAC,CAAC;IACxB,IAAI,EAAE,GAAG,CAAC,CAAC;IACX,IAAI,EAAE,GAAG,IAAI,CAAC,EAAE,CAAC;IACjB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,GAAG,IAAI,EAAE,GAAG,EAAE,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;QAC5C,MAAM,GAAG,GAAG,CAAC,EAAE,GAAG,EAAE,CAAC,GAAG,CAAC,CAAC;QAC1B,IAAI,GAAG,KAAK,EAAE,IAAI,GAAG,KAAK,EAAE;YAAE,MAAM;QACpC,IAAI,qBAAqB,CAAC,EAAE,MAAM,EAAE,GAAG,EAAE,CAAC,GAAG,GAAG;YAAE,EAAE,GAAG,GAAG,CAAC;;YACtD,EAAE,GAAG,GAAG,CAAC;IAChB,CAAC;IACD,OAAO,CAAC,EAAE,GAAG,EAAE,CAAC,GAAG,CAAC,CAAC;AACvB,CAAC;AAED;;;;GAIG;AACH,MAAM,UAAU,YAAY,CAAC,MAAc;IACzC,MAAM,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC;IACzC,IAAI,CAAC,KAAK,CAAC;QAAE,OAAO,CAAC,CAAC;IACtB,OAAO,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,IAAI,CAAC,EAAE,GAAG,GAAG,CAAC,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;AACzE,CAAC","sourcesContent":["/**\n * The pendulum and its linear oscillator, integrated and read off, for the\n * witnesses of `ab-pendulum-linear`'s declared translations. Time is in units\n * of T0 (ω0 = 2π); both motions are released from rest at θ0.\n *\n * Also the closed forms those translations are built from: θ0 recovered from\n * the bound's ε, and the nome q of the pendulum's elliptic modulus.\n *\n * @module atlas/oscillators/pendulum-motion\n * @internal\n */\n\nimport { solveODESystem } from '@danielsimonjr/mathts-functions';\nimport { agm, pendulumPeriodErrorAt } from './bridges-limits.js';\n\n/** ω0² at T0 = 1. @internal */\nexport const W2 = (2 * Math.PI) ** 2;\n/** @internal */\nexport const pendulumAccel = (x: number): number => -W2 * Math.sin(x);\n/** @internal */\nexport const linearAccel = (x: number): number => -W2 * x;\n\n/**\n * One fixed step of x″ = accel(x). MathTS `solveODESystem` with `dt` is the\n * classical RK4 step.\n *\n * @internal\n */\nexport function rk4Step(x: number, v: number, h: number, accel: (x: number) => number): [number, number] {\n  const sol = solveODESystem((_t, y) => [y[1]!, accel(y[0]!)], [x, v], [0, h], { dt: h });\n  const end = sol.y[sol.y.length - 1] ?? [x, v];\n  return [end[0]!, end[1]!];\n}\n\n/**\n * The root in (0, h] of x(τ), the position one partial RK4 step of length τ\n * from (x, v), by the secant method from the bracket (0, x) and (h, nx).\n *\n * The crossing is then located to the integrator's own truncation error,\n * rather than to the O(h³) error of a straight line between the two steps.\n * That truncation error is nearly the same for the pendulum and the linear\n * oscillator at small θ0, so it cancels in their phase difference; the\n * interpolation error does not.\n */\nfunction refineCrossing(x: number, v: number, nx: number, h: number, accel: (x: number) => number): number {\n  let a = 0;\n  let fa = x;\n  let b = h;\n  let fb = nx;\n  for (let k = 0; k < 50 && fb !== 0 && fb !== fa; k++) {\n    const c = b - (fb * (b - a)) / (fb - fa);\n    if (!(c >= -h && c <= 2 * h) || c === b) break;\n    a = b;\n    fa = fb;\n    b = c;\n    fb = rk4Step(x, v, c, accel)[0];\n  }\n  return b;\n}\n\n/** Zero-crossing times of a release from rest at `theta0`, over `periods` of T0 = 1. @internal */\nexport function crossingTimes(accel: (x: number) => number, theta0: number, stepsPerPeriod: number, periods: number): number[] {\n  const h = 1 / stepsPerPeriod;\n  const out: number[] = [];\n  let x = theta0;\n  let v = 0;\n  for (let i = 0; i < stepsPerPeriod * periods; i++) {\n    const [nx, nv] = rk4Step(x, v, h, accel);\n    if (nx === 0 || x > 0 !== nx > 0) out.push(i * h + refineCrossing(x, v, nx, h, accel));\n    x = nx;\n    v = nv;\n  }\n  return out;\n}\n\n/** The phase at `t`, interpolated through the anchors (0, 0) and (c_k, π/2 + kπ). @internal */\nexport function phaseAt(crossings: readonly number[], t: number): number {\n  let prevT = 0;\n  let prevPhase = 0;\n  for (let k = 0; k < crossings.length; k++) {\n    const c = crossings[k]!;\n    const phase = Math.PI / 2 + k * Math.PI;\n    if (t <= c) return prevPhase + ((phase - prevPhase) * (t - prevT)) / (c - prevT);\n    prevT = c;\n    prevPhase = phase;\n  }\n  return Number.NaN;\n}\n\n/**\n * θ0 in [0, π) whose exact relative period error is `eps`, by bisection on the\n * strictly increasing `pendulumPeriodErrorAt`.\n *\n * @returns NaN for a negative or non-finite `eps`.\n * @internal\n */\nexport function theta0OfPeriodError(eps: number): number {\n  if (!Number.isFinite(eps) || eps < 0) return Number.NaN;\n  if (eps === 0) return 0;\n  let lo = 0;\n  let hi = Math.PI;\n  for (let i = 0; i < 200 && hi - lo > 0; i++) {\n    const mid = (lo + hi) / 2;\n    if (mid === lo || mid === hi) break;\n    if (pendulumPeriodErrorAt({ theta0: mid }) < eps) lo = mid;\n    else hi = mid;\n  }\n  return (lo + hi) / 2;\n}\n\n/**\n * The nome q = exp(−π K′/K) of the modulus k = sin(θ0/2), with\n * K′/K = AGM(1, cos(θ0/2)) / AGM(1, sin(θ0/2)).\n * @internal\n */\nexport function pendulumNome(theta0: number): number {\n  const k = Math.sin(Math.abs(theta0) / 2);\n  if (k === 0) return 0;\n  return Math.exp((-Math.PI * agm(1, Math.cos(theta0 / 2))) / agm(1, k));\n}\n"]}
+{"version":3,"file":"pendulum-motion.js","sourceRoot":"","sources":["../../../src/atlas/oscillators/pendulum-motion.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;GAUG;AAEH,OAAO,EAAE,SAAS,EAAE,QAAQ,EAAE,QAAQ,EAAE,MAAM,iCAAiC,CAAC;AAChF,OAAO,EAAE,UAAU,EAAE,MAAM,2BAA2B,CAAC;AACvD,OAAO,EAAE,YAAY,EAAE,MAAM,wCAAwC,CAAC;AACtE,OAAO,EAAE,qBAAqB,EAAE,MAAM,qBAAqB,CAAC;AAE5D,iFAAiF;AACjF,MAAM,CAAC,MAAM,MAAM,GAAG,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC;AAClC,+BAA+B;AAC/B,MAAM,CAAC,MAAM,EAAE,GAAG,MAAM,IAAI,CAAC,CAAC;AAC9B,gBAAgB;AAChB,MAAM,CAAC,MAAM,aAAa,GAAG,CAAC,CAAS,EAAU,EAAE,CAAC,CAAC,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;AACtE,gBAAgB;AAChB,MAAM,CAAC,MAAM,WAAW,GAAG,CAAC,CAAS,EAAU,EAAE,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC;AAE1D;;;;;GAKG;AACH,MAAM,UAAU,OAAO,CAAC,CAAS,EAAE,CAAS,EAAE,CAAS,EAAE,KAA4B;IACnF,MAAM,GAAG,GAAG,YAAY,CAAC,CAAC,EAAE,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAE,EAAE,KAAK,CAAC,CAAC,CAAC,CAAC,CAAE,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;IAC5E,OAAO,CAAC,GAAG,CAAC,CAAC,CAAE,EAAE,GAAG,CAAC,CAAC,CAAE,CAAC,CAAC;AAC5B,CAAC;AAED;;;;;;;;;;;;GAYG;AACH,SAAS,cAAc,CAAC,CAAS,EAAE,CAAS,EAAE,EAAU,EAAE,CAAS,EAAE,KAA4B;IAC/F,IAAI,EAAE,KAAK,CAAC;QAAE,OAAO,CAAC,CAAC;IACvB,MAAM,QAAQ,GAAG,CAAC,GAAW,EAAU,EAAE,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,OAAO,CAAC,CAAC,EAAE,CAAC,EAAE,GAAG,EAAE,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IAC1G,OAAO,QAAQ,CAAC,QAAQ,EAAE,CAAC,EAAE,CAAC,EAAE,EAAE,GAAG,EAAE,KAAK,EAAE,OAAO,EAAE,GAAG,EAAE,CAAC,CAAC;AAChE,CAAC;AAED,kGAAkG;AAClG,MAAM,UAAU,aAAa,CAAC,KAA4B,EAAE,MAAc,EAAE,cAAsB,EAAE,OAAe;IACjH,MAAM,CAAC,GAAG,CAAC,GAAG,cAAc,CAAC;IAC7B,MAAM,GAAG,GAAa,EAAE,CAAC;IACzB,IAAI,CAAC,GAAG,MAAM,CAAC;IACf,IAAI,CAAC,GAAG,CAAC,CAAC;IACV,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,cAAc,GAAG,OAAO,EAAE,CAAC,EAAE,EAAE,CAAC;QAClD,MAAM,CAAC,EAAE,EAAE,EAAE,CAAC,GAAG,OAAO,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,KAAK,CAAC,CAAC;QACzC,IAAI,EAAE,KAAK,CAAC,IAAI,CAAC,GAAG,CAAC,KAAK,EAAE,GAAG,CAAC;YAAE,GAAG,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,GAAG,cAAc,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,EAAE,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC;QACvF,CAAC,GAAG,EAAE,CAAC;QACP,CAAC,GAAG,EAAE,CAAC;IACT,CAAC;IACD,OAAO,GAAG,CAAC;AACb,CAAC;AAED,+FAA+F;AAC/F,MAAM,UAAU,OAAO,CAAC,SAA4B,EAAE,CAAS;IAC7D,IAAI,KAAK,GAAG,CAAC,CAAC;IACd,IAAI,SAAS,GAAG,CAAC,CAAC;IAClB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,SAAS,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;QAC1C,MAAM,CAAC,GAAG,SAAS,CAAC,CAAC,CAAE,CAAC;QACxB,MAAM,KAAK,GAAG,IAAI,CAAC,EAAE,GAAG,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC;QACxC,IAAI,CAAC,IAAI,CAAC;YAAE,OAAO,SAAS,GAAG,CAAC,CAAC,KAAK,GAAG,SAAS,CAAC,GAAG,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC;QACjF,KAAK,GAAG,CAAC,CAAC;QACV,SAAS,GAAG,KAAK,CAAC;IACpB,CAAC;IACD,OAAO,MAAM,CAAC,GAAG,CAAC;AACpB,CAAC;AAED;;;;;;;;;;;;GAYG;AACH,MAAM,UAAU,mBAAmB,CAAC,GAAW;IAC7C,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,GAAG,CAAC,IAAI,GAAG,GAAG,CAAC;QAAE,OAAO,MAAM,CAAC,GAAG,CAAC;IACxD,IAAI,GAAG,KAAK,CAAC;QAAE,OAAO,CAAC,CAAC;IACxB,MAAM,EAAE,GAAG,IAAI,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,MAAM,CAAC,OAAO,CAAC,CAAC;IAC1C,IAAI,qBAAqB,CAAC,EAAE,MAAM,EAAE,EAAE,EAAE,CAAC,IAAI,GAAG;QAAE,OAAO,EAAE,CAAC;IAC5D,OAAO,QAAQ,CAAC,CAAC,MAAM,EAAE,EAAE,CAAC,qBAAqB,CAAC,EAAE,MAAM,EAAE,CAAC,GAAG,GAAG,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,EAAE;QAC9E,GAAG,EAAE,CAAC,GAAG,MAAM,CAAC,OAAO;QACvB,OAAO,EAAE,GAAG;KACb,CAAC,CAAC;AACL,CAAC;AAED;;;;;GAKG;AACH,MAAM,UAAU,YAAY,CAAC,MAAc;IACzC,MAAM,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC;IACzC,IAAI,CAAC,KAAK,CAAC;QAAE,OAAO,CAAC,CAAC;IACtB,MAAM,CAAC,GAAG,UAAU,CAAC,WAAW,EAAE,SAAS,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;IACpD,MAAM,EAAE,GAAG,UAAU,CAAC,WAAW,EAAE,SAAS,CAAC,IAAI,CAAC,GAAG,CAAC,MAAM,GAAG,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;IACzE,OAAO,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,IAAI,CAAC,EAAE,GAAG,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC;AACvC,CAAC;AAED;;;;;;;;;;;;GAYG;AACH,MAAM,UAAU,kBAAkB,CAAC,MAAc,EAAE,CAAS;IAC1D,MAAM,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC;IAC/B,IAAI,CAAC,KAAK,CAAC;QAAE,OAAO,CAAC,CAAC;IACtB,MAAM,CAAC,GAAG,CAAC,IAAI,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,qBAAqB,CAAC,EAAE,MAAM,EAAE,CAAC,CAAC,CAAC;IAClE,OAAO,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,GAAG,QAAQ,CAAC,CAAC,GAAG,MAAM,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;AAC5D,CAAC","sourcesContent":["/**\n * The pendulum and its linear oscillator, integrated and read off, for the\n * witnesses of `ab-pendulum-linear`'s declared translations. Time is in units\n * of T0 (ω0 = 2π); both motions are released from rest at θ0.\n *\n * Also the closed forms those translations are built from: θ0 recovered from\n * the bound's ε, and the nome q of the pendulum's elliptic modulus.\n *\n * @module atlas/oscillators/pendulum-motion\n * @internal\n */\n\nimport { ellipticK, findRoot, jacobiSN } from '@danielsimonjr/mathts-functions';\nimport { numberFrom } from '../../core/mathts-read.js';\nimport { integrateRK4 } from '../../numerical/null-ray-integrator.js';\nimport { pendulumPeriodErrorAt } from './bridges-limits.js';\n\n/** 2π: one period of the angle variable, the T0 = 1 time unit's ω0. @internal */\nexport const TWO_PI = 2 * Math.PI;\n/** ω0² at T0 = 1. @internal */\nexport const W2 = TWO_PI ** 2;\n/** @internal */\nexport const pendulumAccel = (x: number): number => -W2 * Math.sin(x);\n/** @internal */\nexport const linearAccel = (x: number): number => -W2 * x;\n\n/**\n * One fixed step of x″ = accel(x): `integrateRK4` (MathTS `solveODESystem`\n * with `dt`, the classical RK4 step) over one step of length `h`.\n *\n * @internal\n */\nexport function rk4Step(x: number, v: number, h: number, accel: (x: number) => number): [number, number] {\n  const end = integrateRK4((_t, y) => [y[1]!, accel(y[0]!)], [x, v], 0, h, 1);\n  return [end[0]!, end[1]!];\n}\n\n/**\n * The root in (0, h] of x(τ), the position one partial RK4 step of length τ\n * from (x, v), by MathTS `findRoot` (Brent) on the bracket (0, x), (h, nx).\n *\n * The crossing is then located to the integrator's own truncation error,\n * rather than to the O(h³) error of a straight line between the two steps.\n * That truncation error is nearly the same for the pendulum and the linear\n * oscillator at small θ0, so it cancels in their phase difference; the\n * interpolation error does not. `findRoot` evaluates both ends, and a step\n * of length 0 is the start itself (`solveODESystem` cannot take `dt = 0`).\n * Its `tol` stops on the residual (an angle) and on the bracket width (a\n * time in T0), both at 1e-12, below the integrator's error either way.\n */\nfunction refineCrossing(x: number, v: number, nx: number, h: number, accel: (x: number) => number): number {\n  if (nx === 0) return h;\n  const position = (tau: number): number => (tau === 0 ? x : tau === h ? nx : rk4Step(x, v, tau, accel)[0]);\n  return findRoot(position, 0, h, { tol: 1e-12, maxIter: 200 });\n}\n\n/** Zero-crossing times of a release from rest at `theta0`, over `periods` of T0 = 1. @internal */\nexport function crossingTimes(accel: (x: number) => number, theta0: number, stepsPerPeriod: number, periods: number): number[] {\n  const h = 1 / stepsPerPeriod;\n  const out: number[] = [];\n  let x = theta0;\n  let v = 0;\n  for (let i = 0; i < stepsPerPeriod * periods; i++) {\n    const [nx, nv] = rk4Step(x, v, h, accel);\n    if (nx === 0 || x > 0 !== nx > 0) out.push(i * h + refineCrossing(x, v, nx, h, accel));\n    x = nx;\n    v = nv;\n  }\n  return out;\n}\n\n/** The phase at `t`, interpolated through the anchors (0, 0) and (c_k, π/2 + kπ). @internal */\nexport function phaseAt(crossings: readonly number[], t: number): number {\n  let prevT = 0;\n  let prevPhase = 0;\n  for (let k = 0; k < crossings.length; k++) {\n    const c = crossings[k]!;\n    const phase = Math.PI / 2 + k * Math.PI;\n    if (t <= c) return prevPhase + ((phase - prevPhase) * (t - prevT)) / (c - prevT);\n    prevT = c;\n    prevPhase = phase;\n  }\n  return Number.NaN;\n}\n\n/**\n * θ0 in [0, π) whose exact relative period error is `eps`, by MathTS\n * `findRoot` (Brent) on the strictly increasing `pendulumPeriodErrorAt`.\n *\n * The residual is solved as `err/eps − 1`, of order 1, so `findRoot`'s one\n * `tol` (it stops on the residual AND on the bracket width) cannot stop on a\n * residual that is small only because `eps` is. The upper end is just below\n * π, where the error is finite; at π it is `Infinity`, which Brent's\n * interpolation cannot take.\n *\n * @returns NaN for a negative or non-finite `eps`.\n * @internal\n */\nexport function theta0OfPeriodError(eps: number): number {\n  if (!Number.isFinite(eps) || eps < 0) return Number.NaN;\n  if (eps === 0) return 0;\n  const hi = Math.PI * (1 - Number.EPSILON);\n  if (pendulumPeriodErrorAt({ theta0: hi }) <= eps) return hi;\n  return findRoot((theta0) => pendulumPeriodErrorAt({ theta0 }) / eps - 1, 0, hi, {\n    tol: 4 * Number.EPSILON,\n    maxIter: 200,\n  });\n}\n\n/**\n * The nome q = exp(−π K′/K) of the modulus k = sin(θ0/2): K = K(m) and\n * K′ = K(1 − m) at the parameter m = k², both MathTS `ellipticK` (whose\n * route is the AGM, so K′/K = AGM(1, cos(θ0/2)) / AGM(1, sin(θ0/2))).\n * @internal\n */\nexport function pendulumNome(theta0: number): number {\n  const k = Math.sin(Math.abs(theta0) / 2);\n  if (k === 0) return 0;\n  const K = numberFrom('ellipticK', ellipticK(k * k));\n  const Kp = numberFrom('ellipticK', ellipticK(Math.cos(theta0 / 2) ** 2));\n  return Math.exp((-Math.PI * Kp) / K);\n}\n\n/**\n * The exact pendulum released from rest at `theta0`, in units of T0 = 1\n * (ω0 = 2π): θ(t) = 2 asin(k sn(K − ω0 t, m)), with k = sin(θ0/2), the\n * parameter m = k², and K = K(m) the quarter period in the sn argument, read\n * from `pendulumPeriodErrorAt` as (π/2)(1 + ε) so one `ellipticK` read\n * serves both. At t = 0, sn K = 1 and θ = θ0; the period is T = 4K/ω0 =\n * (1 + ε) T0. MathTS `jacobiSN` is a plain number function of (u, m). A\n * closed form, independent of every integrator: `tests/atlas/pendulum-exact.test.ts`\n * checks it against the cn·dn identities (θ″ = −2ω0² k sn u dn u, sin θ =\n * 2k sn u dn u) and against RK4.\n *\n * @internal\n */\nexport function exactPendulumAngle(theta0: number, t: number): number {\n  const k = Math.sin(theta0 / 2);\n  if (k === 0) return 0;\n  const K = (Math.PI / 2) * (1 + pendulumPeriodErrorAt({ theta0 }));\n  return 2 * Math.asin(k * jacobiSN(K - TWO_PI * t, k * k));\n}\n"]}