From 8a3a5b8cbd1f0f36bb7662795177dd6924064ba5 Mon Sep 17 00:00:00 2001 From: wolf Date: Fri, 19 Nov 1999 14:53:12 +0000 Subject: [PATCH] Correct fix of before and document why. git-svn-id: https://svn.dealii.org/trunk@1902 0785d39b-7218-0410-832d-ea1e28bc413d --- deal.II/deal.II/source/multigrid/mg_base.cc | 7 +++++++ deal.II/lac/source/mgbase.cc | 7 +++++++ 2 files changed, 14 insertions(+) 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; -- 2.39.5