Use tolerance when rounding

Author: BenChampion

highs/lp_data/HighsInterface.cpp

@@ -3478,8 +3478,9 @@ bool Highs::infeasibleBoundsOk() {
       if (perform_inward_integer_rounding &&
           lp.integrality_[iCol] == HighsVarType::kInteger) {
         // Assess bounds after inward integer rounding
-        double integer_lower = std::ceil(lower);
-        double integer_upper = std::floor(upper);
+        double feastol = this->options_.mip_feasibility_tolerance;
+        double integer_lower = std::ceil(lower - feastol);
+        double integer_upper = std::floor(upper + feastol);
         assert(integer_lower >= lower);
         assert(integer_upper <= upper);
         performed_inward_integer_rounding =

highs/mip/HighsFeasibilityJump.cpp

@@ -46,8 +46,8 @@ HighsModelStatus HighsMipSolverData::feasibilityJump() {
       fjVarType = external_feasibilityjump::VarType::Continuous;
     } else {
       fjVarType = external_feasibilityjump::VarType::Integer;
-      lower = std::ceil(lower);
-      upper = std::floor(upper);
+      lower = std::ceil(lower - feastol);
+      upper = std::floor(upper + feastol);
     }
 
     const bool legal_bounds = lower <= upper && lower < kHighsInf &&