From: Maximilian Bergbauer <maximilian.bergbauer@tum.de>
Date: Fri, 1 Dec 2023 22:57:13 +0000 (+0100)
Subject: Change assert
X-Git-Tag: relicensing~264^2
X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=refs%2Fpull%2F16318%2Fhead;p=dealii.git

Change assert
---

diff --git a/include/deal.II/matrix_free/matrix_free.templates.h b/include/deal.II/matrix_free/matrix_free.templates.h
index 2477979744..c7bf9a972e 100644
--- a/include/deal.II/matrix_free/matrix_free.templates.h
+++ b/include/deal.II/matrix_free/matrix_free.templates.h
@@ -308,9 +308,6 @@ MatrixFree<dim, Number, VectorizedArrayType>::get_face_iterator(
   const unsigned int fe_component) const
 {
   AssertIndexRange(fe_component, dof_handlers.size());
-  AssertIndexRange(face_batch_index,
-                   n_inner_face_batches() +
-                     (interior ? n_boundary_face_batches() : 0));
 
   AssertIndexRange(lane_index,
                    n_active_entries_per_face_batch(face_batch_index));
@@ -323,6 +320,11 @@ MatrixFree<dim, Number, VectorizedArrayType>::get_face_iterator(
                                     face2cell_info.cells_interior[lane_index] :
                                     face2cell_info.cells_exterior[lane_index];
 
+  Assert(cell_index != numbers::invalid_unsigned_int,
+         ExcMessage(
+           "Invalid cell index for requested face. Exterior cells cannot be"
+           " accessed on boundaries."));
+
   std::pair<unsigned int, unsigned int> index = cell_level_index[cell_index];
 
   return {typename DoFHandler<dim>::cell_iterator(