From: peterrum <peterrmuench@aol.com>
Date: Sat, 19 Jun 2021 13:30:25 +0000 (+0200)
Subject: Fix fpe code
X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=14f31fbf2afa7a7a5a76e6786ad85e9c01417428;p=release-papers.git

Fix fpe code
---

diff --git a/9.3/paper.tex b/9.3/paper.tex
index 1b75055..e99c8d5 100644
--- a/9.3/paper.tex
+++ b/9.3/paper.tex
@@ -935,12 +935,17 @@ in general, located differently within $K$ than for any other cell:
 \end{align*}
 In \dealii{}, this can now be conveniently written as
 \begin{c++}
+phi_curvature.reinit(cell, reference_points);
 phi_normal.reinit(cell, reference_points);
 phi_force.reinit(cell, reference_points);
+
+phi_curvature.evaluate(curvature_values, EvaluationFlags::values);
 phi_normal.evaluate(normal_values, EvaluationFlags::values);
+
 for (unsigned int q = 0; q < n_points; ++q)
-   phi_force.submit_value(phi_curvature.get_value(q) *
-                          phi_normal.get_value(q) * JxW[q], q);
+  phi_force.submit_value(phi_curvature.get_value(q) *
+                         phi_normal.get_value(q) * JxW[q], q);
+                          
 phi_force.integrate(force_values, EvaluationFlags::values);
 \end{c++}
 The quadrature points (at reference positions