debug

Author: filikat

highs/ipm/hipo/factorhighs/HybridSolveHandler.cpp

@@ -80,10 +80,11 @@ void HybridSolveHandler::forwardSolve(std::vector<double>& x) const {
       // temporary space for gemv
       const Int gemv_space = ldSn - nb * j - jb;
       std::vector<double> y(gemv_space);
-
-      callAndTime_dgemv('T', jb, gemv_space, 1.0, &sn_columns_[sn][SnCol_ind],
-                        jb, &x[x_start], 1, 0.0, y.data(), 1, data_);
-      SnCol_ind += jb * gemv_space;
+      if (gemv_space > 0) {
+        callAndTime_dgemv('T', jb, gemv_space, 1.0, &sn_columns_[sn][SnCol_ind],
+                          jb, &x[x_start], 1, 0.0, y.data(), 1, data_);
+        SnCol_ind += jb * gemv_space;
+      }
 #if HIPO_TIMING_LEVEL >= 2
       if (data_) data_->sumTime(kTimeSolveSolve_dense, clock.stop());
 #endif
@@ -174,24 +175,27 @@ void HybridSolveHandler::backwardSolve(std::vector<double>& x) const {
       const Int gemv_space = ldSn - nb * j - jb;
       std::vector<double> y(gemv_space);
 
+      if (gemv_space > 0) {
 #if HIPO_TIMING_LEVEL >= 2
-      clock.start();
+        clock.start();
 #endif
-      // scatter entries into y
-      for (Int i = 0; i < gemv_space; ++i) {
-        const Int row = S_.rows(start_row + nb * j + jb + i);
-        y[i] = x[row];
-      }
+        // scatter entries into y
+        for (Int i = 0; i < gemv_space; ++i) {
+          const Int row = S_.rows(start_row + nb * j + jb + i);
+          y[i] = x[row];
+        }
 #if HIPO_TIMING_LEVEL >= 2
-      if (data_) data_->sumTime(kTimeSolveSolve_sparse, clock.stop());
+        if (data_) data_->sumTime(kTimeSolveSolve_sparse, clock.stop());
 #endif
 
 #if HIPO_TIMING_LEVEL >= 2
-      clock.start();
+        clock.start();
 #endif
-      SnCol_ind -= jb * gemv_space;
-      callAndTime_dgemv('N', jb, gemv_space, -1.0, &sn_columns_[sn][SnCol_ind],
-                        jb, y.data(), 1, 1.0, &x[x_start], 1, data_);
+        SnCol_ind -= jb * gemv_space;
+        callAndTime_dgemv('N', jb, gemv_space, -1.0,
+                          &sn_columns_[sn][SnCol_ind], jb, y.data(), 1, 1.0,
+                          &x[x_start], 1, data_);
+      }
 
       SnCol_ind -= diag_entries;
       callAndTime_dtrsv('U', 'N', 'U', jb, &sn_columns_[sn][SnCol_ind], jb,