save gurobi's Irreducible Infeasible Sets (IIS) for easier debug of i…

…nfeasible optimization problems
joschu · May 23, 2017 · 8da5bb0 · 8da5bb0
1 parent bf823b1
commit 8da5bb0
Show file tree

Hide file tree

Showing 2 changed files with 6 additions and 2 deletions.
diff --git a/src/sco/gurobi_interface.cpp b/src/sco/gurobi_interface.cpp
@@ -188,7 +188,10 @@ CvxOptStatus GurobiModel::optimize(){
     LOG_DEBUG("solver objective value: %.3e", objval);
     return CVX_SOLVED;
   }
-  else if (status == GRB_INFEASIBLE) return CVX_INFEASIBLE;
+  else if (status == GRB_INFEASIBLE) {
+    GRBcomputeIIS(m_model);
+    return CVX_INFEASIBLE;
+  }
   else return CVX_FAILED;
 }
 CvxOptStatus GurobiModel::optimizeFeasRelax(){

diff --git a/src/sco/optimizers.cpp b/src/sco/optimizers.cpp
@@ -289,8 +289,9 @@ OptStatus BasicTrustRegionSQP::optimize() {
         CvxOptStatus status = model_->optimize();
         ++results_.n_qp_solves;
         if (status != CVX_SOLVED) {
-          LOG_ERROR("convex solver failed! set TRAJOPT_LOG_THRESH=DEBUG to see solver output. saving model to /tmp/fail.lp");
+          LOG_ERROR("convex solver failed! set TRAJOPT_LOG_THRESH=DEBUG to see solver output. saving model to /tmp/fail.lp and IIS to /tmp/fail.ilp");
           model_->writeToFile("/tmp/fail.lp");
+          model_->writeToFile("/tmp/fail.ilp");
           retval = OPT_FAILED;
           goto cleanup;
         }