From 9d690e02d67213843b253e00642b2d0d4b81acb8 Mon Sep 17 00:00:00 2001 From: Wolfgang Bangerth Date: Mon, 2 Dec 2019 14:55:53 -0700 Subject: [PATCH] Augment step-61 documentation. --- examples/step-61/step-61.cc | 27 +++++++++++++++------------ 1 file changed, 15 insertions(+), 12 deletions(-) diff --git a/examples/step-61/step-61.cc b/examples/step-61/step-61.cc index 5c77011b15..5b983318f6 100644 --- a/examples/step-61/step-61.cc +++ b/examples/step-61/step-61.cc @@ -19,7 +19,10 @@ // @sect3{Include files} // This program is based on step-7, step-20 and step-51, // so most of the following header files are familiar. We -// need the following: +// need the following, of which only the one that +// imports the FE_DGRaviartThomas class (namely, `deal.II/fe/fe_dg_vector.h`) +// is really new; the FE_DGRaviartThomas implements the "broken" Raviart-Thomas +// space discussed in the introduction: #include #include #include @@ -49,13 +52,9 @@ #include #include #include -#include #include #include #include -#include -#include -#include #include #include @@ -76,7 +75,14 @@ namespace Step61 // // The structure of the class is not fundamentally different from that of // previous tutorial programs, so there is little need to comment on the - // details. + // details with one exception: The class has a member variable `fe_dgrt` + // that corresponds to the "broken" Raviart-Thomas space mentioned in the + // introduction. There is a matching `dof_handler_dgrt` that represents a + // global enumeration of a finite element field created from this element, and + // a vector `darcy_velocity` that holds nodal values for this field. We will + // use these three variables after solving for the pressure to compute a + // postprocessed velocity field for which we can then evaluate the error + // and which we can output for visualization. template class WGDarcyEquation { @@ -98,9 +104,6 @@ namespace Step61 FESystem fe; DoFHandler dof_handler; - FE_DGRaviartThomas fe_dgrt; - DoFHandler dof_handler_dgrt; - AffineConstraints constraints; SparsityPattern sparsity_pattern; @@ -109,7 +112,9 @@ namespace Step61 Vector solution; Vector system_rhs; - Vector darcy_velocity; + FE_DGRaviartThomas fe_dgrt; + DoFHandler dof_handler_dgrt; + Vector darcy_velocity; }; @@ -263,7 +268,6 @@ namespace Step61 , dof_handler(triangulation) , fe_dgrt(degree) , dof_handler_dgrt(triangulation) - {} @@ -1026,7 +1030,6 @@ int main() { try { - dealii::deallog.depth_console(2); Step61::WGDarcyEquation<2> wg_darcy(0); wg_darcy.run(); } -- 2.39.5