drop zeros when reconstructing the solution

scripts/G₂_Adj.jl

@@ -82,61 +82,17 @@ unit = Δ
     show_progress = true,
 )
 
-warm = nothing
+solve_in_loop(
-
+    model,
-let status = JuMP.OPTIMIZE_NOT_CALLED, warm = warm, eps = 1e-9
+    wd,
-    certified, λ = false, 0.0
+    varP;
-    while status ≠ JuMP.OPTIMAL
+    logdir = "./log/G2/r=$HALFRADIUS/Adj-$(UPPER_BOUND)Δ",
-        @time status, warm = PropertyT.solve(
+    optimizer = scs_optimizer(;
-            model,
+        linear_solver = SCS.MKLDirectSolver,
-            scs_optimizer(;
+        eps = 1e-9,
-                linear_solver = SCS.MKLDirectSolver,
+        max_iters = 100_000,
-                eps = eps,
+        accel = 50,
-                max_iters = 100_000,
+        alpha = 1.95,
-                accel = 50,
+    ),
-                alpha = 1.95,
+    data = (elt = elt, unit = unit, halfradius = HALFRADIUS),
-            ),
+)
-            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),
-# )

scripts/utils.jl

@@ -11,11 +11,17 @@ function get_solution(model)
     return solution
 end
 
-function get_solution(model, wd, varP)
+function get_solution(model, wd, varP, eps = 1e-10)
     λ = JuMP.value(model[:λ])
 
-    Qs = [real.(sqrt(JuMP.value.(P))) for P in varP]
+    @info "reconstructing the solution"
-    Q = PropertyT.reconstruct(Qs, wd)
+    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
+
     solution = Dict(:λ => λ, :Q => Q)
 
     return solution