[{"type":"text","text":"Full smoke test at n=300 with validation against E-015."},{"type":"tool_use","id":"Bash_65","name":"Bash","input":{"command":"cd /d/Python/Quantum && ./.venv/Scripts/python.exe - <<'EOF'\nsrc = open(\"research/numerics/scripts/E-021-sinf-sf-curve.py\").read()\nsrc = src.replace('if __name__ == \"__main__\":', 'if False:')\nns = {\"__file__\": \"research/numerics/scripts/E-021-sinf-sf-curve.py\"}\nexec(compile(src, \"e21\", \"exec\"), ns)\nimport numpy as np\nzs = np.load(\"research/numerics/data/E-021-zeros-cache.npy\")\nzeta_str, fd = ns[\"get_zeta_strings\"](300)\neta_str = ns[\"compute_eta\"](zeta_str, fd, zs)\nS_inf = ns[\"compute_S_inf\"](300, zeta_str, fd)\nS_f = ns[\"compute_Sf\"](300, eta_str)\nimport mpmath as mp\nmp.mp.dps = 80\ne015 = {}\nfor line in open(\"research/numerics/data/E-015-li-coefficients-n3000.csv\"):\n    p = line.strip().split(\",\")\n    if p[0].isdigit(): e015[int(p[0])] = mp.mpf(p[1])\nworst = 0; wn = 0\nfor n in sorted(e015):\n    if n > 300: break\n    lam = S_inf[n] - S_f[n] + 1\n    dev = abs(lam - e015[n]) / abs(e015[n])\n    if float(dev) > worst: worst = float(dev); wn = n\nprint(\"worst rel dev vs E-015 (n<=300):\", worst, \"at n =\", wn)\nprint(\"S_f(200) =\", mp.nstr(S_f[200], 15), \" (expect -2.14037668...)\")\nprint(\"S_f(300) =\", mp.nstr(S_f[300], 15))\nEOF","description":"Smoke test v3 eta at n=300","timeout":600000},"extra_content":null}]