From: Martin Kronbichler Date: Thu, 30 Jan 2025 09:28:14 +0000 (+0100) Subject: Update a comment X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=536b07eb6d3249a93c90daa665a2eebbf3617355;p=dealii.git Update a comment --- 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; }