From 3d52c30485cf1e136ca1afe8b844a46aa7e0cdeb Mon Sep 17 00:00:00 2001
From: John Forrest <jjhforrest@gmail.com>
Date: Tue, 9 Apr 2024 13:45:51 +0100
Subject: [PATCH] correct ray when primal says infeasible

---
 src/ClpSimplexPrimal.cpp | 8 +++++++-
 1 file changed, 7 insertions(+), 1 deletion(-)

diff --git a/src/ClpSimplexPrimal.cpp b/src/ClpSimplexPrimal.cpp
index f8bc53b4..65b33d13 100644
--- a/src/ClpSimplexPrimal.cpp
+++ b/src/ClpSimplexPrimal.cpp
@@ -1484,9 +1484,15 @@ void ClpSimplexPrimal::statusOfProblemInPrimal(int &lastCleaned, int type,
               double *saveObjective = CoinCopyOfArray(objective(),
                 numberColumns_);
               memset(objective(), 0, numberColumns_ * sizeof(double));
-              infeasibilityCost_ = 1.0;
+              infeasibilityCost_ = 1.0e200;
               createRim(4);
               nonLinearCost_->checkInfeasibilities(primalTolerance_);
+	      for (int i=0;i<numberColumns_+numberRows_;i++) {
+		if (fabs(cost_[i])!=1.0e200)
+		  cost_[i] = 0.0;
+		else
+		  cost_[i] /= 1.0e200;
+	      }
               gutsOfSolution(NULL, NULL, false);
               memcpy(objective(), saveObjective, numberColumns_ * sizeof(double));
               delete[] saveObjective;