From 61fc51f67187623aeefec29288fd507c8b7ca7d7 Mon Sep 17 00:00:00 2001 From: Maximilian Bergbauer Date: Tue, 20 May 2025 16:24:14 +0200 Subject: [PATCH] Allow adaptive meshes in MappingInfo for empty quadratures --- include/deal.II/non_matching/mapping_info.h | 8 +++++++- 1 file changed, 7 insertions(+), 1 deletion(-) diff --git a/include/deal.II/non_matching/mapping_info.h b/include/deal.II/non_matching/mapping_info.h index 946843b920..d7022f31c1 100644 --- a/include/deal.II/non_matching/mapping_info.h +++ b/include/deal.II/non_matching/mapping_info.h @@ -1532,7 +1532,13 @@ namespace NonMatching const auto &cell_p = cell_m->at_boundary(f_m) ? cell_m : cell_m->neighbor(f_m); const auto f_p = - cell_m->at_boundary(f_m) ? f_m : cell_m->neighbor_of_neighbor(f_m); + cell_m->at_boundary(f_m) ? f_m : cell_m->neighbor_face_no(f_m); + + Assert( + empty || (cell_m->level() == cell_p->level()), + ExcMessage( + "Intersected faces with quadrature points need to have the same " + "refinement level!")); face_number.emplace_back(f_m, f_p); -- 2.39.5