universal-physics-tensor 0.47.1 → 1.0.2

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

@@ -0,0 +1,122 @@
+/**
+ * The pendulum and its linear oscillator, integrated and read off, for the
+ * witnesses of `ab-pendulum-linear`'s declared translations. Time is in units
+ * of T0 (ω0 = 2π); both motions are released from rest at θ0.
+ *
+ * Also the closed forms those translations are built from: θ0 recovered from
+ * the bound's ε, and the nome q of the pendulum's elliptic modulus.
+ *
+ * @module atlas/oscillators/pendulum-motion
+ * @internal
+ */
+import { agm, pendulumPeriodErrorAt } from './bridges-limits.js';
+/** ω0² at T0 = 1. @internal */
+export const W2 = (2 * Math.PI) ** 2;
+/** @internal */
+export const pendulumAccel = (x) => -W2 * Math.sin(x);
+/** @internal */
+export const linearAccel = (x) => -W2 * x;
+/** One classical RK4 step of x″ = accel(x). @internal */
+export function rk4Step(x, v, h, accel) {
+    const k1x = v;
+    const k1v = accel(x);
+    const k2x = v + 0.5 * h * k1v;
+    const k2v = accel(x + 0.5 * h * k1x);
+    const k3x = v + 0.5 * h * k2v;
+    const k3v = accel(x + 0.5 * h * k2x);
+    const k4x = v + h * k3v;
+    const k4v = accel(x + h * k3x);
+    return [x + (h / 6) * (k1x + 2 * k2x + 2 * k3x + k4x), v + (h / 6) * (k1v + 2 * k2v + 2 * k3v + k4v)];
+}
+/**
+ * The root in (0, h] of x(τ), the position one partial RK4 step of length τ
+ * from (x, v), by the secant method from the bracket (0, x) and (h, nx).
+ *
+ * The crossing is then located to the integrator's own truncation error,
+ * rather than to the O(h³) error of a straight line between the two steps.
+ * That truncation error is nearly the same for the pendulum and the linear
+ * oscillator at small θ0, so it cancels in their phase difference; the
+ * interpolation error does not.
+ */
+function refineCrossing(x, v, nx, h, accel) {
+    let a = 0;
+    let fa = x;
+    let b = h;
+    let fb = nx;
+    for (let k = 0; k < 50 && fb !== 0 && fb !== fa; k++) {
+        const c = b - (fb * (b - a)) / (fb - fa);
+        if (!(c >= -h && c <= 2 * h) || c === b)
+            break;
+        a = b;
+        fa = fb;
+        b = c;
+        fb = rk4Step(x, v, c, accel)[0];
+    }
+    return b;
+}
+/** Zero-crossing times of a release from rest at `theta0`, over `periods` of T0 = 1. @internal */
+export function crossingTimes(accel, theta0, stepsPerPeriod, periods) {
+    const h = 1 / stepsPerPeriod;
+    const out = [];
+    let x = theta0;
+    let v = 0;
+    for (let i = 0; i < stepsPerPeriod * periods; i++) {
+        const [nx, nv] = rk4Step(x, v, h, accel);
+        if (nx === 0 || x > 0 !== nx > 0)
+            out.push(i * h + refineCrossing(x, v, nx, h, accel));
+        x = nx;
+        v = nv;
+    }
+    return out;
+}
+/** The phase at `t`, interpolated through the anchors (0, 0) and (c_k, π/2 + kπ). @internal */
+export function phaseAt(crossings, t) {
+    let prevT = 0;
+    let prevPhase = 0;
+    for (let k = 0; k < crossings.length; k++) {
+        const c = crossings[k];
+        const phase = Math.PI / 2 + k * Math.PI;
+        if (t <= c)
+            return prevPhase + ((phase - prevPhase) * (t - prevT)) / (c - prevT);
+        prevT = c;
+        prevPhase = phase;
+    }
+    return Number.NaN;
+}
+/**
+ * θ0 in [0, π) whose exact relative period error is `eps`, by bisection on the
+ * strictly increasing `pendulumPeriodErrorAt`.
+ *
+ * @returns NaN for a negative or non-finite `eps`.
+ * @internal
+ */
+export function theta0OfPeriodError(eps) {
+    if (!Number.isFinite(eps) || eps < 0)
+        return Number.NaN;
+    if (eps === 0)
+        return 0;
+    let lo = 0;
+    let hi = Math.PI;
+    for (let i = 0; i < 200 && hi - lo > 0; i++) {
+        const mid = (lo + hi) / 2;
+        if (mid === lo || mid === hi)
+            break;
+        if (pendulumPeriodErrorAt({ theta0: mid }) < eps)
+            lo = mid;
+        else
+            hi = mid;
+    }
+    return (lo + hi) / 2;
+}
+/**
+ * The nome q = exp(−π K′/K) of the modulus k = sin(θ0/2), with
+ * K′/K = AGM(1, cos(θ0/2)) / AGM(1, sin(θ0/2)).
+ * @internal
+ */
+export function pendulumNome(theta0) {
+    const k = Math.sin(Math.abs(theta0) / 2);
+    if (k === 0)
+        return 0;
+    return Math.exp((-Math.PI * agm(1, Math.cos(theta0 / 2))) / agm(1, k));
+}
+//# sourceMappingURL=pendulum-motion.js.map

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

@@ -0,0 +1 @@
+{"version":3,"file":"pendulum-motion.js","sourceRoot":"","sources":["../../../src/atlas/oscillators/pendulum-motion.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;GAUG;AAEH,OAAO,EAAE,GAAG,EAAE,qBAAqB,EAAE,MAAM,qBAAqB,CAAC;AAEjE,+BAA+B;AAC/B,MAAM,CAAC,MAAM,EAAE,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC,IAAI,CAAC,CAAC;AACrC,gBAAgB;AAChB,MAAM,CAAC,MAAM,aAAa,GAAG,CAAC,CAAS,EAAU,EAAE,CAAC,CAAC,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;AACtE,gBAAgB;AAChB,MAAM,CAAC,MAAM,WAAW,GAAG,CAAC,CAAS,EAAU,EAAE,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC;AAE1D,yDAAyD;AACzD,MAAM,UAAU,OAAO,CAAC,CAAS,EAAE,CAAS,EAAE,CAAS,EAAE,KAA4B;IACnF,MAAM,GAAG,GAAG,CAAC,CAAC;IACd,MAAM,GAAG,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC;IACrB,MAAM,GAAG,GAAG,CAAC,GAAG,GAAG,GAAG,CAAC,GAAG,GAAG,CAAC;IAC9B,MAAM,GAAG,GAAG,KAAK,CAAC,CAAC,GAAG,GAAG,GAAG,CAAC,GAAG,GAAG,CAAC,CAAC;IACrC,MAAM,GAAG,GAAG,CAAC,GAAG,GAAG,GAAG,CAAC,GAAG,GAAG,CAAC;IAC9B,MAAM,GAAG,GAAG,KAAK,CAAC,CAAC,GAAG,GAAG,GAAG,CAAC,GAAG,GAAG,CAAC,CAAC;IACrC,MAAM,GAAG,GAAG,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC;IACxB,MAAM,GAAG,GAAG,KAAK,CAAC,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC,CAAC;IAC/B,OAAO,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC,GAAG,GAAG,GAAG,CAAC,GAAG,GAAG,GAAG,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC,GAAG,GAAG,GAAG,CAAC,GAAG,GAAG,GAAG,GAAG,CAAC,CAAC,CAAC;AACxG,CAAC;AAED;;;;;;;;;GASG;AACH,SAAS,cAAc,CAAC,CAAS,EAAE,CAAS,EAAE,EAAU,EAAE,CAAS,EAAE,KAA4B;IAC/F,IAAI,CAAC,GAAG,CAAC,CAAC;IACV,IAAI,EAAE,GAAG,CAAC,CAAC;IACX,IAAI,CAAC,GAAG,CAAC,CAAC;IACV,IAAI,EAAE,GAAG,EAAE,CAAC;IACZ,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,EAAE,IAAI,EAAE,KAAK,CAAC,IAAI,EAAE,KAAK,EAAE,EAAE,CAAC,EAAE,EAAE,CAAC;QACrD,MAAM,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,EAAE,GAAG,EAAE,CAAC,CAAC;QACzC,IAAI,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,IAAI,CAAC,KAAK,CAAC;YAAE,MAAM;QAC/C,CAAC,GAAG,CAAC,CAAC;QACN,EAAE,GAAG,EAAE,CAAC;QACR,CAAC,GAAG,CAAC,CAAC;QACN,EAAE,GAAG,OAAO,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC;IAClC,CAAC;IACD,OAAO,CAAC,CAAC;AACX,CAAC;AAED,kGAAkG;AAClG,MAAM,UAAU,aAAa,CAAC,KAA4B,EAAE,MAAc,EAAE,cAAsB,EAAE,OAAe;IACjH,MAAM,CAAC,GAAG,CAAC,GAAG,cAAc,CAAC;IAC7B,MAAM,GAAG,GAAa,EAAE,CAAC;IACzB,IAAI,CAAC,GAAG,MAAM,CAAC;IACf,IAAI,CAAC,GAAG,CAAC,CAAC;IACV,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,cAAc,GAAG,OAAO,EAAE,CAAC,EAAE,EAAE,CAAC;QAClD,MAAM,CAAC,EAAE,EAAE,EAAE,CAAC,GAAG,OAAO,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,KAAK,CAAC,CAAC;QACzC,IAAI,EAAE,KAAK,CAAC,IAAI,CAAC,GAAG,CAAC,KAAK,EAAE,GAAG,CAAC;YAAE,GAAG,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,GAAG,cAAc,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,EAAE,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC;QACvF,CAAC,GAAG,EAAE,CAAC;QACP,CAAC,GAAG,EAAE,CAAC;IACT,CAAC;IACD,OAAO,GAAG,CAAC;AACb,CAAC;AAED,+FAA+F;AAC/F,MAAM,UAAU,OAAO,CAAC,SAA4B,EAAE,CAAS;IAC7D,IAAI,KAAK,GAAG,CAAC,CAAC;IACd,IAAI,SAAS,GAAG,CAAC,CAAC;IAClB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,SAAS,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;QAC1C,MAAM,CAAC,GAAG,SAAS,CAAC,CAAC,CAAE,CAAC;QACxB,MAAM,KAAK,GAAG,IAAI,CAAC,EAAE,GAAG,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC;QACxC,IAAI,CAAC,IAAI,CAAC;YAAE,OAAO,SAAS,GAAG,CAAC,CAAC,KAAK,GAAG,SAAS,CAAC,GAAG,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC;QACjF,KAAK,GAAG,CAAC,CAAC;QACV,SAAS,GAAG,KAAK,CAAC;IACpB,CAAC;IACD,OAAO,MAAM,CAAC,GAAG,CAAC;AACpB,CAAC;AAED;;;;;;GAMG;AACH,MAAM,UAAU,mBAAmB,CAAC,GAAW;IAC7C,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,GAAG,CAAC,IAAI,GAAG,GAAG,CAAC;QAAE,OAAO,MAAM,CAAC,GAAG,CAAC;IACxD,IAAI,GAAG,KAAK,CAAC;QAAE,OAAO,CAAC,CAAC;IACxB,IAAI,EAAE,GAAG,CAAC,CAAC;IACX,IAAI,EAAE,GAAG,IAAI,CAAC,EAAE,CAAC;IACjB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,GAAG,IAAI,EAAE,GAAG,EAAE,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;QAC5C,MAAM,GAAG,GAAG,CAAC,EAAE,GAAG,EAAE,CAAC,GAAG,CAAC,CAAC;QAC1B,IAAI,GAAG,KAAK,EAAE,IAAI,GAAG,KAAK,EAAE;YAAE,MAAM;QACpC,IAAI,qBAAqB,CAAC,EAAE,MAAM,EAAE,GAAG,EAAE,CAAC,GAAG,GAAG;YAAE,EAAE,GAAG,GAAG,CAAC;;YACtD,EAAE,GAAG,GAAG,CAAC;IAChB,CAAC;IACD,OAAO,CAAC,EAAE,GAAG,EAAE,CAAC,GAAG,CAAC,CAAC;AACvB,CAAC;AAED;;;;GAIG;AACH,MAAM,UAAU,YAAY,CAAC,MAAc;IACzC,MAAM,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC;IACzC,IAAI,CAAC,KAAK,CAAC;QAAE,OAAO,CAAC,CAAC;IACtB,OAAO,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,IAAI,CAAC,EAAE,GAAG,GAAG,CAAC,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;AACzE,CAAC"}

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

@@ -0,0 +1,38 @@
+/**
+ * `ab-spring-lc`'s declaration that its map carries the PHASE unchanged, so
+ * `ab-pendulum-linear`'s phase translation holds on the path
+ * pendulum → spring → LC.
+ *
+ * The bridge's composed map is q(t_LC) = (q0/x0)·x(t) with t_LC = t·ω_s/ω_LC.
+ * The charge is a fixed nonzero multiple of the displacement at the mapped
+ * time, so the two share every zero crossing and turning point, and the angle
+ * variable (0 at release from rest, 2π per period, uniform in time) of the
+ * circuit at t_LC equals the spring's at t. A phase error against the spring
+ * at t is therefore the same phase error against the circuit at t_LC: nothing
+ * is added and nothing truncated.
+ *
+ * This is a claim about the phase only. It says nothing about the position
+ * error — x and q are different quantities in different units, and a bound in
+ * rad on the angle does not become one on the charge — so `ab-spring-lc`
+ * declares no carriage for `position`.
+ *
+ * The witness W1φ integrates the two oscillators of W1a (m = 2, k = 8,
+ * x0 = 0.37 and L = 0.5, C = 0.25, q0 = −1.9: ω_s = 2, ω_LC = 2√2, and a
+ * NEGATIVE ratio q0/x0) by RK4 at the same step in physical time, and
+ * compares the circuit's phase at the mapped time with the spring's at every
+ * spring zero crossing. The angle variable of each motion is measured from its
+ * own release, so the sign of q0/x0 does not enter.
+ *
+ * @module atlas/oscillators/phase-carriage
+ * @internal
+ */
+import type { ObservableCarriage } from '../translation.js';
+/**
+ * max over the spring's zero crossings c_k of |φ_LC(timeScale · c_k) − (π/2 + kπ)|,
+ * both integrated at `stepsPerSpringPeriod` steps of the spring's period.
+ * @internal
+ */
+export declare function measureCarriedPhaseError(stepsPerSpringPeriod: number, timeScale: number): number;
+/** @internal */
+export declare const SPRING_LC_PHASE_CARRIAGE: ObservableCarriage;
+//# sourceMappingURL=phase-carriage.d.ts.map

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

@@ -0,0 +1 @@
+{"version":3,"file":"phase-carriage.d.ts","sourceRoot":"","sources":["../../../src/atlas/oscillators/phase-carriage.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;;;;;;;;;GA2BG;AAEH,OAAO,KAAK,EAAE,kBAAkB,EAAE,MAAM,mBAAmB,CAAC;AAwB5D;;;;GAIG;AACH,wBAAgB,wBAAwB,CAAC,oBAAoB,EAAE,MAAM,EAAE,SAAS,EAAE,MAAM,GAAG,MAAM,CAUhG;AAWD,gBAAgB;AAChB,eAAO,MAAM,wBAAwB,EAAE,kBAsBtC,CAAC"}

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

@@ -0,0 +1,96 @@
+/**
+ * `ab-spring-lc`'s declaration that its map carries the PHASE unchanged, so
+ * `ab-pendulum-linear`'s phase translation holds on the path
+ * pendulum → spring → LC.
+ *
+ * The bridge's composed map is q(t_LC) = (q0/x0)·x(t) with t_LC = t·ω_s/ω_LC.
+ * The charge is a fixed nonzero multiple of the displacement at the mapped
+ * time, so the two share every zero crossing and turning point, and the angle
+ * variable (0 at release from rest, 2π per period, uniform in time) of the
+ * circuit at t_LC equals the spring's at t. A phase error against the spring
+ * at t is therefore the same phase error against the circuit at t_LC: nothing
+ * is added and nothing truncated.
+ *
+ * This is a claim about the phase only. It says nothing about the position
+ * error — x and q are different quantities in different units, and a bound in
+ * rad on the angle does not become one on the charge — so `ab-spring-lc`
+ * declares no carriage for `position`.
+ *
+ * The witness W1φ integrates the two oscillators of W1a (m = 2, k = 8,
+ * x0 = 0.37 and L = 0.5, C = 0.25, q0 = −1.9: ω_s = 2, ω_LC = 2√2, and a
+ * NEGATIVE ratio q0/x0) by RK4 at the same step in physical time, and
+ * compares the circuit's phase at the mapped time with the spring's at every
+ * spring zero crossing. The angle variable of each motion is measured from its
+ * own release, so the sign of q0/x0 does not enter.
+ *
+ * @module atlas/oscillators/phase-carriage
+ * @internal
+ */
+import { phaseAt, rk4Step } from './pendulum-motion.js';
+const [M, K, L, C, X0, Q0] = [2, 8, 0.5, 0.25, 0.37, -1.9];
+const OMEGA_S = Math.sqrt(K / M);
+const OMEGA_LC = Math.sqrt(1 / (L * C));
+const SPRING_PERIODS = 20;
+/** Zero crossings of u″ = −ω²u from rest at `u0`, stepping `h` in the oscillator's own time. */
+function crossings(omega, u0, h, tEnd) {
+    const out = [];
+    let x = u0;
+    let v = 0;
+    const w2 = omega * omega;
+    for (let i = 0; i * h < tEnd; i++) {
+        const [nx, nv] = rk4Step(x, v, h, (s) => -w2 * s);
+        if (nx === 0 || x > 0 !== nx > 0)
+            out.push(i * h + (h * x) / (x - nx));
+        x = nx;
+        v = nv;
+    }
+    return out;
+}
+/**
+ * max over the spring's zero crossings c_k of |φ_LC(timeScale · c_k) − (π/2 + kπ)|,
+ * both integrated at `stepsPerSpringPeriod` steps of the spring's period.
+ * @internal
+ */
+export function measureCarriedPhaseError(stepsPerSpringPeriod, timeScale) {
+    const tSpring = SPRING_PERIODS * ((2 * Math.PI) / OMEGA_S);
+    const h = (2 * Math.PI) / OMEGA_S / stepsPerSpringPeriod;
+    const spring = crossings(OMEGA_S, X0, h, tSpring);
+    const lc = crossings(OMEGA_LC, Q0, h, tSpring * Math.max(1, timeScale) + 1);
+    let worst = 0;
+    for (let k = 0; k < spring.length; k++) {
+        worst = Math.max(worst, Math.abs(phaseAt(lc, timeScale * spring[k]) - (Math.PI / 2 + k * Math.PI)));
+    }
+    return worst;
+}
+const W1PHI = {
+    id: 'W1φ',
+    evaluate: (steps) => measureCarriedPhaseError(steps, OMEGA_S / OMEGA_LC),
+    target: 0,
+    coarseResolution: 100,
+    fineResolution: 400,
+    tolerance: 1e-6,
+};
+/** @internal */
+export const SPRING_LC_PHASE_CARRIAGE = {
+    bridgeId: 'ab-spring-lc',
+    observable: 'phase',
+    map: 'the circuit\'s angle variable at t_LC = t·ω_s/ω_LC is the spring\'s at t, because q(t_LC) = (q0/x0)·x(t)',
+    derivation: 'a fixed nonzero multiple of the displacement at the mapped time shares its zero crossings and turning points, so the ' +
+        'angle variable is the same function of the mapped time; a phase error against the spring at t is the same ' +
+        'error against the circuit at t_LC, exactly',
+    premises: [
+        'the circuit is released from rest: q(0) = q0, i(0) = 0 (either sign of q0/x0: each phase is measured from its own release)',
+        'lossless and unforced (the bridge\'s own side conditions)',
+    ],
+    timeMap: 't is read on the spring\'s clock (the unit of T0); the circuit is read at t·ω_s/ω_LC, the same t when LC = m/k',
+    witnesses: [
+        {
+            id: 'W1φ',
+            kind: 'numeric',
+            test: 'tests/atlas/spring-lc-phase-carriage.test.ts',
+            tolerance: `RK4 circuit phase at the mapped time within 1e-6 rad of the spring's at every crossing over ${SPRING_PERIODS} spring periods (W1a's m = 2, k = 8, x0 = 0.37; L = 0.5, C = 0.25, q0 = −1.9)`,
+        },
+    ],
+    checks: [W1PHI],
+};
+//# sourceMappingURL=phase-carriage.js.map

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

@@ -0,0 +1 @@
+{"version":3,"file":"phase-carriage.js","sourceRoot":"","sources":["../../../src/atlas/oscillators/phase-carriage.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;;;;;;;;;GA2BG;AAIH,OAAO,EAAE,OAAO,EAAE,OAAO,EAAE,MAAM,sBAAsB,CAAC;AAExD,MAAM,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,EAAE,EAAE,EAAE,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,GAAG,EAAE,IAAI,EAAE,IAAI,EAAE,CAAC,GAAG,CAAC,CAAC;AAC3D,MAAM,OAAO,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;AACjC,MAAM,QAAQ,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;AACxC,MAAM,cAAc,GAAG,EAAE,CAAC;AAE1B,gGAAgG;AAChG,SAAS,SAAS,CAAC,KAAa,EAAE,EAAU,EAAE,CAAS,EAAE,IAAY;IACnE,MAAM,GAAG,GAAa,EAAE,CAAC;IACzB,IAAI,CAAC,GAAG,EAAE,CAAC;IACX,IAAI,CAAC,GAAG,CAAC,CAAC;IACV,MAAM,EAAE,GAAG,KAAK,GAAG,KAAK,CAAC;IACzB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,GAAG,IAAI,EAAE,CAAC,EAAE,EAAE,CAAC;QAClC,MAAM,CAAC,EAAE,EAAE,EAAE,CAAC,GAAG,OAAO,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC;QAClD,IAAI,EAAE,KAAK,CAAC,IAAI,CAAC,GAAG,CAAC,KAAK,EAAE,GAAG,CAAC;YAAE,GAAG,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC;QACvE,CAAC,GAAG,EAAE,CAAC;QACP,CAAC,GAAG,EAAE,CAAC;IACT,CAAC;IACD,OAAO,GAAG,CAAC;AACb,CAAC;AAED;;;;GAIG;AACH,MAAM,UAAU,wBAAwB,CAAC,oBAA4B,EAAE,SAAiB;IACtF,MAAM,OAAO,GAAG,cAAc,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC,GAAG,OAAO,CAAC,CAAC;IAC3D,MAAM,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC,GAAG,OAAO,GAAG,oBAAoB,CAAC;IACzD,MAAM,MAAM,GAAG,SAAS,CAAC,OAAO,EAAE,EAAE,EAAE,CAAC,EAAE,OAAO,CAAC,CAAC;IAClD,MAAM,EAAE,GAAG,SAAS,CAAC,QAAQ,EAAE,EAAE,EAAE,CAAC,EAAE,OAAO,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,SAAS,CAAC,GAAG,CAAC,CAAC,CAAC;IAC5E,IAAI,KAAK,GAAG,CAAC,CAAC;IACd,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,MAAM,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;QACvC,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,EAAE,IAAI,CAAC,GAAG,CAAC,OAAO,CAAC,EAAE,EAAE,SAAS,GAAG,MAAM,CAAC,CAAC,CAAE,CAAC,GAAG,CAAC,IAAI,CAAC,EAAE,GAAG,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;IACvG,CAAC;IACD,OAAO,KAAK,CAAC;AACf,CAAC;AAED,MAAM,KAAK,GAAuB;IAChC,EAAE,EAAE,KAAK;IACT,QAAQ,EAAE,CAAC,KAAK,EAAE,EAAE,CAAC,wBAAwB,CAAC,KAAK,EAAE,OAAO,GAAG,QAAQ,CAAC;IACxE,MAAM,EAAE,CAAC;IACT,gBAAgB,EAAE,GAAG;IACrB,cAAc,EAAE,GAAG;IACnB,SAAS,EAAE,IAAI;CAChB,CAAC;AAEF,gBAAgB;AAChB,MAAM,CAAC,MAAM,wBAAwB,GAAuB;IAC1D,QAAQ,EAAE,cAAc;IACxB,UAAU,EAAE,OAAO;IACnB,GAAG,EAAE,0GAA0G;IAC/G,UAAU,EACR,uHAAuH;QACvH,4GAA4G;QAC5G,4CAA4C;IAC9C,QAAQ,EAAE;QACR,4HAA4H;QAC5H,2DAA2D;KAC5D;IACD,OAAO,EAAE,gHAAgH;IACzH,SAAS,EAAE;QACT;YACE,EAAE,EAAE,KAAK;YACT,IAAI,EAAE,SAAS;YACf,IAAI,EAAE,8CAA8C;YACpD,SAAS,EAAE,+FAA+F,cAAc,+EAA+E;SACxM;KACF;IACD,MAAM,EAAE,CAAC,KAAK,CAAC;CAChB,CAAC"}

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

@@ -0,0 +1,46 @@
+/**
+ * `ab-pendulum-linear`'s declared translation from its bound's quantity, the
+ * relative period error, into the accumulated PHASE error.
+ *
+ * Both motions start from rest at θ0 at t = 0. The phase is the angle
+ * variable: 0 at release and advancing by 2π per period, uniformly in time,
+ * so the pendulum's is 2π t/T and the linear oscillator's 2π t/T0. With the
+ * bound's `ε = T/T0 − 1` (normalised by the reduced model's period),
+ *
+ *   Δφ(t) = 2π t/T0 − 2π t/T = 2π (t/T0) · ε/(1 + ε),
+ *
+ * which reaches a tolerance Δ at `t* = Δ · T0 · (1 + ε)/(2π ε)`. Nothing is
+ * truncated: the map is exact in ε, and ε itself is the bound's closed form.
+ *
+ * It is NOT a bound on the pointwise position error |θ(t) − θ_lin(t)|: the
+ * pendulum's waveform carries a third harmonic of amplitude ≈ θ0³/192, which
+ * is not a phase, and the position error at a given phase error depends on
+ * where in the cycle the observation falls.
+ *
+ * The witness W7p integrates BOTH equations of motion and reads each phase off
+ * its zero crossings — crossing k is at phase π/2 + kπ by the symmetry of a
+ * release from rest, and the angle variable is linear in time between them —
+ * so it measures the drift without using ε or the elliptic integral. It runs
+ * at the fixture θ0 = 0.2 as the regression anchor, and again at the caller's
+ * point, bounded to 512 T0 of integration (`phasePointCheck`).
+ *
+ * @module atlas/oscillators/phase-translation
+ * @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 */
+export declare const PENDULUM_PHASE_TRANSLATION: ObservableTranslation;
+//# sourceMappingURL=phase-translation.d.ts.map

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

@@ -0,0 +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"}

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

@@ -0,0 +1,178 @@
+/**
+ * `ab-pendulum-linear`'s declared translation from its bound's quantity, the
+ * relative period error, into the accumulated PHASE error.
+ *
+ * Both motions start from rest at θ0 at t = 0. The phase is the angle
+ * variable: 0 at release and advancing by 2π per period, uniformly in time,
+ * so the pendulum's is 2π t/T and the linear oscillator's 2π t/T0. With the
+ * bound's `ε = T/T0 − 1` (normalised by the reduced model's period),
+ *
+ *   Δφ(t) = 2π t/T0 − 2π t/T = 2π (t/T0) · ε/(1 + ε),
+ *
+ * which reaches a tolerance Δ at `t* = Δ · T0 · (1 + ε)/(2π ε)`. Nothing is
+ * truncated: the map is exact in ε, and ε itself is the bound's closed form.
+ *
+ * It is NOT a bound on the pointwise position error |θ(t) − θ_lin(t)|: the
+ * pendulum's waveform carries a third harmonic of amplitude ≈ θ0³/192, which
+ * is not a phase, and the position error at a given phase error depends on
+ * where in the cycle the observation falls.
+ *
+ * The witness W7p integrates BOTH equations of motion and reads each phase off
+ * its zero crossings — crossing k is at phase π/2 + kπ by the symmetry of a
+ * release from rest, and the angle variable is linear in time between them —
+ * so it measures the drift without using ε or the elliptic integral. It runs
+ * at the fixture θ0 = 0.2 as the regression anchor, and again at the caller's
+ * point, bounded to 512 T0 of integration (`phasePointCheck`).
+ *
+ * @module atlas/oscillators/phase-translation
+ * @internal
+ */
+import { pendulumPeriodErrorAt } from './bridges-limits.js';
+import { crossingTimes, linearAccel, pendulumAccel, phaseAt, theta0OfPeriodError } from './pendulum-motion.js';
+const TWO_PI = 2 * Math.PI;
+/**
+ * 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 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))
+            return t;
+    }
+    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) {
+    const pend = crossingTimes(pendulumAccel, theta0, stepsPerPeriod, periods);
+    const lin = crossingTimes(linearAccel, theta0, stepsPerPeriod, periods);
+    let prevT = 0;
+    let prevDrift = 0;
+    for (let k = 0; k < pend.length; k++) {
+        const c = pend[k];
+        const drift = phaseAt(lin, c) - (Math.PI / 2 + k * Math.PI);
+        if (!Number.isFinite(drift))
+            break;
+        if (drift >= tolerance)
+            return prevT + ((tolerance - prevDrift) * (c - prevT)) / (drift - prevDrift);
+        prevT = c;
+        prevDrift = drift;
+    }
+    return Number.NaN;
+}
+/** `t* = Δ T0 (1+ε)/(2π ε)`; `Infinity` at ε = 0, where no phase accumulates. */
+function phaseHorizon(eps, tolerance, T0) {
+    if (T0 === undefined || !Number.isFinite(T0) || T0 <= 0 || !Number.isFinite(eps) || eps < 0)
+        return Number.NaN;
+    return eps === 0 ? Infinity : (tolerance * T0 * (1 + eps)) / (TWO_PI * eps);
+}
+const W7P_THETA0 = 0.2;
+const W7P_TOLERANCE = 1;
+/** The point witness integrates at most this many periods of T0, so its cost is bounded. */
+export const PHASE_POINT_MAX_PERIODS = 512;
+/**
+ * The point witness measures the time to this much drift, or to the drift
+ * reached at the cap if less: the map is linear in t, so its slope is what is
+ * checked, whatever tolerance the caller asked about.
+ */
+const PHASE_POINT_DRIFT = 1;
+/**
+ * Below this measured drift the rounding noise of the two motions' crossing
+ * times comes within a factor of ten of the witness tolerance. With the
+ * crossings located by `crossingTimes`' secant refinement, the fine error was
+ * measured at 800 steps per T0 as 1.5e-8 relative at 1.1e-5 rad of drift,
+ * 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
+ * 2e-2 rad up. With the straight-line crossings it replaced, the same error
+ * was 9.5e-6 relative at 2e-4 rad and 1.2e-7 at 2e-2 rad.
+ */
+const PHASE_POINT_MIN_DRIFT = 1e-4;
+/**
+ * The (1+ε) control differs from the target by ε/(1+ε) relative, so it can be
+ * refuted only where ε exceeds this tolerance (θ0 above about 9e-4).
+ */
+const PHASE_POINT_RELATIVE_TOLERANCE = 5e-8;
+/**
+ * W7p at the caller's θ0 (recovered from ε) instead of the fixture: the time
+ * the measured drift reaches min(1 rad, the drift at the cap), within
+ * PHASE_POINT_RELATIVE_TOLERANCE of the closed form. Its control is the same measurement against
+ * the horizon without the (1+ε) factor, which differs from the target by
+ * ε/(1+ε) relative — so the control can fail only where that exceeds the
+ * witness tolerance, and the caller is told when it does not.
+ */
+function phasePointCheck(eps) {
+    if (!Number.isFinite(eps) || eps <= 0)
+        return { unavailable: `no phase drift accumulates at ε = ${eps}` };
+    const theta0 = theta0OfPeriodError(eps);
+    const measured = Math.min(PHASE_POINT_DRIFT, (TWO_PI * PHASE_POINT_MAX_PERIODS * eps) / (1 + eps));
+    if (measured < PHASE_POINT_MIN_DRIFT) {
+        return {
+            unavailable: `the drift within ${PHASE_POINT_MAX_PERIODS} T0 at θ0 = ${Number(theta0.toPrecision(8))} is ${measured} rad, below the ` +
+                `${PHASE_POINT_MIN_DRIFT} rad this witness resolves; not run`,
+        };
+    }
+    const target = phaseHorizon(eps, measured, 1);
+    const periods = Math.ceil(target) + 2;
+    const evaluate = (stepsPerPeriod) => measurePhaseHorizonWithin(theta0, measured, stepsPerPeriod, periods);
+    const spec = { evaluate, coarseResolution: 200, fineResolution: 800, tolerance: PHASE_POINT_RELATIVE_TOLERANCE * target };
+    const capped = measured < PHASE_POINT_DRIFT ? ` (the drift reached in ${PHASE_POINT_MAX_PERIODS} T0)` : '';
+    return {
+        claim: `RK4-measured time to ${measured} rad${capped} of phase drift at θ0 = ${Number(theta0.toPrecision(8))} (recovered from ε) ` +
+            `within ${PHASE_POINT_RELATIVE_TOLERANCE} relative of t* = ${target} T0; 200 and 800 steps per T0, ${periods} T0 integrated`,
+        check: { id: 'W7p@point', ...spec, target },
+        control: { id: 'W7p@point-control', ...spec, target: measured / (TWO_PI * eps) },
+        controlClaim: 'the same measurement against the horizon without the (1+ε) factor, Δ T0/(2π ε)',
+    };
+}
+/** @internal */
+export const PENDULUM_PHASE_TRANSLATION = {
+    bridgeId: 'ab-pendulum-linear',
+    observable: 'phase',
+    unit: 'rad',
+    fromNorm: 'relative period error, normalized by the value of the reduced model',
+    errorKind: 'exact',
+    definition: 'the angle variable of each motion (0 at release from rest at θ0, 2π per period); the error is the ' +
+        "linear oscillator's lead over the pendulum, Δφ(t) = 2πt/T0 − 2πt/T",
+    derivation: 'with ε = T/T0 − 1, Δφ(t) = 2π (t/T0) ε/(1+ε), exact in ε; it reaches Δ at t* = Δ T0 (1+ε)/(2π ε). ' +
+        'ε/(1+ε) increases with ε, so the domain supremum of ε gives the shortest horizon over the domain',
+    premises: [
+        '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',
+    ],
+    parameters: ['T0'],
+    timeUnit: 'the unit of T0',
+    errorAt: (eps, t, { T0 }) => {
+        if (T0 === undefined || !Number.isFinite(T0) || T0 <= 0 || !Number.isFinite(eps) || !Number.isFinite(t))
+            return Number.NaN;
+        return (TWO_PI * (t / T0) * eps) / (1 + eps);
+    },
+    horizonFor: (eps, tolerance, { T0 }) => phaseHorizon(eps, tolerance, T0),
+    boundary: 'inclusive: adequate iff t ≤ t*, i.e. Δφ(t) ≤ the tolerance',
+    notCovered: [
+        'the pointwise position error |θ(t) − θ_lin(t)| (the waveform carries a third harmonic ≈ θ0³/192 that is not a phase)',
+        'the amplitude',
+    ],
+    witnesses: [
+        {
+            id: 'W7p',
+            kind: 'numeric',
+            test: 'tests/atlas/pendulum-phase-translation.test.ts',
+            tolerance: `RK4-measured time to ${W7P_TOLERANCE} rad of phase drift at θ0 = ${W7P_THETA0} within 1e-4 T0 of t*`,
+        },
+    ],
+    checks: [
+        {
+            id: 'W7p',
+            evaluate: (stepsPerPeriod) => measurePhaseHorizon(W7P_THETA0, W7P_TOLERANCE, stepsPerPeriod),
+            target: phaseHorizon(pendulumPeriodErrorAt({ theta0: W7P_THETA0 }), W7P_TOLERANCE, 1),
+            coarseResolution: 50,
+            fineResolution: 200,
+            tolerance: 1e-4,
+        },
+    ],
+    pointCheck: phasePointCheck,
+};
+//# sourceMappingURL=phase-translation.js.map

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

@@ -0,0 +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"}

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

@@ -0,0 +1,59 @@
+/**
+ * `ab-pendulum-linear`'s declared translation from its bound's quantity, the
+ * relative period error ε, into an UPPER BOUND on the pointwise position error
+ * |θ(t) − θ_lin(t)|, in rad.
+ *
+ * Both motions start from rest at θ0. The pendulum is θ(t) = F(φ), φ = 2πt/T,
+ * with the classical Fourier series of the elliptic solution
+ *
+ *   F(φ) = Σ_{n≥0} (−1)ⁿ aₙ cos((2n+1)φ),  aₙ = 8 q^{n+½} / ((2n+1)(1 + q^{2n+1})),
+ *
+ * q the nome of k = sin(θ0/2); F(0) = θ0 = Σ (−1)ⁿ aₙ. The linear motion is
+ * θ0 cos φ_lin, φ_lin = 2πt/T0. Splitting at θ0 cos φ,
+ *
+ *   |θ − θ_lin| ≤ |F(φ) − θ0 cos φ| + θ0 |cos φ − cos φ_lin|
+ *              ≤ 2 Σ_{n≥1} aₙ + 2θ0 |sin(Δφ/2)|,
+ *
+ * because F(φ) − θ0 cos φ = Σ_{n≥1} (−1)ⁿ aₙ (cos((2n+1)φ) − cos φ). With
+ * aₙ ≤ (8/3) q^{n+½} for n ≥ 1 the first term is at most
+ * W̄ = (16/3) q^{3/2}/(1 − q) ≈ θ0³/96, and Δφ = 2π(t/T0) ε/(1+ε) is the
+ * exact phase lead of `phase-translation.ts`. So
+ *
+ *   B(t) = W̄ + 2θ0 sin(min(Δφ(t), π)/2)
+ *
+ * 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.
+ *
+ * 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
+ * motions are in antiphase and the error reaches 2θ0), and the waveform term
+ * cannot be dropped: without W̄, RK4 finds the error above 2θ0 sin(Δφ/2) by
+ * about 3% (witness control).
+ *
+ * @module atlas/oscillators/position-translation
+ * @internal
+ */
+import type { ObservableTranslation } from '../translation.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 θ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[];
+/**
+ * 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.
+ * @internal
+ */
+export declare function fourierSeriesCheck(id: string, variant?: 'plus' | 'minus'): NumericWitnessSpec;
+/** @internal */
+export declare const PENDULUM_POSITION_TRANSLATION: ObservableTranslation;
+//# sourceMappingURL=position-translation.d.ts.map

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

@@ -0,0 +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"}