2022-11-15 18:51:43 +01:00
|
|
|
## Low-level solve
|
|
|
|
|
|
|
|
setwarmstart!(model::JuMP.Model, ::Nothing) = model
|
|
|
|
|
|
|
|
function setwarmstart!(model::JuMP.Model, warmstart)
|
|
|
|
constraint_map = Dict(
|
|
|
|
ct => JuMP.all_constraints(model, ct...) for
|
|
|
|
ct in JuMP.list_of_constraint_types(model)
|
|
|
|
)
|
2023-03-19 20:33:28 +01:00
|
|
|
try
|
|
|
|
JuMP.set_start_value.(JuMP.all_variables(model), warmstart.primal)
|
2022-11-15 18:51:43 +01:00
|
|
|
|
2023-03-19 20:33:28 +01:00
|
|
|
for (ct, idx) in pairs(constraint_map)
|
|
|
|
JuMP.set_start_value.(idx, warmstart.slack[ct])
|
|
|
|
JuMP.set_dual_start_value.(idx, warmstart.dual[ct])
|
|
|
|
end
|
|
|
|
catch e
|
|
|
|
@warn "Warmstarting was not succesfull!" e
|
2022-11-15 18:51:43 +01:00
|
|
|
end
|
|
|
|
return model
|
|
|
|
end
|
|
|
|
|
|
|
|
function getwarmstart(model::JuMP.Model)
|
|
|
|
constraint_map = Dict(
|
|
|
|
ct => JuMP.all_constraints(model, ct...) for
|
|
|
|
ct in JuMP.list_of_constraint_types(model)
|
|
|
|
)
|
|
|
|
|
|
|
|
primal = value.(JuMP.all_variables(model))
|
|
|
|
|
|
|
|
slack = Dict(k => value.(v) for (k, v) in constraint_map)
|
|
|
|
duals = Dict(k => JuMP.dual.(v) for (k, v) in constraint_map)
|
|
|
|
|
2023-03-19 20:33:28 +01:00
|
|
|
return (primal = primal, dual = duals, slack = slack)
|
2022-11-15 18:51:43 +01:00
|
|
|
end
|
|
|
|
|
2023-03-19 20:33:28 +01:00
|
|
|
function solve(m::JuMP.Model, optimizer, warmstart = nothing)
|
2022-11-15 18:51:43 +01:00
|
|
|
JuMP.set_optimizer(m, optimizer)
|
|
|
|
MOIU.attach_optimizer(m)
|
|
|
|
|
|
|
|
m = setwarmstart!(m, warmstart)
|
|
|
|
|
|
|
|
JuMP.optimize!(m)
|
|
|
|
|
|
|
|
status = JuMP.termination_status(m)
|
|
|
|
|
|
|
|
return status, getwarmstart(m)
|
|
|
|
end
|
|
|
|
|
2023-03-19 20:33:28 +01:00
|
|
|
function solve(solverlog::String, m::JuMP.Model, optimizer, warmstart = nothing)
|
2022-11-15 18:51:43 +01:00
|
|
|
isdir(dirname(solverlog)) || mkpath(dirname(solverlog))
|
|
|
|
|
|
|
|
Base.flush(Base.stdout)
|
|
|
|
Base.Libc.flush_cstdio()
|
|
|
|
status, warmstart = open(solverlog, "a+") do logfile
|
|
|
|
redirect_stdout(logfile) do
|
|
|
|
status, warmstart = solve(m, optimizer, warmstart)
|
2022-11-17 11:22:17 +01:00
|
|
|
Base.Libc.flush_cstdio()
|
2023-03-19 20:33:28 +01:00
|
|
|
return status, warmstart
|
2022-11-15 18:51:43 +01:00
|
|
|
end
|
|
|
|
end
|
|
|
|
|
|
|
|
return status, warmstart
|
|
|
|
end
|