From 5aa48dcf8047cbd53dafb0c377e260629675ae9e Mon Sep 17 00:00:00 2001
From: Emily Boudreaux <thomas@boudreauxmail.com>
Date: Fri, 18 Sep 2026 10:07:46 -0400
Subject: [PATCH] Keep coupled evidence audit linear in state size and repeat
 raw timing

---
 experiments/check_coupled_trial.py            |  3 +-
 .../coupled_trial_raw_repeat_2026-09-18.json  | 47 +++++++++++++++++++
 2 files changed, 49 insertions(+), 1 deletion(-)
 create mode 100644 experiments/coupled_trial_raw_repeat_2026-09-18.json

diff --git a/experiments/check_coupled_trial.py b/experiments/check_coupled_trial.py
index 1194394..66ac2cc 100644
--- a/experiments/check_coupled_trial.py
+++ b/experiments/check_coupled_trial.py
@@ -79,6 +79,7 @@ def inspect_case_artifacts(folder):
     check('base_physical_normalization',near_vector(x,[v*a for v,a in zip(physical,sf)],norm(x)))
     layout=rows(folder/(P+'layout.csv'));vl=layout_for(layout,'value',size);rl=layout_for(layout,'residual',size)
     gi=[i for name,(start,n) in vl.items() if name in GRAVITY_VALUES for i in range(start,start+n)]
+    gravity_value_indices=set(gi)
     gr=[i for name,(start,n) in rl.items() if name in GRAVITY_ROWS for i in range(start,start+n)]
     f0,pn,jn=norm(f),norm(p),norm(j)
     check('finite_positive_base_residual',math.isfinite(f0) and f0>0)
@@ -159,7 +160,7 @@ def inspect_case_artifacts(folder):
             gc=get('gravity_correction'+suffix,save_gravity)
             if save_gravity:
                 check(label+'gravity_partition_present',bool(gi) and bool(gr))
-                check(label+'gravity_only_correction',len(gc)==size and all(v==0 for k,v in enumerate(gc) if k not in set(gi)))
+                check(label+'gravity_only_correction',len(gc)==size and all(v==0 for k,v in enumerate(gc) if k not in gravity_value_indices))
                 if flags['before_prepared']:
                     check(label+'gravity_displacement',near_vector(vectors['before'][0],[v+w for v,w in zip(raw_state,gc)],norm(raw_state)+norm(gc)))
                     target=norm([vectors['before'][1][k]-prediction[k] for k in gr])
diff --git a/experiments/coupled_trial_raw_repeat_2026-09-18.json b/experiments/coupled_trial_raw_repeat_2026-09-18.json
new file mode 100644
index 0000000..f075eca
--- /dev/null
+++ b/experiments/coupled_trial_raw_repeat_2026-09-18.json
@@ -0,0 +1,47 @@
+[
+  {
+    "id": "late_raw_coupled",
+    "arguments": [
+      "--index",
+      "1.5",
+      "--angular-momentum",
+      "0",
+      "--no-central",
+      "--max-newton",
+      "0",
+      "--coupled-trial-probe",
+      "--strategy",
+      "coupled",
+      "--probe-state",
+      "/Users/tboudreaux/.codex/.chatgpt-projects/g-p-6a6b3f2194ec8191bebcd59075897f1d/diagnostics/nonlinear_step_2026-09-18/late_v1/n1p5_late_eta003/nonlinear_step_physical_base.txt",
+      "--probe-direction",
+      "/Users/tboudreaux/.codex/.chatgpt-projects/g-p-6a6b3f2194ec8191bebcd59075897f1d/diagnostics/nonlinear_step_2026-09-18/late_v1/n1p5_late_eta003/nonlinear_step_direction_full.txt",
+      "--linear-tol",
+      "0.03",
+      "--restart",
+      "160",
+      "--preconditioner-variant",
+      "gravity-schur",
+      "--schur-sweeps",
+      "1",
+      "--schur-relaxation",
+      "0.5",
+      "--gravity-tol",
+      "1e-10",
+      "--gravity-iterations",
+      "300",
+      "--gravity-restart",
+      "80",
+      "--coupled-tol",
+      "0.03",
+      "--coupled-iterations",
+      "160",
+      "--coupled-restart",
+      "160",
+      "--coupled-max-ratio",
+      "1",
+      "--coupled-attempts",
+      "1"
+    ]
+  }
+]