From 66848f4d6233e1d14d7882ab1a3e07cc880e687b Mon Sep 17 00:00:00 2001 From: Lukas Riedel Date: Mon, 20 Feb 2017 15:28:15 +0100 Subject: [PATCH] adapt_grid integration order is now the same as in DG operator --- dune/dorie/interface/adaptivity.hh | 9 ++++++++- 1 file changed, 8 insertions(+), 1 deletion(-) diff --git a/dune/dorie/interface/adaptivity.hh b/dune/dorie/interface/adaptivity.hh index ccf61645..33de7b1b 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(); -- GitLab