From: Maximilian Bergbauer <maximilian.bergbauer@tum.de>
Date: Wed, 31 May 2023 11:58:42 +0000 (+0200)
Subject: Suggestion from review
X-Git-Tag: v9.5.0-rc1~178^2~1
X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=ddab6298660c2cce811012535e52304942b3f5db;p=dealii.git

Suggestion from review
---

diff --git a/include/deal.II/matrix_free/fe_point_evaluation.h b/include/deal.II/matrix_free/fe_point_evaluation.h
index 25bb4fc7b3..b96ace8a09 100644
--- a/include/deal.II/matrix_free/fe_point_evaluation.h
+++ b/include/deal.II/matrix_free/fe_point_evaluation.h
@@ -1693,20 +1693,22 @@ FEPointEvaluation<n_components_, dim, spacedim, Number>::prepare_evaluate_fast(
     {
       if (use_face_path)
         {
-          ScalarNumber *input  = scratch_data_scalar.begin();
-          ScalarNumber *output = input + dofs_per_component;
-
+          const ScalarNumber *input;
           if (renumber.empty())
-            input = const_cast<ScalarNumber *>(solution_values.data());
+            input = solution_values.data();
           else
             {
               const unsigned int *renumber_ptr =
                 renumber.data() +
                 (component_in_base_element + comp) * dofs_per_component;
               for (unsigned int i = 0; i < dofs_per_component; ++i)
-                input[i] = solution_values[renumber_ptr[i]];
+                scratch_data_scalar[i] = solution_values[renumber_ptr[i]];
+              input = scratch_data_scalar.data();
             }
 
+          ScalarNumber *output =
+            scratch_data_scalar.begin() + dofs_per_component;
+
           internal::FEFaceNormalEvaluationImpl<dim, -1, ScalarNumber>::
             template interpolate<true, false>(1,
                                               evaluation_flags,