diff --git a/src/one-body-sample.cc b/src/one-body-sample.cc
index ea70eb427c96e23b3d565a95e6a62f66f813d240..446866e9c1491aab8c11725a2c5c5b1b37cfc4e9 100644
--- a/src/one-body-sample.cc
+++ b/src/one-body-sample.cc
@@ -384,11 +384,11 @@ int main(int argc, char *argv[]) {
       initialDisplacementProblemSolver.solve();
     }
     VectorType v_initial(finestSize);
+    v_initial = 0.0;
     {
       // Prescribe a homogeneous velocity field in the x-direction
       // This way, the initial condition for v and the Dirichlet
       // condition match up at t=0
-      v_initial = 0.0;
       double v_initial_const;
       velocityDirichletFunction.evaluate(0.0, v_initial_const);
       for (size_t i = 0; i < v_initial.size(); ++i)