Commit 3e456b86 authored by Lukas Riedel's avatar Lukas Riedel

adapt_grid() intorder now depends on FEorder #13

parent f147acad
......@@ -20,6 +20,7 @@ class AdaptivityController
const double alpha; //!< Refinement factor
const double beta; //!< Coarsening factor
const float adaptivityThreshold; //!< Global error threshold below which no refinement is applied
const int intorder; //!< Polynomial order of the local basis functions
int verbose; //!< Verbosity level
const int dim; //!< Spatial dimensions
......@@ -39,6 +40,7 @@ class AdaptivityController
alpha(inifile_.get<float>("adaptivity.refinementFraction")),
beta(inifile_.get<float>("adaptivity.coarseningFraction")),
adaptivityThreshold(inifile_.get<float>("adaptivity.threshold")),
intorder(inifile_.get<int>("grid.FEorder")),
verbose(inifile_.get<int>("output.verbose")),
dim(dim_), param(param_), inifile(inifile_), gv(gv_), fem(fem_)
{
......@@ -134,7 +136,7 @@ class AdaptivityController
t_mark = timer3.elapsed();
timer3.reset();
Dune::PDELab::adapt_grid( *grid, gfs, u, 2 );
Dune::PDELab::adapt_grid( *grid, gfs, u, intorder*2 );
t_adapt = timer3.elapsed();
timer3.reset();
......
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