Add missing milp changes

Author: Opt-Mucca

highs/mip/HighsMipWorker.cpp

@@ -237,6 +237,11 @@ double HighsMipWorker::transformNewIntegerFeasibleSolution(
     }
   }
 
+  const double transformed_solobj =
+    static_cast<double>(static_cast<HighsInt>(mipsolver_.orig_model_->sense_) *
+                            mipsolver_quad_objective_value -
+                        mipsolver_.model_->offset_);
+
   // todo:
   // Possible MIP solution callback
   // if (!mipsolver.submip && feasible && mipsolver.callback_->user_callback &&
@@ -249,6 +254,12 @@ double HighsMipWorker::transformNewIntegerFeasibleSolution(
   //   assert(!interrupt);
   // }
 
+  // Catch the case where the repaired solution now has worse objective
+  // than the current stored solution
+  if (transformed_solobj >= upper_bound && !sol.empty()) {
+    return transformed_solobj;
+  }
+
   if (possibly_store_as_new_incumbent) {
     // Store the solution as incumbent in the original space if there
     // is no solution or if it is feasible
@@ -294,8 +305,5 @@ double HighsMipWorker::transformNewIntegerFeasibleSolution(
     }
   }
   // return the objective value in the transformed space
-  if (mipsolver_.orig_model_->sense_ == ObjSense::kMaximize)
-    return -double(mipsolver_quad_objective_value + mipsolver_.model_->offset_);
-
-  return double(mipsolver_quad_objective_value - mipsolver_.model_->offset_);
+  return transformed_solobj;
 }

highs/mip/HighsSearch.cpp

@@ -678,7 +678,7 @@ HighsInt HighsSearch::selectBranchingCandidate(int64_t maxSbIters,
         if (orbitalFixing)
           nodestack.back().stabilizerOrbits->orbitalFixing(localdom);
         else
-          getSymmetries();
+          getSymmetries().propagateOrbitopes(localdom);
       }
 
       inferences += localdom.getDomainChangeStack().size();