Reproduce

One script recomputes every number on this site from the raw data files. Python 3 standard library only; no network, no random numbers, no clock. Two runs give byte-identical output.

Steps

Download the site folder (or the files under data/ and reproduce/).

In that folder run:

python3 -I reproduce/reproduce.py

For the machine-readable form:

python3 -I reproduce/reproduce.py --json

Compare with the expected output below. Its SHA-256:

255595810f738b36c2c227c2eb9f4860acbf5155061b56e73745163af60c9b9f

Check with:

python3 -I reproduce/reproduce.py | shasum -a 256

Script: reproduce/reproduce.py, sha256 74c1c62a20f7e927138053b58211660d81f0deb24754003dcfce0c48c03f82f9. Data hashes: Raw data. Constants: CODATA 2018 and the IAU solar mass parameter, written out at the top of the script.

Scope: the script recomputes the published numbers with standard formulas. The lab's own working scripts are not published, because they name the lab's private relations; every number they produced for this site is recomputed here from the same data.

Expected output

REPRODUCE: every number on the site, recomputed from data/
E1  alpha = 0 (clock rate follows the potential)                   residuals [-0.13, 0.0, 0.08]           chi2 0.02318    -> agrees
E2  gamma - 1 = 0 (full bending, general relativity)               residuals [0.91, -0.67]                chi2 1.278      -> agrees
E2  gamma - 1 = -1 (clock-rate half only, Soldner's 0.876 arcsec)  residuals [43479.17, 8332.67]          chi2 1.96e+09   -> contradicted
E3  Hafele and Keating's own prediction from the flight logs       residuals [-0.76, -0.09]               chi2 0.5821     -> agrees
E4  gamma x tau0 (time dilation)                                   residuals [-0.14, -0.94]               chi2 1.26       -> agrees
E5  S = 2 sqrt 2 at every apparatus (ideal analysers)              residuals [-8.76, -4.92, -2.04, -5.09] chi2 131.1      -> contradicted
E5  S at most 2 (each end holds its own preset value)              residuals [46.47, 36.5, 2.1, 4.11]     chi2 3513       -> contradicted
E5  Aspect 1982 against its own apparatus prediction 2.70 +/- 0.05 residuals [-0.06]                      chi2 0.003303   -> agrees
E6  return at 720 degrees (two-valued wave)                        residuals [-0.84]                      chi2 0.7091     -> agrees
E6  return at 360 degrees (ordinary vector, rotating body)         residuals [93.89]                      chi2 8816       -> contradicted
E7  general relativity, 6 pi GM / (c^2 a (1 - e^2)) per orbit      residuals [-0.02]                      chi2 0.0002854  -> agrees
E7  weak-field clock-rate map alone (no advance)                   residuals [1074.5]                     chi2 1.155e+06  -> contradicted
E8  Bohr ladder with reduced mass                                  residuals [2298950011.7]               chi2 5.285e+18  -> contradicted
E8  Dirac levels with reduced mass (no QED)                        residuals [-710596536.45]              chi2 5.049e+17  -> contradicted
E9  antimatter falls at 1 g                                        residuals [-1.21]                      chi2 1.473      -> agrees
E9  antimatter rises at -1 g                                       residuals [8.5]                        chi2 72.17      -> contradicted
E5  bound 2 < S <= 2 sqrt 2 -> agrees
E3-part-A  scored 29, agree 16, contradicted 12, open 1, recalled 1
  T1   +0.314    agrees (data recalled from a secondary source, not counted)
  T2   -0.478    agrees
  T3   +0.154    agrees
  T4   +0.909    agrees
  T5   +0.81     agrees
  T6   -         agrees (inside the published bound)
  T7   -0.107    agrees
  T8a  +5.17     contradicted
  T8b  -0.497    agrees
  T9a  +19.9     contradicted
  T9b  -0.12     agrees
  T10  -1        agrees
  T11  -0.556    agrees
  T12  -0.587    agrees
  C1   -         agrees (inside the published bound)
  C2   +0        agrees
  C3   -         open
  S1   -3.99e+07 contradicted
  S2   -225      contradicted
  S3   -1.35e+03 contradicted
  S4   -3.94e+03 contradicted
  K1   +2.45e+03 contradicted
  K1d  -451      contradicted
  K2   +1.26e+03 contradicted
  K2d  -205      contradicted
  K3   +0        agrees
  K4   +0        agrees
  K5   -         contradicted (exact integers differ)
  K6   -         contradicted (exact integers differ)
  N1   -         agrees (inside the published bound)
WAVE  S: two ends, probabilistic 1.4142; two ends, sign 2.0000; read once for both ends 2.8284; ceiling 2.8284
P01  z from 720: -0.84; z from 360: +93.9  [10.1007/BF01321193]
P02  46.5 sigma above the preset-answer limit 2  [10.1103/PhysRevLett.49.91]
P03  2.1 sigma above 2  [10.1038/nature15759]
P04  4.1 sigma above 2  [10.1126/science.aan3211]
P05  a three-direction inverse square predicts no deviation: inside the bound  [10.1103/PhysRevLett.124.101101]
P06  no excitation seen: a fixed mass inside the range is excluded  [10.1103/PhysRevD.98.092001]
P07  combined +0.07 sigma from 0; chi2 = 0.0232, p = 0.999  [10.1103/PhysRevLett.121.231101]
P08  full bending +0.77 sigma; half bending +44270 sigma  [10.1038/nature01997]
P09  against 64.4309 us: residuals [-0.14, -0.94]  [10.1038/268301a0]
P10  residuals against the flight-log predictions [-0.76, -0.09]  [10.1126/science.177.4044.168]
P11  falls at 1 g: -1.21 sigma; rises at -1 g: +8.5 sigma  [10.1038/s41586-023-06527-1]
P12  no rotation term predicts 0: +5.17 sigma  [10.1103/PhysRevLett.106.221101]
P13  beta = 1: -0.69 sigma  [10.3847/1538-3881/aa5be2]
P14  Bohr ladder misses by 2.299e+10 Hz; Dirac levels by -7.106e+09 Hz  [10.1103/PhysRevLett.107.203001]
P15  universal free fall predicts 0: -0.56 sigma  [10.1103/PhysRevLett.129.121102]
P16  ratio 1: -0.59 sigma  [10.1103/PhysRevX.11.041050]
P17  against 1.09e-19 per mm: -0.48 sigma  [10.1038/s41586-021-04349-7]
P18  against 633.1 m/s: +0.81 sigma  [10.1051/0004-6361/202038937]