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 92a782a12d2d29b32ab5cdf10d244c5a92f6d3b8..182fa44f74c736679a3872118b592d6e5dab2a15 100644 (file)
--- 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 f47d71f4d3c493dbb1165a60aa7d0f414eddd783..afeefa5c637ada4eb1d135d7933e941acc836f29 100644 (file)
--- 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));