Merge pull request #2672 from ERGO-Code/hipo-fix

Hipo fix

Author: jajhall

highs/ipm/hipo/factorhighs/FormatHandler.h

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

highs/ipm/hipo/factorhighs/HybridHybridFormatHandler.cpp

@@ -25,7 +25,9 @@ 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;
   frontal_.resize(frontal_size);
   std::memset(frontal_.data(), 0, frontal_size * sizeof(double));
 

highs/ipm/hipo/factorhighs/HybridSolveHandler.cpp

@@ -146,7 +146,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_frontal;
 
     // go through blocks of columns for this supernode in reverse order
     for (Int j = n_blocks - 1; j >= 0; --j) {