From: Martin Kronbichler <martin.kronbichler@uni-a.de>
Date: Tue, 7 Jun 2022 21:02:22 +0000 (+0200)
Subject: Use shortcase when determining isotropic refinement cases
X-Git-Tag: v9.5.0-rc1~1171^2
X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=d4ae87c8a9bbe0a0f683755a6c75992bd9a7c266;p=dealii.git

Use shortcase when determining isotropic refinement cases
---

diff --git a/include/deal.II/base/geometry_info.h b/include/deal.II/base/geometry_info.h
index 4c778352d0..ab7ac139e4 100644
--- a/include/deal.II/base/geometry_info.h
+++ b/include/deal.II/base/geometry_info.h
@@ -3675,6 +3675,10 @@ GeometryInfo<2>::face_refinement_case(
                    RefinementCase<dim>::isotropic_refinement + 1);
   AssertIndexRange(face_no, GeometryInfo<dim>::faces_per_cell);
 
+  // simple special case
+  if (cell_refinement_case == RefinementCase<dim>::cut_xy)
+    return RefinementCase<1>::cut_x;
+
   const RefinementCase<dim - 1>
     ref_cases[RefinementCase<dim>::isotropic_refinement +
               1][GeometryInfo<dim>::faces_per_cell / 2] = {
@@ -3706,6 +3710,10 @@ GeometryInfo<3>::face_refinement_case(
                    RefinementCase<dim>::isotropic_refinement + 1);
   AssertIndexRange(face_no, GeometryInfo<dim>::faces_per_cell);
 
+  // simple special case
+  if (cell_refinement_case == RefinementCase<dim>::cut_xyz)
+    return RefinementCase<dim - 1>::cut_xy;
+
   const RefinementCase<dim - 1>
     ref_cases[RefinementCase<dim>::isotropic_refinement + 1]
              [GeometryInfo<dim>::faces_per_cell / 2] = {
@@ -3816,6 +3824,10 @@ GeometryInfo<3>::line_refinement_case(
                    RefinementCase<dim>::isotropic_refinement + 1);
   AssertIndexRange(line_no, GeometryInfo<dim>::lines_per_cell);
 
+  // simple special case
+  if (cell_refinement_case == RefinementCase<dim>::cut_xyz)
+    return RefinementCase<1>::cut_x;
+
   // array indicating, which simple refine
   // case cuts a line in direction x, y or
   // z. For example, cut_y and everything