From: Peter Munch <peterrmuench@gmail.com>
Date: Sun, 17 Oct 2021 19:37:51 +0000 (+0200)
Subject: Fix MatrixCreator::create_mass/laplace_matrix for hp
X-Git-Tag: v9.4.0-rc1~931^2
X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=refs%2Fpull%2F12840%2Fhead;p=dealii.git

Fix MatrixCreator::create_mass/laplace_matrix for hp
---

diff --git a/include/deal.II/numerics/matrix_creator.templates.h b/include/deal.II/numerics/matrix_creator.templates.h
index 586ca63beb..ff9f2fac5c 100644
--- a/include/deal.II/numerics/matrix_creator.templates.h
+++ b/include/deal.II/numerics/matrix_creator.templates.h
@@ -640,8 +640,8 @@ namespace MatrixCreator
     Assert(matrix.n() == dof.n_dofs(),
            ExcDimensionMismatch(matrix.n(), dof.n_dofs()));
 
-    hp::FECollection<dim, spacedim>      fe_collection(dof.get_fe());
-    hp::QCollection<dim>                 q_collection(q);
+    const auto &         fe_collection = dof.get_fe_collection();
+    hp::QCollection<dim> q_collection(q);
     hp::MappingCollection<dim, spacedim> mapping_collection(mapping);
     MatrixCreator::internal::AssemblerData::Scratch<dim, spacedim, number>
       assembler_data(fe_collection,
@@ -715,8 +715,8 @@ namespace MatrixCreator
     Assert(matrix.n() == dof.n_dofs(),
            ExcDimensionMismatch(matrix.n(), dof.n_dofs()));
 
-    hp::FECollection<dim, spacedim>      fe_collection(dof.get_fe());
-    hp::QCollection<dim>                 q_collection(q);
+    const auto &         fe_collection = dof.get_fe_collection();
+    hp::QCollection<dim> q_collection(q);
     hp::MappingCollection<dim, spacedim> mapping_collection(mapping);
     MatrixCreator::internal::AssemblerData::Scratch<dim, spacedim, number>
       assembler_data(fe_collection,
@@ -1871,8 +1871,8 @@ namespace MatrixCreator
     Assert(matrix.n() == dof.n_dofs(),
            ExcDimensionMismatch(matrix.n(), dof.n_dofs()));
 
-    hp::FECollection<dim, spacedim>      fe_collection(dof.get_fe());
-    hp::QCollection<dim>                 q_collection(q);
+    const auto &         fe_collection = dof.get_fe_collection();
+    hp::QCollection<dim> q_collection(q);
     hp::MappingCollection<dim, spacedim> mapping_collection(mapping);
     MatrixCreator::internal::AssemblerData::Scratch<dim, spacedim, double>
       assembler_data(fe_collection,
@@ -1944,8 +1944,8 @@ namespace MatrixCreator
     Assert(matrix.n() == dof.n_dofs(),
            ExcDimensionMismatch(matrix.n(), dof.n_dofs()));
 
-    hp::FECollection<dim, spacedim>      fe_collection(dof.get_fe());
-    hp::QCollection<dim>                 q_collection(q);
+    const auto &         fe_collection = dof.get_fe_collection();
+    hp::QCollection<dim> q_collection(q);
     hp::MappingCollection<dim, spacedim> mapping_collection(mapping);
     MatrixCreator::internal::AssemblerData::Scratch<dim, spacedim, double>
       assembler_data(fe_collection,