Do not default initialize projection matrix in GMRES

diff --git a/include/deal.II/lac/solver_gmres.h b/include/deal.II/lac/solver_gmres.h
index b4166385a24ede1c3f7ddc6bc4d110ca4c6604a3..9b21890a241ecdcf09dbd436a68f2dacd88a7f60 100644 (file)
--- a/include/deal.II/lac/solver_gmres.h
+++ b/include/deal.II/lac/solver_gmres.h
@@ -879,7 +879,7 @@ SolverGMRES<VectorType>::solve(const MatrixType &        A,
     H_orig.reinit(n_tmp_vectors, n_tmp_vectors - 1);
 
   // matrix used for the orthogonalization process later
-  H.reinit(n_tmp_vectors, n_tmp_vectors - 1);
+  H.reinit(n_tmp_vectors, n_tmp_vectors - 1, /* omit_initialization */ true);
 
   // some additional vectors, also used in the orthogonalization
   gamma.reinit(n_tmp_vectors);