Add warm start

.gitignore

@@ -0,0 +1,2 @@
+deps/lib/
+deps/deps.jl

src/wrappers.jl

@@ -224,7 +224,24 @@ mutable struct Results
     Results() = new(Float64[], Float64[], QPALM.Info(), Float64[], Float64[])
 end
 
-function solve!(model::QPALM.Model, results::QPALM.Results=Results())
+function warm_start!(model::QPALM.Model;
+                     x_warm_start::Maybe{Vector{Float64}} = nothing,
+                     y_warm_start::Maybe{Vector{Float64}} = nothing)
+    if model.workspace == C_NULL
+        error("QPALM Model not setup")
+    end
+
+    ccall(
+        (:qpalm_warm_start, LIBQPALM_PATH),
+        Cvoid,
+        (Ptr{QPALM.Workspace}, Ptr{Float64}, Ptr{Float64}),
+        model.workspace,
+        x_warm_start != nothing ? pointer(x_warm_start) : C_NULL,
+        y_warm_start != nothing ? pointer(y_warm_start) : C_NULL,
+    )
+end
+
+function solve!(model::QPALM.Model, results::QPALM.Results=Results())::QPALM.Results
     ccall(
         (:qpalm_solve, LIBQPALM_PATH),
         Cvoid,
@@ -272,11 +289,12 @@ function solve!(model::QPALM.Model, results::QPALM.Results=Results())
         end
     end
 
-    if results.info.status == :Non_convex
-        results.info.obj_val = NaN
-    end
+    # An objective value is not calculated at this moment.
+    # if results.info.status == :Non_convex
+    #     results.info.obj_val = NaN
+    # end
 
-    results
+    return results
 end
 
 function cleanup!(model::QPALM.Model)

test/runtests.jl

@@ -1,4 +1,5 @@
 include("feasible.jl")
 include("infeasible.jl")
+include("warm_start.jl")
 include("linear_programs.jl")
 include("update.jl")

test/warm_start.jl

@@ -0,0 +1,49 @@
+using Test
+
+@testset "Warm start" begin
+
+    using QPALM
+    using LinearAlgebra
+    using SparseArrays
+
+    n, m = 2, 3
+    act = 3
+
+    Q = sparse(I, n, n)
+    A = sparse(ones(m, n))
+    q = ones(n)
+    bmin = -ones(m)
+    bmax = ones(m)
+
+
+    x0 = ones(n)
+    y0 = ones(m)
+
+    @testset "Not setup" begin
+        model = QPALM.Model()
+
+        @test_throws ErrorException QPALM.warm_start!(model, x_warm_start=x0)
+    end
+
+    @testset "Normal usage (1)" begin
+        model = QPALM.Model()
+
+        QPALM.setup!(model, Q=Q, q=q, A=A, bmin=bmin, bmax=bmax)
+        QPALM.warm_start!(model, x_warm_start=x0, y_warm_start=y0)
+
+        results = QPALM.solve!(model)
+        @test results.info.status == :Solved
+    end
+
+    @testset "Normal usage (2)" begin
+        model = QPALM.Model()
+
+        QPALM.setup!(model, Q=Q, q=q, A=A, bmin=bmin, bmax=bmax)
+        QPALM.warm_start!(model, y_warm_start=y0)
+        QPALM.warm_start!(model, x_warm_start=x0)
+
+        results = QPALM.solve!(model)
+        @test results.info.status == :Solved
+    end
+
+end