universal-physics-tensor 0.47.1 → 1.0.1

package/dist/atlas/oscillators/norm-transport-witness.d.ts

@@ -0,0 +1,37 @@
+/**
+ * The measurement behind W1τ, `ab-spring-lc`'s relative-period norm transport
+ * (`./norm-transport.ts`).
+ *
+ * The pendulum is pushed through the bridge's map q(t_LC) = (q0/θ0)·θ(t),
+ * t_LC = t·ω_s/ω_LC, and integrated as its OWN equation in circuit time,
+ * q'' = −ω_LC² (q0/θ0) sin(θ0 q/q0). Its period is measured by RK4 and compared
+ * with the LC circuit's period 2π√(LC). Nothing here multiplies a pendulum
+ * period by a factor: the image is solved where the circuit lives.
+ *
+ * @module atlas/oscillators/norm-transport-witness
+ * @internal
+ */
+/** The circuit and release a measurement uses. @internal */
+export interface TransportCircuit {
+    readonly L: number;
+    readonly C: number;
+    /** Release charge, > 0. */
+    readonly q0: number;
+}
+/**
+ * First downward zero crossing of the release from rest at `q0 > 0`, by cubic
+ * Hermite interpolation of the bracketing RK4 step (values and derivatives at
+ * both ends), refined by Newton. Released from rest, it is a quarter period.
+ *
+ * @internal
+ */
+export declare function quarterPeriod(accel: (q: number) => number, q0: number, h: number, maxSteps: number): number;
+/**
+ * T_image / T_LC − 1: the relative period error of the pendulum's image against
+ * the circuit, with `stepsPerPeriod` RK4 steps per LC period.
+ *
+ * @returns NaN when `theta0` is not in (0, π) or no crossing is found.
+ * @internal
+ */
+export declare function measureTransportedPeriodError(theta0: number, stepsPerPeriod: number, circuit?: TransportCircuit): number;
+//# sourceMappingURL=norm-transport-witness.d.ts.map

package/dist/atlas/oscillators/norm-transport-witness.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"norm-transport-witness.d.ts","sourceRoot":"","sources":["../../../src/atlas/oscillators/norm-transport-witness.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;GAYG;AAKH,4DAA4D;AAC5D,MAAM,WAAW,gBAAgB;IAC/B,QAAQ,CAAC,CAAC,EAAE,MAAM,CAAC;IACnB,QAAQ,CAAC,CAAC,EAAE,MAAM,CAAC;IACnB,2BAA2B;IAC3B,QAAQ,CAAC,EAAE,EAAE,MAAM,CAAC;CACrB;AAED;;;;;;GAMG;AACH,wBAAgB,aAAa,CAAC,KAAK,EAAE,CAAC,CAAC,EAAE,MAAM,KAAK,MAAM,EAAE,EAAE,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,QAAQ,EAAE,MAAM,GAAG,MAAM,CAmB3G;AAED;;;;;;GAMG;AACH,wBAAgB,6BAA6B,CAC3C,MAAM,EAAE,MAAM,EACd,cAAc,EAAE,MAAM,EACtB,OAAO,GAAE,gBAAgC,GACxC,MAAM,CAOR"}

package/dist/atlas/oscillators/norm-transport-witness.js

@@ -0,0 +1,58 @@
+/**
+ * The measurement behind W1τ, `ab-spring-lc`'s relative-period norm transport
+ * (`./norm-transport.ts`).
+ *
+ * The pendulum is pushed through the bridge's map q(t_LC) = (q0/θ0)·θ(t),
+ * t_LC = t·ω_s/ω_LC, and integrated as its OWN equation in circuit time,
+ * q'' = −ω_LC² (q0/θ0) sin(θ0 q/q0). Its period is measured by RK4 and compared
+ * with the LC circuit's period 2π√(LC). Nothing here multiplies a pendulum
+ * period by a factor: the image is solved where the circuit lives.
+ *
+ * @module atlas/oscillators/norm-transport-witness
+ * @internal
+ */
+import { rk4Step } from './pendulum-motion.js';
+import { W1TAU_FIXTURE } from './norm-transport.js';
+/**
+ * First downward zero crossing of the release from rest at `q0 > 0`, by cubic
+ * Hermite interpolation of the bracketing RK4 step (values and derivatives at
+ * both ends), refined by Newton. Released from rest, it is a quarter period.
+ *
+ * @internal
+ */
+export function quarterPeriod(accel, q0, h, maxSteps) {
+    let x = q0;
+    let v = 0;
+    for (let i = 0; i < maxSteps; i++) {
+        const [nx, nv] = rk4Step(x, v, h, accel);
+        if (x > 0 && nx <= 0) {
+            const [p0, m0, p1, m1] = [x, v * h, nx, nv * h];
+            const p = (s) => (2 * s ** 3 - 3 * s ** 2 + 1) * p0 + (s ** 3 - 2 * s ** 2 + s) * m0 + (-2 * s ** 3 + 3 * s ** 2) * p1 + (s ** 3 - s ** 2) * m1;
+            const dp = (s) => (6 * s ** 2 - 6 * s) * p0 + (3 * s ** 2 - 4 * s + 1) * m0 + (-6 * s ** 2 + 6 * s) * p1 + (3 * s ** 2 - 2 * s) * m1;
+            let s = p0 / (p0 - p1);
+            for (let k = 0; k < 8; k++)
+                s -= p(s) / dp(s);
+            return (i + s) * h;
+        }
+        x = nx;
+        v = nv;
+    }
+    return Number.NaN;
+}
+/**
+ * T_image / T_LC − 1: the relative period error of the pendulum's image against
+ * the circuit, with `stepsPerPeriod` RK4 steps per LC period.
+ *
+ * @returns NaN when `theta0` is not in (0, π) or no crossing is found.
+ * @internal
+ */
+export function measureTransportedPeriodError(theta0, stepsPerPeriod, circuit = W1TAU_FIXTURE) {
+    if (!(theta0 > 0 && theta0 < Math.PI) || !(circuit.q0 > 0))
+        return Number.NaN;
+    const omegaLC = 1 / Math.sqrt(circuit.L * circuit.C);
+    const tLC = (2 * Math.PI) / omegaLC;
+    const accel = (q) => -(omegaLC ** 2) * (circuit.q0 / theta0) * Math.sin((theta0 * q) / circuit.q0);
+    const quarter = quarterPeriod(accel, circuit.q0, tLC / stepsPerPeriod, 2 * stepsPerPeriod);
+    return (4 * quarter) / tLC - 1;
+}
+//# sourceMappingURL=norm-transport-witness.js.map

package/dist/atlas/oscillators/norm-transport-witness.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"norm-transport-witness.js","sourceRoot":"","sources":["../../../src/atlas/oscillators/norm-transport-witness.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;GAYG;AAEH,OAAO,EAAE,OAAO,EAAE,MAAM,sBAAsB,CAAC;AAC/C,OAAO,EAAE,aAAa,EAAE,MAAM,qBAAqB,CAAC;AAUpD;;;;;;GAMG;AACH,MAAM,UAAU,aAAa,CAAC,KAA4B,EAAE,EAAU,EAAE,CAAS,EAAE,QAAgB;IACjG,IAAI,CAAC,GAAG,EAAE,CAAC;IACX,IAAI,CAAC,GAAG,CAAC,CAAC;IACV,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,QAAQ,EAAE,CAAC,EAAE,EAAE,CAAC;QAClC,MAAM,CAAC,EAAE,EAAE,EAAE,CAAC,GAAG,OAAO,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,KAAK,CAAC,CAAC;QACzC,IAAI,CAAC,GAAG,CAAC,IAAI,EAAE,IAAI,CAAC,EAAE,CAAC;YACrB,MAAM,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,EAAE,EAAE,EAAE,GAAG,CAAC,CAAC,CAAC;YAChD,MAAM,CAAC,GAAG,CAAC,CAAS,EAAU,EAAE,CAC9B,CAAC,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC,GAAG,EAAE,CAAC;YACjI,MAAM,EAAE,GAAG,CAAC,CAAS,EAAU,EAAE,CAC/B,CAAC,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC;YACrH,IAAI,CAAC,GAAG,EAAE,GAAG,CAAC,EAAE,GAAG,EAAE,CAAC,CAAC;YACvB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE;gBAAE,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC;YAC9C,OAAO,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC;QACrB,CAAC;QACD,CAAC,GAAG,EAAE,CAAC;QACP,CAAC,GAAG,EAAE,CAAC;IACT,CAAC;IACD,OAAO,MAAM,CAAC,GAAG,CAAC;AACpB,CAAC;AAED;;;;;;GAMG;AACH,MAAM,UAAU,6BAA6B,CAC3C,MAAc,EACd,cAAsB,EACtB,OAAO,GAAqB,aAAa;IAEzC,IAAI,CAAC,CAAC,MAAM,GAAG,CAAC,IAAI,MAAM,GAAG,IAAI,CAAC,EAAE,CAAC,IAAI,CAAC,CAAC,OAAO,CAAC,EAAE,GAAG,CAAC,CAAC;QAAE,OAAO,MAAM,CAAC,GAAG,CAAC;IAC9E,MAAM,OAAO,GAAG,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,OAAO,CAAC,CAAC,GAAG,OAAO,CAAC,CAAC,CAAC,CAAC;IACrD,MAAM,GAAG,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC,GAAG,OAAO,CAAC;IACpC,MAAM,KAAK,GAAG,CAAC,CAAS,EAAU,EAAE,CAAC,CAAC,CAAC,OAAO,IAAI,CAAC,CAAC,GAAG,CAAC,OAAO,CAAC,EAAE,GAAG,MAAM,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,MAAM,GAAG,CAAC,CAAC,GAAG,OAAO,CAAC,EAAE,CAAC,CAAC;IACnH,MAAM,OAAO,GAAG,aAAa,CAAC,KAAK,EAAE,OAAO,CAAC,EAAE,EAAE,GAAG,GAAG,cAAc,EAAE,CAAC,GAAG,cAAc,CAAC,CAAC;IAC3F,OAAO,CAAC,CAAC,GAAG,OAAO,CAAC,GAAG,GAAG,GAAG,CAAC,CAAC;AACjC,CAAC"}

package/dist/atlas/oscillators/norm-transport.d.ts

@@ -0,0 +1,44 @@
+/**
+ * `ab-spring-lc`'s declared norm transport: the spring → LC map carries a bound
+ * in RELATIVE PERIOD ERROR unchanged. Design:
+ * `docs/planning/ADR-transported-norm-composition.md`.
+ *
+ * The bridge's composed map is q(t_LC) = (q0/x0)·x(t) with t_LC = t·ω_s/ω_LC:
+ * an amplitude scale and a uniform time scale. A period is a time, so every
+ * period on the spring side is multiplied by the same ω_s/ω_LC on the circuit
+ * side, and the ratio T_full/T_reduced − 1 is unchanged. K = 1, with no
+ * additive term.
+ *
+ * Declared in ONE direction and ONE norm:
+ *
+ * - LC → spring is not declared. The same algebra gives K = 1 there too, but no
+ *   witness checks it and no route in the atlas uses it, and a declaration
+ *   without its own witness is what the ADR forbids.
+ * - The absolute period error (factor ω_s/ω_LC) and the trajectory error
+ *   (factor q0/x0, metres to coulombs) are not declared: the map is not an
+ *   isometry of either, and a route in those norms stays refused.
+ *
+ * @module atlas/oscillators/norm-transport
+ * @internal
+ */
+import type { NormTransport } from '../types.js';
+/** @internal */
+export declare const SPRING_LC_TRANSPORT_TEST = "tests/atlas/spring-lc-norm-transport.test.ts";
+/**
+ * The circuit W1τ integrates in: L = 0.5, C = 0.25 (W1a's circuit), so
+ * ω_LC = 2√2 while the spring side's ω_s = 2 — the time map is not the identity.
+ * q0 > 0 is the release charge, and θ0 = 0.5 is the edge of ab-pendulum-linear's
+ * domain, where its error is largest.
+ *
+ * @internal
+ */
+export declare const W1TAU_FIXTURE: {
+    readonly L: 0.5;
+    readonly C: 0.25;
+    readonly q0: 0.001;
+    readonly theta0: 0.5;
+    readonly omegaS: 2;
+};
+/** @internal */
+export declare const SPRING_LC_RELATIVE_PERIOD_TRANSPORT: NormTransport;
+//# sourceMappingURL=norm-transport.d.ts.map

package/dist/atlas/oscillators/norm-transport.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"norm-transport.d.ts","sourceRoot":"","sources":["../../../src/atlas/oscillators/norm-transport.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;;;;GAsBG;AAEH,OAAO,KAAK,EAAE,aAAa,EAAE,MAAM,aAAa,CAAC;AAGjD,gBAAgB;AAChB,eAAO,MAAM,wBAAwB,iDAAiD,CAAC;AAEvF;;;;;;;GAOG;AACH,eAAO,MAAM,aAAa;aAAK,CAAC,EAAE,GAAG;aAAE,CAAC,EAAE,IAAI;aAAE,EAAE,EAAE,KAAI;aAAE,MAAM,EAAE,GAAG;aAAE,MAAM,EAAE,CAAC;CAAW,CAAC;AAE5F,gBAAgB;AAChB,eAAO,MAAM,mCAAmC,EAAE,aA+BjD,CAAC"}

package/dist/atlas/oscillators/norm-transport.js

@@ -0,0 +1,66 @@
+/**
+ * `ab-spring-lc`'s declared norm transport: the spring → LC map carries a bound
+ * in RELATIVE PERIOD ERROR unchanged. Design:
+ * `docs/planning/ADR-transported-norm-composition.md`.
+ *
+ * The bridge's composed map is q(t_LC) = (q0/x0)·x(t) with t_LC = t·ω_s/ω_LC:
+ * an amplitude scale and a uniform time scale. A period is a time, so every
+ * period on the spring side is multiplied by the same ω_s/ω_LC on the circuit
+ * side, and the ratio T_full/T_reduced − 1 is unchanged. K = 1, with no
+ * additive term.
+ *
+ * Declared in ONE direction and ONE norm:
+ *
+ * - LC → spring is not declared. The same algebra gives K = 1 there too, but no
+ *   witness checks it and no route in the atlas uses it, and a declaration
+ *   without its own witness is what the ADR forbids.
+ * - The absolute period error (factor ω_s/ω_LC) and the trajectory error
+ *   (factor q0/x0, metres to coulombs) are not declared: the map is not an
+ *   isometry of either, and a route in those norms stays refused.
+ *
+ * @module atlas/oscillators/norm-transport
+ * @internal
+ */
+import { RELATIVE_PERIOD_NORM } from './norms.js';
+/** @internal */
+export const SPRING_LC_TRANSPORT_TEST = 'tests/atlas/spring-lc-norm-transport.test.ts';
+/**
+ * The circuit W1τ integrates in: L = 0.5, C = 0.25 (W1a's circuit), so
+ * ω_LC = 2√2 while the spring side's ω_s = 2 — the time map is not the identity.
+ * q0 > 0 is the release charge, and θ0 = 0.5 is the edge of ab-pendulum-linear's
+ * domain, where its error is largest.
+ *
+ * @internal
+ */
+export const W1TAU_FIXTURE = { L: 0.5, C: 0.25, q0: 1e-3, theta0: 0.5, omegaS: 2 };
+/** @internal */
+export const SPRING_LC_RELATIVE_PERIOD_TRANSPORT = {
+    id: 'nt-spring-lc-relative-period',
+    fromModel: 'model-spring',
+    toModel: 'model-lc',
+    from: RELATIVE_PERIOD_NORM,
+    to: RELATIVE_PERIOD_NORM,
+    K: 1,
+    KAt: () => 1,
+    domain: "every m, k, L, C > 0 and x0, q0 ≠ 0 (the bridge's side conditions); K is the same constant throughout",
+    derivation: 'q(t_LC) = (q0/x0)·x(t) with t_LC = t·ω_s/ω_LC scales every period by the same ω_s/ω_LC, so a ratio of two ' +
+        'periods, T_full/T_reduced − 1, is unchanged: K = 1 and nothing is added',
+    timeMap: {
+        map: 't_LC = t·ω_s/ω_LC, one constant factor for every t',
+        uniform: true,
+        horizon: 'a horizon that reads time only through t/T0 holds unchanged with t read on the circuit clock and ' +
+            'T0 = 2π√(LC): a uniform time scale leaves t/T0 invariant',
+        restateHorizon: (holds) => holds,
+    },
+    uniformity: "preserved: K is one constant over the whole domain, so the incoming bound's uniformity holds unchanged",
+    witness: {
+        id: 'W1τ',
+        kind: 'numeric',
+        test: SPRING_LC_TRANSPORT_TEST,
+        tolerance: "RK4 relative period error of the pendulum's image q'' = −ω_LC²(q0/θ0) sin(θ0 q/q0) against the LC period, " +
+            'in circuit time (L = 0.5, C = 0.25, ω_LC = 2√2 ≠ ω_s = 2), at 640 steps per LC period within 1e-9 of K ' +
+            'times the closed-form elliptic error at θ0 = 0.5; 160 steps is the coarse resolution',
+    },
+    basis: 'numerically-supported',
+};
+//# sourceMappingURL=norm-transport.js.map

package/dist/atlas/oscillators/norm-transport.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"norm-transport.js","sourceRoot":"","sources":["../../../src/atlas/oscillators/norm-transport.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;;;;GAsBG;AAGH,OAAO,EAAE,oBAAoB,EAAE,MAAM,YAAY,CAAC;AAElD,gBAAgB;AAChB,MAAM,CAAC,MAAM,wBAAwB,GAAG,8CAA8C,CAAC;AAEvF;;;;;;;GAOG;AACH,MAAM,CAAC,MAAM,aAAa,GAAG,EAAE,CAAC,EAAE,GAAG,EAAE,CAAC,EAAE,IAAI,EAAE,EAAE,EAAE,IAAI,EAAE,MAAM,EAAE,GAAG,EAAE,MAAM,EAAE,CAAC,EAAW,CAAC;AAE5F,gBAAgB;AAChB,MAAM,CAAC,MAAM,mCAAmC,GAAkB;IAChE,EAAE,EAAE,8BAA8B;IAClC,SAAS,EAAE,cAAc;IACzB,OAAO,EAAE,UAAU;IACnB,IAAI,EAAE,oBAAoB;IAC1B,EAAE,EAAE,oBAAoB;IACxB,CAAC,EAAE,CAAC;IACJ,GAAG,EAAE,GAAG,EAAE,CAAC,CAAC;IACZ,MAAM,EAAE,uGAAuG;IAC/G,UAAU,EACR,4GAA4G;QAC5G,yEAAyE;IAC3E,OAAO,EAAE;QACP,GAAG,EAAE,oDAAoD;QACzD,OAAO,EAAE,IAAI;QACb,OAAO,EACL,mGAAmG;YACnG,0DAA0D;QAC5D,cAAc,EAAE,CAAC,KAAK,EAAE,EAAE,CAAC,KAAK;KACjC;IACD,UAAU,EAAE,wGAAwG;IACpH,OAAO,EAAE;QACP,EAAE,EAAE,KAAK;QACT,IAAI,EAAE,SAAS;QACf,IAAI,EAAE,wBAAwB;QAC9B,SAAS,EACP,4GAA4G;YAC5G,yGAAyG;YACzG,sFAAsF;KACzF;IACD,KAAK,EAAE,uBAAuB;CAC/B,CAAC"}

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

@@ -0,0 +1,11 @@
+/**
+ * Norm names the oscillator records share, so a bound and the transport that
+ * carries it name the same norm by the same string. `boundPath` compares norms
+ * verbatim, and a norm spelled twice is a norm that can drift.
+ *
+ * @module atlas/oscillators/norms
+ * @internal
+ */
+/** The norm `ab-pendulum-linear` states its bound in. @internal */
+export declare const RELATIVE_PERIOD_NORM = "relative period error, normalized by the value of the reduced model";
+//# sourceMappingURL=norms.d.ts.map

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

@@ -0,0 +1 @@
+{"version":3,"file":"norms.d.ts","sourceRoot":"","sources":["../../../src/atlas/oscillators/norms.ts"],"names":[],"mappings":"AAAA;;;;;;;GAOG;AAEH,mEAAmE;AACnE,eAAO,MAAM,oBAAoB,wEAAwE,CAAC"}

package/dist/atlas/oscillators/norms.js

@@ -0,0 +1,11 @@
+/**
+ * Norm names the oscillator records share, so a bound and the transport that
+ * carries it name the same norm by the same string. `boundPath` compares norms
+ * verbatim, and a norm spelled twice is a norm that can drift.
+ *
+ * @module atlas/oscillators/norms
+ * @internal
+ */
+/** The norm `ab-pendulum-linear` states its bound in. @internal */
+export const RELATIVE_PERIOD_NORM = 'relative period error, normalized by the value of the reduced model';
+//# sourceMappingURL=norms.js.map

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

@@ -0,0 +1 @@
+{"version":3,"file":"norms.js","sourceRoot":"","sources":["../../../src/atlas/oscillators/norms.ts"],"names":[],"mappings":"AAAA;;;;;;;GAOG;AAEH,mEAAmE;AACnE,MAAM,CAAC,MAAM,oBAAoB,GAAG,qEAAqE,CAAC"}

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

@@ -0,0 +1,38 @@
+/**
+ * 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
+ */
+/** ω0² at T0 = 1. @internal */
+export declare const W2: number;
+/** @internal */
+export declare const pendulumAccel: (x: number) => number;
+/** @internal */
+export declare const linearAccel: (x: number) => number;
+/** One classical RK4 step of x″ = accel(x). @internal */
+export declare function rk4Step(x: number, v: number, h: number, accel: (x: number) => number): [number, number];
+/** Zero-crossing times of a release from rest at `theta0`, over `periods` of T0 = 1. @internal */
+export declare function crossingTimes(accel: (x: number) => number, theta0: number, stepsPerPeriod: number, periods: number): number[];
+/** The phase at `t`, interpolated through the anchors (0, 0) and (c_k, π/2 + kπ). @internal */
+export declare function phaseAt(crossings: readonly number[], t: number): number;
+/**
+ * θ0 in [0, π) whose exact relative period error is `eps`, by bisection on the
+ * strictly increasing `pendulumPeriodErrorAt`.
+ *
+ * @returns NaN for a negative or non-finite `eps`.
+ * @internal
+ */
+export declare function theta0OfPeriodError(eps: number): number;
+/**
+ * The nome q = exp(−π K′/K) of the modulus k = sin(θ0/2), with
+ * K′/K = AGM(1, cos(θ0/2)) / AGM(1, sin(θ0/2)).
+ * @internal
+ */
+export declare function pendulumNome(theta0: number): number;
+//# sourceMappingURL=pendulum-motion.d.ts.map

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

@@ -0,0 +1 @@
+{"version":3,"file":"pendulum-motion.d.ts","sourceRoot":"","sources":["../../../src/atlas/oscillators/pendulum-motion.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;GAUG;AAIH,+BAA+B;AAC/B,eAAO,MAAM,EAAE,QAAqB,CAAC;AACrC,gBAAgB;AAChB,eAAO,MAAM,aAAa,MAAO,MAAM,KAAG,MAA2B,CAAC;AACtE,gBAAgB;AAChB,eAAO,MAAM,WAAW,MAAO,MAAM,KAAG,MAAiB,CAAC;AAE1D,yDAAyD;AACzD,wBAAgB,OAAO,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,KAAK,EAAE,CAAC,CAAC,EAAE,MAAM,KAAK,MAAM,GAAG,CAAC,MAAM,EAAE,MAAM,CAAC,CAUvG;AA4BD,kGAAkG;AAClG,wBAAgB,aAAa,CAAC,KAAK,EAAE,CAAC,CAAC,EAAE,MAAM,KAAK,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,cAAc,EAAE,MAAM,EAAE,OAAO,EAAE,MAAM,GAAG,MAAM,EAAE,CAY7H;AAED,+FAA+F;AAC/F,wBAAgB,OAAO,CAAC,SAAS,EAAE,SAAS,MAAM,EAAE,EAAE,CAAC,EAAE,MAAM,GAAG,MAAM,CAWvE;AAED;;;;;;GAMG;AACH,wBAAgB,mBAAmB,CAAC,GAAG,EAAE,MAAM,GAAG,MAAM,CAYvD;AAED;;;;GAIG;AACH,wBAAgB,YAAY,CAAC,MAAM,EAAE,MAAM,GAAG,MAAM,CAInD"}

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"}