Merge pull request #13688 from peterrum/fe_eval_user_define_batches_g…

…ather_evaluate Fix FEEvaluation::gather_evaluate() for user-defined batches
dealii · May 8, 2022 · 782b4a5 · 782b4a5
2 parents 3e39d6d + 718a3f6
commit 782b4a5
Showing 1 changed file with 4 additions and 0 deletions.
diff --git a/include/deal.II/matrix_free/fe_evaluation.h b/include/deal.II/matrix_free/fe_evaluation.h
@@ -7269,6 +7269,10 @@ namespace internal
   VectorizedArrayType *
   check_vector_access_inplace(const EvaluatorType &fe_eval, VectorType &vector)
   {
+    // for user-defined cell batches this functionality is not supported
+    if (fe_eval.get_current_cell_index() == numbers::invalid_unsigned_int)
+      return nullptr;
+
     const unsigned int cell     = fe_eval.get_cell_or_face_batch_id();
     const auto &       dof_info = fe_eval.get_dof_info();