From 536b07eb6d3249a93c90daa665a2eebbf3617355 Mon Sep 17 00:00:00 2001 From: Martin Kronbichler Date: Thu, 30 Jan 2025 10:28:14 +0100 Subject: [PATCH] Update a comment --- include/deal.II/matrix_free/face_setup_internal.h | 7 ++++--- 1 file changed, 4 insertions(+), 3 deletions(-) diff --git a/include/deal.II/matrix_free/face_setup_internal.h b/include/deal.II/matrix_free/face_setup_internal.h index 80eacba4a6..ff07874873 100644 --- a/include/deal.II/matrix_free/face_setup_internal.h +++ b/include/deal.II/matrix_free/face_setup_internal.h @@ -1048,14 +1048,15 @@ namespace internal info.subface_index = GeometryInfo::max_children_per_cell; Assert(neighbor->level() <= cell->level(), ExcInternalError()); - if (cell->level() > neighbor->level()) + + // for dim > 1 and hanging faces we must set a subface index + if (dim > 1 && cell->level() > neighbor->level()) { if (cell->has_periodic_neighbor(face_no)) info.subface_index = cell->periodic_neighbor_of_coarser_periodic_neighbor(face_no) .second; - // else if must be dim > 1 because dim == 1 adm creates no subfaces - else if (dim > 1) + else info.subface_index = cell->neighbor_of_coarser_neighbor(face_no).second; } -- 2.39.5