Commit 66848f4d authored by Lukas Riedel's avatar Lukas Riedel

adapt_grid integration order is now the same as in DG operator

parent d2b43918
......@@ -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;
/// Initialize members from config file parameters.
......@@ -174,7 +177,11 @@ public:
t_mark = timer3.elapsed();
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();
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment