From 263af398eb99b08bac90f4e9c276df8bfa707f98 Mon Sep 17 00:00:00 2001
From: Marek Kaluba <kalmar@amu.edu.pl>
Date: Tue, 11 Apr 2023 10:31:54 +0200
Subject: [PATCH] =?UTF-8?q?G=E2=82=82=5FAdj:=20bring=20back=20solve=5Fin?=
 =?UTF-8?q?=5Floop?=
MIME-Version: 1.0
Content-Type: text/plain; charset=UTF-8
Content-Transfer-Encoding: 8bit

---
 scripts/G₂_Adj.jl | 72 +++++++++------------------------------------
 1 file changed, 14 insertions(+), 58 deletions(-)

diff --git a/scripts/G₂_Adj.jl b/scripts/G₂_Adj.jl
index e9acfc9..8a499d0 100644
--- a/scripts/G₂_Adj.jl
+++ b/scripts/G₂_Adj.jl
@@ -82,61 +82,17 @@ unit = Δ
     show_progress = true,
 )
 
-warm = nothing
-
-let status = JuMP.OPTIMIZE_NOT_CALLED, warm = warm, eps = 1e-9
-    certified, λ = false, 0.0
-    while status ≠ JuMP.OPTIMAL
-        @time status, warm = PropertyT.solve(
-            model,
-            scs_optimizer(;
-                linear_solver = SCS.MKLDirectSolver,
-                eps = eps,
-                max_iters = 100_000,
-                accel = 50,
-                alpha = 1.95,
-            ),
-            warm,
-        )
-
-        @info "reconstructing the solution"
-        Q = @time let wd = wd, Ps = [JuMP.value.(P) for P in varP], eps = eps
-            PropertyT.__droptol!.(Ps, 100eps)
-            Qs = real.(sqrt.(Ps))
-            PropertyT.__droptol!.(Qs, eps)
-
-            PropertyT.reconstruct(Qs, wd)
-        end
-
-        @info "certifying the solution"
-        @time certified, λ = PropertyT.certify_solution(
-            elt,
-            unit,
-            JuMP.objective_value(model),
-            Q;
-            halfradius = HALFRADIUS,
-            augmented = true,
-        )
-    end
-
-    if certified && λ > 0
-        Κ(λ, S) = round(sqrt(2λ / length(S)), Base.RoundDown; digits = 5)
-        @info "Certified result: $G has property (T):" N λ Κ(λ, S)
-    else
-        @info "Could NOT certify the result:" certified λ
-    end
-end
-
-# solve_in_loop(
-#     model,
-#     wd,
-#     varP;
-#     logdir = "./log/G2/r=$HALFRADIUS/Adj-InfΔ",
-#     optimizer = scs_optimizer(;
-#         eps = 1e-10,
-#         max_iters = 50_000,
-#         accel = 50,
-#         alpha = 1.95,
-#     ),
-#     data = (elt = elt, unit = unit, halfradius = HALFRADIUS),
-# )
+solve_in_loop(
+    model,
+    wd,
+    varP;
+    logdir = "./log/G2/r=$HALFRADIUS/Adj-$(UPPER_BOUND)Δ",
+    optimizer = scs_optimizer(;
+        linear_solver = SCS.MKLDirectSolver,
+        eps = 1e-9,
+        max_iters = 100_000,
+        accel = 50,
+        alpha = 1.95,
+    ),
+    data = (elt = elt, unit = unit, halfradius = HALFRADIUS),
+)