From 8ce64c0f17e89653ee48e0f3015308a48d4b4686 Mon Sep 17 00:00:00 2001
From: Maximilian Bergbauer <maximilian.bergbauer@tum.de>
Date: Tue, 21 Mar 2023 11:25:06 +0100
Subject: [PATCH] Use reinit() function internal

---
 .../deal.II/matrix_free/fe_point_evaluation.h | 20 ++-----------------
 1 file changed, 2 insertions(+), 18 deletions(-)

diff --git a/include/deal.II/matrix_free/fe_point_evaluation.h b/include/deal.II/matrix_free/fe_point_evaluation.h
index 1f55b1487a..958ce2e29a 100644
--- a/include/deal.II/matrix_free/fe_point_evaluation.h
+++ b/include/deal.II/matrix_free/fe_point_evaluation.h
@@ -970,15 +970,7 @@ FEPointEvaluation<n_components, dim, spacedim, Number>::evaluate(
   const EvaluationFlags::EvaluationFlags &evaluation_flag)
 {
   if (!is_reinitialized)
-    {
-      const_cast<unsigned int &>(n_q_points) =
-        mapping_info->get_unit_points(current_cell_index, current_face_number)
-          .size();
-      if (update_flags & update_values)
-        values.resize(n_q_points, numbers::signaling_nan<value_type>());
-      if (update_flags & update_gradients)
-        gradients.resize(n_q_points, numbers::signaling_nan<gradient_type>());
-    }
+    reinit(numbers::invalid_unsigned_int);
 
   if (n_q_points == 0)
     return;
@@ -1121,15 +1113,7 @@ FEPointEvaluation<n_components, dim, spacedim, Number>::integrate(
   const EvaluationFlags::EvaluationFlags &integration_flags)
 {
   if (!is_reinitialized)
-    {
-      const_cast<unsigned int &>(n_q_points) =
-        mapping_info->get_unit_points(current_cell_index, current_face_number)
-          .size();
-      if (update_flags & update_values)
-        values.resize(n_q_points, numbers::signaling_nan<value_type>());
-      if (update_flags & update_gradients)
-        gradients.resize(n_q_points, numbers::signaling_nan<gradient_type>());
-    }
+    reinit(numbers::invalid_unsigned_int);
 
   if (n_q_points == 0) // no evaluation points provided
     {
-- 
2.39.5