From 5b66f4a5252e766ddb55d252acb6bdc85fec59c5 Mon Sep 17 00:00:00 2001
From: Elias Pipping <elias.pipping@fu-berlin.de>
Date: Tue, 15 Nov 2011 20:56:31 +0100
Subject: [PATCH] Get levels from refinements
---
src/one-body-sample.cc | 4 +++-
1 file changed, 3 insertions(+), 1 deletion(-)
diff --git a/src/one-body-sample.cc b/src/one-body-sample.cc
index 6eaaaf13..85d849f9 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();
--
GitLab