]> https://gitweb.dealii.org/ - release-papers.git/commitdiff
Fix fpe code
authorpeterrum <peterrmuench@aol.com>
Sat, 19 Jun 2021 13:30:25 +0000 (15:30 +0200)
committerpeterrum <peterrmuench@aol.com>
Sat, 19 Jun 2021 13:30:25 +0000 (15:30 +0200)
9.3/paper.tex

index 1b7505578c988715c7827e14bd02a074d3cffa04..e99c8d5762d7ed35f768a1aa7b5fd7598f44b439 100644 (file)
@@ -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

In the beginning the Universe was created. This has made a lot of people very angry and has been widely regarded as a bad move.

Douglas Adams


Typeset in Trocchi and Trocchi Bold Sans Serif.