diff --git a/src/one-body-sample.cc b/src/one-body-sample.cc
index a3dccca8b2cdb75b748f91c8bba35ff52f20bc15..562c4cb7a96d3ac5dbe85616c88412fe9455fb9a 100644
--- a/src/one-body-sample.cc
+++ b/src/one-body-sample.cc
@@ -151,16 +151,15 @@ int main() {
     size_t const solver_maxIterations = 100;
     double const solver_tolerance = 1e-4;
 
+    // Set up grid
     typedef Dune::YaspGrid<dim> GridType;
-
     Dune::FieldVector<double, dim> const end_points(
         1); // nth dimension (zero-indexed) goes from 0 to end_points[n]
-    Dune::FieldVector<int, dim> const elements(
-        2); // number of elements in each direction
-    Dune::FieldVector<bool, dim> const periodic(false);
-
-    GridType grid(end_points, elements, periodic, 0);
-
+    GridType grid(
+        end_points,
+        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(refinements);
 
     typedef P1NodalBasis<GridType::LeafGridView, double> P1Basis;