diff --git a/src/one-body-sample.cc b/src/one-body-sample.cc
index 3501cda8b9ccb140a40ac611ec2b7466a998c34d..6597700be83dc29ad328014b542ebafbbed250f6 100644
--- a/src/one-body-sample.cc
+++ b/src/one-body-sample.cc
@@ -185,16 +185,6 @@ int main() {
     u1 = 0;
     VectorType u2 = u1;
 
-    // TODO: populate on S_F
-    std::vector<double> coefficientOfFriction;
-    coefficientOfFriction.resize(grid.size(grid.maxLevel(), dim));
-    std::fill(coefficientOfFriction.begin(), coefficientOfFriction.end(), 0.0);
-
-    // TODO: populate on S_F
-    std::vector<double> normalStress;
-    normalStress.resize(grid.size(grid.maxLevel(), dim));
-    std::fill(normalStress.begin(), normalStress.end(), 0.0);
-
     VectorType f(grid.size(grid.maxLevel(), dim));
     assemble_neumann<GridType, SmallVector, P1Basis>(grid, p1Basis,
                                                      neumannNodes, f);
@@ -215,6 +205,17 @@ int main() {
           frictionalBoundaryAssembler, nodalIntegrals,
           true); // resize the output vector and zero all of its entries
 
+      // TODO: populate on S_F
+      std::vector<double> normalStress;
+      normalStress.resize(grid.size(grid.maxLevel(), dim));
+      std::fill(normalStress.begin(), normalStress.end(), 0.0);
+
+      // TODO: populate on S_F
+      std::vector<double> coefficientOfFriction;
+      coefficientOfFriction.resize(grid.size(grid.maxLevel(), dim));
+      std::fill(coefficientOfFriction.begin(), coefficientOfFriction.end(),
+                0.0);
+
       Dune::GlobalNonlinearity<dim, Dune::LinearFunction> myGlobalNonlinearity(
           coefficientOfFriction, normalStress, nodalIntegrals);