Merge pull request #14142 from peterrum/fpe_grad_type

Fix type in FPE
dealii · Jul 18, 2022 · 1967d5b · 1967d5b
2 parents 2163360 + 1d271bf
commit 1967d5b
Showing 1 changed file with 2 additions and 2 deletions.
diff --git a/include/deal.II/matrix_free/fe_point_evaluation.h b/include/deal.II/matrix_free/fe_point_evaluation.h
@@ -951,7 +951,7 @@ FEPointEvaluation<n_components, dim, spacedim, Number>::evaluate(
                     static_cast<typename internal::FEPointEvaluation::
                                   EvaluatorTypeTraits<dim,
                                                       n_components,
-                                                      double>::gradient_type>(
+                                                      Number>::gradient_type>(
                       apply_transformation(mapping_info->get_mapping_data()
                                              .inverse_jacobians[i + j]
                                              .transpose(),
@@ -1093,7 +1093,7 @@ FEPointEvaluation<n_components, dim, spacedim, Number>::integrate(
               {
                 gradients[i + j] =
                   static_cast<typename internal::FEPointEvaluation::
-                                EvaluatorTypeTraits<dim, n_components, double>::
+                                EvaluatorTypeTraits<dim, n_components, Number>::
                                   gradient_type>(apply_transformation(
                     mapping_info->get_mapping_data().inverse_jacobians[i + j],
                     gradients[i + j]));