universal-physics-tensor 10.0.0 → 12.0.0

package/dist/atlas/oscillators/phase-translation.d.ts

@@ -28,17 +28,6 @@
  * @internal
  */
 import type { ObservableTranslation } from '../translation.js';
-/**
- * The time, in units of T0, at which the MEASURED phase lead of the linear
- * oscillator over the pendulum first reaches `tolerance` rad, both integrated
- * by RK4 at `stepsPerPeriod` steps of T0.
- *
- * @returns NaN if the drift does not reach `tolerance` within 2^16 periods.
- * @internal
- */
-export declare function measurePhaseHorizon(theta0: number, tolerance: number, stepsPerPeriod: number): number;
-/** As `measurePhaseHorizon`, over exactly `periods` of T0; NaN if the drift does not reach `tolerance` in them. @internal */
-export declare function measurePhaseHorizonWithin(theta0: number, tolerance: number, stepsPerPeriod: number, periods: number): number;
 /** The point witness integrates at most this many periods of T0, so its cost is bounded. */
 export declare const PHASE_POINT_MAX_PERIODS = 512;
 /** @internal */

package/dist/atlas/oscillators/phase-translation.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"phase-translation.d.ts","sourceRoot":"","sources":["../../../src/atlas/oscillators/phase-translation.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;;;;;;;;;;GA4BG;AAEH,OAAO,KAAK,EAAE,qBAAqB,EAAE,MAAM,mBAAmB,CAAC;AAM/D;;;;;;;GAOG;AACH,wBAAgB,mBAAmB,CAAC,MAAM,EAAE,MAAM,EAAE,SAAS,EAAE,MAAM,EAAE,cAAc,EAAE,MAAM,GAAG,MAAM,CAMrG;AAED,6HAA6H;AAC7H,wBAAgB,yBAAyB,CAAC,MAAM,EAAE,MAAM,EAAE,SAAS,EAAE,MAAM,EAAE,cAAc,EAAE,MAAM,EAAE,OAAO,EAAE,MAAM,GAAG,MAAM,CAc5H;AAWD,4FAA4F;AAC5F,eAAO,MAAM,uBAAuB,MAAM,CAAC;AAyD3C,gBAAgB;AAChB,eAAO,MAAM,0BAA0B,EAAE,qBA+CxC,CAAC"}
+{"version":3,"file":"phase-translation.d.ts","sourceRoot":"","sources":["../../../src/atlas/oscillators/phase-translation.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;;;;;;;;;;GA4BG;AAEH,OAAO,KAAK,EAAE,qBAAqB,EAAE,MAAM,mBAAmB,CAAC;AA8C/D,4FAA4F;AAC5F,eAAO,MAAM,uBAAuB,MAAM,CAAC;AAyD3C,gBAAgB;AAChB,eAAO,MAAM,0BAA0B,EAAE,qBA+CxC,CAAC"}

package/dist/atlas/oscillators/phase-translation.js

@@ -28,8 +28,7 @@
  * @internal
  */
 import { pendulumPeriodErrorAt } from './bridges-limits.js';
-import { crossingTimes, linearAccel, pendulumAccel, phaseAt, theta0OfPeriodError } from './pendulum-motion.js';
-const TWO_PI = 2 * Math.PI;
+import { crossingTimes, linearAccel, pendulumAccel, phaseAt, theta0OfPeriodError, TWO_PI } from './pendulum-motion.js';
 /**
  * The time, in units of T0, at which the MEASURED phase lead of the linear
  * oscillator over the pendulum first reaches `tolerance` rad, both integrated
@@ -38,7 +37,7 @@ const TWO_PI = 2 * Math.PI;
  * @returns NaN if the drift does not reach `tolerance` within 2^16 periods.
  * @internal
  */
-export function measurePhaseHorizon(theta0, tolerance, stepsPerPeriod) {
+function measurePhaseHorizon(theta0, tolerance, stepsPerPeriod) {
     for (let periods = 16; periods <= 2 ** 16; periods *= 2) {
         const t = measurePhaseHorizonWithin(theta0, tolerance, stepsPerPeriod, periods);
         if (!Number.isNaN(t))
@@ -47,7 +46,7 @@ export function measurePhaseHorizon(theta0, tolerance, stepsPerPeriod) {
     return Number.NaN;
 }
 /** As `measurePhaseHorizon`, over exactly `periods` of T0; NaN if the drift does not reach `tolerance` in them. @internal */
-export function measurePhaseHorizonWithin(theta0, tolerance, stepsPerPeriod, periods) {
+function measurePhaseHorizonWithin(theta0, tolerance, stepsPerPeriod, periods) {
     const pend = crossingTimes(pendulumAccel, theta0, stepsPerPeriod, periods);
     const lin = crossingTimes(linearAccel, theta0, stepsPerPeriod, periods);
     let prevT = 0;

package/dist/atlas/oscillators/phase-translation.js.map

@@ -1 +1 @@
-{"version":3,"file":"phase-translation.js","sourceRoot":"","sources":["../../../src/atlas/oscillators/phase-translation.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;;;;;;;;;;GA4BG;AAGH,OAAO,EAAE,qBAAqB,EAAE,MAAM,qBAAqB,CAAC;AAC5D,OAAO,EAAE,aAAa,EAAE,WAAW,EAAE,aAAa,EAAE,OAAO,EAAE,mBAAmB,EAAE,MAAM,sBAAsB,CAAC;AAE/G,MAAM,MAAM,GAAG,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC;AAE3B;;;;;;;GAOG;AACH,MAAM,UAAU,mBAAmB,CAAC,MAAc,EAAE,SAAiB,EAAE,cAAsB;IAC3F,KAAK,IAAI,OAAO,GAAG,EAAE,EAAE,OAAO,IAAI,CAAC,IAAI,EAAE,EAAE,OAAO,IAAI,CAAC,EAAE,CAAC;QACxD,MAAM,CAAC,GAAG,yBAAyB,CAAC,MAAM,EAAE,SAAS,EAAE,cAAc,EAAE,OAAO,CAAC,CAAC;QAChF,IAAI,CAAC,MAAM,CAAC,KAAK,CAAC,CAAC,CAAC;YAAE,OAAO,CAAC,CAAC;IACjC,CAAC;IACD,OAAO,MAAM,CAAC,GAAG,CAAC;AACpB,CAAC;AAED,6HAA6H;AAC7H,MAAM,UAAU,yBAAyB,CAAC,MAAc,EAAE,SAAiB,EAAE,cAAsB,EAAE,OAAe;IAClH,MAAM,IAAI,GAAG,aAAa,CAAC,aAAa,EAAE,MAAM,EAAE,cAAc,EAAE,OAAO,CAAC,CAAC;IAC3E,MAAM,GAAG,GAAG,aAAa,CAAC,WAAW,EAAE,MAAM,EAAE,cAAc,EAAE,OAAO,CAAC,CAAC;IACxE,IAAI,KAAK,GAAG,CAAC,CAAC;IACd,IAAI,SAAS,GAAG,CAAC,CAAC;IAClB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;QACrC,MAAM,CAAC,GAAG,IAAI,CAAC,CAAC,CAAE,CAAC;QACnB,MAAM,KAAK,GAAG,OAAO,CAAC,GAAG,EAAE,CAAC,CAAC,GAAG,CAAC,IAAI,CAAC,EAAE,GAAG,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC,CAAC;QAC5D,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,KAAK,CAAC;YAAE,MAAM;QACnC,IAAI,KAAK,IAAI,SAAS;YAAE,OAAO,KAAK,GAAG,CAAC,CAAC,SAAS,GAAG,SAAS,CAAC,GAAG,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,GAAG,CAAC,KAAK,GAAG,SAAS,CAAC,CAAC;QACrG,KAAK,GAAG,CAAC,CAAC;QACV,SAAS,GAAG,KAAK,CAAC;IACpB,CAAC;IACD,OAAO,MAAM,CAAC,GAAG,CAAC;AACpB,CAAC;AAED,iFAAiF;AACjF,SAAS,YAAY,CAAC,GAAW,EAAE,SAAiB,EAAE,EAAsB;IAC1E,IAAI,EAAE,KAAK,SAAS,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,EAAE,CAAC,IAAI,EAAE,IAAI,CAAC,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,GAAG,CAAC,IAAI,GAAG,GAAG,CAAC;QAAE,OAAO,MAAM,CAAC,GAAG,CAAC;IAC/G,OAAO,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,SAAS,GAAG,EAAE,GAAG,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,GAAG,CAAC,MAAM,GAAG,GAAG,CAAC,CAAC;AAC9E,CAAC;AAED,MAAM,UAAU,GAAG,GAAG,CAAC;AACvB,MAAM,aAAa,GAAG,CAAC,CAAC;AAExB,4FAA4F;AAC5F,MAAM,CAAC,MAAM,uBAAuB,GAAG,GAAG,CAAC;AAC3C;;;;GAIG;AACH,MAAM,iBAAiB,GAAG,CAAC,CAAC;AAC5B;;;;;;;;GAQG;AACH,MAAM,qBAAqB,GAAG,IAAI,CAAC;AACnC;;;GAGG;AACH,MAAM,8BAA8B,GAAG,IAAI,CAAC;AAE5C;;;;;;;GAOG;AACH,SAAS,eAAe,CAAC,GAAW;IAClC,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,GAAG,CAAC,IAAI,GAAG,IAAI,CAAC;QAAE,OAAO,EAAE,WAAW,EAAE,qCAAqC,GAAG,EAAE,EAAE,CAAC;IAC1G,MAAM,MAAM,GAAG,mBAAmB,CAAC,GAAG,CAAC,CAAC;IACxC,MAAM,QAAQ,GAAG,IAAI,CAAC,GAAG,CAAC,iBAAiB,EAAE,CAAC,MAAM,GAAG,uBAAuB,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC;IACnG,IAAI,QAAQ,GAAG,qBAAqB,EAAE,CAAC;QACrC,OAAO;YACL,WAAW,EACT,oBAAoB,uBAAuB,eAAe,MAAM,CAAC,MAAM,CAAC,WAAW,CAAC,CAAC,CAAC,CAAC,OAAO,QAAQ,kBAAkB;gBACxH,GAAG,qBAAqB,qCAAqC;SAChE,CAAC;IACJ,CAAC;IACD,MAAM,MAAM,GAAG,YAAY,CAAC,GAAG,EAAE,QAAQ,EAAE,CAAC,CAAC,CAAC;IAC9C,MAAM,OAAO,GAAG,IAAI,CAAC,IAAI,CAAC,MAAM,CAAC,GAAG,CAAC,CAAC;IACtC,MAAM,QAAQ,GAAG,CAAC,cAAsB,EAAU,EAAE,CAAC,yBAAyB,CAAC,MAAM,EAAE,QAAQ,EAAE,cAAc,EAAE,OAAO,CAAC,CAAC;IAC1H,MAAM,IAAI,GAAG,EAAE,QAAQ,EAAE,gBAAgB,EAAE,GAAG,EAAE,cAAc,EAAE,GAAG,EAAE,SAAS,EAAE,8BAA8B,GAAG,MAAM,EAAE,CAAC;IAC1H,MAAM,MAAM,GAAG,QAAQ,GAAG,iBAAiB,CAAC,CAAC,CAAC,0BAA0B,uBAAuB,MAAM,CAAC,CAAC,CAAC,EAAE,CAAC;IAC3G,OAAO;QACL,KAAK,EACH,wBAAwB,QAAQ,OAAO,MAAM,2BAA2B,MAAM,CAAC,MAAM,CAAC,WAAW,CAAC,CAAC,CAAC,CAAC,sBAAsB;YAC3H,UAAU,8BAA8B,qBAAqB,MAAM,kCAAkC,OAAO,gBAAgB;QAC9H,KAAK,EAAE,EAAE,EAAE,EAAE,WAAW,EAAE,GAAG,IAAI,EAAE,MAAM,EAAE;QAC3C,OAAO,EAAE,EAAE,EAAE,EAAE,mBAAmB,EAAE,GAAG,IAAI,EAAE,MAAM,EAAE,QAAQ,GAAG,CAAC,MAAM,GAAG,GAAG,CAAC,EAAE;QAChF,YAAY,EAAE,gFAAgF;KAC/F,CAAC;AACJ,CAAC;AAED,gBAAgB;AAChB,MAAM,CAAC,MAAM,0BAA0B,GAA0B;IAC/D,QAAQ,EAAE,oBAAoB;IAC9B,UAAU,EAAE,OAAO;IACnB,IAAI,EAAE,KAAK;IACX,QAAQ,EAAE,qEAAqE;IAC/E,SAAS,EAAE,OAAO;IAClB,UAAU,EACR,oGAAoG;QACpG,oEAAoE;IACtE,UAAU,EACR,oGAAoG;QACpG,kGAAkG;IACpG,QAAQ,EAAE;QACR,6DAA6D;QAC7D,gEAAgE;KACjE;IACD,UAAU,EAAE,CAAC,IAAI,CAAC;IAClB,QAAQ,EAAE,gBAAgB;IAC1B,OAAO,EAAE,CAAC,GAAG,EAAE,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE;QAC1B,IAAI,EAAE,KAAK,SAAS,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,EAAE,CAAC,IAAI,EAAE,IAAI,CAAC,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,GAAG,CAAC,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC;YAAE,OAAO,MAAM,CAAC,GAAG,CAAC;QAC3H,OAAO,CAAC,MAAM,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC;IAC/C,CAAC;IACD,UAAU,EAAE,CAAC,GAAG,EAAE,SAAS,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,YAAY,CAAC,GAAG,EAAE,SAAS,EAAE,EAAE,CAAC;IACxE,QAAQ,EAAE,4DAA4D;IACtE,UAAU,EAAE;QACV,sHAAsH;QACtH,eAAe;KAChB;IACD,SAAS,EAAE;QACT;YACE,EAAE,EAAE,KAAK;YACT,IAAI,EAAE,SAAS;YACf,IAAI,EAAE,gDAAgD;YACtD,SAAS,EAAE,wBAAwB,aAAa,+BAA+B,UAAU,uBAAuB;SACjH;KACF;IACD,MAAM,EAAE;QACN;YACE,EAAE,EAAE,KAAK;YACT,QAAQ,EAAE,CAAC,cAAc,EAAE,EAAE,CAAC,mBAAmB,CAAC,UAAU,EAAE,aAAa,EAAE,cAAc,CAAC;YAC5F,MAAM,EAAE,YAAY,CAAC,qBAAqB,CAAC,EAAE,MAAM,EAAE,UAAU,EAAE,CAAC,EAAE,aAAa,EAAE,CAAC,CAAC;YACrF,gBAAgB,EAAE,EAAE;YACpB,cAAc,EAAE,GAAG;YACnB,SAAS,EAAE,IAAI;SAChB;KACF;IACD,UAAU,EAAE,eAAe;CAC5B,CAAC","sourcesContent":["/**\n * `ab-pendulum-linear`'s declared translation from its bound's quantity, the\n * relative period error, into the accumulated PHASE error.\n *\n * Both motions start from rest at θ0 at t = 0. The phase is the angle\n * variable: 0 at release and advancing by 2π per period, uniformly in time,\n * so the pendulum's is 2π t/T and the linear oscillator's 2π t/T0. With the\n * bound's `ε = T/T0 − 1` (normalised by the reduced model's period),\n *\n *   Δφ(t) = 2π t/T0 − 2π t/T = 2π (t/T0) · ε/(1 + ε),\n *\n * which reaches a tolerance Δ at `t* = Δ · T0 · (1 + ε)/(2π ε)`. Nothing is\n * truncated: the map is exact in ε, and ε itself is the bound's closed form.\n *\n * It is NOT a bound on the pointwise position error |θ(t) − θ_lin(t)|: the\n * pendulum's waveform carries a third harmonic of amplitude ≈ θ0³/192, which\n * is not a phase, and the position error at a given phase error depends on\n * where in the cycle the observation falls.\n *\n * The witness W7p integrates BOTH equations of motion and reads each phase off\n * its zero crossings — crossing k is at phase π/2 + kπ by the symmetry of a\n * release from rest, and the angle variable is linear in time between them —\n * so it measures the drift without using ε or the elliptic integral. It runs\n * at the fixture θ0 = 0.2 as the regression anchor, and again at the caller's\n * point, bounded to 512 T0 of integration (`phasePointCheck`).\n *\n * @module atlas/oscillators/phase-translation\n * @internal\n */\n\nimport type { ObservableTranslation } from '../translation.js';\nimport { pendulumPeriodErrorAt } from './bridges-limits.js';\nimport { crossingTimes, linearAccel, pendulumAccel, phaseAt, theta0OfPeriodError } from './pendulum-motion.js';\n\nconst TWO_PI = 2 * Math.PI;\n\n/**\n * The time, in units of T0, at which the MEASURED phase lead of the linear\n * oscillator over the pendulum first reaches `tolerance` rad, both integrated\n * by RK4 at `stepsPerPeriod` steps of T0.\n *\n * @returns NaN if the drift does not reach `tolerance` within 2^16 periods.\n * @internal\n */\nexport function measurePhaseHorizon(theta0: number, tolerance: number, stepsPerPeriod: number): number {\n  for (let periods = 16; periods <= 2 ** 16; periods *= 2) {\n    const t = measurePhaseHorizonWithin(theta0, tolerance, stepsPerPeriod, periods);\n    if (!Number.isNaN(t)) return t;\n  }\n  return Number.NaN;\n}\n\n/** As `measurePhaseHorizon`, over exactly `periods` of T0; NaN if the drift does not reach `tolerance` in them. @internal */\nexport function measurePhaseHorizonWithin(theta0: number, tolerance: number, stepsPerPeriod: number, periods: number): number {\n  const pend = crossingTimes(pendulumAccel, theta0, stepsPerPeriod, periods);\n  const lin = crossingTimes(linearAccel, theta0, stepsPerPeriod, periods);\n  let prevT = 0;\n  let prevDrift = 0;\n  for (let k = 0; k < pend.length; k++) {\n    const c = pend[k]!;\n    const drift = phaseAt(lin, c) - (Math.PI / 2 + k * Math.PI);\n    if (!Number.isFinite(drift)) break;\n    if (drift >= tolerance) return prevT + ((tolerance - prevDrift) * (c - prevT)) / (drift - prevDrift);\n    prevT = c;\n    prevDrift = drift;\n  }\n  return Number.NaN;\n}\n\n/** `t* = Δ T0 (1+ε)/(2π ε)`; `Infinity` at ε = 0, where no phase accumulates. */\nfunction phaseHorizon(eps: number, tolerance: number, T0: number | undefined): number {\n  if (T0 === undefined || !Number.isFinite(T0) || T0 <= 0 || !Number.isFinite(eps) || eps < 0) return Number.NaN;\n  return eps === 0 ? Infinity : (tolerance * T0 * (1 + eps)) / (TWO_PI * eps);\n}\n\nconst W7P_THETA0 = 0.2;\nconst W7P_TOLERANCE = 1;\n\n/** The point witness integrates at most this many periods of T0, so its cost is bounded. */\nexport const PHASE_POINT_MAX_PERIODS = 512;\n/**\n * The point witness measures the time to this much drift, or to the drift\n * reached at the cap if less: the map is linear in t, so its slope is what is\n * checked, whatever tolerance the caller asked about.\n */\nconst PHASE_POINT_DRIFT = 1;\n/**\n * Below this measured drift the rounding noise of the two motions' crossing\n * times comes within a factor of ten of the witness tolerance. With the\n * crossings located by `crossingTimes`' secant refinement, the fine error was\n * measured at 800 steps per T0 as 1.5e-8 relative at 1.1e-5 rad of drift,\n * 7.4e-9 at 7.2e-5 rad, 2.5e-9 at 2.0e-4 rad and 2.4e-10 to 9.2e-10 from\n * 2e-2 rad up. With the straight-line crossings it replaced, the same error\n * was 9.5e-6 relative at 2e-4 rad and 1.2e-7 at 2e-2 rad.\n */\nconst PHASE_POINT_MIN_DRIFT = 1e-4;\n/**\n * The (1+ε) control differs from the target by ε/(1+ε) relative, so it can be\n * refuted only where ε exceeds this tolerance (θ0 above about 9e-4).\n */\nconst PHASE_POINT_RELATIVE_TOLERANCE = 5e-8;\n\n/**\n * W7p at the caller's θ0 (recovered from ε) instead of the fixture: the time\n * the measured drift reaches min(1 rad, the drift at the cap), within\n * PHASE_POINT_RELATIVE_TOLERANCE of the closed form. Its control is the same measurement against\n * the horizon without the (1+ε) factor, which differs from the target by\n * ε/(1+ε) relative — so the control can fail only where that exceeds the\n * witness tolerance, and the caller is told when it does not.\n */\nfunction phasePointCheck(eps: number): ReturnType<ObservableTranslation['pointCheck']> {\n  if (!Number.isFinite(eps) || eps <= 0) return { unavailable: `no phase drift accumulates at ε = ${eps}` };\n  const theta0 = theta0OfPeriodError(eps);\n  const measured = Math.min(PHASE_POINT_DRIFT, (TWO_PI * PHASE_POINT_MAX_PERIODS * eps) / (1 + eps));\n  if (measured < PHASE_POINT_MIN_DRIFT) {\n    return {\n      unavailable:\n        `the drift within ${PHASE_POINT_MAX_PERIODS} T0 at θ0 = ${Number(theta0.toPrecision(8))} is ${measured} rad, below the ` +\n        `${PHASE_POINT_MIN_DRIFT} rad this witness resolves; not run`,\n    };\n  }\n  const target = phaseHorizon(eps, measured, 1);\n  const periods = Math.ceil(target) + 2;\n  const evaluate = (stepsPerPeriod: number): number => measurePhaseHorizonWithin(theta0, measured, stepsPerPeriod, periods);\n  const spec = { evaluate, coarseResolution: 200, fineResolution: 800, tolerance: PHASE_POINT_RELATIVE_TOLERANCE * target };\n  const capped = measured < PHASE_POINT_DRIFT ? ` (the drift reached in ${PHASE_POINT_MAX_PERIODS} T0)` : '';\n  return {\n    claim:\n      `RK4-measured time to ${measured} rad${capped} of phase drift at θ0 = ${Number(theta0.toPrecision(8))} (recovered from ε) ` +\n      `within ${PHASE_POINT_RELATIVE_TOLERANCE} relative of t* = ${target} T0; 200 and 800 steps per T0, ${periods} T0 integrated`,\n    check: { id: 'W7p@point', ...spec, target },\n    control: { id: 'W7p@point-control', ...spec, target: measured / (TWO_PI * eps) },\n    controlClaim: 'the same measurement against the horizon without the (1+ε) factor, Δ T0/(2π ε)',\n  };\n}\n\n/** @internal */\nexport const PENDULUM_PHASE_TRANSLATION: ObservableTranslation = {\n  bridgeId: 'ab-pendulum-linear',\n  observable: 'phase',\n  unit: 'rad',\n  fromNorm: 'relative period error, normalized by the value of the reduced model',\n  errorKind: 'exact',\n  definition:\n    'the angle variable of each motion (0 at release from rest at θ0, 2π per period); the error is the ' +\n    \"linear oscillator's lead over the pendulum, Δφ(t) = 2πt/T0 − 2πt/T\",\n  derivation:\n    'with ε = T/T0 − 1, Δφ(t) = 2π (t/T0) ε/(1+ε), exact in ε; it reaches Δ at t* = Δ T0 (1+ε)/(2π ε). ' +\n    'ε/(1+ε) increases with ε, so the domain supremum of ε gives the shortest horizon over the domain',\n  premises: [\n    'both motions are released from rest at the same θ0 at t = 0',\n    'T0 = 2π √(ℓ/g) is the linear period, in the unit t is given in',\n  ],\n  parameters: ['T0'],\n  timeUnit: 'the unit of T0',\n  errorAt: (eps, t, { T0 }) => {\n    if (T0 === undefined || !Number.isFinite(T0) || T0 <= 0 || !Number.isFinite(eps) || !Number.isFinite(t)) return Number.NaN;\n    return (TWO_PI * (t / T0) * eps) / (1 + eps);\n  },\n  horizonFor: (eps, tolerance, { T0 }) => phaseHorizon(eps, tolerance, T0),\n  boundary: 'inclusive: adequate iff t ≤ t*, i.e. Δφ(t) ≤ the tolerance',\n  notCovered: [\n    'the pointwise position error |θ(t) − θ_lin(t)| (the waveform carries a third harmonic ≈ θ0³/192 that is not a phase)',\n    'the amplitude',\n  ],\n  witnesses: [\n    {\n      id: 'W7p',\n      kind: 'numeric',\n      test: 'tests/atlas/pendulum-phase-translation.test.ts',\n      tolerance: `RK4-measured time to ${W7P_TOLERANCE} rad of phase drift at θ0 = ${W7P_THETA0} within 1e-4 T0 of t*`,\n    },\n  ],\n  checks: [\n    {\n      id: 'W7p',\n      evaluate: (stepsPerPeriod) => measurePhaseHorizon(W7P_THETA0, W7P_TOLERANCE, stepsPerPeriod),\n      target: phaseHorizon(pendulumPeriodErrorAt({ theta0: W7P_THETA0 }), W7P_TOLERANCE, 1),\n      coarseResolution: 50,\n      fineResolution: 200,\n      tolerance: 1e-4,\n    },\n  ],\n  pointCheck: phasePointCheck,\n};\n"]}
+{"version":3,"file":"phase-translation.js","sourceRoot":"","sources":["../../../src/atlas/oscillators/phase-translation.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;;;;;;;;;;GA4BG;AAGH,OAAO,EAAE,qBAAqB,EAAE,MAAM,qBAAqB,CAAC;AAC5D,OAAO,EAAE,aAAa,EAAE,WAAW,EAAE,aAAa,EAAE,OAAO,EAAE,mBAAmB,EAAE,MAAM,EAAE,MAAM,sBAAsB,CAAC;AAEvH;;;;;;;GAOG;AACH,SAAS,mBAAmB,CAAC,MAAc,EAAE,SAAiB,EAAE,cAAsB;IACpF,KAAK,IAAI,OAAO,GAAG,EAAE,EAAE,OAAO,IAAI,CAAC,IAAI,EAAE,EAAE,OAAO,IAAI,CAAC,EAAE,CAAC;QACxD,MAAM,CAAC,GAAG,yBAAyB,CAAC,MAAM,EAAE,SAAS,EAAE,cAAc,EAAE,OAAO,CAAC,CAAC;QAChF,IAAI,CAAC,MAAM,CAAC,KAAK,CAAC,CAAC,CAAC;YAAE,OAAO,CAAC,CAAC;IACjC,CAAC;IACD,OAAO,MAAM,CAAC,GAAG,CAAC;AACpB,CAAC;AAED,6HAA6H;AAC7H,SAAS,yBAAyB,CAAC,MAAc,EAAE,SAAiB,EAAE,cAAsB,EAAE,OAAe;IAC3G,MAAM,IAAI,GAAG,aAAa,CAAC,aAAa,EAAE,MAAM,EAAE,cAAc,EAAE,OAAO,CAAC,CAAC;IAC3E,MAAM,GAAG,GAAG,aAAa,CAAC,WAAW,EAAE,MAAM,EAAE,cAAc,EAAE,OAAO,CAAC,CAAC;IACxE,IAAI,KAAK,GAAG,CAAC,CAAC;IACd,IAAI,SAAS,GAAG,CAAC,CAAC;IAClB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;QACrC,MAAM,CAAC,GAAG,IAAI,CAAC,CAAC,CAAE,CAAC;QACnB,MAAM,KAAK,GAAG,OAAO,CAAC,GAAG,EAAE,CAAC,CAAC,GAAG,CAAC,IAAI,CAAC,EAAE,GAAG,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC,CAAC;QAC5D,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,KAAK,CAAC;YAAE,MAAM;QACnC,IAAI,KAAK,IAAI,SAAS;YAAE,OAAO,KAAK,GAAG,CAAC,CAAC,SAAS,GAAG,SAAS,CAAC,GAAG,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,GAAG,CAAC,KAAK,GAAG,SAAS,CAAC,CAAC;QACrG,KAAK,GAAG,CAAC,CAAC;QACV,SAAS,GAAG,KAAK,CAAC;IACpB,CAAC;IACD,OAAO,MAAM,CAAC,GAAG,CAAC;AACpB,CAAC;AAED,iFAAiF;AACjF,SAAS,YAAY,CAAC,GAAW,EAAE,SAAiB,EAAE,EAAsB;IAC1E,IAAI,EAAE,KAAK,SAAS,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,EAAE,CAAC,IAAI,EAAE,IAAI,CAAC,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,GAAG,CAAC,IAAI,GAAG,GAAG,CAAC;QAAE,OAAO,MAAM,CAAC,GAAG,CAAC;IAC/G,OAAO,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,SAAS,GAAG,EAAE,GAAG,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,GAAG,CAAC,MAAM,GAAG,GAAG,CAAC,CAAC;AAC9E,CAAC;AAED,MAAM,UAAU,GAAG,GAAG,CAAC;AACvB,MAAM,aAAa,GAAG,CAAC,CAAC;AAExB,4FAA4F;AAC5F,MAAM,CAAC,MAAM,uBAAuB,GAAG,GAAG,CAAC;AAC3C;;;;GAIG;AACH,MAAM,iBAAiB,GAAG,CAAC,CAAC;AAC5B;;;;;;;;GAQG;AACH,MAAM,qBAAqB,GAAG,IAAI,CAAC;AACnC;;;GAGG;AACH,MAAM,8BAA8B,GAAG,IAAI,CAAC;AAE5C;;;;;;;GAOG;AACH,SAAS,eAAe,CAAC,GAAW;IAClC,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,GAAG,CAAC,IAAI,GAAG,IAAI,CAAC;QAAE,OAAO,EAAE,WAAW,EAAE,qCAAqC,GAAG,EAAE,EAAE,CAAC;IAC1G,MAAM,MAAM,GAAG,mBAAmB,CAAC,GAAG,CAAC,CAAC;IACxC,MAAM,QAAQ,GAAG,IAAI,CAAC,GAAG,CAAC,iBAAiB,EAAE,CAAC,MAAM,GAAG,uBAAuB,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC;IACnG,IAAI,QAAQ,GAAG,qBAAqB,EAAE,CAAC;QACrC,OAAO;YACL,WAAW,EACT,oBAAoB,uBAAuB,eAAe,MAAM,CAAC,MAAM,CAAC,WAAW,CAAC,CAAC,CAAC,CAAC,OAAO,QAAQ,kBAAkB;gBACxH,GAAG,qBAAqB,qCAAqC;SAChE,CAAC;IACJ,CAAC;IACD,MAAM,MAAM,GAAG,YAAY,CAAC,GAAG,EAAE,QAAQ,EAAE,CAAC,CAAC,CAAC;IAC9C,MAAM,OAAO,GAAG,IAAI,CAAC,IAAI,CAAC,MAAM,CAAC,GAAG,CAAC,CAAC;IACtC,MAAM,QAAQ,GAAG,CAAC,cAAsB,EAAU,EAAE,CAAC,yBAAyB,CAAC,MAAM,EAAE,QAAQ,EAAE,cAAc,EAAE,OAAO,CAAC,CAAC;IAC1H,MAAM,IAAI,GAAG,EAAE,QAAQ,EAAE,gBAAgB,EAAE,GAAG,EAAE,cAAc,EAAE,GAAG,EAAE,SAAS,EAAE,8BAA8B,GAAG,MAAM,EAAE,CAAC;IAC1H,MAAM,MAAM,GAAG,QAAQ,GAAG,iBAAiB,CAAC,CAAC,CAAC,0BAA0B,uBAAuB,MAAM,CAAC,CAAC,CAAC,EAAE,CAAC;IAC3G,OAAO;QACL,KAAK,EACH,wBAAwB,QAAQ,OAAO,MAAM,2BAA2B,MAAM,CAAC,MAAM,CAAC,WAAW,CAAC,CAAC,CAAC,CAAC,sBAAsB;YAC3H,UAAU,8BAA8B,qBAAqB,MAAM,kCAAkC,OAAO,gBAAgB;QAC9H,KAAK,EAAE,EAAE,EAAE,EAAE,WAAW,EAAE,GAAG,IAAI,EAAE,MAAM,EAAE;QAC3C,OAAO,EAAE,EAAE,EAAE,EAAE,mBAAmB,EAAE,GAAG,IAAI,EAAE,MAAM,EAAE,QAAQ,GAAG,CAAC,MAAM,GAAG,GAAG,CAAC,EAAE;QAChF,YAAY,EAAE,gFAAgF;KAC/F,CAAC;AACJ,CAAC;AAED,gBAAgB;AAChB,MAAM,CAAC,MAAM,0BAA0B,GAA0B;IAC/D,QAAQ,EAAE,oBAAoB;IAC9B,UAAU,EAAE,OAAO;IACnB,IAAI,EAAE,KAAK;IACX,QAAQ,EAAE,qEAAqE;IAC/E,SAAS,EAAE,OAAO;IAClB,UAAU,EACR,oGAAoG;QACpG,oEAAoE;IACtE,UAAU,EACR,oGAAoG;QACpG,kGAAkG;IACpG,QAAQ,EAAE;QACR,6DAA6D;QAC7D,gEAAgE;KACjE;IACD,UAAU,EAAE,CAAC,IAAI,CAAC;IAClB,QAAQ,EAAE,gBAAgB;IAC1B,OAAO,EAAE,CAAC,GAAG,EAAE,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE;QAC1B,IAAI,EAAE,KAAK,SAAS,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,EAAE,CAAC,IAAI,EAAE,IAAI,CAAC,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,GAAG,CAAC,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC;YAAE,OAAO,MAAM,CAAC,GAAG,CAAC;QAC3H,OAAO,CAAC,MAAM,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC;IAC/C,CAAC;IACD,UAAU,EAAE,CAAC,GAAG,EAAE,SAAS,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,YAAY,CAAC,GAAG,EAAE,SAAS,EAAE,EAAE,CAAC;IACxE,QAAQ,EAAE,4DAA4D;IACtE,UAAU,EAAE;QACV,sHAAsH;QACtH,eAAe;KAChB;IACD,SAAS,EAAE;QACT;YACE,EAAE,EAAE,KAAK;YACT,IAAI,EAAE,SAAS;YACf,IAAI,EAAE,gDAAgD;YACtD,SAAS,EAAE,wBAAwB,aAAa,+BAA+B,UAAU,uBAAuB;SACjH;KACF;IACD,MAAM,EAAE;QACN;YACE,EAAE,EAAE,KAAK;YACT,QAAQ,EAAE,CAAC,cAAc,EAAE,EAAE,CAAC,mBAAmB,CAAC,UAAU,EAAE,aAAa,EAAE,cAAc,CAAC;YAC5F,MAAM,EAAE,YAAY,CAAC,qBAAqB,CAAC,EAAE,MAAM,EAAE,UAAU,EAAE,CAAC,EAAE,aAAa,EAAE,CAAC,CAAC;YACrF,gBAAgB,EAAE,EAAE;YACpB,cAAc,EAAE,GAAG;YACnB,SAAS,EAAE,IAAI;SAChB;KACF;IACD,UAAU,EAAE,eAAe;CAC5B,CAAC","sourcesContent":["/**\n * `ab-pendulum-linear`'s declared translation from its bound's quantity, the\n * relative period error, into the accumulated PHASE error.\n *\n * Both motions start from rest at θ0 at t = 0. The phase is the angle\n * variable: 0 at release and advancing by 2π per period, uniformly in time,\n * so the pendulum's is 2π t/T and the linear oscillator's 2π t/T0. With the\n * bound's `ε = T/T0 − 1` (normalised by the reduced model's period),\n *\n *   Δφ(t) = 2π t/T0 − 2π t/T = 2π (t/T0) · ε/(1 + ε),\n *\n * which reaches a tolerance Δ at `t* = Δ · T0 · (1 + ε)/(2π ε)`. Nothing is\n * truncated: the map is exact in ε, and ε itself is the bound's closed form.\n *\n * It is NOT a bound on the pointwise position error |θ(t) − θ_lin(t)|: the\n * pendulum's waveform carries a third harmonic of amplitude ≈ θ0³/192, which\n * is not a phase, and the position error at a given phase error depends on\n * where in the cycle the observation falls.\n *\n * The witness W7p integrates BOTH equations of motion and reads each phase off\n * its zero crossings — crossing k is at phase π/2 + kπ by the symmetry of a\n * release from rest, and the angle variable is linear in time between them —\n * so it measures the drift without using ε or the elliptic integral. It runs\n * at the fixture θ0 = 0.2 as the regression anchor, and again at the caller's\n * point, bounded to 512 T0 of integration (`phasePointCheck`).\n *\n * @module atlas/oscillators/phase-translation\n * @internal\n */\n\nimport type { ObservableTranslation } from '../translation.js';\nimport { pendulumPeriodErrorAt } from './bridges-limits.js';\nimport { crossingTimes, linearAccel, pendulumAccel, phaseAt, theta0OfPeriodError, TWO_PI } from './pendulum-motion.js';\n\n/**\n * The time, in units of T0, at which the MEASURED phase lead of the linear\n * oscillator over the pendulum first reaches `tolerance` rad, both integrated\n * by RK4 at `stepsPerPeriod` steps of T0.\n *\n * @returns NaN if the drift does not reach `tolerance` within 2^16 periods.\n * @internal\n */\nfunction measurePhaseHorizon(theta0: number, tolerance: number, stepsPerPeriod: number): number {\n  for (let periods = 16; periods <= 2 ** 16; periods *= 2) {\n    const t = measurePhaseHorizonWithin(theta0, tolerance, stepsPerPeriod, periods);\n    if (!Number.isNaN(t)) return t;\n  }\n  return Number.NaN;\n}\n\n/** As `measurePhaseHorizon`, over exactly `periods` of T0; NaN if the drift does not reach `tolerance` in them. @internal */\nfunction measurePhaseHorizonWithin(theta0: number, tolerance: number, stepsPerPeriod: number, periods: number): number {\n  const pend = crossingTimes(pendulumAccel, theta0, stepsPerPeriod, periods);\n  const lin = crossingTimes(linearAccel, theta0, stepsPerPeriod, periods);\n  let prevT = 0;\n  let prevDrift = 0;\n  for (let k = 0; k < pend.length; k++) {\n    const c = pend[k]!;\n    const drift = phaseAt(lin, c) - (Math.PI / 2 + k * Math.PI);\n    if (!Number.isFinite(drift)) break;\n    if (drift >= tolerance) return prevT + ((tolerance - prevDrift) * (c - prevT)) / (drift - prevDrift);\n    prevT = c;\n    prevDrift = drift;\n  }\n  return Number.NaN;\n}\n\n/** `t* = Δ T0 (1+ε)/(2π ε)`; `Infinity` at ε = 0, where no phase accumulates. */\nfunction phaseHorizon(eps: number, tolerance: number, T0: number | undefined): number {\n  if (T0 === undefined || !Number.isFinite(T0) || T0 <= 0 || !Number.isFinite(eps) || eps < 0) return Number.NaN;\n  return eps === 0 ? Infinity : (tolerance * T0 * (1 + eps)) / (TWO_PI * eps);\n}\n\nconst W7P_THETA0 = 0.2;\nconst W7P_TOLERANCE = 1;\n\n/** The point witness integrates at most this many periods of T0, so its cost is bounded. */\nexport const PHASE_POINT_MAX_PERIODS = 512;\n/**\n * The point witness measures the time to this much drift, or to the drift\n * reached at the cap if less: the map is linear in t, so its slope is what is\n * checked, whatever tolerance the caller asked about.\n */\nconst PHASE_POINT_DRIFT = 1;\n/**\n * Below this measured drift the rounding noise of the two motions' crossing\n * times comes within a factor of ten of the witness tolerance. With the\n * crossings located by `crossingTimes`' secant refinement, the fine error was\n * measured at 800 steps per T0 as 1.5e-8 relative at 1.1e-5 rad of drift,\n * 7.4e-9 at 7.2e-5 rad, 2.5e-9 at 2.0e-4 rad and 2.4e-10 to 9.2e-10 from\n * 2e-2 rad up. With the straight-line crossings it replaced, the same error\n * was 9.5e-6 relative at 2e-4 rad and 1.2e-7 at 2e-2 rad.\n */\nconst PHASE_POINT_MIN_DRIFT = 1e-4;\n/**\n * The (1+ε) control differs from the target by ε/(1+ε) relative, so it can be\n * refuted only where ε exceeds this tolerance (θ0 above about 9e-4).\n */\nconst PHASE_POINT_RELATIVE_TOLERANCE = 5e-8;\n\n/**\n * W7p at the caller's θ0 (recovered from ε) instead of the fixture: the time\n * the measured drift reaches min(1 rad, the drift at the cap), within\n * PHASE_POINT_RELATIVE_TOLERANCE of the closed form. Its control is the same measurement against\n * the horizon without the (1+ε) factor, which differs from the target by\n * ε/(1+ε) relative — so the control can fail only where that exceeds the\n * witness tolerance, and the caller is told when it does not.\n */\nfunction phasePointCheck(eps: number): ReturnType<ObservableTranslation['pointCheck']> {\n  if (!Number.isFinite(eps) || eps <= 0) return { unavailable: `no phase drift accumulates at ε = ${eps}` };\n  const theta0 = theta0OfPeriodError(eps);\n  const measured = Math.min(PHASE_POINT_DRIFT, (TWO_PI * PHASE_POINT_MAX_PERIODS * eps) / (1 + eps));\n  if (measured < PHASE_POINT_MIN_DRIFT) {\n    return {\n      unavailable:\n        `the drift within ${PHASE_POINT_MAX_PERIODS} T0 at θ0 = ${Number(theta0.toPrecision(8))} is ${measured} rad, below the ` +\n        `${PHASE_POINT_MIN_DRIFT} rad this witness resolves; not run`,\n    };\n  }\n  const target = phaseHorizon(eps, measured, 1);\n  const periods = Math.ceil(target) + 2;\n  const evaluate = (stepsPerPeriod: number): number => measurePhaseHorizonWithin(theta0, measured, stepsPerPeriod, periods);\n  const spec = { evaluate, coarseResolution: 200, fineResolution: 800, tolerance: PHASE_POINT_RELATIVE_TOLERANCE * target };\n  const capped = measured < PHASE_POINT_DRIFT ? ` (the drift reached in ${PHASE_POINT_MAX_PERIODS} T0)` : '';\n  return {\n    claim:\n      `RK4-measured time to ${measured} rad${capped} of phase drift at θ0 = ${Number(theta0.toPrecision(8))} (recovered from ε) ` +\n      `within ${PHASE_POINT_RELATIVE_TOLERANCE} relative of t* = ${target} T0; 200 and 800 steps per T0, ${periods} T0 integrated`,\n    check: { id: 'W7p@point', ...spec, target },\n    control: { id: 'W7p@point-control', ...spec, target: measured / (TWO_PI * eps) },\n    controlClaim: 'the same measurement against the horizon without the (1+ε) factor, Δ T0/(2π ε)',\n  };\n}\n\n/** @internal */\nexport const PENDULUM_PHASE_TRANSLATION: ObservableTranslation = {\n  bridgeId: 'ab-pendulum-linear',\n  observable: 'phase',\n  unit: 'rad',\n  fromNorm: 'relative period error, normalized by the value of the reduced model',\n  errorKind: 'exact',\n  definition:\n    'the angle variable of each motion (0 at release from rest at θ0, 2π per period); the error is the ' +\n    \"linear oscillator's lead over the pendulum, Δφ(t) = 2πt/T0 − 2πt/T\",\n  derivation:\n    'with ε = T/T0 − 1, Δφ(t) = 2π (t/T0) ε/(1+ε), exact in ε; it reaches Δ at t* = Δ T0 (1+ε)/(2π ε). ' +\n    'ε/(1+ε) increases with ε, so the domain supremum of ε gives the shortest horizon over the domain',\n  premises: [\n    'both motions are released from rest at the same θ0 at t = 0',\n    'T0 = 2π √(ℓ/g) is the linear period, in the unit t is given in',\n  ],\n  parameters: ['T0'],\n  timeUnit: 'the unit of T0',\n  errorAt: (eps, t, { T0 }) => {\n    if (T0 === undefined || !Number.isFinite(T0) || T0 <= 0 || !Number.isFinite(eps) || !Number.isFinite(t)) return Number.NaN;\n    return (TWO_PI * (t / T0) * eps) / (1 + eps);\n  },\n  horizonFor: (eps, tolerance, { T0 }) => phaseHorizon(eps, tolerance, T0),\n  boundary: 'inclusive: adequate iff t ≤ t*, i.e. Δφ(t) ≤ the tolerance',\n  notCovered: [\n    'the pointwise position error |θ(t) − θ_lin(t)| (the waveform carries a third harmonic ≈ θ0³/192 that is not a phase)',\n    'the amplitude',\n  ],\n  witnesses: [\n    {\n      id: 'W7p',\n      kind: 'numeric',\n      test: 'tests/atlas/pendulum-phase-translation.test.ts',\n      tolerance: `RK4-measured time to ${W7P_TOLERANCE} rad of phase drift at θ0 = ${W7P_THETA0} within 1e-4 T0 of t*`,\n    },\n  ],\n  checks: [\n    {\n      id: 'W7p',\n      evaluate: (stepsPerPeriod) => measurePhaseHorizon(W7P_THETA0, W7P_TOLERANCE, stepsPerPeriod),\n      target: phaseHorizon(pendulumPeriodErrorAt({ theta0: W7P_THETA0 }), W7P_TOLERANCE, 1),\n      coarseResolution: 50,\n      fineResolution: 200,\n      tolerance: 1e-4,\n    },\n  ],\n  pointCheck: phasePointCheck,\n};\n"]}

package/dist/atlas/oscillators/position-translation.d.ts

@@ -23,9 +23,16 @@
  *
  * bounds the error at every t; the clamp at π keeps B non-decreasing in t and
  * is still a bound, since |cos a − cos b| ≤ 2. It is rigorous GIVEN the
- * Fourier series, which is quoted, not derived here. Witness W7xs checks the
- * series against RK4 over one period at each θ0 of a grid across the domain;
- * that is a check at those points, not a derivation.
+ * Fourier series, which is quoted, not derived here. The exact check
+ * (`fourierSeriesExactDeviation`, run by `tests/atlas/pendulum-exact.test.ts`;
+ * not a registered witness, because it has no resolution for a runner to
+ * refine) compares the series with the EXACT solution θ = 2 asin(k sn(K − ω0t, k²))
+ * (`exactPendulumAngle`) over a dense grid of amplitudes and phases at
+ * 1e-12 rad (1e-13 above θ0 = 0.002; the floor below that is the nome's, see
+ * `FOURIER_SERIES_EXACT_TOLERANCE`), and W7xs against RK4 over one period at
+ * six amplitudes at 1e-9;
+ * both are checks at grid points, not a derivation. The coefficients also
+ * have a second route, the FFT of one exact period (`fourierCoefficientsByFft`).
  *
  * B is an upper bound, not the error: past its horizon the BOUND exceeds the
  * tolerance, and the error may not. It is near-sharp where Δφ → π (the
@@ -37,17 +44,71 @@
  * @internal
  */
 import type { ObservableTranslation } from '../translation.js';
+import type { DominanceWitnessSpec } from '../witness-dominance.js';
 import type { NumericWitnessSpec } from '../witness-numeric.js';
 /** W̄ = (16/3) q^{3/2}/(1 − q): the bound on the waveform distortion |F(φ) − θ0 cos φ|. @internal */
 export declare function waveformBound(theta0: number): number;
 /** The pendulum's fundamental Fourier coefficient a₀ = 8√q/(1 + q). @internal */
 export declare function fundamentalCoefficient(theta0: number): number;
-/** The point witness integrates at most this many periods of T0. */
-export declare const POSITION_POINT_MAX_PERIODS = 512;
-/** a₀ measured by projecting one RK4 pendulum period onto cos φ (trapezoid, exact for a periodic integrand). */
-export declare function measureFundamental(theta0: number, stepsPerPeriod: number): number;
+/**
+ * The exact route to W7x's measured value: |θ_exact(t) − θ0 cos(2πt)| from the
+ * closed forms, at the same 200 samples per T0 as {@link dominanceSamples},
+ * against B(t) at θ0 recovered from `eps`. The resolution is ignored, so the
+ * runner's coarse-to-fine uncertainty u is exactly 0 here and its rule
+ * `bound − measured > u` is `bound > measured`: u is a property of the RK4
+ * route, not of the claim. The registered W7x keeps the RK4 route.
+ *
+ * @internal
+ */
+export declare function exactDominanceSpec(id: string, eps: number): DominanceWitnessSpec;
 /** The θ0 at which W7xs checks the Fourier series: a grid across the bridge's domain θ0 ≤ 0.5. @internal */
 export declare const FOURIER_SERIES_THETA0S: readonly number[];
+/** The pendulum's Fourier coefficient aₙ at amplitude θ0 (the closed form, from the nome). @internal */
+export declare function fourierCoefficient(theta0: number, n: number, variant?: 'plus' | 'minus'): number;
+/** The grid the exact check runs on: amplitudes θ0 = 0.5 j/500, j = 1..500, and phases φ = 2π i/256. @internal */
+export declare const FOURIER_SERIES_EXACT_GRID: {
+    readonly amplitudes: 500;
+    readonly phases: 256;
+};
+/**
+ * The exact check's tolerance over the whole grid, in rad. The measured deviation is
+ * 2.3e-13, all of it at the smallest amplitude θ0 = 0.001 and on the SERIES
+ * side: the nome's K′ = K(cos²(θ0/2)) is read from a complementary parameter
+ * that rounds at 1e-16 absolute, 4e-10 relative at m = 2.5e-7, so a₀ is low
+ * by 2.3e-10 of θ0 there. The exact solution agrees with an independent AGM
+ * sn to the bit at that point.
+ * @internal
+ */
+export declare const FOURIER_SERIES_EXACT_TOLERANCE = 1e-12;
+/** Above θ0 = 0.002 the same check holds at 1e-13 rad (measured 3.8e-14 at θ0 = 0.002, below 1.5e-14 from θ0 = 0.01). @internal */
+export declare const FOURIER_SERIES_EXACT_TIGHT: {
+    readonly theta0Min: 0.002;
+    readonly tolerance: 1e-13;
+};
+/**
+ * The exact check, or its control: max over the grid of |θ_exact(t) − F(φ)|, t = φT/2π,
+ * with θ_exact the elliptic solution (`exactPendulumAngle`) and F the 20-term
+ * series, over the amplitudes at or above `theta0Min`. A closed form against
+ * a closed form, so no resolution and no convergence ratio: the number is the
+ * deviation itself, held to {@link FOURIER_SERIES_EXACT_TOLERANCE} over the
+ * whole grid and to {@link FOURIER_SERIES_EXACT_TIGHT} above its floor by the
+ * test.
+ *
+ * @internal
+ */
+export declare function fourierSeriesExactDeviation(variant?: 'plus' | 'minus', grid?: {
+    readonly amplitudes: number;
+    readonly phases: number;
+}, theta0Min?: number): number;
+/**
+ * The first `count` coefficients aₙ from the FFT of one exact period sampled
+ * at `samples` points, MathTS `rfft`: aₙ = (−1)ⁿ 2 re[2n+1]/N. A second
+ * route to the series, independent of the nome. `samples` must exceed
+ * 2·count + 1.
+ *
+ * @internal
+ */
+export declare function fourierCoefficientsByFft(theta0: number, samples: number, count: number): number[];
 /**
  * W7xs, or its control: the Fourier series against RK4 of the pendulum, at 200
  * and 800 steps per period, within 1e-9 rad at every step and every grid θ0.

package/dist/atlas/oscillators/position-translation.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"position-translation.d.ts","sourceRoot":"","sources":["../../../src/atlas/oscillators/position-translation.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GAqCG;AAEH,OAAO,KAAK,EAAE,qBAAqB,EAAE,MAAM,mBAAmB,CAAC;AAE/D,OAAO,KAAK,EAAE,kBAAkB,EAAE,MAAM,uBAAuB,CAAC;AAMhE,qGAAqG;AACrG,wBAAgB,aAAa,CAAC,MAAM,EAAE,MAAM,GAAG,MAAM,CAGpD;AAED,iFAAiF;AACjF,wBAAgB,sBAAsB,CAAC,MAAM,EAAE,MAAM,GAAG,MAAM,CAG7D;AAsBD,oEAAoE;AACpE,eAAO,MAAM,0BAA0B,MAAM,CAAC;AA8C9C,gHAAgH;AAChH,wBAAgB,kBAAkB,CAAC,MAAM,EAAE,MAAM,EAAE,cAAc,EAAE,MAAM,GAAG,MAAM,CAWjF;AAED,4GAA4G;AAC5G,eAAO,MAAM,sBAAsB,EAAE,SAAS,MAAM,EAAoC,CAAC;AAmCzF;;;;GAIG;AACH,wBAAgB,kBAAkB,CAAC,EAAE,EAAE,MAAM,EAAE,OAAO,GAAE,MAAM,GAAG,OAAgB,GAAG,kBAAkB,CASrG;AA6BD,gBAAgB;AAChB,eAAO,MAAM,6BAA6B,EAAE,qBAqD3C,CAAC"}
+{"version":3,"file":"position-translation.d.ts","sourceRoot":"","sources":["../../../src/atlas/oscillators/position-translation.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GA4CG;AAGH,OAAO,KAAK,EAAE,qBAAqB,EAAE,MAAM,mBAAmB,CAAC;AAC/D,OAAO,KAAK,EAAE,oBAAoB,EAAE,MAAM,yBAAyB,CAAC;AACpE,OAAO,KAAK,EAAE,kBAAkB,EAAE,MAAM,uBAAuB,CAAC;AAYhE,qGAAqG;AACrG,wBAAgB,aAAa,CAAC,MAAM,EAAE,MAAM,GAAG,MAAM,CAGpD;AAED,iFAAiF;AACjF,wBAAgB,sBAAsB,CAAC,MAAM,EAAE,MAAM,GAAG,MAAM,CAG7D;AA2DD;;;;;;;;;GASG;AACH,wBAAgB,kBAAkB,CAAC,EAAE,EAAE,MAAM,EAAE,GAAG,EAAE,MAAM,GAAG,oBAAoB,CAqBhF;AA4BD,4GAA4G;AAC5G,eAAO,MAAM,sBAAsB,EAAE,SAAS,MAAM,EAAoC,CAAC;AAUzF,wGAAwG;AACxG,wBAAgB,kBAAkB,CAAC,MAAM,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,OAAO,GAAE,MAAM,GAAG,OAAgB,GAAG,MAAM,CAExG;AAgBD,kHAAkH;AAClH,eAAO,MAAM,yBAAyB;aAAK,UAAU,EAAE,GAAG;aAAE,MAAM,EAAE,GAAG;CAAW,CAAC;AACnF;;;;;;;;GAQG;AACH,eAAO,MAAM,8BAA8B,QAAQ,CAAC;AACpD,mIAAmI;AACnI,eAAO,MAAM,0BAA0B;aAAK,SAAS,EAAE,KAAK;aAAE,SAAS,EAAE,KAAK;CAAW,CAAC;AAE1F;;;;;;;;;;GAUG;AACH,wBAAgB,2BAA2B,CACzC,OAAO,GAAE,MAAM,GAAG,OAAgB,EAClC,IAAI,GAAE;IAAE,QAAQ,CAAC,UAAU,EAAE,MAAM,CAAC;IAAC,QAAQ,CAAC,MAAM,EAAE,MAAM,CAAA;CAA8B,EAC1F,SAAS,SAAI,GACZ,MAAM,CAgBR;AAED;;;;;;;GAOG;AACH,wBAAgB,wBAAwB,CAAC,MAAM,EAAE,MAAM,EAAE,OAAO,EAAE,MAAM,EAAE,KAAK,EAAE,MAAM,GAAG,MAAM,EAAE,CAOjG;AAkBD;;;;GAIG;AACH,wBAAgB,kBAAkB,CAAC,EAAE,EAAE,MAAM,EAAE,OAAO,GAAE,MAAM,GAAG,OAAgB,GAAG,kBAAkB,CASrG;AA6BD,gBAAgB;AAChB,eAAO,MAAM,6BAA6B,EAAE,qBAuD3C,CAAC"}

package/dist/atlas/oscillators/position-translation.js

@@ -23,9 +23,16 @@
  *
  * bounds the error at every t; the clamp at π keeps B non-decreasing in t and
  * is still a bound, since |cos a − cos b| ≤ 2. It is rigorous GIVEN the
- * Fourier series, which is quoted, not derived here. Witness W7xs checks the
- * series against RK4 over one period at each θ0 of a grid across the domain;
- * that is a check at those points, not a derivation.
+ * Fourier series, which is quoted, not derived here. The exact check
+ * (`fourierSeriesExactDeviation`, run by `tests/atlas/pendulum-exact.test.ts`;
+ * not a registered witness, because it has no resolution for a runner to
+ * refine) compares the series with the EXACT solution θ = 2 asin(k sn(K − ω0t, k²))
+ * (`exactPendulumAngle`) over a dense grid of amplitudes and phases at
+ * 1e-12 rad (1e-13 above θ0 = 0.002; the floor below that is the nome's, see
+ * `FOURIER_SERIES_EXACT_TOLERANCE`), and W7xs against RK4 over one period at
+ * six amplitudes at 1e-9;
+ * both are checks at grid points, not a derivation. The coefficients also
+ * have a second route, the FFT of one exact period (`fourierCoefficientsByFft`).
  *
  * B is an upper bound, not the error: past its horizon the BOUND exceeds the
  * tolerance, and the error may not. It is near-sharp where Δφ → π (the
@@ -36,9 +43,9 @@
  * @module atlas/oscillators/position-translation
  * @internal
  */
+import { rfft } from '@danielsimonjr/mathts-functions';
 import { pendulumPeriodErrorAt } from './bridges-limits.js';
-import { linearAccel, pendulumAccel, pendulumNome, rk4Step, theta0OfPeriodError } from './pendulum-motion.js';
-const TWO_PI = 2 * Math.PI;
+import { exactPendulumAngle, linearAccel, pendulumAccel, pendulumNome, rk4Step, theta0OfPeriodError, TWO_PI, } from './pendulum-motion.js';
 /** W̄ = (16/3) q^{3/2}/(1 − q): the bound on the waveform distortion |F(φ) − θ0 cos φ|. @internal */
 export function waveformBound(theta0) {
     const q = pendulumNome(theta0);
@@ -71,12 +78,14 @@ function positionHorizon(eps, tolerance, T0) {
     return (2 * Math.asin(r) * T0 * (1 + eps)) / (TWO_PI * eps);
 }
 /** The point witness integrates at most this many periods of T0. */
-export const POSITION_POINT_MAX_PERIODS = 512;
+const POSITION_POINT_MAX_PERIODS = 512;
 /** Samples per T0 both resolutions report at. */
 const SAMPLES_PER_PERIOD = 200;
 /**
  * |θ(t) − θ_lin(t)| from RK4 of both motions, and B(t) at the same samples,
- * over `periods` of T0; `stepsPerPeriod` must be a multiple of 200.
+ * over `periods` of T0; `stepsPerPeriod` must be a multiple of 200. The
+ * exact route is {@link exactDominanceSpec}; this one is what W7x registers,
+ * so the integrator stays part of what the witness exercises.
  */
 function dominanceSamples(theta0, periods, stepsPerPeriod, bound) {
     const stride = stepsPerPeriod / SAMPLES_PER_PERIOD;
@@ -102,6 +111,38 @@ function dominanceSamples(theta0, periods, stepsPerPeriod, bound) {
 function dominanceWindow(eps) {
     return Math.min(POSITION_POINT_MAX_PERIODS, (1 + eps) / (2 * eps));
 }
+/**
+ * The exact route to W7x's measured value: |θ_exact(t) − θ0 cos(2πt)| from the
+ * closed forms, at the same 200 samples per T0 as {@link dominanceSamples},
+ * against B(t) at θ0 recovered from `eps`. The resolution is ignored, so the
+ * runner's coarse-to-fine uncertainty u is exactly 0 here and its rule
+ * `bound − measured > u` is `bound > measured`: u is a property of the RK4
+ * route, not of the claim. The registered W7x keeps the RK4 route.
+ *
+ * @internal
+ */
+export function exactDominanceSpec(id, eps) {
+    const theta0 = theta0OfPeriodError(eps);
+    const periods = dominanceWindow(eps);
+    const bound = positionBound(eps);
+    const samples = Math.floor(periods * SAMPLES_PER_PERIOD);
+    return {
+        id,
+        kind: 'dominance',
+        evaluate: () => {
+            const measured = [];
+            const b = [];
+            for (let i = 1; i <= samples; i++) {
+                const t = i / SAMPLES_PER_PERIOD;
+                measured.push(Math.abs(exactPendulumAngle(theta0, t) - theta0 * Math.cos(TWO_PI * t)));
+                b.push(bound(t));
+            }
+            return { measured, bound: b };
+        },
+        coarseResolution: 200,
+        fineResolution: 800,
+    };
+}
 function dominanceSpec(id, eps, drift, waveform) {
     const theta0 = theta0OfPeriodError(eps);
     const periods = dominanceWindow(eps);
@@ -114,7 +155,7 @@ function dominanceSpec(id, eps, drift, waveform) {
     };
 }
 /** a₀ measured by projecting one RK4 pendulum period onto cos φ (trapezoid, exact for a periodic integrand). */
-export function measureFundamental(theta0, stepsPerPeriod) {
+function measureFundamental(theta0, stepsPerPeriod) {
     const T = 1 + pendulumPeriodErrorAt({ theta0 });
     const h = T / stepsPerPeriod;
     let x = theta0;
@@ -130,6 +171,15 @@ export function measureFundamental(theta0, stepsPerPeriod) {
 export const FOURIER_SERIES_THETA0S = [0.01, 0.1, 0.2, 0.3, 0.4, 0.5];
 const FOURIER_SERIES_TERMS = 20;
 const FOURIER_SERIES_TOLERANCE = 1e-9;
+/** aₙ = 8 q^{n+½}/((2n+1)(1 + q^{2n+1})) from the nome; `'minus'` is the (1 − q^{2n+1}) transcription error the controls use. */
+function coefficientFromNome(q, n, variant) {
+    const odd = 2 * n + 1;
+    return (8 * q ** (n + 0.5)) / (odd * (variant === 'plus' ? 1 + q ** odd : 1 - q ** odd));
+}
+/** The pendulum's Fourier coefficient aₙ at amplitude θ0 (the closed form, from the nome). @internal */
+export function fourierCoefficient(theta0, n, variant = 'plus') {
+    return coefficientFromNome(pendulumNome(theta0), n, variant);
+}
 /**
  * F(φ) = Σ_{n<20} (−1)ⁿ aₙ cos((2n+1)φ). `'minus'` puts (1 − q^{2n+1}) in the
  * denominator of aₙ, a transcription error the control uses.
@@ -138,12 +188,71 @@ function fourierSeries(theta0, phi, variant) {
     const q = pendulumNome(theta0);
     let s = 0;
     for (let n = 0; n < FOURIER_SERIES_TERMS; n++) {
-        const odd = 2 * n + 1;
-        const a = (8 * q ** (n + 0.5)) / (odd * (variant === 'plus' ? 1 + q ** odd : 1 - q ** odd));
-        s += (n % 2 === 0 ? a : -a) * Math.cos(odd * phi);
+        const a = coefficientFromNome(q, n, variant);
+        s += (n % 2 === 0 ? a : -a) * Math.cos((2 * n + 1) * phi);
     }
     return s;
 }
+/** The grid the exact check runs on: amplitudes θ0 = 0.5 j/500, j = 1..500, and phases φ = 2π i/256. @internal */
+export const FOURIER_SERIES_EXACT_GRID = { amplitudes: 500, phases: 256 };
+/**
+ * The exact check's tolerance over the whole grid, in rad. The measured deviation is
+ * 2.3e-13, all of it at the smallest amplitude θ0 = 0.001 and on the SERIES
+ * side: the nome's K′ = K(cos²(θ0/2)) is read from a complementary parameter
+ * that rounds at 1e-16 absolute, 4e-10 relative at m = 2.5e-7, so a₀ is low
+ * by 2.3e-10 of θ0 there. The exact solution agrees with an independent AGM
+ * sn to the bit at that point.
+ * @internal
+ */
+export const FOURIER_SERIES_EXACT_TOLERANCE = 1e-12;
+/** Above θ0 = 0.002 the same check holds at 1e-13 rad (measured 3.8e-14 at θ0 = 0.002, below 1.5e-14 from θ0 = 0.01). @internal */
+export const FOURIER_SERIES_EXACT_TIGHT = { theta0Min: 0.002, tolerance: 1e-13 };
+/**
+ * The exact check, or its control: max over the grid of |θ_exact(t) − F(φ)|, t = φT/2π,
+ * with θ_exact the elliptic solution (`exactPendulumAngle`) and F the 20-term
+ * series, over the amplitudes at or above `theta0Min`. A closed form against
+ * a closed form, so no resolution and no convergence ratio: the number is the
+ * deviation itself, held to {@link FOURIER_SERIES_EXACT_TOLERANCE} over the
+ * whole grid and to {@link FOURIER_SERIES_EXACT_TIGHT} above its floor by the
+ * test.
+ *
+ * @internal
+ */
+export function fourierSeriesExactDeviation(variant = 'plus', grid = FOURIER_SERIES_EXACT_GRID, theta0Min = 0) {
+    let worst = 0;
+    for (let j = 1; j <= grid.amplitudes; j++) {
+        const theta0 = (0.5 * j) / grid.amplitudes;
+        if (theta0 < theta0Min)
+            continue;
+        const T = 1 + pendulumPeriodErrorAt({ theta0 });
+        const q = pendulumNome(theta0);
+        const a = Array.from({ length: FOURIER_SERIES_TERMS }, (_, n) => coefficientFromNome(q, n, variant));
+        for (let i = 0; i < grid.phases; i++) {
+            const phi = (TWO_PI * i) / grid.phases;
+            let s = 0;
+            for (let n = 0; n < FOURIER_SERIES_TERMS; n++)
+                s += (n % 2 === 0 ? a[n] : -a[n]) * Math.cos((2 * n + 1) * phi);
+            worst = Math.max(worst, Math.abs(exactPendulumAngle(theta0, (i * T) / grid.phases) - s));
+        }
+    }
+    return worst;
+}
+/**
+ * The first `count` coefficients aₙ from the FFT of one exact period sampled
+ * at `samples` points, MathTS `rfft`: aₙ = (−1)ⁿ 2 re[2n+1]/N. A second
+ * route to the series, independent of the nome. `samples` must exceed
+ * 2·count + 1.
+ *
+ * @internal
+ */
+export function fourierCoefficientsByFft(theta0, samples, count) {
+    if (!Number.isInteger(samples) || samples <= 2 * count + 1) {
+        throw new RangeError(`fourierCoefficientsByFft: ${samples} samples cannot resolve ${count} odd harmonics`);
+    }
+    const T = 1 + pendulumPeriodErrorAt({ theta0 });
+    const { re } = rfft(Array.from({ length: samples }, (_, i) => exactPendulumAngle(theta0, (i * T) / samples)));
+    return Array.from({ length: count }, (_, n) => ((n % 2 === 0 ? 2 : -2) * re[2 * n + 1]) / samples);
+}
 /** max over the grid, and over one period T = T0(1+ε) sampled at every RK4 step, of |θ_RK4(t) − F(2πt/T)|. */
 function fourierSeriesDeviation(stepsPerPeriod, variant) {
     let worst = 0;
@@ -214,7 +323,9 @@ export const PENDULUM_POSITION_TRANSLATION = {
         'both motions are released from rest at the same θ0 at t = 0',
         'T0 = 2π √(ℓ/g) is the linear period, in the unit t is given in',
         'the classical Fourier series of the pendulum (all aₙ > 0, Σ (−1)ⁿ aₙ = θ0) — quoted, not derived here; ' +
-            `machine-checked against RK4 at θ0 ∈ {${FOURIER_SERIES_THETA0S.join(', ')}} (W7xs), not between those points`,
+            `machine-checked against the elliptic solution θ = 2 asin(k sn(K − ω0t, k²)) over ${FOURIER_SERIES_EXACT_GRID.amplitudes} amplitudes in (0, 0.5] ` +
+            `and ${FOURIER_SERIES_EXACT_GRID.phases} phases within ${FOURIER_SERIES_EXACT_TOLERANCE} rad, and within ${FOURIER_SERIES_EXACT_TIGHT.tolerance} rad ` +
+            `above θ0 = ${FOURIER_SERIES_EXACT_TIGHT.theta0Min} (tests/atlas/pendulum-exact.test.ts), and against RK4 at θ0 ∈ {${FOURIER_SERIES_THETA0S.join(', ')}} (W7xs)`,
     ],
     parameters: ['T0'],
     timeUnit: 'the unit of T0',

package/dist/atlas/oscillators/position-translation.js.map

@@ -1 +1 @@
-{"version":3,"file":"position-translation.js","sourceRoot":"","sources":["../../../src/atlas/oscillators/position-translation.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GAqCG;AAKH,OAAO,EAAE,qBAAqB,EAAE,MAAM,qBAAqB,CAAC;AAC5D,OAAO,EAAE,WAAW,EAAE,aAAa,EAAE,YAAY,EAAE,OAAO,EAAE,mBAAmB,EAAE,MAAM,sBAAsB,CAAC;AAE9G,MAAM,MAAM,GAAG,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC;AAE3B,qGAAqG;AACrG,MAAM,UAAU,aAAa,CAAC,MAAc;IAC1C,MAAM,CAAC,GAAG,YAAY,CAAC,MAAM,CAAC,CAAC;IAC/B,OAAO,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,IAAI,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;AACzC,CAAC;AAED,iFAAiF;AACjF,MAAM,UAAU,sBAAsB,CAAC,MAAc;IACnD,MAAM,CAAC,GAAG,YAAY,CAAC,MAAM,CAAC,CAAC;IAC/B,OAAO,CAAC,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;AACtC,CAAC;AAED,MAAM,IAAI,GAAG,CAAC,GAAW,EAAE,OAAe,EAAU,EAAE,CAAC,CAAC,MAAM,GAAG,OAAO,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC;AAE5F,qGAAqG;AACrG,SAAS,aAAa,CAAC,GAAW,EAAE,KAAK,GAAG,CAAC,EAAE,QAAQ,GAAG,CAAC;IACzD,MAAM,MAAM,GAAG,mBAAmB,CAAC,GAAG,CAAC,CAAC;IACxC,MAAM,KAAK,GAAG,QAAQ,GAAG,aAAa,CAAC,MAAM,CAAC,CAAC;IAC/C,OAAO,CAAC,OAAO,EAAE,EAAE,CAAC,KAAK,GAAG,KAAK,GAAG,CAAC,GAAG,MAAM,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,EAAE,OAAO,CAAC,EAAE,IAAI,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC;AACvG,CAAC;AAED,SAAS,eAAe,CAAC,GAAW,EAAE,SAAiB,EAAE,EAAsB;IAC7E,IAAI,EAAE,KAAK,SAAS,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,EAAE,CAAC,IAAI,EAAE,IAAI,CAAC,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,GAAG,CAAC,IAAI,GAAG,GAAG,CAAC;QAAE,OAAO,MAAM,CAAC,GAAG,CAAC;IAC/G,IAAI,GAAG,KAAK,CAAC;QAAE,OAAO,QAAQ,CAAC;IAC/B,MAAM,MAAM,GAAG,mBAAmB,CAAC,GAAG,CAAC,CAAC;IACxC,MAAM,KAAK,GAAG,aAAa,CAAC,MAAM,CAAC,CAAC;IACpC,IAAI,SAAS,GAAG,KAAK;QAAE,OAAO,MAAM,CAAC,GAAG,CAAC;IACzC,MAAM,CAAC,GAAG,CAAC,SAAS,GAAG,KAAK,CAAC,GAAG,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC;IAC7C,IAAI,CAAC,IAAI,CAAC;QAAE,OAAO,QAAQ,CAAC;IAC5B,OAAO,CAAC,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,GAAG,CAAC,MAAM,GAAG,GAAG,CAAC,CAAC;AAC9D,CAAC;AAED,oEAAoE;AACpE,MAAM,CAAC,MAAM,0BAA0B,GAAG,GAAG,CAAC;AAC9C,iDAAiD;AACjD,MAAM,kBAAkB,GAAG,GAAG,CAAC;AAE/B;;;GAGG;AACH,SAAS,gBAAgB,CAAC,MAAc,EAAE,OAAe,EAAE,cAAsB,EAAE,KAA4B;IAC7G,MAAM,MAAM,GAAG,cAAc,GAAG,kBAAkB,CAAC;IACnD,MAAM,CAAC,GAAG,CAAC,GAAG,cAAc,CAAC;IAC7B,MAAM,QAAQ,GAAa,EAAE,CAAC;IAC9B,MAAM,CAAC,GAAa,EAAE,CAAC;IACvB,IAAI,CAAC,GAAG,MAAM,CAAC;IACf,IAAI,CAAC,GAAG,CAAC,CAAC;IACV,IAAI,CAAC,GAAG,MAAM,CAAC;IACf,IAAI,CAAC,GAAG,CAAC,CAAC;IACV,MAAM,KAAK,GAAG,IAAI,CAAC,KAAK,CAAC,OAAO,GAAG,kBAAkB,CAAC,GAAG,MAAM,CAAC;IAChE,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,KAAK,EAAE,CAAC,EAAE,EAAE,CAAC;QAChC,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,OAAO,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,aAAa,CAAC,CAAC;QACzC,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,OAAO,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,WAAW,CAAC,CAAC;QACvC,IAAI,CAAC,GAAG,MAAM,KAAK,CAAC,EAAE,CAAC;YACrB,QAAQ,CAAC,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;YAC/B,CAAC,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;QACvB,CAAC;IACH,CAAC;IACD,OAAO,EAAE,QAAQ,EAAE,KAAK,EAAE,CAAC,EAAE,CAAC;AAChC,CAAC;AAED,mHAAmH;AACnH,SAAS,eAAe,CAAC,GAAW;IAClC,OAAO,IAAI,CAAC,GAAG,CAAC,0BAA0B,EAAE,CAAC,CAAC,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC;AACrE,CAAC;AAED,SAAS,aAAa,CAAC,EAAU,EAAE,GAAW,EAAE,KAAa,EAAE,QAAgB;IAC7E,MAAM,MAAM,GAAG,mBAAmB,CAAC,GAAG,CAAC,CAAC;IACxC,MAAM,OAAO,GAAG,eAAe,CAAC,GAAG,CAAC,CAAC;IACrC,OAAO;QACL,EAAE;QACF,IAAI,EAAE,WAAW;QACjB,QAAQ,EAAE,CAAC,cAAc,EAAE,EAAE,CAAC,gBAAgB,CAAC,MAAM,EAAE,OAAO,EAAE,cAAc,EAAE,aAAa,CAAC,GAAG,EAAE,KAAK,EAAE,QAAQ,CAAC,CAAC;QACpH,gBAAgB,EAAE,GAAG;QACrB,cAAc,EAAE,GAAG;KACpB,CAAC;AACJ,CAAC;AAED,gHAAgH;AAChH,MAAM,UAAU,kBAAkB,CAAC,MAAc,EAAE,cAAsB;IACvE,MAAM,CAAC,GAAG,CAAC,GAAG,qBAAqB,CAAC,EAAE,MAAM,EAAE,CAAC,CAAC;IAChD,MAAM,CAAC,GAAG,CAAC,GAAG,cAAc,CAAC;IAC7B,IAAI,CAAC,GAAG,MAAM,CAAC;IACf,IAAI,CAAC,GAAG,CAAC,CAAC;IACV,IAAI,GAAG,GAAG,MAAM,CAAC;IACjB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,cAAc,EAAE,CAAC,EAAE,EAAE,CAAC;QACxC,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,OAAO,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,aAAa,CAAC,CAAC;QACzC,GAAG,IAAI,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,MAAM,GAAG,CAAC,CAAC,GAAG,cAAc,CAAC,CAAC;IACrD,CAAC;IACD,OAAO,CAAC,CAAC,GAAG,GAAG,CAAC,GAAG,cAAc,CAAC;AACpC,CAAC;AAED,4GAA4G;AAC5G,MAAM,CAAC,MAAM,sBAAsB,GAAsB,CAAC,IAAI,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,CAAC,CAAC;AACzF,MAAM,oBAAoB,GAAG,EAAE,CAAC;AAChC,MAAM,wBAAwB,GAAG,IAAI,CAAC;AAEtC;;;GAGG;AACH,SAAS,aAAa,CAAC,MAAc,EAAE,GAAW,EAAE,OAAyB;IAC3E,MAAM,CAAC,GAAG,YAAY,CAAC,MAAM,CAAC,CAAC;IAC/B,IAAI,CAAC,GAAG,CAAC,CAAC;IACV,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,oBAAoB,EAAE,CAAC,EAAE,EAAE,CAAC;QAC9C,MAAM,GAAG,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;QACtB,MAAM,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC,OAAO,KAAK,MAAM,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,IAAI,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,IAAI,GAAG,CAAC,CAAC,CAAC;QAC5F,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,GAAG,GAAG,CAAC,CAAC;IACpD,CAAC;IACD,OAAO,CAAC,CAAC;AACX,CAAC;AAED,8GAA8G;AAC9G,SAAS,sBAAsB,CAAC,cAAsB,EAAE,OAAyB;IAC/E,IAAI,KAAK,GAAG,CAAC,CAAC;IACd,KAAK,MAAM,MAAM,IAAI,sBAAsB,EAAE,CAAC;QAC5C,MAAM,CAAC,GAAG,CAAC,CAAC,GAAG,qBAAqB,CAAC,EAAE,MAAM,EAAE,CAAC,CAAC,GAAG,cAAc,CAAC;QACnE,IAAI,CAAC,GAAG,MAAM,CAAC;QACf,IAAI,CAAC,GAAG,CAAC,CAAC;QACV,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,EAAE,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,aAAa,CAAC,MAAM,EAAE,CAAC,EAAE,OAAO,CAAC,CAAC,CAAC,CAAC;QACzE,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,cAAc,EAAE,CAAC,EAAE,EAAE,CAAC;YACzC,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,OAAO,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,aAAa,CAAC,CAAC;YACzC,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,EAAE,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,aAAa,CAAC,MAAM,EAAE,CAAC,MAAM,GAAG,CAAC,CAAC,GAAG,cAAc,EAAE,OAAO,CAAC,CAAC,CAAC,CAAC;QACvG,CAAC;IACH,CAAC;IACD,OAAO,KAAK,CAAC;AACf,CAAC;AAED;;;;GAIG;AACH,MAAM,UAAU,kBAAkB,CAAC,EAAU,EAAE,OAAO,GAAqB,MAAM;IAC/E,OAAO;QACL,EAAE;QACF,QAAQ,EAAE,CAAC,cAAc,EAAE,EAAE,CAAC,sBAAsB,CAAC,cAAc,EAAE,OAAO,CAAC;QAC7E,MAAM,EAAE,CAAC;QACT,gBAAgB,EAAE,GAAG;QACrB,cAAc,EAAE,GAAG;QACnB,SAAS,EAAE,wBAAwB;KACpC,CAAC;AACJ,CAAC;AAED,MAAM,UAAU,GAAG,GAAG,CAAC;AACvB,MAAM,OAAO,GAAG,qBAAqB,CAAC,EAAE,MAAM,EAAE,UAAU,EAAE,CAAC,CAAC;AAE9D,MAAM,IAAI,GAAuB;IAC/B,EAAE,EAAE,MAAM;IACV,QAAQ,EAAE,CAAC,cAAc,EAAE,EAAE,CAAC,kBAAkB,CAAC,UAAU,EAAE,cAAc,CAAC;IAC5E,MAAM,EAAE,sBAAsB,CAAC,UAAU,CAAC;IAC1C,gBAAgB,EAAE,GAAG;IACrB,cAAc,EAAE,GAAG;IACnB,SAAS,EAAE,KAAK;CACjB,CAAC;AAEF,SAAS,kBAAkB,CAAC,GAAW;IACrC,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,GAAG,CAAC,IAAI,GAAG,IAAI,CAAC;QAAE,OAAO,EAAE,WAAW,EAAE,iBAAiB,GAAG,+CAA+C,EAAE,CAAC;IACnI,MAAM,MAAM,GAAG,mBAAmB,CAAC,GAAG,CAAC,CAAC;IACxC,MAAM,OAAO,GAAG,eAAe,CAAC,GAAG,CAAC,CAAC;IACrC,MAAM,MAAM,GAAG,OAAO,KAAK,0BAA0B,CAAC,CAAC,CAAC,eAAe,0BAA0B,wBAAwB,CAAC,CAAC,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,CAAC,MAAM,CAAC,CAAC,CAAC,yBAAyB,CAAC;IACjL,OAAO;QACL,KAAK,EACH,gDAAgD,MAAM,CAAC,MAAM,CAAC,WAAW,CAAC,CAAC,CAAC,CAAC,iCAAiC,OAAO,MAAM,MAAM,IAAI;YACrI,gGAAgG;QAClG,KAAK,EAAE,aAAa,CAAC,WAAW,EAAE,GAAG,EAAE,CAAC,EAAE,CAAC,CAAC;QAC5C,OAAO,EAAE,aAAa,CAAC,mBAAmB,EAAE,GAAG,EAAE,CAAC,EAAE,CAAC,CAAC;QACtD,YAAY,EAAE,kEAAkE;KACjF,CAAC;AACJ,CAAC;AAED,gBAAgB;AAChB,MAAM,CAAC,MAAM,6BAA6B,GAA0B;IAClE,QAAQ,EAAE,oBAAoB;IAC9B,UAAU,EAAE,UAAU;IACtB,IAAI,EAAE,KAAK;IACX,QAAQ,EAAE,qEAAqE;IAC/E,SAAS,EAAE,aAAa;IACxB,UAAU,EAAE,yFAAyF;IACrG,UAAU,EACR,uGAAuG;QACvG,0GAA0G;QAC1G,uIAAuI;QACvI,4FAA4F;IAC9F,QAAQ,EAAE;QACR,6DAA6D;QAC7D,gEAAgE;QAChE,yGAAyG;YACvG,wCAAwC,sBAAsB,CAAC,IAAI,CAAC,IAAI,CAAC,oCAAoC;KAChH;IACD,UAAU,EAAE,CAAC,IAAI,CAAC;IAClB,QAAQ,EAAE,gBAAgB;IAC1B,OAAO,EAAE,CAAC,GAAG,EAAE,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE;QAC1B,IAAI,EAAE,KAAK,SAAS,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,EAAE,CAAC,IAAI,EAAE,IAAI,CAAC,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,GAAG,CAAC,IAAI,GAAG,GAAG,CAAC,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC;YAAE,OAAO,MAAM,CAAC,GAAG,CAAC;QACtI,OAAO,aAAa,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC;IACpC,CAAC;IACD,UAAU,EAAE,CAAC,GAAG,EAAE,SAAS,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,eAAe,CAAC,GAAG,EAAE,SAAS,EAAE,EAAE,CAAC;IAC3E,KAAK,EAAE,CAAC,GAAG,EAAE,EAAE,CAAC,aAAa,CAAC,mBAAmB,CAAC,GAAG,CAAC,CAAC;IACvD,QAAQ,EACN,+IAA+I;IACjJ,UAAU,EAAE,CAAC,wDAAwD,EAAE,sBAAsB,EAAE,eAAe,CAAC;IAC/G,SAAS,EAAE;QACT;YACE,EAAE,EAAE,KAAK;YACT,IAAI,EAAE,SAAS;YACf,IAAI,EAAE,mDAAmD;YACzD,SAAS,EAAE,6DAA6D,UAAU,8CAA8C;SACjI;QACD;YACE,EAAE,EAAE,MAAM;YACV,IAAI,EAAE,SAAS;YACf,IAAI,EAAE,mDAAmD;YACzD,SAAS,EAAE,iDAAiD,UAAU,4BAA4B;SACnG;QACD;YACE,EAAE,EAAE,MAAM;YACV,IAAI,EAAE,SAAS;YACf,IAAI,EAAE,mDAAmD;YACzD,SAAS,EACP,uBAAuB,oBAAoB,kBAAkB,wBAAwB,2CAA2C;gBAChI,YAAY,sBAAsB,CAAC,IAAI,CAAC,IAAI,CAAC,GAAG;SACnD;KACF;IACD,MAAM,EAAE,CAAC,aAAa,CAAC,KAAK,EAAE,OAAO,EAAE,CAAC,EAAE,CAAC,CAAC,EAAE,IAAI,EAAE,kBAAkB,CAAC,MAAM,CAAC,CAAC;IAC/E,UAAU,EAAE,kBAAkB;CAC/B,CAAC","sourcesContent":["/**\n * `ab-pendulum-linear`'s declared translation from its bound's quantity, the\n * relative period error ε, into an UPPER BOUND on the pointwise position error\n * |θ(t) − θ_lin(t)|, in rad.\n *\n * Both motions start from rest at θ0. The pendulum is θ(t) = F(φ), φ = 2πt/T,\n * with the classical Fourier series of the elliptic solution\n *\n *   F(φ) = Σ_{n≥0} (−1)ⁿ aₙ cos((2n+1)φ),  aₙ = 8 q^{n+½} / ((2n+1)(1 + q^{2n+1})),\n *\n * q the nome of k = sin(θ0/2); F(0) = θ0 = Σ (−1)ⁿ aₙ. The linear motion is\n * θ0 cos φ_lin, φ_lin = 2πt/T0. Splitting at θ0 cos φ,\n *\n *   |θ − θ_lin| ≤ |F(φ) − θ0 cos φ| + θ0 |cos φ − cos φ_lin|\n *              ≤ 2 Σ_{n≥1} aₙ + 2θ0 |sin(Δφ/2)|,\n *\n * because F(φ) − θ0 cos φ = Σ_{n≥1} (−1)ⁿ aₙ (cos((2n+1)φ) − cos φ). With\n * aₙ ≤ (8/3) q^{n+½} for n ≥ 1 the first term is at most\n * W̄ = (16/3) q^{3/2}/(1 − q) ≈ θ0³/96, and Δφ = 2π(t/T0) ε/(1+ε) is the\n * exact phase lead of `phase-translation.ts`. So\n *\n *   B(t) = W̄ + 2θ0 sin(min(Δφ(t), π)/2)\n *\n * bounds the error at every t; the clamp at π keeps B non-decreasing in t and\n * is still a bound, since |cos a − cos b| ≤ 2. It is rigorous GIVEN the\n * Fourier series, which is quoted, not derived here. Witness W7xs checks the\n * series against RK4 over one period at each θ0 of a grid across the domain;\n * that is a check at those points, not a derivation.\n *\n * B is an upper bound, not the error: past its horizon the BOUND exceeds the\n * tolerance, and the error may not. It is near-sharp where Δφ → π (the\n * motions are in antiphase and the error reaches 2θ0), and the waveform term\n * cannot be dropped: without W̄, RK4 finds the error above 2θ0 sin(Δφ/2) by\n * about 3% (witness control).\n *\n * @module atlas/oscillators/position-translation\n * @internal\n */\n\nimport type { ObservableTranslation } from '../translation.js';\nimport type { DominanceWitnessSpec } from '../witness-dominance.js';\nimport type { NumericWitnessSpec } from '../witness-numeric.js';\nimport { pendulumPeriodErrorAt } from './bridges-limits.js';\nimport { linearAccel, pendulumAccel, pendulumNome, rk4Step, theta0OfPeriodError } from './pendulum-motion.js';\n\nconst TWO_PI = 2 * Math.PI;\n\n/** W̄ = (16/3) q^{3/2}/(1 − q): the bound on the waveform distortion |F(φ) − θ0 cos φ|. @internal */\nexport function waveformBound(theta0: number): number {\n  const q = pendulumNome(theta0);\n  return ((16 / 3) * q ** 1.5) / (1 - q);\n}\n\n/** The pendulum's fundamental Fourier coefficient a₀ = 8√q/(1 + q). @internal */\nexport function fundamentalCoefficient(theta0: number): number {\n  const q = pendulumNome(theta0);\n  return (8 * Math.sqrt(q)) / (1 + q);\n}\n\nconst lead = (eps: number, tOverT0: number): number => (TWO_PI * tOverT0 * eps) / (1 + eps);\n\n/** B(t/T0) at θ0 recovered from ε; `drift` and `waveform` scale its two terms (1, 1 is B itself). */\nfunction positionBound(eps: number, drift = 1, waveform = 1): (tOverT0: number) => number {\n  const theta0 = theta0OfPeriodError(eps);\n  const floor = waveform * waveformBound(theta0);\n  return (tOverT0) => floor + drift * 2 * theta0 * Math.sin(Math.min(lead(eps, tOverT0), Math.PI) / 2);\n}\n\nfunction positionHorizon(eps: number, tolerance: number, T0: number | undefined): number {\n  if (T0 === undefined || !Number.isFinite(T0) || T0 <= 0 || !Number.isFinite(eps) || eps < 0) return Number.NaN;\n  if (eps === 0) return Infinity;\n  const theta0 = theta0OfPeriodError(eps);\n  const floor = waveformBound(theta0);\n  if (tolerance < floor) return Number.NaN;\n  const r = (tolerance - floor) / (2 * theta0);\n  if (r >= 1) return Infinity;\n  return (2 * Math.asin(r) * T0 * (1 + eps)) / (TWO_PI * eps);\n}\n\n/** The point witness integrates at most this many periods of T0. */\nexport const POSITION_POINT_MAX_PERIODS = 512;\n/** Samples per T0 both resolutions report at. */\nconst SAMPLES_PER_PERIOD = 200;\n\n/**\n * |θ(t) − θ_lin(t)| from RK4 of both motions, and B(t) at the same samples,\n * over `periods` of T0; `stepsPerPeriod` must be a multiple of 200.\n */\nfunction dominanceSamples(theta0: number, periods: number, stepsPerPeriod: number, bound: (t: number) => number) {\n  const stride = stepsPerPeriod / SAMPLES_PER_PERIOD;\n  const h = 1 / stepsPerPeriod;\n  const measured: number[] = [];\n  const b: number[] = [];\n  let x = theta0;\n  let v = 0;\n  let y = theta0;\n  let u = 0;\n  const steps = Math.floor(periods * SAMPLES_PER_PERIOD) * stride;\n  for (let i = 1; i <= steps; i++) {\n    [x, v] = rk4Step(x, v, h, pendulumAccel);\n    [y, u] = rk4Step(y, u, h, linearAccel);\n    if (i % stride === 0) {\n      measured.push(Math.abs(x - y));\n      b.push(bound(i * h));\n    }\n  }\n  return { measured, bound: b };\n}\n\n/** The window a dominance witness covers: up to antiphase (Δφ = π), past which B ≥ 2θ0 holds trivially, capped. */\nfunction dominanceWindow(eps: number): number {\n  return Math.min(POSITION_POINT_MAX_PERIODS, (1 + eps) / (2 * eps));\n}\n\nfunction dominanceSpec(id: string, eps: number, drift: number, waveform: number): DominanceWitnessSpec {\n  const theta0 = theta0OfPeriodError(eps);\n  const periods = dominanceWindow(eps);\n  return {\n    id,\n    kind: 'dominance',\n    evaluate: (stepsPerPeriod) => dominanceSamples(theta0, periods, stepsPerPeriod, positionBound(eps, drift, waveform)),\n    coarseResolution: 200,\n    fineResolution: 800,\n  };\n}\n\n/** a₀ measured by projecting one RK4 pendulum period onto cos φ (trapezoid, exact for a periodic integrand). */\nexport function measureFundamental(theta0: number, stepsPerPeriod: number): number {\n  const T = 1 + pendulumPeriodErrorAt({ theta0 });\n  const h = T / stepsPerPeriod;\n  let x = theta0;\n  let v = 0;\n  let sum = theta0;\n  for (let i = 1; i < stepsPerPeriod; i++) {\n    [x, v] = rk4Step(x, v, h, pendulumAccel);\n    sum += x * Math.cos((TWO_PI * i) / stepsPerPeriod);\n  }\n  return (2 * sum) / stepsPerPeriod;\n}\n\n/** The θ0 at which W7xs checks the Fourier series: a grid across the bridge's domain θ0 ≤ 0.5. @internal */\nexport const FOURIER_SERIES_THETA0S: readonly number[] = [0.01, 0.1, 0.2, 0.3, 0.4, 0.5];\nconst FOURIER_SERIES_TERMS = 20;\nconst FOURIER_SERIES_TOLERANCE = 1e-9;\n\n/**\n * F(φ) = Σ_{n<20} (−1)ⁿ aₙ cos((2n+1)φ). `'minus'` puts (1 − q^{2n+1}) in the\n * denominator of aₙ, a transcription error the control uses.\n */\nfunction fourierSeries(theta0: number, phi: number, variant: 'plus' | 'minus'): number {\n  const q = pendulumNome(theta0);\n  let s = 0;\n  for (let n = 0; n < FOURIER_SERIES_TERMS; n++) {\n    const odd = 2 * n + 1;\n    const a = (8 * q ** (n + 0.5)) / (odd * (variant === 'plus' ? 1 + q ** odd : 1 - q ** odd));\n    s += (n % 2 === 0 ? a : -a) * Math.cos(odd * phi);\n  }\n  return s;\n}\n\n/** max over the grid, and over one period T = T0(1+ε) sampled at every RK4 step, of |θ_RK4(t) − F(2πt/T)|. */\nfunction fourierSeriesDeviation(stepsPerPeriod: number, variant: 'plus' | 'minus'): number {\n  let worst = 0;\n  for (const theta0 of FOURIER_SERIES_THETA0S) {\n    const h = (1 + pendulumPeriodErrorAt({ theta0 })) / stepsPerPeriod;\n    let x = theta0;\n    let v = 0;\n    worst = Math.max(worst, Math.abs(x - fourierSeries(theta0, 0, variant)));\n    for (let i = 1; i <= stepsPerPeriod; i++) {\n      [x, v] = rk4Step(x, v, h, pendulumAccel);\n      worst = Math.max(worst, Math.abs(x - fourierSeries(theta0, (TWO_PI * i) / stepsPerPeriod, variant)));\n    }\n  }\n  return worst;\n}\n\n/**\n * W7xs, or its control: the Fourier series against RK4 of the pendulum, at 200\n * and 800 steps per period, within 1e-9 rad at every step and every grid θ0.\n * @internal\n */\nexport function fourierSeriesCheck(id: string, variant: 'plus' | 'minus' = 'plus'): NumericWitnessSpec {\n  return {\n    id,\n    evaluate: (stepsPerPeriod) => fourierSeriesDeviation(stepsPerPeriod, variant),\n    target: 0,\n    coarseResolution: 200,\n    fineResolution: 800,\n    tolerance: FOURIER_SERIES_TOLERANCE,\n  };\n}\n\nconst W7X_THETA0 = 0.3;\nconst W7X_EPS = pendulumPeriodErrorAt({ theta0: W7X_THETA0 });\n\nconst W7XA: NumericWitnessSpec = {\n  id: 'W7xa',\n  evaluate: (stepsPerPeriod) => measureFundamental(W7X_THETA0, stepsPerPeriod),\n  target: fundamentalCoefficient(W7X_THETA0),\n  coarseResolution: 200,\n  fineResolution: 800,\n  tolerance: 1e-10,\n};\n\nfunction positionPointCheck(eps: number): ReturnType<ObservableTranslation['pointCheck']> {\n  if (!Number.isFinite(eps) || eps <= 0) return { unavailable: `θ0 = 0 at ε = ${eps}: both motions are at rest and the error is 0` };\n  const theta0 = theta0OfPeriodError(eps);\n  const periods = dominanceWindow(eps);\n  const capped = periods === POSITION_POINT_MAX_PERIODS ? `, capped at ${POSITION_POINT_MAX_PERIODS} T0 (antiphase is at ${(1 + eps) / (2 * eps)} T0)` : ' (to antiphase, Δφ = π)';\n  return {\n    claim:\n      `RK4 |θ − θ_lin| ≤ B(t) at every sample, θ0 = ${Number(theta0.toPrecision(8))} (recovered from ε), over t ≤ ${periods} T0${capped}; ` +\n      `200 and 800 steps per T0, sampled 200 per T0; the coarse-to-fine difference is the uncertainty`,\n    check: dominanceSpec('W7x@point', eps, 1, 1),\n    control: dominanceSpec('W7x@point-control', eps, 1, 0),\n    controlClaim: 'the same bound without its waveform term W̄, 2θ0 sin(Δφ/2) alone',\n  };\n}\n\n/** @internal */\nexport const PENDULUM_POSITION_TRANSLATION: ObservableTranslation = {\n  bridgeId: 'ab-pendulum-linear',\n  observable: 'position',\n  unit: 'rad',\n  fromNorm: 'relative period error, normalized by the value of the reduced model',\n  errorKind: 'upper-bound',\n  definition: 'the pointwise angle error |θ(t) − θ_lin(t)| of the two motions released from rest at θ0',\n  derivation:\n    'an UPPER BOUND B(t) = W̄ + 2θ0 sin(min(Δφ(t), π)/2), Δφ = 2π (t/T0) ε/(1+ε) the exact phase lead and ' +\n    'W̄ = (16/3) q^{3/2}/(1−q) ≈ θ0³/96 the waveform term (q the nome of sin(θ0/2); θ0 is recovered from ε). ' +\n    'From the Fourier series θ = Σ (−1)ⁿ aₙ cos((2n+1)φ): |θ − θ0 cos φ| ≤ 2 Σ_{n≥1} aₙ ≤ W̄, and θ0|cos φ − cos φ_lin| ≤ 2θ0|sin(Δφ/2)|. ' +\n    'B increases with ε, so the domain supremum of ε gives the shortest horizon over the domain',\n  premises: [\n    'both motions are released from rest at the same θ0 at t = 0',\n    'T0 = 2π √(ℓ/g) is the linear period, in the unit t is given in',\n    'the classical Fourier series of the pendulum (all aₙ > 0, Σ (−1)ⁿ aₙ = θ0) — quoted, not derived here; ' +\n      `machine-checked against RK4 at θ0 ∈ {${FOURIER_SERIES_THETA0S.join(', ')}} (W7xs), not between those points`,\n  ],\n  parameters: ['T0'],\n  timeUnit: 'the unit of T0',\n  errorAt: (eps, t, { T0 }) => {\n    if (T0 === undefined || !Number.isFinite(T0) || T0 <= 0 || !Number.isFinite(eps) || eps < 0 || !Number.isFinite(t)) return Number.NaN;\n    return positionBound(eps)(t / T0);\n  },\n  horizonFor: (eps, tolerance, { T0 }) => positionHorizon(eps, tolerance, T0),\n  floor: (eps) => waveformBound(theta0OfPeriodError(eps)),\n  boundary:\n    'inclusive: certified iff t ≤ t*, i.e. B(t) ≤ the tolerance; past t* the BOUND exceeds it and the error may not (undetermined, not inadequate)',\n  notCovered: ['the phase (declared separately: --tolerance=phase:EPS)', 'the angular velocity', 'the amplitude'],\n  witnesses: [\n    {\n      id: 'W7x',\n      kind: 'numeric',\n      test: 'tests/atlas/pendulum-position-translation.test.ts',\n      tolerance: `RK4 |θ − θ_lin| ≤ B(t) at every sample to antiphase, θ0 = ${W7X_THETA0}, by more than the coarse-to-fine difference`,\n    },\n    {\n      id: 'W7xa',\n      kind: 'numeric',\n      test: 'tests/atlas/pendulum-position-translation.test.ts',\n      tolerance: `RK4-projected fundamental coefficient at θ0 = ${W7X_THETA0} within 1e-10 of 8√q/(1+q)`,\n    },\n    {\n      id: 'W7xs',\n      kind: 'numeric',\n      test: 'tests/atlas/pendulum-position-translation.test.ts',\n      tolerance:\n        `the Fourier series (${FOURIER_SERIES_TERMS} terms) within ${FOURIER_SERIES_TOLERANCE} rad of RK4 at every step of one period, ` +\n        `at θ0 ∈ {${FOURIER_SERIES_THETA0S.join(', ')}}`,\n    },\n  ],\n  checks: [dominanceSpec('W7x', W7X_EPS, 1, 1), W7XA, fourierSeriesCheck('W7xs')],\n  pointCheck: positionPointCheck,\n};\n"]}
+{"version":3,"file":"position-translation.js","sourceRoot":"","sources":["../../../src/atlas/oscillators/position-translation.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GA4CG;AAEH,OAAO,EAAE,IAAI,EAAE,MAAM,iCAAiC,CAAC;AAIvD,OAAO,EAAE,qBAAqB,EAAE,MAAM,qBAAqB,CAAC;AAC5D,OAAO,EACL,kBAAkB,EAClB,WAAW,EACX,aAAa,EACb,YAAY,EACZ,OAAO,EACP,mBAAmB,EACnB,MAAM,GACP,MAAM,sBAAsB,CAAC;AAE9B,qGAAqG;AACrG,MAAM,UAAU,aAAa,CAAC,MAAc;IAC1C,MAAM,CAAC,GAAG,YAAY,CAAC,MAAM,CAAC,CAAC;IAC/B,OAAO,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,IAAI,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;AACzC,CAAC;AAED,iFAAiF;AACjF,MAAM,UAAU,sBAAsB,CAAC,MAAc;IACnD,MAAM,CAAC,GAAG,YAAY,CAAC,MAAM,CAAC,CAAC;IAC/B,OAAO,CAAC,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;AACtC,CAAC;AAED,MAAM,IAAI,GAAG,CAAC,GAAW,EAAE,OAAe,EAAU,EAAE,CAAC,CAAC,MAAM,GAAG,OAAO,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC;AAE5F,qGAAqG;AACrG,SAAS,aAAa,CAAC,GAAW,EAAE,KAAK,GAAG,CAAC,EAAE,QAAQ,GAAG,CAAC;IACzD,MAAM,MAAM,GAAG,mBAAmB,CAAC,GAAG,CAAC,CAAC;IACxC,MAAM,KAAK,GAAG,QAAQ,GAAG,aAAa,CAAC,MAAM,CAAC,CAAC;IAC/C,OAAO,CAAC,OAAO,EAAE,EAAE,CAAC,KAAK,GAAG,KAAK,GAAG,CAAC,GAAG,MAAM,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,EAAE,OAAO,CAAC,EAAE,IAAI,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC;AACvG,CAAC;AAED,SAAS,eAAe,CAAC,GAAW,EAAE,SAAiB,EAAE,EAAsB;IAC7E,IAAI,EAAE,KAAK,SAAS,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,EAAE,CAAC,IAAI,EAAE,IAAI,CAAC,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,GAAG,CAAC,IAAI,GAAG,GAAG,CAAC;QAAE,OAAO,MAAM,CAAC,GAAG,CAAC;IAC/G,IAAI,GAAG,KAAK,CAAC;QAAE,OAAO,QAAQ,CAAC;IAC/B,MAAM,MAAM,GAAG,mBAAmB,CAAC,GAAG,CAAC,CAAC;IACxC,MAAM,KAAK,GAAG,aAAa,CAAC,MAAM,CAAC,CAAC;IACpC,IAAI,SAAS,GAAG,KAAK;QAAE,OAAO,MAAM,CAAC,GAAG,CAAC;IACzC,MAAM,CAAC,GAAG,CAAC,SAAS,GAAG,KAAK,CAAC,GAAG,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC;IAC7C,IAAI,CAAC,IAAI,CAAC;QAAE,OAAO,QAAQ,CAAC;IAC5B,OAAO,CAAC,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,GAAG,CAAC,MAAM,GAAG,GAAG,CAAC,CAAC;AAC9D,CAAC;AAED,oEAAoE;AACpE,MAAM,0BAA0B,GAAG,GAAG,CAAC;AACvC,iDAAiD;AACjD,MAAM,kBAAkB,GAAG,GAAG,CAAC;AAE/B;;;;;GAKG;AACH,SAAS,gBAAgB,CAAC,MAAc,EAAE,OAAe,EAAE,cAAsB,EAAE,KAA4B;IAC7G,MAAM,MAAM,GAAG,cAAc,GAAG,kBAAkB,CAAC;IACnD,MAAM,CAAC,GAAG,CAAC,GAAG,cAAc,CAAC;IAC7B,MAAM,QAAQ,GAAa,EAAE,CAAC;IAC9B,MAAM,CAAC,GAAa,EAAE,CAAC;IACvB,IAAI,CAAC,GAAG,MAAM,CAAC;IACf,IAAI,CAAC,GAAG,CAAC,CAAC;IACV,IAAI,CAAC,GAAG,MAAM,CAAC;IACf,IAAI,CAAC,GAAG,CAAC,CAAC;IACV,MAAM,KAAK,GAAG,IAAI,CAAC,KAAK,CAAC,OAAO,GAAG,kBAAkB,CAAC,GAAG,MAAM,CAAC;IAChE,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,KAAK,EAAE,CAAC,EAAE,EAAE,CAAC;QAChC,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,OAAO,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,aAAa,CAAC,CAAC;QACzC,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,OAAO,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,WAAW,CAAC,CAAC;QACvC,IAAI,CAAC,GAAG,MAAM,KAAK,CAAC,EAAE,CAAC;YACrB,QAAQ,CAAC,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;YAC/B,CAAC,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;QACvB,CAAC;IACH,CAAC;IACD,OAAO,EAAE,QAAQ,EAAE,KAAK,EAAE,CAAC,EAAE,CAAC;AAChC,CAAC;AAED,mHAAmH;AACnH,SAAS,eAAe,CAAC,GAAW;IAClC,OAAO,IAAI,CAAC,GAAG,CAAC,0BAA0B,EAAE,CAAC,CAAC,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC;AACrE,CAAC;AAED;;;;;;;;;GASG;AACH,MAAM,UAAU,kBAAkB,CAAC,EAAU,EAAE,GAAW;IACxD,MAAM,MAAM,GAAG,mBAAmB,CAAC,GAAG,CAAC,CAAC;IACxC,MAAM,OAAO,GAAG,eAAe,CAAC,GAAG,CAAC,CAAC;IACrC,MAAM,KAAK,GAAG,aAAa,CAAC,GAAG,CAAC,CAAC;IACjC,MAAM,OAAO,GAAG,IAAI,CAAC,KAAK,CAAC,OAAO,GAAG,kBAAkB,CAAC,CAAC;IACzD,OAAO;QACL,EAAE;QACF,IAAI,EAAE,WAAW;QACjB,QAAQ,EAAE,GAAG,EAAE;YACb,MAAM,QAAQ,GAAa,EAAE,CAAC;YAC9B,MAAM,CAAC,GAAa,EAAE,CAAC;YACvB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,OAAO,EAAE,CAAC,EAAE,EAAE,CAAC;gBAClC,MAAM,CAAC,GAAG,CAAC,GAAG,kBAAkB,CAAC;gBACjC,QAAQ,CAAC,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,kBAAkB,CAAC,MAAM,EAAE,CAAC,CAAC,GAAG,MAAM,GAAG,IAAI,CAAC,GAAG,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;gBACvF,CAAC,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC;YACnB,CAAC;YACD,OAAO,EAAE,QAAQ,EAAE,KAAK,EAAE,CAAC,EAAE,CAAC;QAChC,CAAC;QACD,gBAAgB,EAAE,GAAG;QACrB,cAAc,EAAE,GAAG;KACpB,CAAC;AACJ,CAAC;AAED,SAAS,aAAa,CAAC,EAAU,EAAE,GAAW,EAAE,KAAa,EAAE,QAAgB;IAC7E,MAAM,MAAM,GAAG,mBAAmB,CAAC,GAAG,CAAC,CAAC;IACxC,MAAM,OAAO,GAAG,eAAe,CAAC,GAAG,CAAC,CAAC;IACrC,OAAO;QACL,EAAE;QACF,IAAI,EAAE,WAAW;QACjB,QAAQ,EAAE,CAAC,cAAc,EAAE,EAAE,CAAC,gBAAgB,CAAC,MAAM,EAAE,OAAO,EAAE,cAAc,EAAE,aAAa,CAAC,GAAG,EAAE,KAAK,EAAE,QAAQ,CAAC,CAAC;QACpH,gBAAgB,EAAE,GAAG;QACrB,cAAc,EAAE,GAAG;KACpB,CAAC;AACJ,CAAC;AAED,gHAAgH;AAChH,SAAS,kBAAkB,CAAC,MAAc,EAAE,cAAsB;IAChE,MAAM,CAAC,GAAG,CAAC,GAAG,qBAAqB,CAAC,EAAE,MAAM,EAAE,CAAC,CAAC;IAChD,MAAM,CAAC,GAAG,CAAC,GAAG,cAAc,CAAC;IAC7B,IAAI,CAAC,GAAG,MAAM,CAAC;IACf,IAAI,CAAC,GAAG,CAAC,CAAC;IACV,IAAI,GAAG,GAAG,MAAM,CAAC;IACjB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,cAAc,EAAE,CAAC,EAAE,EAAE,CAAC;QACxC,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,OAAO,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,aAAa,CAAC,CAAC;QACzC,GAAG,IAAI,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,MAAM,GAAG,CAAC,CAAC,GAAG,cAAc,CAAC,CAAC;IACrD,CAAC;IACD,OAAO,CAAC,CAAC,GAAG,GAAG,CAAC,GAAG,cAAc,CAAC;AACpC,CAAC;AAED,4GAA4G;AAC5G,MAAM,CAAC,MAAM,sBAAsB,GAAsB,CAAC,IAAI,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,CAAC,CAAC;AACzF,MAAM,oBAAoB,GAAG,EAAE,CAAC;AAChC,MAAM,wBAAwB,GAAG,IAAI,CAAC;AAEtC,iIAAiI;AACjI,SAAS,mBAAmB,CAAC,CAAS,EAAE,CAAS,EAAE,OAAyB;IAC1E,MAAM,GAAG,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;IACtB,OAAO,CAAC,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC,OAAO,KAAK,MAAM,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,IAAI,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,IAAI,GAAG,CAAC,CAAC,CAAC;AAC3F,CAAC;AAED,wGAAwG;AACxG,MAAM,UAAU,kBAAkB,CAAC,MAAc,EAAE,CAAS,EAAE,OAAO,GAAqB,MAAM;IAC9F,OAAO,mBAAmB,CAAC,YAAY,CAAC,MAAM,CAAC,EAAE,CAAC,EAAE,OAAO,CAAC,CAAC;AAC/D,CAAC;AAED;;;GAGG;AACH,SAAS,aAAa,CAAC,MAAc,EAAE,GAAW,EAAE,OAAyB;IAC3E,MAAM,CAAC,GAAG,YAAY,CAAC,MAAM,CAAC,CAAC;IAC/B,IAAI,CAAC,GAAG,CAAC,CAAC;IACV,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,oBAAoB,EAAE,CAAC,EAAE,EAAE,CAAC;QAC9C,MAAM,CAAC,GAAG,mBAAmB,CAAC,CAAC,EAAE,CAAC,EAAE,OAAO,CAAC,CAAC;QAC7C,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC;IAC5D,CAAC;IACD,OAAO,CAAC,CAAC;AACX,CAAC;AAED,kHAAkH;AAClH,MAAM,CAAC,MAAM,yBAAyB,GAAG,EAAE,UAAU,EAAE,GAAG,EAAE,MAAM,EAAE,GAAG,EAAW,CAAC;AACnF;;;;;;;;GAQG;AACH,MAAM,CAAC,MAAM,8BAA8B,GAAG,KAAK,CAAC;AACpD,mIAAmI;AACnI,MAAM,CAAC,MAAM,0BAA0B,GAAG,EAAE,SAAS,EAAE,KAAK,EAAE,SAAS,EAAE,KAAK,EAAW,CAAC;AAE1F;;;;;;;;;;GAUG;AACH,MAAM,UAAU,2BAA2B,CACzC,OAAO,GAAqB,MAAM,EAClC,IAAI,GAA6D,yBAAyB,EAC1F,SAAS,GAAG,CAAC;IAEb,IAAI,KAAK,GAAG,CAAC,CAAC;IACd,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,IAAI,CAAC,UAAU,EAAE,CAAC,EAAE,EAAE,CAAC;QAC1C,MAAM,MAAM,GAAG,CAAC,GAAG,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,UAAU,CAAC;QAC3C,IAAI,MAAM,GAAG,SAAS;YAAE,SAAS;QACjC,MAAM,CAAC,GAAG,CAAC,GAAG,qBAAqB,CAAC,EAAE,MAAM,EAAE,CAAC,CAAC;QAChD,MAAM,CAAC,GAAG,YAAY,CAAC,MAAM,CAAC,CAAC;QAC/B,MAAM,CAAC,GAAG,KAAK,CAAC,IAAI,CAAC,EAAE,MAAM,EAAE,oBAAoB,EAAE,EAAE,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,mBAAmB,CAAC,CAAC,EAAE,CAAC,EAAE,OAAO,CAAC,CAAC,CAAC;QACrG,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;YACrC,MAAM,GAAG,GAAG,CAAC,MAAM,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,MAAM,CAAC;YACvC,IAAI,CAAC,GAAG,CAAC,CAAC;YACV,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,oBAAoB,EAAE,CAAC,EAAE;gBAAE,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAE,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC;YACjH,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,EAAE,IAAI,CAAC,GAAG,CAAC,kBAAkB,CAAC,MAAM,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;QAC3F,CAAC;IACH,CAAC;IACD,OAAO,KAAK,CAAC;AACf,CAAC;AAED;;;;;;;GAOG;AACH,MAAM,UAAU,wBAAwB,CAAC,MAAc,EAAE,OAAe,EAAE,KAAa;IACrF,IAAI,CAAC,MAAM,CAAC,SAAS,CAAC,OAAO,CAAC,IAAI,OAAO,IAAI,CAAC,GAAG,KAAK,GAAG,CAAC,EAAE,CAAC;QAC3D,MAAM,IAAI,UAAU,CAAC,6BAA6B,OAAO,2BAA2B,KAAK,gBAAgB,CAAC,CAAC;IAC7G,CAAC;IACD,MAAM,CAAC,GAAG,CAAC,GAAG,qBAAqB,CAAC,EAAE,MAAM,EAAE,CAAC,CAAC;IAChD,MAAM,EAAE,EAAE,EAAE,GAAG,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,EAAE,MAAM,EAAE,OAAO,EAAE,EAAE,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,kBAAkB,CAAC,MAAM,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,OAAO,CAAC,CAAC,CAAC,CAAC;IAC9G,OAAO,KAAK,CAAC,IAAI,CAAC,EAAE,MAAM,EAAE,KAAK,EAAE,EAAE,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAE,CAAC,GAAG,OAAO,CAAC,CAAC;AACtG,CAAC;AAED,8GAA8G;AAC9G,SAAS,sBAAsB,CAAC,cAAsB,EAAE,OAAyB;IAC/E,IAAI,KAAK,GAAG,CAAC,CAAC;IACd,KAAK,MAAM,MAAM,IAAI,sBAAsB,EAAE,CAAC;QAC5C,MAAM,CAAC,GAAG,CAAC,CAAC,GAAG,qBAAqB,CAAC,EAAE,MAAM,EAAE,CAAC,CAAC,GAAG,cAAc,CAAC;QACnE,IAAI,CAAC,GAAG,MAAM,CAAC;QACf,IAAI,CAAC,GAAG,CAAC,CAAC;QACV,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,EAAE,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,aAAa,CAAC,MAAM,EAAE,CAAC,EAAE,OAAO,CAAC,CAAC,CAAC,CAAC;QACzE,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,cAAc,EAAE,CAAC,EAAE,EAAE,CAAC;YACzC,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,OAAO,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,aAAa,CAAC,CAAC;YACzC,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,EAAE,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,aAAa,CAAC,MAAM,EAAE,CAAC,MAAM,GAAG,CAAC,CAAC,GAAG,cAAc,EAAE,OAAO,CAAC,CAAC,CAAC,CAAC;QACvG,CAAC;IACH,CAAC;IACD,OAAO,KAAK,CAAC;AACf,CAAC;AAED;;;;GAIG;AACH,MAAM,UAAU,kBAAkB,CAAC,EAAU,EAAE,OAAO,GAAqB,MAAM;IAC/E,OAAO;QACL,EAAE;QACF,QAAQ,EAAE,CAAC,cAAc,EAAE,EAAE,CAAC,sBAAsB,CAAC,cAAc,EAAE,OAAO,CAAC;QAC7E,MAAM,EAAE,CAAC;QACT,gBAAgB,EAAE,GAAG;QACrB,cAAc,EAAE,GAAG;QACnB,SAAS,EAAE,wBAAwB;KACpC,CAAC;AACJ,CAAC;AAED,MAAM,UAAU,GAAG,GAAG,CAAC;AACvB,MAAM,OAAO,GAAG,qBAAqB,CAAC,EAAE,MAAM,EAAE,UAAU,EAAE,CAAC,CAAC;AAE9D,MAAM,IAAI,GAAuB;IAC/B,EAAE,EAAE,MAAM;IACV,QAAQ,EAAE,CAAC,cAAc,EAAE,EAAE,CAAC,kBAAkB,CAAC,UAAU,EAAE,cAAc,CAAC;IAC5E,MAAM,EAAE,sBAAsB,CAAC,UAAU,CAAC;IAC1C,gBAAgB,EAAE,GAAG;IACrB,cAAc,EAAE,GAAG;IACnB,SAAS,EAAE,KAAK;CACjB,CAAC;AAEF,SAAS,kBAAkB,CAAC,GAAW;IACrC,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,GAAG,CAAC,IAAI,GAAG,IAAI,CAAC;QAAE,OAAO,EAAE,WAAW,EAAE,iBAAiB,GAAG,+CAA+C,EAAE,CAAC;IACnI,MAAM,MAAM,GAAG,mBAAmB,CAAC,GAAG,CAAC,CAAC;IACxC,MAAM,OAAO,GAAG,eAAe,CAAC,GAAG,CAAC,CAAC;IACrC,MAAM,MAAM,GAAG,OAAO,KAAK,0BAA0B,CAAC,CAAC,CAAC,eAAe,0BAA0B,wBAAwB,CAAC,CAAC,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,CAAC,MAAM,CAAC,CAAC,CAAC,yBAAyB,CAAC;IACjL,OAAO;QACL,KAAK,EACH,gDAAgD,MAAM,CAAC,MAAM,CAAC,WAAW,CAAC,CAAC,CAAC,CAAC,iCAAiC,OAAO,MAAM,MAAM,IAAI;YACrI,gGAAgG;QAClG,KAAK,EAAE,aAAa,CAAC,WAAW,EAAE,GAAG,EAAE,CAAC,EAAE,CAAC,CAAC;QAC5C,OAAO,EAAE,aAAa,CAAC,mBAAmB,EAAE,GAAG,EAAE,CAAC,EAAE,CAAC,CAAC;QACtD,YAAY,EAAE,kEAAkE;KACjF,CAAC;AACJ,CAAC;AAED,gBAAgB;AAChB,MAAM,CAAC,MAAM,6BAA6B,GAA0B;IAClE,QAAQ,EAAE,oBAAoB;IAC9B,UAAU,EAAE,UAAU;IACtB,IAAI,EAAE,KAAK;IACX,QAAQ,EAAE,qEAAqE;IAC/E,SAAS,EAAE,aAAa;IACxB,UAAU,EAAE,yFAAyF;IACrG,UAAU,EACR,uGAAuG;QACvG,0GAA0G;QAC1G,uIAAuI;QACvI,4FAA4F;IAC9F,QAAQ,EAAE;QACR,6DAA6D;QAC7D,gEAAgE;QAChE,yGAAyG;YACvG,oFAAoF,yBAAyB,CAAC,UAAU,0BAA0B;YAClJ,OAAO,yBAAyB,CAAC,MAAM,kBAAkB,8BAA8B,oBAAoB,0BAA0B,CAAC,SAAS,OAAO;YACtJ,cAAc,0BAA0B,CAAC,SAAS,mEAAmE,sBAAsB,CAAC,IAAI,CAAC,IAAI,CAAC,UAAU;KACnK;IACD,UAAU,EAAE,CAAC,IAAI,CAAC;IAClB,QAAQ,EAAE,gBAAgB;IAC1B,OAAO,EAAE,CAAC,GAAG,EAAE,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE;QAC1B,IAAI,EAAE,KAAK,SAAS,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,EAAE,CAAC,IAAI,EAAE,IAAI,CAAC,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,GAAG,CAAC,IAAI,GAAG,GAAG,CAAC,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC;YAAE,OAAO,MAAM,CAAC,GAAG,CAAC;QACtI,OAAO,aAAa,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC;IACpC,CAAC;IACD,UAAU,EAAE,CAAC,GAAG,EAAE,SAAS,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,eAAe,CAAC,GAAG,EAAE,SAAS,EAAE,EAAE,CAAC;IAC3E,KAAK,EAAE,CAAC,GAAG,EAAE,EAAE,CAAC,aAAa,CAAC,mBAAmB,CAAC,GAAG,CAAC,CAAC;IACvD,QAAQ,EACN,+IAA+I;IACjJ,UAAU,EAAE,CAAC,wDAAwD,EAAE,sBAAsB,EAAE,eAAe,CAAC;IAC/G,SAAS,EAAE;QACT;YACE,EAAE,EAAE,KAAK;YACT,IAAI,EAAE,SAAS;YACf,IAAI,EAAE,mDAAmD;YACzD,SAAS,EAAE,6DAA6D,UAAU,8CAA8C;SACjI;QACD;YACE,EAAE,EAAE,MAAM;YACV,IAAI,EAAE,SAAS;YACf,IAAI,EAAE,mDAAmD;YACzD,SAAS,EAAE,iDAAiD,UAAU,4BAA4B;SACnG;QACD;YACE,EAAE,EAAE,MAAM;YACV,IAAI,EAAE,SAAS;YACf,IAAI,EAAE,mDAAmD;YACzD,SAAS,EACP,uBAAuB,oBAAoB,kBAAkB,wBAAwB,2CAA2C;gBAChI,YAAY,sBAAsB,CAAC,IAAI,CAAC,IAAI,CAAC,GAAG;SACnD;KACF;IACD,MAAM,EAAE,CAAC,aAAa,CAAC,KAAK,EAAE,OAAO,EAAE,CAAC,EAAE,CAAC,CAAC,EAAE,IAAI,EAAE,kBAAkB,CAAC,MAAM,CAAC,CAAC;IAC/E,UAAU,EAAE,kBAAkB;CAC/B,CAAC","sourcesContent":["/**\n * `ab-pendulum-linear`'s declared translation from its bound's quantity, the\n * relative period error ε, into an UPPER BOUND on the pointwise position error\n * |θ(t) − θ_lin(t)|, in rad.\n *\n * Both motions start from rest at θ0. The pendulum is θ(t) = F(φ), φ = 2πt/T,\n * with the classical Fourier series of the elliptic solution\n *\n *   F(φ) = Σ_{n≥0} (−1)ⁿ aₙ cos((2n+1)φ),  aₙ = 8 q^{n+½} / ((2n+1)(1 + q^{2n+1})),\n *\n * q the nome of k = sin(θ0/2); F(0) = θ0 = Σ (−1)ⁿ aₙ. The linear motion is\n * θ0 cos φ_lin, φ_lin = 2πt/T0. Splitting at θ0 cos φ,\n *\n *   |θ − θ_lin| ≤ |F(φ) − θ0 cos φ| + θ0 |cos φ − cos φ_lin|\n *              ≤ 2 Σ_{n≥1} aₙ + 2θ0 |sin(Δφ/2)|,\n *\n * because F(φ) − θ0 cos φ = Σ_{n≥1} (−1)ⁿ aₙ (cos((2n+1)φ) − cos φ). With\n * aₙ ≤ (8/3) q^{n+½} for n ≥ 1 the first term is at most\n * W̄ = (16/3) q^{3/2}/(1 − q) ≈ θ0³/96, and Δφ = 2π(t/T0) ε/(1+ε) is the\n * exact phase lead of `phase-translation.ts`. So\n *\n *   B(t) = W̄ + 2θ0 sin(min(Δφ(t), π)/2)\n *\n * bounds the error at every t; the clamp at π keeps B non-decreasing in t and\n * is still a bound, since |cos a − cos b| ≤ 2. It is rigorous GIVEN the\n * Fourier series, which is quoted, not derived here. The exact check\n * (`fourierSeriesExactDeviation`, run by `tests/atlas/pendulum-exact.test.ts`;\n * not a registered witness, because it has no resolution for a runner to\n * refine) compares the series with the EXACT solution θ = 2 asin(k sn(K − ω0t, k²))\n * (`exactPendulumAngle`) over a dense grid of amplitudes and phases at\n * 1e-12 rad (1e-13 above θ0 = 0.002; the floor below that is the nome's, see\n * `FOURIER_SERIES_EXACT_TOLERANCE`), and W7xs against RK4 over one period at\n * six amplitudes at 1e-9;\n * both are checks at grid points, not a derivation. The coefficients also\n * have a second route, the FFT of one exact period (`fourierCoefficientsByFft`).\n *\n * B is an upper bound, not the error: past its horizon the BOUND exceeds the\n * tolerance, and the error may not. It is near-sharp where Δφ → π (the\n * motions are in antiphase and the error reaches 2θ0), and the waveform term\n * cannot be dropped: without W̄, RK4 finds the error above 2θ0 sin(Δφ/2) by\n * about 3% (witness control).\n *\n * @module atlas/oscillators/position-translation\n * @internal\n */\n\nimport { rfft } from '@danielsimonjr/mathts-functions';\nimport type { ObservableTranslation } from '../translation.js';\nimport type { DominanceWitnessSpec } from '../witness-dominance.js';\nimport type { NumericWitnessSpec } from '../witness-numeric.js';\nimport { pendulumPeriodErrorAt } from './bridges-limits.js';\nimport {\n  exactPendulumAngle,\n  linearAccel,\n  pendulumAccel,\n  pendulumNome,\n  rk4Step,\n  theta0OfPeriodError,\n  TWO_PI,\n} from './pendulum-motion.js';\n\n/** W̄ = (16/3) q^{3/2}/(1 − q): the bound on the waveform distortion |F(φ) − θ0 cos φ|. @internal */\nexport function waveformBound(theta0: number): number {\n  const q = pendulumNome(theta0);\n  return ((16 / 3) * q ** 1.5) / (1 - q);\n}\n\n/** The pendulum's fundamental Fourier coefficient a₀ = 8√q/(1 + q). @internal */\nexport function fundamentalCoefficient(theta0: number): number {\n  const q = pendulumNome(theta0);\n  return (8 * Math.sqrt(q)) / (1 + q);\n}\n\nconst lead = (eps: number, tOverT0: number): number => (TWO_PI * tOverT0 * eps) / (1 + eps);\n\n/** B(t/T0) at θ0 recovered from ε; `drift` and `waveform` scale its two terms (1, 1 is B itself). */\nfunction positionBound(eps: number, drift = 1, waveform = 1): (tOverT0: number) => number {\n  const theta0 = theta0OfPeriodError(eps);\n  const floor = waveform * waveformBound(theta0);\n  return (tOverT0) => floor + drift * 2 * theta0 * Math.sin(Math.min(lead(eps, tOverT0), Math.PI) / 2);\n}\n\nfunction positionHorizon(eps: number, tolerance: number, T0: number | undefined): number {\n  if (T0 === undefined || !Number.isFinite(T0) || T0 <= 0 || !Number.isFinite(eps) || eps < 0) return Number.NaN;\n  if (eps === 0) return Infinity;\n  const theta0 = theta0OfPeriodError(eps);\n  const floor = waveformBound(theta0);\n  if (tolerance < floor) return Number.NaN;\n  const r = (tolerance - floor) / (2 * theta0);\n  if (r >= 1) return Infinity;\n  return (2 * Math.asin(r) * T0 * (1 + eps)) / (TWO_PI * eps);\n}\n\n/** The point witness integrates at most this many periods of T0. */\nconst POSITION_POINT_MAX_PERIODS = 512;\n/** Samples per T0 both resolutions report at. */\nconst SAMPLES_PER_PERIOD = 200;\n\n/**\n * |θ(t) − θ_lin(t)| from RK4 of both motions, and B(t) at the same samples,\n * over `periods` of T0; `stepsPerPeriod` must be a multiple of 200. The\n * exact route is {@link exactDominanceSpec}; this one is what W7x registers,\n * so the integrator stays part of what the witness exercises.\n */\nfunction dominanceSamples(theta0: number, periods: number, stepsPerPeriod: number, bound: (t: number) => number) {\n  const stride = stepsPerPeriod / SAMPLES_PER_PERIOD;\n  const h = 1 / stepsPerPeriod;\n  const measured: number[] = [];\n  const b: number[] = [];\n  let x = theta0;\n  let v = 0;\n  let y = theta0;\n  let u = 0;\n  const steps = Math.floor(periods * SAMPLES_PER_PERIOD) * stride;\n  for (let i = 1; i <= steps; i++) {\n    [x, v] = rk4Step(x, v, h, pendulumAccel);\n    [y, u] = rk4Step(y, u, h, linearAccel);\n    if (i % stride === 0) {\n      measured.push(Math.abs(x - y));\n      b.push(bound(i * h));\n    }\n  }\n  return { measured, bound: b };\n}\n\n/** The window a dominance witness covers: up to antiphase (Δφ = π), past which B ≥ 2θ0 holds trivially, capped. */\nfunction dominanceWindow(eps: number): number {\n  return Math.min(POSITION_POINT_MAX_PERIODS, (1 + eps) / (2 * eps));\n}\n\n/**\n * The exact route to W7x's measured value: |θ_exact(t) − θ0 cos(2πt)| from the\n * closed forms, at the same 200 samples per T0 as {@link dominanceSamples},\n * against B(t) at θ0 recovered from `eps`. The resolution is ignored, so the\n * runner's coarse-to-fine uncertainty u is exactly 0 here and its rule\n * `bound − measured > u` is `bound > measured`: u is a property of the RK4\n * route, not of the claim. The registered W7x keeps the RK4 route.\n *\n * @internal\n */\nexport function exactDominanceSpec(id: string, eps: number): DominanceWitnessSpec {\n  const theta0 = theta0OfPeriodError(eps);\n  const periods = dominanceWindow(eps);\n  const bound = positionBound(eps);\n  const samples = Math.floor(periods * SAMPLES_PER_PERIOD);\n  return {\n    id,\n    kind: 'dominance',\n    evaluate: () => {\n      const measured: number[] = [];\n      const b: number[] = [];\n      for (let i = 1; i <= samples; i++) {\n        const t = i / SAMPLES_PER_PERIOD;\n        measured.push(Math.abs(exactPendulumAngle(theta0, t) - theta0 * Math.cos(TWO_PI * t)));\n        b.push(bound(t));\n      }\n      return { measured, bound: b };\n    },\n    coarseResolution: 200,\n    fineResolution: 800,\n  };\n}\n\nfunction dominanceSpec(id: string, eps: number, drift: number, waveform: number): DominanceWitnessSpec {\n  const theta0 = theta0OfPeriodError(eps);\n  const periods = dominanceWindow(eps);\n  return {\n    id,\n    kind: 'dominance',\n    evaluate: (stepsPerPeriod) => dominanceSamples(theta0, periods, stepsPerPeriod, positionBound(eps, drift, waveform)),\n    coarseResolution: 200,\n    fineResolution: 800,\n  };\n}\n\n/** a₀ measured by projecting one RK4 pendulum period onto cos φ (trapezoid, exact for a periodic integrand). */\nfunction measureFundamental(theta0: number, stepsPerPeriod: number): number {\n  const T = 1 + pendulumPeriodErrorAt({ theta0 });\n  const h = T / stepsPerPeriod;\n  let x = theta0;\n  let v = 0;\n  let sum = theta0;\n  for (let i = 1; i < stepsPerPeriod; i++) {\n    [x, v] = rk4Step(x, v, h, pendulumAccel);\n    sum += x * Math.cos((TWO_PI * i) / stepsPerPeriod);\n  }\n  return (2 * sum) / stepsPerPeriod;\n}\n\n/** The θ0 at which W7xs checks the Fourier series: a grid across the bridge's domain θ0 ≤ 0.5. @internal */\nexport const FOURIER_SERIES_THETA0S: readonly number[] = [0.01, 0.1, 0.2, 0.3, 0.4, 0.5];\nconst FOURIER_SERIES_TERMS = 20;\nconst FOURIER_SERIES_TOLERANCE = 1e-9;\n\n/** aₙ = 8 q^{n+½}/((2n+1)(1 + q^{2n+1})) from the nome; `'minus'` is the (1 − q^{2n+1}) transcription error the controls use. */\nfunction coefficientFromNome(q: number, n: number, variant: 'plus' | 'minus'): number {\n  const odd = 2 * n + 1;\n  return (8 * q ** (n + 0.5)) / (odd * (variant === 'plus' ? 1 + q ** odd : 1 - q ** odd));\n}\n\n/** The pendulum's Fourier coefficient aₙ at amplitude θ0 (the closed form, from the nome). @internal */\nexport function fourierCoefficient(theta0: number, n: number, variant: 'plus' | 'minus' = 'plus'): number {\n  return coefficientFromNome(pendulumNome(theta0), n, variant);\n}\n\n/**\n * F(φ) = Σ_{n<20} (−1)ⁿ aₙ cos((2n+1)φ). `'minus'` puts (1 − q^{2n+1}) in the\n * denominator of aₙ, a transcription error the control uses.\n */\nfunction fourierSeries(theta0: number, phi: number, variant: 'plus' | 'minus'): number {\n  const q = pendulumNome(theta0);\n  let s = 0;\n  for (let n = 0; n < FOURIER_SERIES_TERMS; n++) {\n    const a = coefficientFromNome(q, n, variant);\n    s += (n % 2 === 0 ? a : -a) * Math.cos((2 * n + 1) * phi);\n  }\n  return s;\n}\n\n/** The grid the exact check runs on: amplitudes θ0 = 0.5 j/500, j = 1..500, and phases φ = 2π i/256. @internal */\nexport const FOURIER_SERIES_EXACT_GRID = { amplitudes: 500, phases: 256 } as const;\n/**\n * The exact check's tolerance over the whole grid, in rad. The measured deviation is\n * 2.3e-13, all of it at the smallest amplitude θ0 = 0.001 and on the SERIES\n * side: the nome's K′ = K(cos²(θ0/2)) is read from a complementary parameter\n * that rounds at 1e-16 absolute, 4e-10 relative at m = 2.5e-7, so a₀ is low\n * by 2.3e-10 of θ0 there. The exact solution agrees with an independent AGM\n * sn to the bit at that point.\n * @internal\n */\nexport const FOURIER_SERIES_EXACT_TOLERANCE = 1e-12;\n/** Above θ0 = 0.002 the same check holds at 1e-13 rad (measured 3.8e-14 at θ0 = 0.002, below 1.5e-14 from θ0 = 0.01). @internal */\nexport const FOURIER_SERIES_EXACT_TIGHT = { theta0Min: 0.002, tolerance: 1e-13 } as const;\n\n/**\n * The exact check, or its control: max over the grid of |θ_exact(t) − F(φ)|, t = φT/2π,\n * with θ_exact the elliptic solution (`exactPendulumAngle`) and F the 20-term\n * series, over the amplitudes at or above `theta0Min`. A closed form against\n * a closed form, so no resolution and no convergence ratio: the number is the\n * deviation itself, held to {@link FOURIER_SERIES_EXACT_TOLERANCE} over the\n * whole grid and to {@link FOURIER_SERIES_EXACT_TIGHT} above its floor by the\n * test.\n *\n * @internal\n */\nexport function fourierSeriesExactDeviation(\n  variant: 'plus' | 'minus' = 'plus',\n  grid: { readonly amplitudes: number; readonly phases: number } = FOURIER_SERIES_EXACT_GRID,\n  theta0Min = 0,\n): number {\n  let worst = 0;\n  for (let j = 1; j <= grid.amplitudes; j++) {\n    const theta0 = (0.5 * j) / grid.amplitudes;\n    if (theta0 < theta0Min) continue;\n    const T = 1 + pendulumPeriodErrorAt({ theta0 });\n    const q = pendulumNome(theta0);\n    const a = Array.from({ length: FOURIER_SERIES_TERMS }, (_, n) => coefficientFromNome(q, n, variant));\n    for (let i = 0; i < grid.phases; i++) {\n      const phi = (TWO_PI * i) / grid.phases;\n      let s = 0;\n      for (let n = 0; n < FOURIER_SERIES_TERMS; n++) s += (n % 2 === 0 ? a[n]! : -a[n]!) * Math.cos((2 * n + 1) * phi);\n      worst = Math.max(worst, Math.abs(exactPendulumAngle(theta0, (i * T) / grid.phases) - s));\n    }\n  }\n  return worst;\n}\n\n/**\n * The first `count` coefficients aₙ from the FFT of one exact period sampled\n * at `samples` points, MathTS `rfft`: aₙ = (−1)ⁿ 2 re[2n+1]/N. A second\n * route to the series, independent of the nome. `samples` must exceed\n * 2·count + 1.\n *\n * @internal\n */\nexport function fourierCoefficientsByFft(theta0: number, samples: number, count: number): number[] {\n  if (!Number.isInteger(samples) || samples <= 2 * count + 1) {\n    throw new RangeError(`fourierCoefficientsByFft: ${samples} samples cannot resolve ${count} odd harmonics`);\n  }\n  const T = 1 + pendulumPeriodErrorAt({ theta0 });\n  const { re } = rfft(Array.from({ length: samples }, (_, i) => exactPendulumAngle(theta0, (i * T) / samples)));\n  return Array.from({ length: count }, (_, n) => ((n % 2 === 0 ? 2 : -2) * re[2 * n + 1]!) / samples);\n}\n\n/** max over the grid, and over one period T = T0(1+ε) sampled at every RK4 step, of |θ_RK4(t) − F(2πt/T)|. */\nfunction fourierSeriesDeviation(stepsPerPeriod: number, variant: 'plus' | 'minus'): number {\n  let worst = 0;\n  for (const theta0 of FOURIER_SERIES_THETA0S) {\n    const h = (1 + pendulumPeriodErrorAt({ theta0 })) / stepsPerPeriod;\n    let x = theta0;\n    let v = 0;\n    worst = Math.max(worst, Math.abs(x - fourierSeries(theta0, 0, variant)));\n    for (let i = 1; i <= stepsPerPeriod; i++) {\n      [x, v] = rk4Step(x, v, h, pendulumAccel);\n      worst = Math.max(worst, Math.abs(x - fourierSeries(theta0, (TWO_PI * i) / stepsPerPeriod, variant)));\n    }\n  }\n  return worst;\n}\n\n/**\n * W7xs, or its control: the Fourier series against RK4 of the pendulum, at 200\n * and 800 steps per period, within 1e-9 rad at every step and every grid θ0.\n * @internal\n */\nexport function fourierSeriesCheck(id: string, variant: 'plus' | 'minus' = 'plus'): NumericWitnessSpec {\n  return {\n    id,\n    evaluate: (stepsPerPeriod) => fourierSeriesDeviation(stepsPerPeriod, variant),\n    target: 0,\n    coarseResolution: 200,\n    fineResolution: 800,\n    tolerance: FOURIER_SERIES_TOLERANCE,\n  };\n}\n\nconst W7X_THETA0 = 0.3;\nconst W7X_EPS = pendulumPeriodErrorAt({ theta0: W7X_THETA0 });\n\nconst W7XA: NumericWitnessSpec = {\n  id: 'W7xa',\n  evaluate: (stepsPerPeriod) => measureFundamental(W7X_THETA0, stepsPerPeriod),\n  target: fundamentalCoefficient(W7X_THETA0),\n  coarseResolution: 200,\n  fineResolution: 800,\n  tolerance: 1e-10,\n};\n\nfunction positionPointCheck(eps: number): ReturnType<ObservableTranslation['pointCheck']> {\n  if (!Number.isFinite(eps) || eps <= 0) return { unavailable: `θ0 = 0 at ε = ${eps}: both motions are at rest and the error is 0` };\n  const theta0 = theta0OfPeriodError(eps);\n  const periods = dominanceWindow(eps);\n  const capped = periods === POSITION_POINT_MAX_PERIODS ? `, capped at ${POSITION_POINT_MAX_PERIODS} T0 (antiphase is at ${(1 + eps) / (2 * eps)} T0)` : ' (to antiphase, Δφ = π)';\n  return {\n    claim:\n      `RK4 |θ − θ_lin| ≤ B(t) at every sample, θ0 = ${Number(theta0.toPrecision(8))} (recovered from ε), over t ≤ ${periods} T0${capped}; ` +\n      `200 and 800 steps per T0, sampled 200 per T0; the coarse-to-fine difference is the uncertainty`,\n    check: dominanceSpec('W7x@point', eps, 1, 1),\n    control: dominanceSpec('W7x@point-control', eps, 1, 0),\n    controlClaim: 'the same bound without its waveform term W̄, 2θ0 sin(Δφ/2) alone',\n  };\n}\n\n/** @internal */\nexport const PENDULUM_POSITION_TRANSLATION: ObservableTranslation = {\n  bridgeId: 'ab-pendulum-linear',\n  observable: 'position',\n  unit: 'rad',\n  fromNorm: 'relative period error, normalized by the value of the reduced model',\n  errorKind: 'upper-bound',\n  definition: 'the pointwise angle error |θ(t) − θ_lin(t)| of the two motions released from rest at θ0',\n  derivation:\n    'an UPPER BOUND B(t) = W̄ + 2θ0 sin(min(Δφ(t), π)/2), Δφ = 2π (t/T0) ε/(1+ε) the exact phase lead and ' +\n    'W̄ = (16/3) q^{3/2}/(1−q) ≈ θ0³/96 the waveform term (q the nome of sin(θ0/2); θ0 is recovered from ε). ' +\n    'From the Fourier series θ = Σ (−1)ⁿ aₙ cos((2n+1)φ): |θ − θ0 cos φ| ≤ 2 Σ_{n≥1} aₙ ≤ W̄, and θ0|cos φ − cos φ_lin| ≤ 2θ0|sin(Δφ/2)|. ' +\n    'B increases with ε, so the domain supremum of ε gives the shortest horizon over the domain',\n  premises: [\n    'both motions are released from rest at the same θ0 at t = 0',\n    'T0 = 2π √(ℓ/g) is the linear period, in the unit t is given in',\n    'the classical Fourier series of the pendulum (all aₙ > 0, Σ (−1)ⁿ aₙ = θ0) — quoted, not derived here; ' +\n      `machine-checked against the elliptic solution θ = 2 asin(k sn(K − ω0t, k²)) over ${FOURIER_SERIES_EXACT_GRID.amplitudes} amplitudes in (0, 0.5] ` +\n      `and ${FOURIER_SERIES_EXACT_GRID.phases} phases within ${FOURIER_SERIES_EXACT_TOLERANCE} rad, and within ${FOURIER_SERIES_EXACT_TIGHT.tolerance} rad ` +\n      `above θ0 = ${FOURIER_SERIES_EXACT_TIGHT.theta0Min} (tests/atlas/pendulum-exact.test.ts), and against RK4 at θ0 ∈ {${FOURIER_SERIES_THETA0S.join(', ')}} (W7xs)`,\n  ],\n  parameters: ['T0'],\n  timeUnit: 'the unit of T0',\n  errorAt: (eps, t, { T0 }) => {\n    if (T0 === undefined || !Number.isFinite(T0) || T0 <= 0 || !Number.isFinite(eps) || eps < 0 || !Number.isFinite(t)) return Number.NaN;\n    return positionBound(eps)(t / T0);\n  },\n  horizonFor: (eps, tolerance, { T0 }) => positionHorizon(eps, tolerance, T0),\n  floor: (eps) => waveformBound(theta0OfPeriodError(eps)),\n  boundary:\n    'inclusive: certified iff t ≤ t*, i.e. B(t) ≤ the tolerance; past t* the BOUND exceeds it and the error may not (undetermined, not inadequate)',\n  notCovered: ['the phase (declared separately: --tolerance=phase:EPS)', 'the angular velocity', 'the amplitude'],\n  witnesses: [\n    {\n      id: 'W7x',\n      kind: 'numeric',\n      test: 'tests/atlas/pendulum-position-translation.test.ts',\n      tolerance: `RK4 |θ − θ_lin| ≤ B(t) at every sample to antiphase, θ0 = ${W7X_THETA0}, by more than the coarse-to-fine difference`,\n    },\n    {\n      id: 'W7xa',\n      kind: 'numeric',\n      test: 'tests/atlas/pendulum-position-translation.test.ts',\n      tolerance: `RK4-projected fundamental coefficient at θ0 = ${W7X_THETA0} within 1e-10 of 8√q/(1+q)`,\n    },\n    {\n      id: 'W7xs',\n      kind: 'numeric',\n      test: 'tests/atlas/pendulum-position-translation.test.ts',\n      tolerance:\n        `the Fourier series (${FOURIER_SERIES_TERMS} terms) within ${FOURIER_SERIES_TOLERANCE} rad of RK4 at every step of one period, ` +\n        `at θ0 ∈ {${FOURIER_SERIES_THETA0S.join(', ')}}`,\n    },\n  ],\n  checks: [dominanceSpec('W7x', W7X_EPS, 1, 1), W7XA, fourierSeriesCheck('W7xs')],\n  pointCheck: positionPointCheck,\n};\n"]}

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

@@ -16,6 +16,4 @@
 import type { AtlasRejection } from '../types.js';
 /** The rejected bridges of the oscillator pilot. @internal */
 export declare const ATLAS_REJECTIONS: readonly AtlasRejection[];
-/** Look up a rejection by id, throwing on an unknown one. @internal */
-export declare function getAtlasRejection(id: string): AtlasRejection;
 //# sourceMappingURL=rejections.d.ts.map

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

@@ -1 +1 @@
-{"version":3,"file":"rejections.d.ts","sourceRoot":"","sources":["../../../src/atlas/oscillators/rejections.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;GAcG;AAEH,OAAO,KAAK,EAAE,cAAc,EAAE,MAAM,aAAa,CAAC;AAElD,8DAA8D;AAC9D,eAAO,MAAM,gBAAgB,EAAE,SAAS,cAAc,EA8BrD,CAAC;AAEF,uEAAuE;AACvE,wBAAgB,iBAAiB,CAAC,EAAE,EAAE,MAAM,GAAG,cAAc,CAI5D"}
+{"version":3,"file":"rejections.d.ts","sourceRoot":"","sources":["../../../src/atlas/oscillators/rejections.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;GAcG;AAEH,OAAO,KAAK,EAAE,cAAc,EAAE,MAAM,aAAa,CAAC;AAElD,8DAA8D;AAC9D,eAAO,MAAM,gBAAgB,EAAE,SAAS,cAAc,EA8BrD,CAAC"}

package/dist/atlas/oscillators/rejections.js

@@ -44,11 +44,4 @@ export const ATLAS_REJECTIONS = [
         ],
     },
 ];
-/** Look up a rejection by id, throwing on an unknown one. @internal */
-export function getAtlasRejection(id) {
-    const found = ATLAS_REJECTIONS.find((r) => r.id === id);
-    if (found === undefined)
-        throw new Error(`unknown atlas rejection id: ${id}`);
-    return found;
-}
 //# sourceMappingURL=rejections.js.map

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

@@ -1 +1 @@
-{"version":3,"file":"rejections.js","sourceRoot":"","sources":["../../../src/atlas/oscillators/rejections.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;GAcG;AAIH,8DAA8D;AAC9D,MAAM,CAAC,MAAM,gBAAgB,GAA8B;IACzD;QACE,EAAE,EAAE,oBAAoB;QACxB,OAAO,EAAE,mBAAmB;QAC5B,QAAQ,EAAE,CAAC,oBAAoB,CAAC;QAChC,UAAU,EAAE,UAAU;QACtB,MAAM,EACJ,6EAA6E;YAC7E,6EAA6E;YAC7E,6EAA6E;YAC7E,2EAA2E;YAC3E,+EAA+E;YAC/E,gFAAgF;YAChF,4BAA4B;QAC9B,cAAc,EAAE,WAAW;QAC3B,SAAS,EAAE;YACT;gBACE,EAAE,EAAE,IAAI;gBACR,IAAI,EAAE,UAAU;gBAChB,IAAI,EAAE,wCAAwC;gBAC9C,SAAS,EAAE,uEAAuE;aACnF;YACD;gBACE,EAAE,EAAE,KAAK;gBACT,IAAI,EAAE,SAAS;gBACf,IAAI,EAAE,wCAAwC;gBAC9C,SAAS,EAAE,8GAA8G;aAC1H;SACF;KACF;CACF,CAAC;AAEF,uEAAuE;AACvE,MAAM,UAAU,iBAAiB,CAAC,EAAU;IAC1C,MAAM,KAAK,GAAG,gBAAgB,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,EAAE,KAAK,EAAE,CAAC,CAAC;IACxD,IAAI,KAAK,KAAK,SAAS;QAAE,MAAM,IAAI,KAAK,CAAC,+BAA+B,EAAE,EAAE,CAAC,CAAC;IAC9E,OAAO,KAAK,CAAC;AACf,CAAC","sourcesContent":["/**\n * The Phase 0 REJECTION record: a claimed bridge the atlas records as refuted,\n * together with the witness that refutes it.\n *\n * `ax-cubic-spring-lc` claims the cubic (Duffing-type) spring\n * `m x″ + k x + β x³ = 0` is exactly equivalent to the LC circuit\n * `L q″ + q/C = 0`. It is not. The cubic spring's variable set\n * {m, k, β, x0} carries exactly ONE dimensionless group, `β x0² / k`, which\n * measures the nonlinearity and has no counterpart on the LC side: the LC\n * variable set is dimensionally independent and admits no group at all. A\n * surviving group with no image cannot be matched by any change of variables,\n * so no exact equivalence exists.\n *\n * @module atlas/oscillators/rejections\n */\n\nimport type { AtlasRejection } from '../types.js';\n\n/** The rejected bridges of the oscillator pilot. @internal */\nexport const ATLAS_REJECTIONS: readonly AtlasRejection[] = [\n  {\n    id: 'ax-cubic-spring-lc',\n    claimed: 'exact-equivalence',\n    premises: ['model-cubic-spring'],\n    conclusion: 'model-lc',\n    reason:\n      'The cubic spring retains one dimensionless group, β x0² / k, that measures ' +\n      'its nonlinearity, and the LC variable set has none (W3). A surviving group ' +\n      'is necessary evidence, not sufficient. The decisive fact is isochrony: the ' +\n      'LC period is the same at every amplitude, the cubic spring period is not ' +\n      '(W3b), and a change of variables that rescales time by a constant multiplies ' +\n      'every period by the same factor. So no such change of variables makes the two ' +\n      'models exactly equivalent.',\n    survivingGroup: 'β x0² / k',\n    witnesses: [\n      {\n        id: 'W3',\n        kind: 'symbolic',\n        test: 'tests/atlas/oscillators-coarse.test.ts',\n        tolerance: 'exact integer π-group exponents; LC verdict dimensionally-independent',\n      },\n      {\n        id: 'W3b',\n        kind: 'numeric',\n        test: 'tests/atlas/oscillators-coarse.test.ts',\n        tolerance: 'linear period 2π within 1e-6 at A = 0.5, 1, 2; cubic-spring periods spread by more than 10% at β x0²/k = 0.1',\n      },\n    ],\n  },\n];\n\n/** Look up a rejection by id, throwing on an unknown one. @internal */\nexport function getAtlasRejection(id: string): AtlasRejection {\n  const found = ATLAS_REJECTIONS.find((r) => r.id === id);\n  if (found === undefined) throw new Error(`unknown atlas rejection id: ${id}`);\n  return found;\n}\n"]}
+{"version":3,"file":"rejections.js","sourceRoot":"","sources":["../../../src/atlas/oscillators/rejections.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;GAcG;AAIH,8DAA8D;AAC9D,MAAM,CAAC,MAAM,gBAAgB,GAA8B;IACzD;QACE,EAAE,EAAE,oBAAoB;QACxB,OAAO,EAAE,mBAAmB;QAC5B,QAAQ,EAAE,CAAC,oBAAoB,CAAC;QAChC,UAAU,EAAE,UAAU;QACtB,MAAM,EACJ,6EAA6E;YAC7E,6EAA6E;YAC7E,6EAA6E;YAC7E,2EAA2E;YAC3E,+EAA+E;YAC/E,gFAAgF;YAChF,4BAA4B;QAC9B,cAAc,EAAE,WAAW;QAC3B,SAAS,EAAE;YACT;gBACE,EAAE,EAAE,IAAI;gBACR,IAAI,EAAE,UAAU;gBAChB,IAAI,EAAE,wCAAwC;gBAC9C,SAAS,EAAE,uEAAuE;aACnF;YACD;gBACE,EAAE,EAAE,KAAK;gBACT,IAAI,EAAE,SAAS;gBACf,IAAI,EAAE,wCAAwC;gBAC9C,SAAS,EAAE,8GAA8G;aAC1H;SACF;KACF;CACF,CAAC","sourcesContent":["/**\n * The Phase 0 REJECTION record: a claimed bridge the atlas records as refuted,\n * together with the witness that refutes it.\n *\n * `ax-cubic-spring-lc` claims the cubic (Duffing-type) spring\n * `m x″ + k x + β x³ = 0` is exactly equivalent to the LC circuit\n * `L q″ + q/C = 0`. It is not. The cubic spring's variable set\n * {m, k, β, x0} carries exactly ONE dimensionless group, `β x0² / k`, which\n * measures the nonlinearity and has no counterpart on the LC side: the LC\n * variable set is dimensionally independent and admits no group at all. A\n * surviving group with no image cannot be matched by any change of variables,\n * so no exact equivalence exists.\n *\n * @module atlas/oscillators/rejections\n */\n\nimport type { AtlasRejection } from '../types.js';\n\n/** The rejected bridges of the oscillator pilot. @internal */\nexport const ATLAS_REJECTIONS: readonly AtlasRejection[] = [\n  {\n    id: 'ax-cubic-spring-lc',\n    claimed: 'exact-equivalence',\n    premises: ['model-cubic-spring'],\n    conclusion: 'model-lc',\n    reason:\n      'The cubic spring retains one dimensionless group, β x0² / k, that measures ' +\n      'its nonlinearity, and the LC variable set has none (W3). A surviving group ' +\n      'is necessary evidence, not sufficient. The decisive fact is isochrony: the ' +\n      'LC period is the same at every amplitude, the cubic spring period is not ' +\n      '(W3b), and a change of variables that rescales time by a constant multiplies ' +\n      'every period by the same factor. So no such change of variables makes the two ' +\n      'models exactly equivalent.',\n    survivingGroup: 'β x0² / k',\n    witnesses: [\n      {\n        id: 'W3',\n        kind: 'symbolic',\n        test: 'tests/atlas/oscillators-coarse.test.ts',\n        tolerance: 'exact integer π-group exponents; LC verdict dimensionally-independent',\n      },\n      {\n        id: 'W3b',\n        kind: 'numeric',\n        test: 'tests/atlas/oscillators-coarse.test.ts',\n        tolerance: 'linear period 2π within 1e-6 at A = 0.5, 1, 2; cubic-spring periods spread by more than 10% at β x0²/k = 0.1',\n      },\n    ],\n  },\n];\n"]}