]> https://gitweb.dealii.org/ - dealii.git/commitdiff
do not use variable steps for smoothing in the multigrid cycle, but rather increase...
authorBaerbel Jannsen <baerbel.janssen@gmail.com>
Wed, 10 Feb 2010 09:03:22 +0000 (09:03 +0000)
committerBaerbel Jannsen <baerbel.janssen@gmail.com>
Wed, 10 Feb 2010 09:03:22 +0000 (09:03 +0000)
git-svn-id: https://svn.dealii.org/trunk@20540 0785d39b-7218-0410-832d-ea1e28bc413d

deal.II/examples/step-16/step-16.cc

index b458ba22e8e66dd42c0f5eba6a3b1ee93f0af946..1326e013e9cc61e4f24fc66e0bae12371a1894b7 100644 (file)
@@ -365,7 +365,7 @@ void LaplaceProblem<dim>::assemble_multigrid ()
   ZeroFunction<dim>                    homogeneous_dirichlet_bc (1);
   dirichlet_boundary[0] = &homogeneous_dirichlet_bc;
 
-  std::vector<std::set<unsigned int> > boundary_indices(triangulation.n_levels());
+  std::vector<std::set<unsigned int> > boundary_indices (triangulation.n_levels());
   MGTools::make_boundary_list (mg_dof_handler, dirichlet_boundary,
                               boundary_indices);
 
@@ -485,7 +485,7 @@ void LaplaceProblem<dim>::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<dim>::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

In the beginning the Universe was created. This has made a lot of people very angry and has been widely regarded as a bad move.

Douglas Adams


Typeset in Trocchi and Trocchi Bold Sans Serif.