From: Maximilian Bergbauer Date: Tue, 21 Mar 2023 10:25:06 +0000 (+0100) Subject: Use reinit() function internal X-Git-Tag: v9.5.0-rc1~413^2 X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=refs%2Fpull%2F14926%2Fhead;p=dealii.git Use reinit() function internal --- 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::evaluate( const EvaluationFlags::EvaluationFlags &evaluation_flag) { if (!is_reinitialized) - { - const_cast(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()); - if (update_flags & update_gradients) - gradients.resize(n_q_points, numbers::signaling_nan()); - } + reinit(numbers::invalid_unsigned_int); if (n_q_points == 0) return; @@ -1121,15 +1113,7 @@ FEPointEvaluation::integrate( const EvaluationFlags::EvaluationFlags &integration_flags) { if (!is_reinitialized) - { - const_cast(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()); - if (update_flags & update_gradients) - gradients.resize(n_q_points, numbers::signaling_nan()); - } + reinit(numbers::invalid_unsigned_int); if (n_q_points == 0) // no evaluation points provided {