From 2ebf1c819cce25f787de63ddc9ff7474d6ea21f6 Mon Sep 17 00:00:00 2001
From: John Forrest <jjhforrest@gmail.com>
Date: Fri, 22 Mar 2024 13:46:34 +0000
Subject: [PATCH] resolve in ClpSolve if primal infeasible

---
 Clp/src/ClpSolve.cpp | 4 +++-
 1 file changed, 3 insertions(+), 1 deletion(-)

diff --git a/Clp/src/ClpSolve.cpp b/Clp/src/ClpSolve.cpp
index 824ab9a7..0e693be2 100644
--- a/Clp/src/ClpSolve.cpp
+++ b/Clp/src/ClpSolve.cpp
@@ -3818,7 +3818,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