From: Maximilian Bergbauer 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::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(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:: template interpolate(1, evaluation_flags,