// `sqrt' and `fabs' functions:
#include <cmath>
+ // The final step in importing
+ // deal.II is this: All deal.II
+ // functions and classes are in a
+ // namespace <code>dealii</code>, to
+ // make sure they don't clash with
+ // symbols from other libraries you
+ // may want to use in conjunction
+ // with deal.II. One could use these
+ // functions and classes by prefixing
+ // every use of these names by
+ // <code>dealii::</code>, but that
+ // would quickly become cumbersome
+ // and annoying. Rather, we simply
+ // import the entire deal.II
+ // namespace for general use:
using namespace dealii;
// @sect3{Creating the first mesh}
#include <fstream>
#include <cmath>
+ // The last step is as in previous
+ // programs:
using namespace dealii;
// Now, as we want to compute the
#include <iomanip>
#include <cmath>
+ // The last step is as in all
+ // previous programs:
using namespace dealii;
// Then we declare a class which
#include <iostream>
#include <fstream>
+ // The last step is as in all
+ // previous programs:
using namespace dealii;
// @sect3{Equation data}
#include <list>
#include <sstream>
+ // The last step is as in all
+ // previous programs:
using namespace dealii;
// @sect3{Evaluation of the solution}
#include <numeric>
#include <sstream>
+ // The last step is as in all
+ // previous programs:
using namespace dealii;
// @sect3{Evaluating the solution}
#include <iostream>
#include <sstream>
+ // The last step is as in all
+ // previous programs:
using namespace dealii;
// The first thing we have here is a helper
#include <fstream>
#include <sstream>
+ // The last step is as in all
+ // previous programs:
using namespace dealii;
// This class is based on the same
#include <iostream>
#include <sstream>
+ // The last step is as in all
+ // previous programs:
using namespace dealii;
// Now, here comes the declaration of the
#include <sstream>
#include <iomanip>
+ // The last step is as in all
+ // previous programs:
using namespace dealii;
// So much for the header files. As a
#include <iostream>
#include <fstream>
+ // As mentioned in the first few
+ // tutorial programs, all names in
+ // deal.II are declared in a
+ // namespace <code>dealii</code>. To
+ // make using these function and
+ // class names simpler, we import the
+ // entire content of that namespace
+ // into the global scope:
using namespace dealii;
// Before we start with the actual program,
// And this is again needed for C++ output:
#include <fstream>
+ // Finally, as in step-1, we import
+ // the deal.II namespace into the
+ // global scope:
using namespace dealii;
// @sect3{Mesh generation}
// offers such functionality:
#include <base/tensor_function.h>
+ // The last step is as in all
+ // previous programs:
using namespace dealii;
// @sect3{The <code>MixedLaplaceProblem</code> class template}
// offers such functionality:
#include <base/tensor_function.h>
+ // The last step is as in all
+ // previous programs:
using namespace dealii;
// @sect3{The ``TwoPhaseFlowProblem'' class template}
// double digit time steps.
#include <base/utilities.h>
+ // The last step is as in all
+ // previous programs:
using namespace dealii;
// diameter.
#include <grid/grid_tools.h>
+ // The last step is as in all
+ // previous programs:
using namespace dealii;
// @sect3{The "forward problem" class template}
#include <iostream>
#include <sstream>
+ // The last step is as in all
+ // previous programs:
using namespace dealii;
// The following global variable is used to determine whether the
#include <numeric>
#include <grid/tria_boundary.h>
+ // The last step is as in all
+ // previous programs:
using namespace dealii;
class PointCloudSurface : public StraightBoundary<3>
#include <lac/compressed_sparsity_pattern.h>
+ // The last step is as in all
+ // previous programs:
using namespace dealii;
#include <fstream>
#include <iostream>
+ // ...and this is to import the
+ // deal.II namespace into the global
+ // scope:
using namespace dealii;
// @sect3{The <code>LaplaceProblem</code> class}
// the main() function below for that):
#include <base/logstream.h>
+ // The final step, as in previous
+ // programs, is to import all the
+ // deal.II class and function names
+ // into the global namespace:
using namespace dealii;
// @sect3{The <code>LaplaceProblem</code> class template}
// <code>ostringstream</code>:
#include <sstream>
+ // Finally, this has been discussed
+ // in previous tutorial programs
+ // before:
using namespace dealii;
// problems.
#include <numerics/error_estimator.h>
+ // Finally, this is as in previous
+ // programs:
using namespace dealii;
#include <fstream>
#include <iostream>
+ // The last step is as in previous
+ // programs:
using namespace dealii;
// @sect3{Equation data}
#include <fstream>
#include <iostream>
+ // The last step is as in previous
+ // programs:
using namespace dealii;
// @sect3{The <code>ElasticProblem</code> class template}
# define M_PI 3.14159265358979323846
#endif
+ // The last step is as in previous
+ // programs:
using namespace dealii;
// @sect3{AdvectionProblem class declaration}