Add back extra space

Author: filikat

highs/ipm/hipo/factorhighs/FormatHandler.h

@@ -88,6 +88,8 @@ class FormatHandler {
   virtual void extremeEntries() {}
 };
 
+const Int64 extra_space = 10;
+
 }  // namespace hipo
 
 #endif

highs/ipm/hipo/factorhighs/HybridHybridFormatHandler.cpp

@@ -21,10 +21,13 @@ HybridHybridFormatHandler::HybridHybridFormatHandler(
 void HybridHybridFormatHandler::initFrontal() {
   const Int n_blocks = (sn_size_ - 1) / nb_ + 1;
   diag_start_.resize(n_blocks);
-  Int64 frontal_size = getDiagStart(ldf_, sn_size_, nb_, n_blocks, diag_start_);
+  Int64 frontal_size =
+      getDiagStart(ldf_, sn_size_, nb_, n_blocks, diag_start_) + extra_space;
   frontal_.resize(frontal_size);
   std::memset(frontal_.data(), 0, frontal_size * sizeof(double));
 
+  // extra_space is not needed, but it avoid some weird issue on windows
+
   // frontal_ is actually allocated just the first time, then the memory is
   // reused from the previous factorisations and just initialised.
 }

highs/ipm/hipo/factorhighs/HybridSolveHandler.cpp

@@ -145,7 +145,7 @@ void HybridSolveHandler::backwardSolve(std::vector<double>& x) const {
 
     // index to access snColumns[sn]
     // initialised with the total number of entries of snColumns[sn]
-    Int64 SnCol_ind = sn_columns_[sn].size();
+    Int64 SnCol_ind = sn_columns_[sn].size() - extra_space;
 
     // go through blocks of columns for this supernode in reverse order
     for (Int j = n_blocks - 1; j >= 0; --j) {