From: Timo Heister 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%2Fheads%2Fsimplex_code;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 fe(degree); MappingFE mapping(FE_SimplexP(1)); -QGaussSimplex quad(degree + 1); +QGaussSimplex quadrature(degree + 1); -DoFHandler dof_handler(tria); +DoFHandler dof_handler(tria); dof_handler.distribute_dofs(fe); -FEValues fe_values(mapping, fe, quad, flags); +FEValues 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 local_rhs; +Vector 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++}