]> https://gitweb.dealii.org/ - dealii.git/commitdiff
Do not default initialize projection matrix in GMRES 14300/head
authorMartin Kronbichler <martin.kronbichler@uni-a.de>
Thu, 22 Sep 2022 07:24:15 +0000 (09:24 +0200)
committerMartin Kronbichler <martin.kronbichler@uni-a.de>
Thu, 22 Sep 2022 07:24:15 +0000 (09:24 +0200)
include/deal.II/lac/solver_gmres.h

index b4166385a24ede1c3f7ddc6bc4d110ca4c6604a3..9b21890a241ecdcf09dbd436a68f2dacd88a7f60 100644 (file)
@@ -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);

In the beginning the Universe was created. This has made a lot of people very angry and has been widely regarded as a bad move.

Douglas Adams


Typeset in Trocchi and Trocchi Bold Sans Serif.