lpqp solvers don't always get starting points

jump-dev · Aug 26, 2016 · 90354e1 · 90354e1
1 parent cb330ce
commit 90354e1
Showing 1 changed file with 7 additions and 0 deletions.
diff --git a/src/nonlinear_algorithm.jl b/src/nonlinear_algorithm.jl
@@ -450,6 +450,12 @@ function MathProgBase.optimize!(m::PajaritoNonlinearModel)
 
     m.nlp_load_timer = 0.0
 
+    # if we haven't gotten a starting point,
+    # (e.g., if acting as MIQP solver) assume zero is okay)
+    if any(isnan,m.solution)
+        m.solution = zeros(length(m.solution))
+    end
+
     populatelinearmatrix(m)
     # solve it
     # MathProgBase.optimize!(nlp_model)
@@ -482,6 +488,7 @@ function MathProgBase.optimize!(m::PajaritoNonlinearModel)
 
     # pass in starting point
     #MathProgBase.setwarmstart!(nlp_model, m.solution)
+    @show m.solution[1:m.numVar]
     MathProgBase.setwarmstart!(ini_nlp_model, m.solution[1:m.numVar])
 
     start_nlp = time()