From 47f9d16324e76a7341c3173a38614b336dfb35af Mon Sep 17 00:00:00 2001 From: Marc Fehling Date: Tue, 12 May 2020 15:49:09 +0200 Subject: [PATCH] No more restriction on isotropic refinement for hp::Refinement::predict_error(). --- include/deal.II/hp/refinement.h | 2 -- source/hp/refinement.cc | 6 ------ 2 files changed, 8 deletions(-) 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::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)); -- 2.39.5