From: Martin Kronbichler <martin.kronbichler@uni-a.de>
Date: Thu, 22 Sep 2022 07:24:15 +0000 (+0200)
Subject: Do not default initialize projection matrix in GMRES
X-Git-Tag: v9.5.0-rc1~933^2
X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=dd14ad789fcb56d0dc5fa9e6b5c9af4b6e650305;p=dealii.git

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 b4166385a2..9b21890a24 100644
--- 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);