]> https://gitweb.dealii.org/ - dealii.git/commitdiff
Allow adaptive meshes in MappingInfo for empty quadratures 18476/head
authorMaximilian Bergbauer <maximilian.bergbauer@tum.de>
Tue, 20 May 2025 14:24:14 +0000 (16:24 +0200)
committerMaximilian Bergbauer <maximilian.bergbauer@tum.de>
Tue, 20 May 2025 14:24:14 +0000 (16:24 +0200)
include/deal.II/non_matching/mapping_info.h

index 946843b920b6f91f4d2ea3da1d1d0041748ee0de..d7022f31c1c52f3ccc2276c146b1ab8ed3bc6c35 100644 (file)
@@ -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);
 

In the beginning the Universe was created. This has made a lot of people very angry and has been widely regarded as a bad move.

Douglas Adams


Typeset in Trocchi and Trocchi Bold Sans Serif.