resolve in ClpSolve if primal infeasible

jhmgoossens · Mar 22, 2024 · 555b17b · 555b17b
1 parent 576b37f
commit 555b17b
Showing 1 changed file with 3 additions and 1 deletion.
diff --git a/src/ClpSolve.cpp b/src/ClpSolve.cpp
@@ -3850,7 +3850,9 @@ int ClpSimplex::initialSolve(ClpSolve &options)
       int savePerturbation = perturbation();
       if (savePerturbation == 50)
         setPerturbation(51); // small
-      if (!finalStatus || finalStatus == 2 || (moreSpecialOptions_ & 2) == 0 || fabs(sumDual) + fabs(sumPrimal) < 1.0e-3) {
+      if ((finalStatus>=0 && finalStatus <= 2) ||
+	  (moreSpecialOptions_ & 2) == 0 ||
+	  fabs(sumDual) + fabs(sumPrimal) < 1.0e-3) {
         if (finalStatus == 2) {
           if (sumDual > 1.0e-4) {
             // unbounded - get feasible first