Search conservative hard-iron boundary

Author: mgrouch

tools/ou_iii_mag_hard_iron_boundary.py

@@ -0,0 +1,105 @@
+#!/usr/bin/env python3
+"""Search the conservative edge of refinement-stage hard-iron shrinkage.
+
+The broad sweep showed that the lowest-information calm JONSWAP record is the
+first failure when every record is allowed to fit.  This follow-up holds the
+information floor at 0.25 (just above the observed minimum 0.244) and searches
+small shrinkage gains where a regularized estimate might help without turning
+that weak record on.
+"""
+from __future__ import annotations
+
+import ou_iii_mag_hard_iron_shrinkage as exp
+
+
+def install_patches():
+    originals = (
+        exp.FUSION.read_text(), exp.MAG_TUNER.read_text(),
+        exp.SIM.read_text(), exp.COMMON_H.read_text(),
+    )
+    f0, m0, s0, h0 = originals
+
+    m = exp.replace_once(m0, exp.MAG_CFG_OLD, exp.MAG_CFG_NEW, "MagAutoTuner config")
+    m = exp.replace_once(m, exp.MAG_SOLVE_OLD, exp.MAG_SOLVE_NEW, "MagAutoTuner shrink")
+    exp.MAG_TUNER.write_text(m)
+
+    f = exp.replace_once(f0, exp.CFG_OLD, exp.CFG_NEW, "OU-III config")
+    f = exp.replace_once(f, exp.BEGIN_OLD, exp.BEGIN_NEW, "initial hard-iron disable")
+    f = exp.replace_once(f, exp.REFINE_CFG_OLD, exp.REFINE_CFG_NEW, "refinement config")
+    f = exp.replace_once(f, exp.REFINE_USE_OLD, exp.REFINE_USE_NEW, "refinement hard-iron use")
+    f = exp.replace_once(f, exp.PUBLIC_OLD, exp.PUBLIC_NEW, "public diagnostics")
+    exp.FUSION.write_text(f)
+
+    sim = exp.replace_once(s0, exp.ENV_OLD, exp.ENV_NEW, "sim env")
+    sim = exp.replace_once(sim, exp.REFINE_DIAG_OLD, exp.REFINE_DIAG_NEW, "refinement diagnostics")
+    exp.SIM.write_text(sim)
+
+    h = h0
+    if "#include <cstdlib>" not in h:
+        h = h.replace("#include <cstdint>\n", "#include <cstdint>\n#include <cstdlib>\n", 1)
+    h = exp.replace_once(h, exp.SIG_OLD, exp.SIG_NEW, "mag sigma")
+    exp.COMMON_H.write_text(h)
+    return originals
+
+
+def restore(originals):
+    f0, m0, s0, h0 = originals
+    exp.FUSION.write_text(f0)
+    exp.MAG_TUNER.write_text(m0)
+    exp.SIM.write_text(s0)
+    exp.COMMON_H.write_text(h0)
+
+
+def main():
+    originals = install_patches()
+    try:
+        exp.build()
+        base, _ = exp.run()
+        ranking = []
+
+        # Search the transition where seven records are identifiable while the
+        # weakest calm-sea fit remains suppressed.  Include nearby thresholds
+        # because the second-weakest record may itself be harmful.
+        candidates = []
+        for info in (0.25, 0.50, 1.0):
+            for shrink in (0.005, 0.010, 0.020, 0.030, 0.040, 0.050, 0.075, 0.100):
+                candidates.append((info, shrink))
+
+        for info, shrink in candidates:
+            cand, hi = exp.run(
+                enable=True, min_info=info, shrink=shrink,
+                mag_scale=1.02, rs_coeff=0.3475,
+            )
+            mx, n, geo, rows = exp.report(
+                f"BOUNDARY min_info={info:.3f} shrink={shrink:.3f}",
+                base, cand, hi,
+            )
+            ranking.append((mx, geo, -n, info, shrink, rows, hi))
+            if mx < 1.0 and n == 56:
+                print(
+                    f"STRICT_WIN min_info={info:.6f} shrink={shrink:.6f} "
+                    f"max={mx:.9f} geo={geo:.9f}", flush=True,
+                )
+                print("STRICT_DOMINANCE=PASS", flush=True)
+                return
+
+        print("=== BOUNDARY RANKING ===", flush=True)
+        for row in sorted(ranking)[:15]:
+            print(
+                f"min_info={row[3]:.3f} shrink={row[4]:.3f} "
+                f"max={row[0]:.9f} geo={row[1]:.9f} improved={-row[2]}/56",
+                flush=True,
+            )
+        best = sorted(ranking)[0]
+        print(
+            f"BOUNDARY_BEST min_info={best[3]:.6f} shrink={best[4]:.6f} "
+            f"max={best[0]:.9f} geo={best[1]:.9f} improved={-best[2]}/56",
+            flush=True,
+        )
+        print("STRICT_DOMINANCE=FAIL", flush=True)
+    finally:
+        restore(originals)
+
+
+if __name__ == "__main__":
+    main()