From: Marc Fehling <marc.fehling@gmx.net>
Date: Tue, 12 May 2020 13:49:09 +0000 (+0200)
Subject: No more restriction on isotropic refinement for hp::Refinement::predict_error().
X-Git-Tag: v9.2.0-rc1~13^2~1
X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=47f9d16324e76a7341c3173a38614b336dfb35af;p=dealii.git

No more restriction on isotropic refinement for hp::Refinement::predict_error().
---

diff --git a/include/deal.II/hp/refinement.h b/include/deal.II/hp/refinement.h
index 92a782a12d..182fa44f74 100644
--- a/include/deal.II/hp/refinement.h
+++ b/include/deal.II/hp/refinement.h
@@ -554,8 +554,6 @@ namespace hp
      * $\gamma_\text{p}^2 = 0.4$, $\gamma_\text{h}^2 = 4$,
      * $\gamma_\text{n}^2 = 1$.
      *
-     * @note This feature is currently only implemented for isotropic refinement.
-     *
      * @note We want to predict the error by how adaptation will actually happen.
      *   Thus, this function needs to be called after
      *   Triangulation::prepare_coarsening_and_refinement().
diff --git a/source/hp/refinement.cc b/source/hp/refinement.cc
index f47d71f4d3..afeefa5c63 100644
--- a/source/hp/refinement.cc
+++ b/source/hp/refinement.cc
@@ -599,12 +599,6 @@ namespace hp
             // step 2: algebraic decay with h-adaptation
             if (cell->refine_flag_set())
               {
-                Assert(
-                  cell->refine_flag_set() ==
-                    RefinementCase<dim>::isotropic_refinement,
-                  ExcMessage(
-                    "Error prediction is only valid for isotropic refinement!"));
-
                 predicted_errors[cell->active_cell_index()] *=
                   (gamma_h * std::pow(.5, future_fe_degree));