velocity. Consequently, we will introduce a member variable of type
MappingQ that
denotes such a mapping (step-10 and step-11 introduce such mappings
-for the first time) and that we will use in all computations on cells
-that are adjacent to the boundary. Since this only affects a
-relatively small fraction of cells, the additional effort is not very
-large and we will take the luxury of using a quartic mapping for these
-cells.
+for the first time) and that we will use in all computations.
<h3> Parallelization on clusters </h3>
//
// - The <code>mapping</code> variable is used to denote a higher-order
// polynomial mapping. As mentioned in the introduction, we use this
- // mapping when forming integrals through quadrature for all cells that
- // are adjacent to either the inner or outer boundaries of our domain
- // where the boundary is curved.
+ // mapping when forming integrals through quadrature for all cells.
//
// - In a bit of naming confusion, you will notice below that some of the
// variables from namespace TrilinosWrappers are taken from namespace
// examples.
//
// The class is constructed in a way to allow for arbitrary order of
- // approximation of both the domain (through high order mapping) and the
+ // approximation of both the domain (through high order mappings) and the
// finite element space. The order of the finite element space and of the
// mapping can be selected in the constructor of the class.
const unsigned int mapping_degree)
: fe(fe_degree)
, dof_handler(tria)
- , mapping(mapping_degree, true)
+ , mapping(mapping_degree)
, wind(dim)
, singular_quadrature_order(5)
, n_cycles(4)