From 1fbc2b5ea7b9134991637c2d0038b267f6771a26 Mon Sep 17 00:00:00 2001 From: wolf Date: Fri, 19 Nov 1999 16:20:28 +0000 Subject: [PATCH] . git-svn-id: https://svn.dealii.org/trunk@1907 0785d39b-7218-0410-832d-ea1e28bc413d --- deal.II/deal.II/source/multigrid/multigrid.cc | 11 ++++++----- deal.II/deal.II/source/numerics/multigrid.cc | 11 ++++++----- 2 files changed, 12 insertions(+), 10 deletions(-) diff --git a/deal.II/deal.II/source/multigrid/multigrid.cc b/deal.II/deal.II/source/multigrid/multigrid.cc index a593df249e..e5d997ccf6 100644 --- a/deal.II/deal.II/source/multigrid/multigrid.cc +++ b/deal.II/deal.II/source/multigrid/multigrid.cc @@ -41,10 +41,10 @@ MG::MG(const MGDoFHandler &mg_dof_handler, template void -MG::level_vmult(const unsigned int level, - Vector &result, - const Vector &u, - const Vector &/* rhs */) +MG::level_vmult (const unsigned int level, + Vector &result, + const Vector &u, + const Vector &/* rhs */) { (*level_matrices)[level].vmult(result,u); result.scale(-1.); @@ -101,7 +101,7 @@ void MGTransferPrebuilt::build_matrices (const MGDoFHandler &mg_dof) // cell prolongation_sparsities.back().reinit (mg_dof.n_dofs(level+1), mg_dof.n_dofs(level), -//TODO: evil hack, must be corrected!!! +//TODO: evil hack, must be corrected! dofs_per_cell+1); for (typename MGDoFHandler::cell_iterator cell=mg_dof.begin(level); @@ -194,6 +194,7 @@ void MGTransferPrebuilt::restrict_and_add (const unsigned int from_level, { Assert ((from_level >= 1) && (from_level<=prolongation_matrices.size()), ExcIndexRange (from_level, 1, prolongation_matrices.size()+1)); + prolongation_matrices[from_level-1].Tvmult_add (dst, src); }; diff --git a/deal.II/deal.II/source/numerics/multigrid.cc b/deal.II/deal.II/source/numerics/multigrid.cc index a593df249e..e5d997ccf6 100644 --- a/deal.II/deal.II/source/numerics/multigrid.cc +++ b/deal.II/deal.II/source/numerics/multigrid.cc @@ -41,10 +41,10 @@ MG::MG(const MGDoFHandler &mg_dof_handler, template void -MG::level_vmult(const unsigned int level, - Vector &result, - const Vector &u, - const Vector &/* rhs */) +MG::level_vmult (const unsigned int level, + Vector &result, + const Vector &u, + const Vector &/* rhs */) { (*level_matrices)[level].vmult(result,u); result.scale(-1.); @@ -101,7 +101,7 @@ void MGTransferPrebuilt::build_matrices (const MGDoFHandler &mg_dof) // cell prolongation_sparsities.back().reinit (mg_dof.n_dofs(level+1), mg_dof.n_dofs(level), -//TODO: evil hack, must be corrected!!! +//TODO: evil hack, must be corrected! dofs_per_cell+1); for (typename MGDoFHandler::cell_iterator cell=mg_dof.begin(level); @@ -194,6 +194,7 @@ void MGTransferPrebuilt::restrict_and_add (const unsigned int from_level, { Assert ((from_level >= 1) && (from_level<=prolongation_matrices.size()), ExcIndexRange (from_level, 1, prolongation_matrices.size()+1)); + prolongation_matrices[from_level-1].Tvmult_add (dst, src); }; -- 2.39.5