]> https://gitweb.dealii.org/ - dealii.git/commitdiff
No more restriction on isotropic refinement for hp::Refinement::predict_error().
authorMarc Fehling <marc.fehling@gmx.net>
Tue, 12 May 2020 13:49:09 +0000 (15:49 +0200)
committerMarc Fehling <marc.fehling@gmx.net>
Tue, 12 May 2020 13:49:18 +0000 (15:49 +0200)
include/deal.II/hp/refinement.h
source/hp/refinement.cc

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

In the beginning the Universe was created. This has made a lot of people very angry and has been widely regarded as a bad move.

Douglas Adams


Typeset in Trocchi and Trocchi Bold Sans Serif.