From: Maximilian Bergbauer Date: Thu, 16 Mar 2023 16:35:55 +0000 (+0100) Subject: Select correct code path for neighbor elements in FEEvaluation X-Git-Tag: v9.5.0-rc1~455^2 X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=refs%2Fpull%2F14903%2Fhead;p=dealii.git Select correct code path for neighbor elements in FEEvaluation --- diff --git a/include/deal.II/matrix_free/fe_evaluation.h b/include/deal.II/matrix_free/fe_evaluation.h index e08b6b7b08..f415c98261 100644 --- a/include/deal.II/matrix_free/fe_evaluation.h +++ b/include/deal.II/matrix_free/fe_evaluation.h @@ -3422,9 +3422,7 @@ FEEvaluationBase:: internal::is_vectorizable::value> vector_selector; - const bool is_neighbor_cells = !is_face && !this->is_interior_face(); - const bool use_vectorized_path = - !(masking_is_active || is_neighbor_cells || has_hn_constraints); + const bool use_vectorized_path = !(masking_is_active || has_hn_constraints); const std::size_t dofs_per_component = this->data->dofs_per_component_on_cell; std::array values_dofs; diff --git a/tests/matrix_free/fe_evaluation_shift.cc b/tests/matrix_free/fe_evaluation_shift.cc index 31b501fdd5..320f69913f 100644 --- a/tests/matrix_free/fe_evaluation_shift.cc +++ b/tests/matrix_free/fe_evaluation_shift.cc @@ -185,7 +185,7 @@ namespace dealii this->cell_type = this->matrix_free->get_mapping_info().cell_type[cell_batch_index]; - this->cell = cell_batch_index; + this->cell = numbers::invalid_unsigned_int; this->face_orientations[0] = 0; this->subface_index = GeometryInfo::max_children_per_cell; this->face_numbers[0] = face_number;