hsdp_descent_wrapper.jl

# First-order steepest descent potential reduction for HSDP

include("hsdp_descent.jl")

m = 8
n = 10
𝜀 = 0.5
max_iter = 100
converged_vec = zeros(max_iter)
reached_epsilon_vec = zeros(max_iter)
steps_vec = zeros(max_iter)
steps_limit_vec = zeros(max_iter)
cond_vec = zeros(max_iter)

for i = 1:max_iter
    println(i)
    converged_vec[i], reached_epsilon_vec[i], steps_vec[i], steps_limit_vec[i], cond_vec[i] = hsdp_descent(m,n,𝜀)
end

writedlm("hsdp_descent_results.txt", [converged_vec reached_epsilon_vec steps_vec steps_limit_vec cond_vec])