From: Martin Kronbichler 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=refs%2Fpull%2F13921%2Fhead;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::isotropic_refinement + 1); AssertIndexRange(face_no, GeometryInfo::faces_per_cell); + // simple special case + if (cell_refinement_case == RefinementCase::cut_xy) + return RefinementCase<1>::cut_x; + const RefinementCase ref_cases[RefinementCase::isotropic_refinement + 1][GeometryInfo::faces_per_cell / 2] = { @@ -3706,6 +3710,10 @@ GeometryInfo<3>::face_refinement_case( RefinementCase::isotropic_refinement + 1); AssertIndexRange(face_no, GeometryInfo::faces_per_cell); + // simple special case + if (cell_refinement_case == RefinementCase::cut_xyz) + return RefinementCase::cut_xy; + const RefinementCase ref_cases[RefinementCase::isotropic_refinement + 1] [GeometryInfo::faces_per_cell / 2] = { @@ -3816,6 +3824,10 @@ GeometryInfo<3>::line_refinement_case( RefinementCase::isotropic_refinement + 1); AssertIndexRange(line_no, GeometryInfo::lines_per_cell); + // simple special case + if (cell_refinement_case == RefinementCase::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