From d33008705d84342912e6ee0282f2622aa1b2828f Mon Sep 17 00:00:00 2001 From: Baerbel Jannsen Date: Wed, 10 Feb 2010 09:03:22 +0000 Subject: [PATCH] do not use variable steps for smoothing in the multigrid cycle, but rather increase the nuber of smoothing steps git-svn-id: https://svn.dealii.org/trunk@20540 0785d39b-7218-0410-832d-ea1e28bc413d --- deal.II/examples/step-16/step-16.cc | 6 ++---- 1 file changed, 2 insertions(+), 4 deletions(-) diff --git a/deal.II/examples/step-16/step-16.cc b/deal.II/examples/step-16/step-16.cc index b458ba22e8..1326e013e9 100644 --- a/deal.II/examples/step-16/step-16.cc +++ b/deal.II/examples/step-16/step-16.cc @@ -365,7 +365,7 @@ void LaplaceProblem::assemble_multigrid () ZeroFunction homogeneous_dirichlet_bc (1); dirichlet_boundary[0] = &homogeneous_dirichlet_bc; - std::vector > boundary_indices(triangulation.n_levels()); + std::vector > boundary_indices (triangulation.n_levels()); MGTools::make_boundary_list (mg_dof_handler, dirichlet_boundary, boundary_indices); @@ -485,7 +485,7 @@ void LaplaceProblem::solve () mg_smoother.initialize(mg_matrices, smoother_data); // Do two smoothing steps per level - mg_smoother.set_steps(1); + mg_smoother.set_steps(2); // Since the SOR method is not // symmetric, but we use conjugate // gradient iteration below, here @@ -494,8 +494,6 @@ void LaplaceProblem::solve () // symmetric operator even for // nonsymmetric smoothers. mg_smoother.set_symmetric(true); - mg_smoother.set_variable(true); - // We must wrap our matrices in an // object having the required -- 2.39.5