minor optimisation

2021-04-30 23:18:20 +10:00
parent 90c4a0d1e1
commit 7f6234517c
1 changed files with 5 additions and 1 deletions
@@ -166,7 +166,11 @@ SymEigsSolver<eT, SelectionRule, OpType>::restart(uword k)
    nnz = ncv - k + i + 1;
    Mat<eT> V(fac_V.memptr(), dim_n, nnz, false);
    Col<eT> q(Q.colptr(i), nnz, false);
-    Vs.col(i) = V * q;
+    // OLD CODE:
+    // Vs.col(i) = V * q;
+    // NEW CODE:
+    Col<eT> v(Vs.colptr(i), dim_n, false, true);
+    v = V * q;
    }
  Vs.col(k) = fac_V * Q.col(k);
  fac_V.head_cols(k + 1) = Vs;