]> https://gitweb.dealii.org/ - dealii-svn.git/commitdiff
speed up FGMRES by initializing the aux vector only once. The values get overwritten...
authorheister <heister@0785d39b-7218-0410-832d-ea1e28bc413d>
Fri, 15 Oct 2010 08:47:45 +0000 (08:47 +0000)
committerheister <heister@0785d39b-7218-0410-832d-ea1e28bc413d>
Fri, 15 Oct 2010 08:47:45 +0000 (08:47 +0000)
git-svn-id: https://svn.dealii.org/trunk@22348 0785d39b-7218-0410-832d-ea1e28bc413d

deal.II/lac/include/lac/solver_gmres.h

index a1c76a1b7b73e6cb3a69cdf205494296f94e9ac1..e43f84667467cdab11e2df8791b15139e735a937 100644 (file)
@@ -892,9 +892,9 @@ SolverFGMRES<VECTOR>::solve (
   // Iteration starts here
 
   VECTOR* aux = this->memory.alloc();
+  aux->reinit(x);
   do
     {
-      aux->reinit(x);
       A.vmult(*aux, x);
       aux->sadd(-1., 1., b);
 

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.