summaryrefslogtreecommitdiff
path: root/ep_run
diff options
context:
space:
mode:
authorYuren Hao <yurenh2@illinois.edu>2026-07-25 14:10:23 -0500
committerYuren Hao <yurenh2@illinois.edu>2026-07-25 14:10:23 -0500
commit26c7c7df7c6aca6dae6ec2c9c13747be9d968ec8 (patch)
tree7dd7dd4ffe34da9c23eb920182951f3885f745df /ep_run
parent1c89d1b5bee0349aa6a2f39f44ac68e751562a94 (diff)
RESULT 61b: 哨戒21连报警=地板噪声假阳性(resK=1.4e-6完美收敛,rho伪读数),R61'合规血统天花板沉到1.5e-3'撤回; b15真实限制因子=未知,押后BP锚点; 探针判决规则修复(rho>1仅在resK离地板时计数)
Co-Authored-By: Claude Fable 5 <noreply@anthropic.com> Claude-Session: https://claude.ai/code/session_014FAPDWQ49M5Ye3NpTndTpn
Diffstat (limited to 'ep_run')
-rw-r--r--ep_run/probe_cx3_rhogrid.py4
1 files changed, 3 insertions, 1 deletions
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} '