From: Peter Munch <peterrmuench@gmail.com>
Date: Sat, 7 May 2022 19:17:42 +0000 (+0200)
Subject: Fix FEEvaluation::gather_evaluate() for user-defined batches
X-Git-Tag: v9.4.0-rc1~265^2
X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=refs%2Fpull%2F13688%2Fhead;p=dealii.git

Fix FEEvaluation::gather_evaluate() for user-defined batches
---

diff --git a/include/deal.II/matrix_free/fe_evaluation.h b/include/deal.II/matrix_free/fe_evaluation.h
index 5c80c2f994..fa8aafadaa 100644
--- 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();