From: Martin Kronbichler 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::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);