tools(ra1n): add guarded margin value setter

Author: inaciovasquez2020

tools/set_ra1n_margin_values.py

@@ -0,0 +1,56 @@
+#!/usr/bin/env python3
+import argparse
+import json
+import math
+from pathlib import Path
+
+p = Path("artifacts/ra1n_margin/ra1n_margin_certificate.json")
+data = json.loads(p.read_text(encoding="utf-8"))
+
+ap = argparse.ArgumentParser()
+ap.add_argument("--M-phi", type=float, required=True)
+ap.add_argument("--m-phi", type=float, required=True)
+ap.add_argument("--c-fin", type=float, required=True)
+ap.add_argument("--epsilon", type=float, required=True)
+ap.add_argument("--eta-supp", type=float, required=True)
+ap.add_argument("--eta-err", type=float, required=True)
+args = ap.parse_args()
+
+delta = math.sqrt(1.0 - args.eta_supp) * args.epsilon - args.eta_err
+
+if not args.m_phi > 0:
+    raise SystemExit("m_phi must be positive")
+if not args.c_fin > 0:
+    raise SystemExit("c_fin must be positive")
+if not args.epsilon > 0:
+    raise SystemExit("epsilon must be positive")
+if not (0 <= args.eta_supp < 1):
+    raise SystemExit("eta_supp must lie in [0,1)")
+if not args.eta_err >= 0:
+    raise SystemExit("eta_err must be nonnegative")
+if not delta > 0:
+    raise SystemExit(f"Delta_RA1n is not positive: {delta}")
+
+data["status"] = "CERTIFIED_RA1N_MARGIN_POSITIVE"
+
+data["cutoff_certificate"]["M_phi"] = args.M_phi
+data["cutoff_certificate"]["m_phi_lower_bound"] = args.m_phi
+data["cutoff_certificate"]["verified"] = True
+
+data["finite_projection_certificate"]["c_fin"] = args.c_fin
+data["finite_projection_certificate"]["verified"] = True
+
+data["angle_gap_certificate"]["epsilon"] = args.epsilon
+data["angle_gap_certificate"]["verified"] = True
+
+data["support_overlap_certificate"]["eta_supp"] = args.eta_supp
+data["support_overlap_certificate"]["verified"] = True
+
+data["error_certificate"]["eta_err"] = args.eta_err
+data["error_certificate"]["verified"] = True
+
+data["unified_margin"]["Delta_RA1n"] = delta
+data["unified_margin"]["verified_positive"] = True
+
+p.write_text(json.dumps(data, indent=2) + "\n", encoding="utf-8")
+print(f"RA1n certified margin set: Delta_RA1n={delta}")