Removed debugZeroBasicDuals and reverted call to lpsolver.run() in HighsLpRelaxation.cpp

Author: jajhall

highs/mip/HighsLpRelaxation.cpp

@@ -1151,7 +1151,9 @@ HighsLpRelaxation::Status HighsLpRelaxation::run(bool resolve_on_error) {
       use_simplex = true;
     }
   }
-  if (use_simplex) callstatus = lpsolver.run();
+  if (use_simplex) {
+    callstatus = lpsolver.run();
+  }
   // Revert the value of lpsolver.options_.solver
   lpsolver.setOptionValue("solver", solver);
   if (!mipsolver.submip) {

highs/simplex/HEkk.cpp

@@ -3632,8 +3632,6 @@ HighsStatus HEkk::returnFromSolve(const HighsStatus return_status) {
   } else {
     return_dual_solution_status_ = kSolutionStatusInfeasible;
   }
-  // Basic duals are nonzero
-  //  assert(debugZeroBasicDuals());
   assert(debugNoShiftsOrPerturbations());
   computePrimalObjectiveValue();
   if (!options_->log_dev_level) {

highs/simplex/HEkk.h

@@ -409,7 +409,6 @@ class HEkk {
   HighsDebugStatus debugSimplexDualInfeasible(const std::string message,
                                               const bool force_report = false);
   HighsDebugStatus debugComputeDual(const bool initialise = false) const;
-  bool debugZeroBasicDuals() const;
   bool debugNoShiftsOrPerturbations() const;
   friend class HEkkPrimal;
   friend class HEkkDual;

highs/simplex/HEkkDebug.cpp

@@ -1694,16 +1694,6 @@ HighsDebugStatus HEkk::debugSimplexDualInfeasible(const std::string message,
   return HighsDebugStatus::kOk;
 }
 
-bool HEkk::debugZeroBasicDuals() const {
-  const std::vector<double> dual = this->info_.workDual_;
-  const std::vector<HighsInt> basis = this->basis_.basicIndex_;
-  for (HighsInt iRow = 0; iRow < this->lp_.num_row_; iRow++)
-    if (dual[basis[iRow]]) {
-      return false;
-    }
-  return true;
-}
-
 bool HEkk::debugNoShiftsOrPerturbations() const {
   switch (model_status_) {
     case HighsModelStatus::kOptimal: {