]> https://gitweb.dealii.org/ - dealii-svn.git/commitdiff
Simplify smoother usage a little bit. Avoid setting elements in vvectors that will...
authorkronbichler <kronbichler@0785d39b-7218-0410-832d-ea1e28bc413d>
Thu, 4 Feb 2010 11:22:51 +0000 (11:22 +0000)
committerkronbichler <kronbichler@0785d39b-7218-0410-832d-ea1e28bc413d>
Thu, 4 Feb 2010 11:22:51 +0000 (11:22 +0000)
git-svn-id: https://svn.dealii.org/trunk@20492 0785d39b-7218-0410-832d-ea1e28bc413d

deal.II/deal.II/include/multigrid/mg_smoother.h

index 7f4ff69a48f753f5ce02c11a8a274e76a038f99d..13fa78ecc6c0cc340a2d3496f9fc3279761678ad 100644 (file)
@@ -954,8 +954,8 @@ MGSmootherPrecondition<MATRIX, RELAX, VECTOR>::smooth(
 
   try
     {
-      r->reinit(u);
-      d->reinit(u);
+      r->reinit(u,true);
+      d->reinit(u,true);
 
       bool T = this->transpose;
       if (this->symmetric && (steps2 % 2 == 0))
@@ -980,7 +980,7 @@ MGSmootherPrecondition<MATRIX, RELAX, VECTOR>::smooth(
            if (this->debug > 0)
              deallog << 'N';
            matrices[level].vmult(*r,u);
-           r->sadd(-1.,1.,rhs);
+           r->sadd(-1.,rhs);
            if (this->debug > 2)
              deallog << ' ' << r->l2_norm() << ' ';
            smoothers[level].vmult(*d, *r);

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.