diff --git a/src/one-body-sample.cc b/src/one-body-sample.cc index 6eaaaf133ca8073e7ade79e05ccb726d2c7c7d4f..85d849f9bc7a92be38c88670e7edc36b83672a4b 100644 --- a/src/one-body-sample.cc +++ b/src/one-body-sample.cc @@ -152,6 +152,8 @@ int main(int argc, char *argv[]) { parset.get<size_t>("solver.maxiterations"); auto const solver_tolerance = parset.get<double>("solver.tolerance"); + auto const levels = parset.get<int>("grid.refinements") + 1; + // {{{ Set up grid typedef Dune::YaspGrid<dim> GridType; Dune::FieldVector<double, dim> const end_points( @@ -161,7 +163,7 @@ int main(int argc, char *argv[]) { Dune::FieldVector<int, dim>(2), // number of elements in each direction Dune::FieldVector<bool, dim>(false), // non-periodic in each direction 0); // zero overlap (whatever that is) - grid.globalRefine(parset.get<int>("grid.refinements")); + grid.globalRefine(levels - 1); typedef GridType::LeafGridView GridView; GridView const leafView = grid.leafView();