From: Timo Heister <timo.heister@gmail.com>
Date: Wed, 23 Jun 2021 17:25:23 +0000 (-0400)
Subject: simplex code suggestions
X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=refs%2Fpull%2F72%2Fhead;p=release-papers.git

simplex code suggestions
---

diff --git a/9.3/paper.tex b/9.3/paper.tex
index d515dff..1c484ce 100644
--- a/9.3/paper.tex
+++ b/9.3/paper.tex
@@ -606,12 +606,12 @@ quadrature class as follows:
 \begin{c++}
 FE_SimplexP<dim, spacedim> fe(degree);
 MappingFE<dim, spacedim> mapping(FE_SimplexP<dim, spacedim>(1));
-QGaussSimplex<dim> quad(degree + 1);
+QGaussSimplex<dim> quadrature(degree + 1);
 
-DoFHandler<dim> dof_handler(tria);
+DoFHandler<dim, spacedim> dof_handler(tria);
 dof_handler.distribute_dofs(fe);
 
-FEValues<dim, spacedim> fe_values(mapping, fe, quad, flags);
+FEValues<dim, spacedim> fe_values(mapping, fe, quadrature, flags);
 \end{c++}
 The list of currently supported finite element classes is provided in Table~\ref{tab:simplex:fe}. Currently,
 only linear and quadratic mappings via the \texttt{MappingFE} and
@@ -685,11 +685,12 @@ Furthermore, for mixed meshes, the number of degrees of freedom will differ betw
 cell-local arrays need to be resized for each
 cell (as has previously already been the case in the $hp$-context):
 \begin{c++}
-std::vector<double> local_rhs;
+Vector<double> local_rhs;
 for (const auto &cell : dof_handler.active_cell_iterators())
  {
    hp_fe_values.reinit(cell);
    local_rhs.resize(cell->get_fe().n_dofs_per_cell());
+   // ...
  }
 \end{c++}