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

dealii · May 7, 2022 · c565b19 · c565b19
1 parent 9fe9ba2
commit c565b19
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 macro cells 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();