diff --git a/dune/dorie/interface/adaptivity.hh b/dune/dorie/interface/adaptivity.hh index ccf61645ad70622d05b8bc19f7dc8ff6fca0d9d2..33de7b1bb2d550bf249cca4e6d395cd5f9b1d809 100644 --- a/dune/dorie/interface/adaptivity.hh +++ b/dune/dorie/interface/adaptivity.hh @@ -58,6 +58,9 @@ private: const float adaptivityThreshold; //!< Global error threshold below which no refinement is applied const int verbose; //!< output verbosity of this object + const int intorderadd = 2; + const int quadrature_factor = 2; + public: /// Initialize members from config file parameters. @@ -174,7 +177,11 @@ public: t_mark = timer3.elapsed(); timer3.reset(); - Dune::PDELab::adapt_grid(grid, gfs, uold, unew, 2*order); + // calculate integration order like in DG operator + const int intorder = intorderadd + quadrature_factor * order; + + // adapt the grid + Dune::PDELab::adapt_grid(grid, gfs, uold, unew, intorder); t_adapt = timer3.elapsed(); timer3.reset();