From 14f31fbf2afa7a7a5a76e6786ad85e9c01417428 Mon Sep 17 00:00:00 2001
From: peterrum <peterrmuench@aol.com>
Date: Sat, 19 Jun 2021 15:30:25 +0200
Subject: [PATCH] Fix fpe code

---
 9.3/paper.tex | 9 +++++++--
 1 file changed, 7 insertions(+), 2 deletions(-)

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
-- 
2.39.5