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 e08b6b7b08f54b4a9a0cb9828444aab013ae6d2a..f415c9826194c0978cb38a9028f7901646470b94 100644 (file)
--- a/include/deal.II/matrix_free/fe_evaluation.h
+++ b/include/deal.II/matrix_free/fe_evaluation.h
@@ -3422,9 +3422,7 @@ FEEvaluationBase<dim, n_components_, Number, is_face, VectorizedArrayType>::
                          internal::is_vectorizable<VectorType, Number>::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<VectorizedArrayType *, n_components> values_dofs;

diff --git a/tests/matrix_free/fe_evaluation_shift.cc b/tests/matrix_free/fe_evaluation_shift.cc
index 31b501fdd5334a2c6acdcf07dc2707c081d3d4d2..320f69913f80eb9dfa27a317caaef556bd7d214f 100644 (file)
--- 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<dim>::max_children_per_cell;
     this->face_numbers[0]      = face_number;