From: wolf Date: Fri, 19 Nov 1999 14:53:12 +0000 (+0000) Subject: Correct fix of before and document why. X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=8a3a5b8cbd1f0f36bb7662795177dd6924064ba5;p=dealii-svn.git Correct fix of before and document why. git-svn-id: https://svn.dealii.org/trunk@1902 0785d39b-7218-0410-832d-ea1e28bc413d --- diff --git a/deal.II/deal.II/source/multigrid/mg_base.cc b/deal.II/deal.II/source/multigrid/mg_base.cc index 8ea78b4a28..25fa926268 100644 --- a/deal.II/deal.II/source/multigrid/mg_base.cc +++ b/deal.II/deal.II/source/multigrid/mg_base.cc @@ -106,6 +106,13 @@ MGBase::level_mgstep(const unsigned int level, // do recursion level_mgstep(level-1, pre_smooth, post_smooth, coarse_grid_solver); + + // reset size of the auxiliary + // vector, since it has been + // resized in the recursive call to + // level_mgstep directly above + t.reinit(solution[level].size()); + // do coarse grid correction transfer->prolongate(level, t, solution[level-1]); solution[level] += t; diff --git a/deal.II/lac/source/mgbase.cc b/deal.II/lac/source/mgbase.cc index 8ea78b4a28..25fa926268 100644 --- a/deal.II/lac/source/mgbase.cc +++ b/deal.II/lac/source/mgbase.cc @@ -106,6 +106,13 @@ MGBase::level_mgstep(const unsigned int level, // do recursion level_mgstep(level-1, pre_smooth, post_smooth, coarse_grid_solver); + + // reset size of the auxiliary + // vector, since it has been + // resized in the recursive call to + // level_mgstep directly above + t.reinit(solution[level].size()); + // do coarse grid correction transfer->prolongate(level, t, solution[level-1]); solution[level] += t;