From 7d78733bedb69892fffb8657d423d0c40dc9764c Mon Sep 17 00:00:00 2001 From: David Wells Date: Mon, 27 Jan 2025 15:41:47 -0500 Subject: [PATCH] MatrixFreeFunctions::HangingNodes: reuse cell->get_fe(). --- .../deal.II/matrix_free/hanging_nodes_internal.h | 16 +++++++--------- 1 file changed, 7 insertions(+), 9 deletions(-) diff --git a/include/deal.II/matrix_free/hanging_nodes_internal.h b/include/deal.II/matrix_free/hanging_nodes_internal.h index 5172d26825..b6cb9b45fa 100644 --- a/include/deal.II/matrix_free/hanging_nodes_internal.h +++ b/include/deal.II/matrix_free/hanging_nodes_internal.h @@ -626,7 +626,6 @@ namespace internal return; const auto &fe = cell->get_fe(); - AssertDimension(fe.n_unique_faces(), 1); std::vector> @@ -640,14 +639,14 @@ namespace internal std::vector idx_offset = {0}; for (unsigned int base_element_index = 0; - base_element_index < cell->get_fe().n_base_elements(); + base_element_index < fe.n_base_elements(); ++base_element_index) for (unsigned int c = 0; - c < cell->get_fe().element_multiplicity(base_element_index); + c < fe.element_multiplicity(base_element_index); ++c) idx_offset.push_back( idx_offset.back() + - cell->get_fe().base_element(base_element_index).n_dofs_per_cell()); + fe.base_element(base_element_index).n_dofs_per_cell()); std::vector neighbor_dofs_all(idx_offset.back()); std::vector neighbor_dofs_all_temp( @@ -712,10 +711,10 @@ namespace internal index = partitioner->global_to_local(index); for (unsigned int base_element_index = 0, comp = 0; - base_element_index < cell->get_fe().n_base_elements(); + base_element_index < fe.n_base_elements(); ++base_element_index) for (unsigned int c = 0; - c < cell->get_fe().element_multiplicity(base_element_index); + c < fe.element_multiplicity(base_element_index); ++c, ++comp) { if (supported_components[cell->active_fe_index()][comp] == @@ -819,11 +818,10 @@ namespace internal neighbor_cell.line_orientation(local_line_neighbor); for (unsigned int base_element_index = 0, comp = 0; - base_element_index < cell->get_fe().n_base_elements(); + base_element_index < fe.n_base_elements(); ++base_element_index) for (unsigned int c = 0; - c < - cell->get_fe().element_multiplicity(base_element_index); + c < fe.element_multiplicity(base_element_index); ++c, ++comp) { if (supported_components[cell->active_fe_index()][comp] == -- 2.39.5