From 26c7c7df7c6aca6dae6ec2c9c13747be9d968ec8 Mon Sep 17 00:00:00 2001 From: Yuren Hao Date: Sat, 25 Jul 2026 14:10:23 -0500 Subject: =?UTF-8?q?RESULT=2061b:=20=E5=93=A8=E6=88=9221=E8=BF=9E=E6=8A=A5?= =?UTF-8?q?=E8=AD=A6=3D=E5=9C=B0=E6=9D=BF=E5=99=AA=E5=A3=B0=E5=81=87?= =?UTF-8?q?=E9=98=B3=E6=80=A7(resK=3D1.4e-6=E5=AE=8C=E7=BE=8E=E6=94=B6?= =?UTF-8?q?=E6=95=9B,rho=E4=BC=AA=E8=AF=BB=E6=95=B0),R61'=E5=90=88?= =?UTF-8?q?=E8=A7=84=E8=A1=80=E7=BB=9F=E5=A4=A9=E8=8A=B1=E6=9D=BF=E6=B2=89?= =?UTF-8?q?=E5=88=B01.5e-3'=E6=92=A4=E5=9B=9E;=20b15=E7=9C=9F=E5=AE=9E?= =?UTF-8?q?=E9=99=90=E5=88=B6=E5=9B=A0=E5=AD=90=3D=E6=9C=AA=E7=9F=A5,?= =?UTF-8?q?=E6=8A=BC=E5=90=8EBP=E9=94=9A=E7=82=B9;=20=E6=8E=A2=E9=92=88?= =?UTF-8?q?=E5=88=A4=E5=86=B3=E8=A7=84=E5=88=99=E4=BF=AE=E5=A4=8D(rho>1?= =?UTF-8?q?=E4=BB=85=E5=9C=A8resK=E7=A6=BB=E5=9C=B0=E6=9D=BF=E6=97=B6?= =?UTF-8?q?=E8=AE=A1=E6=95=B0)?= MIME-Version: 1.0 Content-Type: text/plain; charset=UTF-8 Content-Transfer-Encoding: 8bit Co-Authored-By: Claude Fable 5 Claude-Session: https://claude.ai/code/session_014FAPDWQ49M5Ye3NpTndTpn --- ep_run/probe_cx3_rhogrid.py | 4 +++- 1 file changed, 3 insertions(+), 1 deletion(-) (limited to 'ep_run') diff --git a/ep_run/probe_cx3_rhogrid.py b/ep_run/probe_cx3_rhogrid.py index 7c27a2d..e154fec 100644 --- a/ep_run/probe_cx3_rhogrid.py +++ b/ep_run/probe_cx3_rhogrid.py @@ -123,7 +123,9 @@ for spec in a.ckpts.split(','): if res_list[i - 1] > 1e-7] tail = ratios[-6:] if len(ratios) >= 6 else ratios rho_tail = float(np.median(tail)) if tail else float('nan') - diverged = (not np.isfinite(res_list[-1])) or res_list[-1] > 1e-2 or rho_tail > 1.0 + # resK-primary verdict: rho~1 AT THE FP FLOOR is noise/noise (documented artifact), + # not divergence — rho only counts when the residual is materially off-floor + diverged = (not np.isfinite(res_list[-1])) or res_list[-1] > 1e-2 or (rho_tail > 1.0 and res_list[-1] > 1e-4) verdict = 'DIVERGE' if diverged else 'converge' if diverged and bstar is None: bstar = beta print(f'{beta:>7.3f} {res_list[0]:>10.2e} {(res_list[1] if len(res_list)>1 else float("nan")):>10.2e} ' -- cgit v1.2.3