From 3883e8e48af5902d086ff5746dbb4f12da706844 Mon Sep 17 00:00:00 2001 From: Maximilian Bergbauer Date: Fri, 1 Dec 2023 23:57:13 +0100 Subject: [PATCH] Change assert --- include/deal.II/matrix_free/matrix_free.templates.h | 8 +++++--- 1 file changed, 5 insertions(+), 3 deletions(-) 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::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::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 index = cell_level_index[cell_index]; return {typename DoFHandler::cell_iterator( -- 2.39.5