From: Maximilian Bergbauer Date: Tue, 20 May 2025 14:24:14 +0000 (+0200) Subject: Allow adaptive meshes in MappingInfo for empty quadratures X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=refs%2Fpull%2F18476%2Fhead;p=dealii.git Allow adaptive meshes in MappingInfo for empty quadratures --- 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);