Merge branch 'main' into cbrissette/cuopt-bindings

Author: ramakrishnap-nv

conda/environments/all_cuda-129_arch-aarch64.yaml

@@ -65,7 +65,7 @@ dependencies:
 - python>=3.11,<3.15
 - pyyaml>=6.0.0
 - rapids-build-backend>=0.4.0,<0.5.0
-- rapids-logger==0.2.*,>=0.0.0a0
+- rapids-logger==0.3.*
 - re2
 - requests
 - rmm==26.10.*,>=0.0.0a0

conda/environments/all_cuda-129_arch-x86_64.yaml

@@ -65,7 +65,7 @@ dependencies:
 - python>=3.11,<3.15
 - pyyaml>=6.0.0
 - rapids-build-backend>=0.4.0,<0.5.0
-- rapids-logger==0.2.*,>=0.0.0a0
+- rapids-logger==0.3.*
 - re2
 - requests
 - rmm==26.10.*,>=0.0.0a0

conda/environments/all_cuda-133_arch-aarch64.yaml

@@ -65,7 +65,7 @@ dependencies:
 - python>=3.11,<3.15
 - pyyaml>=6.0.0
 - rapids-build-backend>=0.4.0,<0.5.0
-- rapids-logger==0.2.*,>=0.0.0a0
+- rapids-logger==0.3.*
 - re2
 - requests
 - rmm==26.10.*,>=0.0.0a0

conda/environments/all_cuda-133_arch-x86_64.yaml

@@ -65,7 +65,7 @@ dependencies:
 - python>=3.11,<3.15
 - pyyaml>=6.0.0
 - rapids-build-backend>=0.4.0,<0.5.0
-- rapids-logger==0.2.*,>=0.0.0a0
+- rapids-logger==0.3.*
 - re2
 - requests
 - rmm==26.10.*,>=0.0.0a0

conda/recipes/libcuopt/recipe.yaml

@@ -82,7 +82,7 @@ cache:
       - cuda-version =${{ cuda_version }}
       - libraft-headers =${{ minor_version }}
       - librmm =${{ minor_version }}
-      - rapids-logger =0.2
+      - rapids-logger =0.3
       - cuda-nvtx-dev
       - libcudss-dev >=0.7,<0.8
       - libcurand-dev
@@ -124,7 +124,7 @@ outputs:
       host:
         - libboost-devel
         - cuda-version =${{ cuda_version }}
-        - rapids-logger =0.2
+        - rapids-logger =0.3
         - librmm =${{ minor_version }}
         - libcublas
         - libcudss-dev >=0.7,<0.8

cpp/CMakeLists.txt

@@ -295,7 +295,7 @@ FetchContent_MakeAvailable(papilo)
 FetchContent_Declare(
         pslp
         GIT_REPOSITORY "https://github.com/dance858/PSLP.git"
-        GIT_TAG "v0.0.8"
+        GIT_TAG "v0.0.10"
         GIT_PROGRESS TRUE
         EXCLUDE_FROM_ALL
         SYSTEM

cpp/src/math_optimization/solver_settings.cu

@@ -94,8 +94,8 @@ solver_settings_t<i_t, f_t>::solver_settings_t() : pdlp_settings(), mip_settings
     {CUOPT_RELATIVE_PRIMAL_TOLERANCE, &pdlp_settings.tolerances.relative_primal_tolerance, f_t(0.0), f_t(1e-1), f_t(1e-4)},
     {CUOPT_ABSOLUTE_GAP_TOLERANCE, &pdlp_settings.tolerances.absolute_gap_tolerance, f_t(0.0), f_t(1e-1), f_t(1e-4)},
     {CUOPT_RELATIVE_GAP_TOLERANCE, &pdlp_settings.tolerances.relative_gap_tolerance, f_t(0.0), f_t(1e-1), f_t(1e-4)},
-    {CUOPT_MIP_ABSOLUTE_TOLERANCE, &mip_settings.tolerances.absolute_tolerance, f_t(0.0), f_t(1e-1), f_t(1e-4)},
-    {CUOPT_MIP_RELATIVE_TOLERANCE, &mip_settings.tolerances.relative_tolerance, f_t(0.0), f_t(1e-1), f_t(1e-4)},
+    {CUOPT_MIP_ABSOLUTE_TOLERANCE, &mip_settings.tolerances.absolute_tolerance, f_t(0.0), f_t(1e-1), f_t(1e-6)},
+    {CUOPT_MIP_RELATIVE_TOLERANCE, &mip_settings.tolerances.relative_tolerance, f_t(0.0), f_t(1e-1), f_t(1e-12)},
     {CUOPT_MIP_INTEGRALITY_TOLERANCE, &mip_settings.tolerances.integrality_tolerance, f_t(0.0), f_t(1e-1), f_t(1e-5)},
     {CUOPT_MIP_ABSOLUTE_GAP, &mip_settings.tolerances.absolute_mip_gap, f_t(0.0), std::numeric_limits<f_t>::infinity(), std::max(f_t(1e-10), std::numeric_limits<f_t>::epsilon())},
     {CUOPT_MIP_RELATIVE_GAP, &mip_settings.tolerances.relative_mip_gap, f_t(0.0), f_t(1e-1), f_t(1e-4)},

cpp/src/mip_heuristics/feasibility_jump/feasibility_jump.cu

@@ -1036,9 +1036,20 @@ void fj_t<i_t, f_t>::resize_vectors(const raft::handle_t* handle_ptr)
   climbers[0]->grid_delta_buf.resize(update_weights_launch_dims.first.x, handle_ptr->get_stream());
 
   // FJ related vars
-  cstr_weights.resize(pb_ptr->n_constraints, handle_ptr->get_stream());
-  cstr_right_weights.resize(pb_ptr->n_constraints, handle_ptr->get_stream());
-  cstr_left_weights.resize(pb_ptr->n_constraints, handle_ptr->get_stream());
+  // the problem can gain constraints between two runs (e.g. the objective cutting plane added by
+  // the feasibility pump), and resize leaves the new elements uninitialized: give them the default
+  // weight
+  auto resize_weights = [&](rmm::device_uvector<f_t>& weights) {
+    const auto old_size = weights.size();
+    weights.resize(pb_ptr->n_constraints, handle_ptr->get_stream());
+    if (old_size < weights.size()) {
+      thrust::uninitialized_fill(
+        handle_ptr->get_thrust_policy(), weights.begin() + old_size, weights.end(), 1.);
+    }
+  };
+  resize_weights(cstr_weights);
+  resize_weights(cstr_right_weights);
+  resize_weights(cstr_left_weights);
   constraint_lower_bounds_csr.resize(pb_ptr->coefficients.size(), handle_ptr->get_stream());
   constraint_upper_bounds_csr.resize(pb_ptr->coefficients.size(), handle_ptr->get_stream());
   cstr_coeff_reciprocal.resize(pb_ptr->coefficients.size(), handle_ptr->get_stream());

cpp/src/mip_heuristics/presolve/third_party_presolve.cpp

@@ -724,7 +724,7 @@ void set_presolve_options(papilo::Presolve<f_t>& presolver,
 {
   presolver.getPresolveOptions().tlim    = time_limit;
   presolver.getPresolveOptions().threads = num_cpu_threads;  //  user setting or  0 (automatic)
-  presolver.getPresolveOptions().feastol = 1e-5;
+  presolver.getPresolveOptions().feastol = absolute_tolerance;
   if (max_rounds > 0) { presolver.getPresolveOptions().maxrounds = max_rounds; }
   if (dual_postsolve) {
     presolver.getPresolveOptions().componentsmaxint = -1;

dependencies.yaml

@@ -565,7 +565,7 @@ dependencies:
     common:
       - output_types: [conda, requirements, pyproject]
         packages:
-          - rapids-logger==0.2.*,>=0.0.0a0
+          - rapids-logger==0.3.*
       - output_types: requirements
         packages:
           # pip recognizes the index as a global option for the requirements.txt file