Skip to content

Commit

Permalink
Merge pull request #15981 from kronbichler/remove_indirection
Browse files Browse the repository at this point in the history
Matrix-free FEEvaluation: Avoid indirection in call to sum factorization
  • Loading branch information
tamiko committed Sep 12, 2023
2 parents ca43d31 + 2ed155d commit bcc7101
Showing 1 changed file with 11 additions and 9 deletions.
20 changes: 11 additions & 9 deletions include/deal.II/matrix_free/fe_evaluation.h
Original file line number Diff line number Diff line change
Expand Up @@ -31,7 +31,6 @@

#include <deal.II/matrix_free/evaluation_flags.h>
#include <deal.II/matrix_free/evaluation_kernels.h>
#include <deal.II/matrix_free/evaluation_selector.h>
#include <deal.II/matrix_free/evaluation_template_factory.h>
#include <deal.II/matrix_free/fe_evaluation_data.h>
#include <deal.II/matrix_free/hanging_nodes_internal.h>
Expand Down Expand Up @@ -7854,8 +7853,11 @@ FEEvaluation<dim,

if constexpr (fe_degree > -1)
{
SelectEvaluator<dim, fe_degree, n_q_points_1d, VectorizedArrayType>::
evaluate(n_components, evaluation_flag_actual, values_array, *this);
internal::FEEvaluationImplSelector<dim, VectorizedArrayType, false>::
template run<fe_degree, n_q_points_1d>(n_components,
evaluation_flag_actual,
values_array,
*this);
}
else
{
Expand Down Expand Up @@ -8084,12 +8086,12 @@ FEEvaluation<dim,

if constexpr (fe_degree > -1)
{
SelectEvaluator<dim, fe_degree, n_q_points_1d, VectorizedArrayType>::
integrate(n_components,
integration_flag_actual,
values_array,
*this,
sum_into_values_array);
internal::FEEvaluationImplSelector<dim, VectorizedArrayType, true>::
template run<fe_degree, n_q_points_1d>(n_components,
integration_flag_actual,
values_array,
*this,
sum_into_values_array);
}
else
{
Expand Down

0 comments on commit bcc7101

Please sign in to comment.