From: Peter Munch 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 fe_collection(dof.get_fe()); - hp::QCollection q_collection(q); + const auto & fe_collection = dof.get_fe_collection(); + hp::QCollection q_collection(q); hp::MappingCollection mapping_collection(mapping); MatrixCreator::internal::AssemblerData::Scratch assembler_data(fe_collection, @@ -715,8 +715,8 @@ namespace MatrixCreator Assert(matrix.n() == dof.n_dofs(), ExcDimensionMismatch(matrix.n(), dof.n_dofs())); - hp::FECollection fe_collection(dof.get_fe()); - hp::QCollection q_collection(q); + const auto & fe_collection = dof.get_fe_collection(); + hp::QCollection q_collection(q); hp::MappingCollection mapping_collection(mapping); MatrixCreator::internal::AssemblerData::Scratch assembler_data(fe_collection, @@ -1871,8 +1871,8 @@ namespace MatrixCreator Assert(matrix.n() == dof.n_dofs(), ExcDimensionMismatch(matrix.n(), dof.n_dofs())); - hp::FECollection fe_collection(dof.get_fe()); - hp::QCollection q_collection(q); + const auto & fe_collection = dof.get_fe_collection(); + hp::QCollection q_collection(q); hp::MappingCollection mapping_collection(mapping); MatrixCreator::internal::AssemblerData::Scratch assembler_data(fe_collection, @@ -1944,8 +1944,8 @@ namespace MatrixCreator Assert(matrix.n() == dof.n_dofs(), ExcDimensionMismatch(matrix.n(), dof.n_dofs())); - hp::FECollection fe_collection(dof.get_fe()); - hp::QCollection q_collection(q); + const auto & fe_collection = dof.get_fe_collection(); + hp::QCollection q_collection(q); hp::MappingCollection mapping_collection(mapping); MatrixCreator::internal::AssemblerData::Scratch assembler_data(fe_collection,