From dd14ad789fcb56d0dc5fa9e6b5c9af4b6e650305 Mon Sep 17 00:00:00 2001 From: Martin Kronbichler Date: Thu, 22 Sep 2022 09:24:15 +0200 Subject: [PATCH] Do not default initialize projection matrix in GMRES --- include/deal.II/lac/solver_gmres.h | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) 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::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); -- 2.39.5