]> https://gitweb.dealii.org/ - dealii.git/commitdiff
Select correct code path for neighbor elements in FEEvaluation 14903/head
authorMaximilian Bergbauer <maximilian.bergbauer@tum.de>
Thu, 16 Mar 2023 16:35:55 +0000 (17:35 +0100)
committerMaximilian Bergbauer <maximilian.bergbauer@tum.de>
Fri, 17 Mar 2023 15:15:59 +0000 (16:15 +0100)
include/deal.II/matrix_free/fe_evaluation.h
tests/matrix_free/fe_evaluation_shift.cc

index e08b6b7b08f54b4a9a0cb9828444aab013ae6d2a..f415c9826194c0978cb38a9028f7901646470b94 100644 (file)
@@ -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;
index 31b501fdd5334a2c6acdcf07dc2707c081d3d4d2..320f69913f80eb9dfa27a317caaef556bd7d214f 100644 (file)
@@ -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;

In the beginning the Universe was created. This has made a lot of people very angry and has been widely regarded as a bad move.

Douglas Adams


Typeset in Trocchi and Trocchi Bold Sans Serif.