From: Baerbel Jannsen Date: Wed, 10 Feb 2010 09:03:22 +0000 (+0000) Subject: do not use variable steps for smoothing in the multigrid cycle, but rather increase... X-Git-Tag: v8.0.0~6512 X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=d33008705d84342912e6ee0282f2622aa1b2828f;p=dealii.git 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 --- 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