From 529738be296bf98f9233ca819d0d06390bc31570 Mon Sep 17 00:00:00 2001 From: Peter Munch Date: Thu, 9 Dec 2021 17:49:54 +0100 Subject: [PATCH] FEEvaluationBase::apply_hanging_node_constraints(): clean up --- include/deal.II/matrix_free/fe_evaluation.h | 12 +++++------- 1 file changed, 5 insertions(+), 7 deletions(-) diff --git a/include/deal.II/matrix_free/fe_evaluation.h b/include/deal.II/matrix_free/fe_evaluation.h index 954113378d..be9a88ce8f 100644 --- a/include/deal.II/matrix_free/fe_evaluation.h +++ b/include/deal.II/matrix_free/fe_evaluation.h @@ -779,9 +779,8 @@ protected: /** * Apply hanging-node constraints. */ - template void - apply_hanging_node_constraints() const; + apply_hanging_node_constraints(const bool transpose) const; /** * This is the general array for all data fields. @@ -4433,16 +4432,15 @@ template -template inline void FEEvaluationBase:: - apply_hanging_node_constraints() const + apply_hanging_node_constraints(const bool transpose) const { if (this->dof_info == nullptr || this->dof_info->hanging_node_constraint_masks.size() == 0) return; // nothing to do with faces - unsigned int n_vectorization_actual = + const unsigned int n_vectorization_actual = this->dof_info ->n_vectorization_lanes_filled[this->dof_access_index][this->cell]; @@ -4526,7 +4524,7 @@ FEEvaluationBase:: internal::VectorReader reader; read_write_operation(reader, src_data.first, src_data.second, mask, true); - apply_hanging_node_constraints(); + apply_hanging_node_constraints(false); # ifdef DEBUG this->dof_values_initialized = true; @@ -4583,7 +4581,7 @@ FEEvaluationBase:: internal::ExcAccessToUninitializedField()); # endif - apply_hanging_node_constraints(); + apply_hanging_node_constraints(true); const auto dst_data = internal::get_vector_data( dst, -- 2.39.5