diff --git a/src/one-body-sample.cc b/src/one-body-sample.cc
index 3dda82d36f3c07ef7920b91f4a1e74e2173d9067..cedecca703ddc349774834be8cb9d8aff846ce51 100644
--- a/src/one-body-sample.cc
+++ b/src/one-body-sample.cc
@@ -343,6 +343,7 @@ int main(int argc, char *argv[]) {
       _ell += gravityFunctional;
     };
     VectorType ell(finestSize);
+    createRHS(0.0, ell);
 
     // {{{ Initial conditions
     SingletonVectorType alpha_initial(finestSize);
@@ -374,7 +375,6 @@ int main(int argc, char *argv[]) {
         myGlobalNonlinearity->updateState(alpha_initial);
         // NOTE: We assume differentiability of Psi at v0 here!
         myGlobalNonlinearity->addGradient(v_initial, problem_rhs_initial);
-        createRHS(0.0, ell);
         problem_rhs_initial -= ell;
         problem_rhs_initial *= -1.0;
       }