[{"type":"text","text":"Quick smoke test at small N_MAX before the full detached run."},{"type":"tool_use","id":"Bash_43","name":"Bash","input":{"command":"cd /d/Python/Quantum && ./.venv/Scripts/python.exe - <<'EOF'\nimport sys\nsys.argv = [\"x\", \"--force\"]\nimport importlib.util\nspec = importlib.util.spec_from_file_location(\n    \"e21\", \"research/numerics/scripts/E-021-sinf-sf-curve.py\")\nm = importlib.util.module_from_spec(spec)\nspec.loader.exec_module.__self__ if False else None\n# load module without running main\nimport types\nsrc = open(\"research/numerics/scripts/E-021-sinf-sf-curve.py\").read()\nsrc = src.replace('if __name__ == \"__main__\":', 'if False:')\nns = {}\nexec(compile(src, \"e21\", \"exec\"), ns)\nns[\"N_MAX\"] = 300\nzeta_str, fd = ns[\"get_zeta_strings\"](300)\neta_str = ns[\"compute_eta\"](zeta_str, fd)\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\n# validate vs E-015\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\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    worst = max(worst, float(dev))\nprint(\"worst rel dev vs E-015 (n<=300):\", worst)\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 fixed P2 at n=300","timeout":600000},"extra_content":null}]