BUG: Gracefully fail for OOM (Out of Memory) in Presolve

src/Highs.h

@@ -23,6 +23,7 @@
 #include "model/HighsModel.h"
 #include "presolve/ICrash.h"
 #include "presolve/PresolveComponent.h"
+#include "util/HighsExceptions.h"
 
 /**
  * @brief Return the version

src/presolve/HPresolve.cpp

@@ -3685,10 +3685,18 @@ HPresolve::Result HPresolve::rowPresolve(HighsPostsolveStack& postsolve_stack,
               //	      assert(non_fractional);
               if (!non_fractional) return Result::kPrimalInfeasible;
             }
-            postsolve_stack.fixedColAtLower(nonzero.index(),
-                                            model->col_lower_[nonzero.index()],
-                                            model->col_cost_[nonzero.index()],
-                                            getColumnVector(nonzero.index()));
+            try {
+              postsolve_stack.fixedColAtLower(
+                  nonzero.index(), model->col_lower_[nonzero.index()],
+                  model->col_cost_[nonzero.index()],
+                  getColumnVector(nonzero.index()));
+            } catch (const DataStackOverflow& e) {
+              highsLogUser(options->log_options, HighsLogType::kError,
+                           "Problem too large for Presolve, try without\n");
+              highsLogUser(options->log_options, HighsLogType::kInfo,
+                           "Specific error raised:\n%s\n", e.what());
+              return Result::kPrimalInfeasible;
+            }
             if (model->col_upper_[nonzero.index()] >
                 model->col_lower_[nonzero.index()])
               changeColUpper(nonzero.index(),

src/util/HighsDataStack.h

@@ -16,9 +16,11 @@
 #define UTIL_HIGHS_DATA_STACK_H_
 
 #include <cstring>
+#include <string>
 #include <type_traits>
 #include <vector>
 
+#include "util/HighsExceptions.h"
 #include "util/HighsInt.h"
 
 #if __GNUG__ && __GNUC__ < 5 && !defined(__clang__)
@@ -37,8 +39,17 @@ class HighsDataStack {
   template <typename T,
             typename std::enable_if<IS_TRIVIALLY_COPYABLE(T), int>::type = 0>
   void push(const T& r) {
-    std::size_t dataSize = data.size();
-    data.resize(dataSize + sizeof(T));
+    HighsInt dataSize = data.size();
+    HighsInt newSize = dataSize + sizeof(T);
+    try {
+      data.resize(newSize);
+    } catch (const std::length_error& e) {
+      throw DataStackOverflow(
+          "Failed to resize the vector. Requested new size: " +
+          std::to_string(newSize) + ". Size to add is " +
+          std::to_string(sizeof(T)) +
+          ". Current size: " + std::to_string(data.size()) + ".");
+    }
     std::memcpy(data.data() + dataSize, &r, sizeof(T));
   }
 

src/util/HighsExceptions.h

@@ -0,0 +1,8 @@
+#pragma once
+#include <stdexcept>
+
+class DataStackOverflow : public std::runtime_error {
+ public:
+  explicit DataStackOverflow(const std::string& msg)
+      : std::runtime_error(msg) {}
+};