From: Peter Munch <peterrmuench@gmail.com>
Date: Mon, 22 Nov 2021 10:32:36 +0000 (+0100)
Subject: Fix CellAccessor::neighbor_child_on_subface()
X-Git-Tag: v9.4.0-rc1~811^2
X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=d9cb20b3042eadef4dc6e1b28a07429fdc341621;p=dealii.git

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 d34ccd09e4..4dfe657f56 100644
--- 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 fd7dac0727..a9a8655462 100644
--- 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);