universal-physics-tensor 0.47.1 → 1.0.1

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

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

@@ -0,0 +1,254 @@
+/**
+ * `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 { pendulumPeriodErrorAt } from './bridges-limits.js';
+import { linearAccel, pendulumAccel, pendulumNome, rk4Step, theta0OfPeriodError } from './pendulum-motion.js';
+const TWO_PI = 2 * Math.PI;
+/** 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);
+    return ((16 / 3) * q ** 1.5) / (1 - q);
+}
+/** The pendulum's fundamental Fourier coefficient a₀ = 8√q/(1 + q). @internal */
+export function fundamentalCoefficient(theta0) {
+    const q = pendulumNome(theta0);
+    return (8 * Math.sqrt(q)) / (1 + q);
+}
+const lead = (eps, tOverT0) => (TWO_PI * tOverT0 * eps) / (1 + eps);
+/** B(t/T0) at θ0 recovered from ε; `drift` and `waveform` scale its two terms (1, 1 is B itself). */
+function positionBound(eps, drift = 1, waveform = 1) {
+    const theta0 = theta0OfPeriodError(eps);
+    const floor = waveform * waveformBound(theta0);
+    return (tOverT0) => floor + drift * 2 * theta0 * Math.sin(Math.min(lead(eps, tOverT0), Math.PI) / 2);
+}
+function positionHorizon(eps, tolerance, T0) {
+    if (T0 === undefined || !Number.isFinite(T0) || T0 <= 0 || !Number.isFinite(eps) || eps < 0)
+        return Number.NaN;
+    if (eps === 0)
+        return Infinity;
+    const theta0 = theta0OfPeriodError(eps);
+    const floor = waveformBound(theta0);
+    if (tolerance < floor)
+        return Number.NaN;
+    const r = (tolerance - floor) / (2 * theta0);
+    if (r >= 1)
+        return Infinity;
+    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;
+/** 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.
+ */
+function dominanceSamples(theta0, periods, stepsPerPeriod, bound) {
+    const stride = stepsPerPeriod / SAMPLES_PER_PERIOD;
+    const h = 1 / stepsPerPeriod;
+    const measured = [];
+    const b = [];
+    let x = theta0;
+    let v = 0;
+    let y = theta0;
+    let u = 0;
+    const steps = Math.floor(periods * SAMPLES_PER_PERIOD) * stride;
+    for (let i = 1; i <= steps; i++) {
+        [x, v] = rk4Step(x, v, h, pendulumAccel);
+        [y, u] = rk4Step(y, u, h, linearAccel);
+        if (i % stride === 0) {
+            measured.push(Math.abs(x - y));
+            b.push(bound(i * h));
+        }
+    }
+    return { measured, bound: b };
+}
+/** The window a dominance witness covers: up to antiphase (Δφ = π), past which B ≥ 2θ0 holds trivially, capped. */
+function dominanceWindow(eps) {
+    return Math.min(POSITION_POINT_MAX_PERIODS, (1 + eps) / (2 * eps));
+}
+function dominanceSpec(id, eps, drift, waveform) {
+    const theta0 = theta0OfPeriodError(eps);
+    const periods = dominanceWindow(eps);
+    return {
+        id,
+        kind: 'dominance',
+        evaluate: (stepsPerPeriod) => dominanceSamples(theta0, periods, stepsPerPeriod, positionBound(eps, drift, waveform)),
+        coarseResolution: 200,
+        fineResolution: 800,
+    };
+}
+/** a₀ measured by projecting one RK4 pendulum period onto cos φ (trapezoid, exact for a periodic integrand). */
+export function measureFundamental(theta0, stepsPerPeriod) {
+    const T = 1 + pendulumPeriodErrorAt({ theta0 });
+    const h = T / stepsPerPeriod;
+    let x = theta0;
+    let v = 0;
+    let sum = theta0;
+    for (let i = 1; i < stepsPerPeriod; i++) {
+        [x, v] = rk4Step(x, v, h, pendulumAccel);
+        sum += x * Math.cos((TWO_PI * i) / stepsPerPeriod);
+    }
+    return (2 * sum) / stepsPerPeriod;
+}
+/** The θ0 at which W7xs checks the Fourier series: a grid across the bridge's domain θ0 ≤ 0.5. @internal */
+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;
+/**
+ * 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.
+ */
+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);
+    }
+    return s;
+}
+/** 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;
+    for (const theta0 of FOURIER_SERIES_THETA0S) {
+        const h = (1 + pendulumPeriodErrorAt({ theta0 })) / stepsPerPeriod;
+        let x = theta0;
+        let v = 0;
+        worst = Math.max(worst, Math.abs(x - fourierSeries(theta0, 0, variant)));
+        for (let i = 1; i <= stepsPerPeriod; i++) {
+            [x, v] = rk4Step(x, v, h, pendulumAccel);
+            worst = Math.max(worst, Math.abs(x - fourierSeries(theta0, (TWO_PI * i) / stepsPerPeriod, variant)));
+        }
+    }
+    return worst;
+}
+/**
+ * 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 function fourierSeriesCheck(id, variant = 'plus') {
+    return {
+        id,
+        evaluate: (stepsPerPeriod) => fourierSeriesDeviation(stepsPerPeriod, variant),
+        target: 0,
+        coarseResolution: 200,
+        fineResolution: 800,
+        tolerance: FOURIER_SERIES_TOLERANCE,
+    };
+}
+const W7X_THETA0 = 0.3;
+const W7X_EPS = pendulumPeriodErrorAt({ theta0: W7X_THETA0 });
+const W7XA = {
+    id: 'W7xa',
+    evaluate: (stepsPerPeriod) => measureFundamental(W7X_THETA0, stepsPerPeriod),
+    target: fundamentalCoefficient(W7X_THETA0),
+    coarseResolution: 200,
+    fineResolution: 800,
+    tolerance: 1e-10,
+};
+function positionPointCheck(eps) {
+    if (!Number.isFinite(eps) || eps <= 0)
+        return { unavailable: `θ0 = 0 at ε = ${eps}: both motions are at rest and the error is 0` };
+    const theta0 = theta0OfPeriodError(eps);
+    const periods = dominanceWindow(eps);
+    const capped = periods === POSITION_POINT_MAX_PERIODS ? `, capped at ${POSITION_POINT_MAX_PERIODS} T0 (antiphase is at ${(1 + eps) / (2 * eps)} T0)` : ' (to antiphase, Δφ = π)';
+    return {
+        claim: `RK4 |θ − θ_lin| ≤ B(t) at every sample, θ0 = ${Number(theta0.toPrecision(8))} (recovered from ε), over t ≤ ${periods} T0${capped}; ` +
+            `200 and 800 steps per T0, sampled 200 per T0; the coarse-to-fine difference is the uncertainty`,
+        check: dominanceSpec('W7x@point', eps, 1, 1),
+        control: dominanceSpec('W7x@point-control', eps, 1, 0),
+        controlClaim: 'the same bound without its waveform term W̄, 2θ0 sin(Δφ/2) alone',
+    };
+}
+/** @internal */
+export const PENDULUM_POSITION_TRANSLATION = {
+    bridgeId: 'ab-pendulum-linear',
+    observable: 'position',
+    unit: 'rad',
+    fromNorm: 'relative period error, normalized by the value of the reduced model',
+    errorKind: 'upper-bound',
+    definition: 'the pointwise angle error |θ(t) − θ_lin(t)| of the two motions released from rest at θ0',
+    derivation: 'an UPPER BOUND B(t) = W̄ + 2θ0 sin(min(Δφ(t), π)/2), Δφ = 2π (t/T0) ε/(1+ε) the exact phase lead and ' +
+        'W̄ = (16/3) q^{3/2}/(1−q) ≈ θ0³/96 the waveform term (q the nome of sin(θ0/2); θ0 is recovered from ε). ' +
+        '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)|. ' +
+        'B 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',
+        '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`,
+    ],
+    parameters: ['T0'],
+    timeUnit: 'the unit of T0',
+    errorAt: (eps, t, { T0 }) => {
+        if (T0 === undefined || !Number.isFinite(T0) || T0 <= 0 || !Number.isFinite(eps) || eps < 0 || !Number.isFinite(t))
+            return Number.NaN;
+        return positionBound(eps)(t / T0);
+    },
+    horizonFor: (eps, tolerance, { T0 }) => positionHorizon(eps, tolerance, T0),
+    floor: (eps) => waveformBound(theta0OfPeriodError(eps)),
+    boundary: '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)',
+    notCovered: ['the phase (declared separately: --tolerance=phase:EPS)', 'the angular velocity', 'the amplitude'],
+    witnesses: [
+        {
+            id: 'W7x',
+            kind: 'numeric',
+            test: 'tests/atlas/pendulum-position-translation.test.ts',
+            tolerance: `RK4 |θ − θ_lin| ≤ B(t) at every sample to antiphase, θ0 = ${W7X_THETA0}, by more than the coarse-to-fine difference`,
+        },
+        {
+            id: 'W7xa',
+            kind: 'numeric',
+            test: 'tests/atlas/pendulum-position-translation.test.ts',
+            tolerance: `RK4-projected fundamental coefficient at θ0 = ${W7X_THETA0} within 1e-10 of 8√q/(1+q)`,
+        },
+        {
+            id: 'W7xs',
+            kind: 'numeric',
+            test: 'tests/atlas/pendulum-position-translation.test.ts',
+            tolerance: `the Fourier series (${FOURIER_SERIES_TERMS} terms) within ${FOURIER_SERIES_TOLERANCE} rad of RK4 at every step of one period, ` +
+                `at θ0 ∈ {${FOURIER_SERIES_THETA0S.join(', ')}}`,
+        },
+    ],
+    checks: [dominanceSpec('W7x', W7X_EPS, 1, 1), W7XA, fourierSeriesCheck('W7xs')],
+    pointCheck: positionPointCheck,
+};
+//# sourceMappingURL=position-translation.js.map

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

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