From: Wolfgang Bangerth <bangerth@math.tamu.edu>
Date: Fri, 12 Feb 2010 13:16:36 +0000 (+0000)
Subject: Document a couple functions.
X-Git-Tag: v8.0.0~6490
X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=32305202575ce97cf9c11888bccd3613db11bb2d;p=dealii.git

Document a couple functions.


git-svn-id: https://svn.dealii.org/trunk@20566 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 9766fd58e8..bd9382e5ba 100644
--- a/deal.II/examples/step-16/step-16.cc
+++ b/deal.II/examples/step-16/step-16.cc
@@ -547,6 +547,24 @@ void LaplaceProblem<dim>::solve ()
 
 
 
+                                 // @sect4{Postprocessing}
+
+				 // The following two functions postprocess a
+				 // solution once it is computed. In
+				 // particular, the first one refines the mesh
+				 // at the beginning of each cycle while the
+				 // second one outputs results at the end of
+				 // each such cycle. The functions are almost
+				 // unchanged from those in step-6, with the
+				 // exception of two minor differences: The
+				 // KellyErrorEstimator::estimate function
+				 // wants an argument of type DoFHandler, not
+				 // MGDoFHandler, and so we have to cast from
+				 // derived to base class; and we generate
+				 // output in VTK format, to use the more
+				 // modern visualization programs available
+				 // today compared to those that were
+				 // available when step-6 was written.
 template <int dim>
 void LaplaceProblem<dim>::refine_grid ()
 {
@@ -564,7 +582,6 @@ void LaplaceProblem<dim>::refine_grid ()
 }
 
 
-
 template <int dim>
 void LaplaceProblem<dim>::output_results (const unsigned int cycle) const
 {
@@ -584,7 +601,16 @@ void LaplaceProblem<dim>::output_results (const unsigned int cycle) const
 }
 
 
+                                 // @sect4{LaplaceProblem::run}
 
+				 // Like several of the functions above, this
+				 // is almost exactly a copy of of the
+				 // corresponding function in step-6. The only
+				 // difference is the call to
+				 // <code>assemble_multigrid</code> that takes
+				 // care of forming the matrices on every
+				 // level that we need in the multigrid
+				 // method.
 template <int dim>
 void LaplaceProblem<dim>::run ()
 {