Add HighsCDouble to obj calc

Author: Opt-Mucca

highs/mip/HighsDebugSol.cpp

@@ -52,8 +52,7 @@ void HighsDebugSol::activate() {
         if (line.find("# Rows") != std::string::npos) {
           break;
         }
-        if (!incolsection)
-          continue;
+        if (!incolsection) continue;
 
         std::istringstream linestream(line);
         linestream >> varname;
@@ -70,8 +69,8 @@ void HighsDebugSol::activate() {
 
       HighsCDouble debugsolobj = 0.0;
       for (HighsInt i = 0; i != mipsolver->orig_model_->num_col_; ++i)
-        debugsolobj +=
-            mipsolver->orig_model_->col_cost_[i] * debugOrigSolution[i];
+        debugsolobj += mipsolver->orig_model_->col_cost_[i] *
+                       HighsCDouble(debugOrigSolution[i]);
 
       debugSolObjective = double(debugsolobj + mipsolver->orig_model_->offset_);
       debugSolActive = true;

highs/mip/HighsMipSolverData.cpp

@@ -694,7 +694,6 @@ void HighsMipSolverData::init() {
 }
 
 void HighsMipSolverData::runPresolve(const HighsInt presolve_reduction_limit) {
-
   mipsolver.timer_.start(mipsolver.timer_.presolve_clock);
   presolve::HPresolve presolve;
   if (!presolve.okSetInput(mipsolver, presolve_reduction_limit)) {
@@ -1018,7 +1017,8 @@ void HighsMipSolverData::runSetup() {
     debugSolution.debugSolObjective = 0;
     HighsCDouble debugsolobj = 0.0;
     for (HighsInt i = 0; i != mipsolver.model_->num_col_; ++i)
-      debugsolobj += mipsolver.colCost(i) * debugSolution.debugSolution[i];
+      debugsolobj +=
+          mipsolver.colCost(i) * HighsCDouble(debugSolution.debugSolution[i]);
     debugSolution.debugSolObjective = static_cast<double>(debugsolobj);
     debugSolution.registerDomain(domain);
     assert(checkSolution(debugSolution.debugSolution));