diff --git a/src/one-body-sample.cc b/src/one-body-sample.cc
index ae780ea130a3db955b9a0629ef9dbcbd9b4a7772..97955995d5e6407e49eebb9522b1f6e92fead977 100644
--- a/src/one-body-sample.cc
+++ b/src/one-body-sample.cc
@@ -178,6 +178,7 @@ int main(int argc, char *argv[]) {
         .assemble(localStiffness, stiffnessMatrix);
     EnergyNorm<OperatorType, VectorType> energyNorm(stiffnessMatrix);
 
+    // Set up the boundary
     Dune::BitSetVector<dim> ignoreNodes(grid.size(grid.maxLevel(), dim), false);
     Dune::BitSetVector<1> neumannNodes(grid.size(grid.maxLevel(), dim), false);
     Dune::BitSetVector<1> frictionalNodes(grid.size(grid.maxLevel(), dim),