Make fe constant

diff --git a/include/deal.II/matrix_free/util.h b/include/deal.II/matrix_free/util.h
index de7ce1723ce792ee4bd5a4fa45dd7180e7bc575e..2a3cc35add7347e31b1fa441000aec008cea7546 100644 (file)
--- a/include/deal.II/matrix_free/util.h
+++ b/include/deal.II/matrix_free/util.h
@@ -61,7 +61,7 @@ namespace internal
 
           for (unsigned int i = 1; i <= 3; ++i)
             {
-              FE_SimplexP<dim> fe(i);
+              const FE_SimplexP<dim> fe(i);
               if (quad == Quadrature<dim>(fe.get_unit_support_points()))
                 return {ReferenceCells::get_simplex<dim>(),
                         dealii::hp::QCollection<dim - 1>(Quadrature<dim - 1>(
@@ -152,7 +152,7 @@ namespace internal
 
           for (unsigned int i = 1; i <= 3; ++i)
             {
-              FE_SimplexP<dim> fe(i);
+              const FE_SimplexP<dim> fe(i);
               if (quad == Quadrature<dim>(fe.get_unit_support_points()))
                 {
                   if (dim == 2)