]> https://gitweb.dealii.org/ - dealii.git/commitdiff
Write more.
authorWolfgang Bangerth <bangerth@math.tamu.edu>
Tue, 9 Feb 2010 23:45:58 +0000 (23:45 +0000)
committerWolfgang Bangerth <bangerth@math.tamu.edu>
Tue, 9 Feb 2010 23:45:58 +0000 (23:45 +0000)
git-svn-id: https://svn.dealii.org/trunk@20530 0785d39b-7218-0410-832d-ea1e28bc413d

deal.II/examples/step-16/doc/intro.dox

index f09995f383740dea951edd42c832a21a8ac652ae..5be8845f81c3e5a20bfdadc5957e07c8a0f763ca 100644 (file)
@@ -49,6 +49,42 @@ However, in comparison to this framework, the implementation in
 deal.II has to take into account the fact that we want to solve linear
 systems on adaptively refined meshes. This leads to the complication
 that it isn't quite as clear any more what exactly a "level" in a
-multilevel hierarchy of a mesh is.
+multilevel hierarchy of a mesh is. The following image shows what we
+consider to be a "level":
 
+<p align="center">
+  @image html "hanging_nodes.png" ""
+</p>
 
+In other words, the fine level in this mesh consists only of the
+degrees of freedom that are defined on the refined cells, but does not
+extend to that part of the domain that is not refined. While this
+guarantees that the overall effort grows as ${\cal O}(N)$ as necessary
+for optimal multigrid complexity, it leads to problems when defining
+where to smooth and what boundary conditions to pose for the operators
+defined on individual levels if the level boundary is not an external
+boundary. These questions are discussed in detail in the
+@ref mg_paper "Multigrid paper by Janssen and Kanschat" that describes
+the implementation in deal.II.
+
+
+
+<h3>The testcase</h3>
+
+The problem we solve here is exactly the same as in
+@ref step_6 "step-6", the only difference being the solver we use
+here. You may want to look there for a definition of what we solve,
+right hand side and boundary conditions. Obviously, the program would
+also work if we changed the geometry and other pieces of data that
+defines this particular problem.
+
+The things that are new are all those parts that concern the
+multigrid. In particular, this includes the following members of the
+main class:
+- <code>LaplaceProblem::mg_dof_handler</code>
+- <code>LaplaceProblem::mg_sparsity</code>
+- <code>LaplaceProblem::mg_matrices</code>
+- <code>LaplaceProblem::mg_interface_matrices_up</code>
+- <code>LaplaceProblem::assemble_multigrid ()</code>
+- <code>LaplaceProblem::solve ()</code>
+Take a look at these functions.

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.