Fix CellAccessor::neighbor_child_on_subface()

diff --git a/include/deal.II/grid/reference_cell.h b/include/deal.II/grid/reference_cell.h
index d34ccd09e49c6414131ce65e1bbc567da02d1301..4dfe657f563fbf176a0045787f6e33b56405fc42 100644 (file)
--- a/include/deal.II/grid/reference_cell.h
+++ b/include/deal.II/grid/reference_cell.h
@@ -1104,7 +1104,7 @@ ReferenceCell::child_cell_on_face(
       const bool face_rotation    = Utilities::get_bit(face_orientation_raw, 1);
 
       return GeometryInfo<2>::child_cell_on_face(
-        RefinementCase<2>(RefinementPossibilities<2>::no_refinement),
+        RefinementCase<2>(RefinementPossibilities<2>::isotropic_refinement),
         face,
         subface,
         face_orientation,
@@ -1130,7 +1130,7 @@ ReferenceCell::child_cell_on_face(
       const bool face_rotation    = Utilities::get_bit(face_orientation_raw, 1);
 
       return GeometryInfo<3>::child_cell_on_face(
-        RefinementCase<3>(RefinementPossibilities<3>::no_refinement),
+        RefinementCase<3>(RefinementPossibilities<3>::isotropic_refinement),
         face,
         subface,
         face_orientation,

diff --git a/source/grid/tria_accessor.cc b/source/grid/tria_accessor.cc
index fd7dac0727ca323b257c9820174dcedc461d0666..a9a8655462ceb8429ddb109010fc30eb2f79d038 100644 (file)
--- a/source/grid/tria_accessor.cc
+++ b/source/grid/tria_accessor.cc
@@ -2894,8 +2894,8 @@ CellAccessor<dim, spacedim>::neighbor_child_on_subface(
                   subface;
 
               const unsigned int neighbor_child_index =
-                ReferenceCells::Triangle.child_cell_on_face(neighbor_face,
-                                                            neighbor_subface);
+                neighbor_cell->reference_cell().child_cell_on_face(
+                  neighbor_face, neighbor_subface);
               const TriaIterator<CellAccessor<dim, spacedim>> sub_neighbor =
                 neighbor_cell->child(neighbor_child_index);