diff --git a/src/one-body-sample.cc b/src/one-body-sample.cc index 60d96e82394bc0aac2635cf162708e6973a6cf13..2f7787defe183a2184c44fea29fef43b43aa3a40 100644 --- a/src/one-body-sample.cc +++ b/src/one-body-sample.cc @@ -334,6 +334,10 @@ int main(int argc, char *argv[]) { VectorType ell(finestSize); // {{{ Initial conditions + SingletonVectorType alpha_initial(finestSize); + alpha_initial = + std::log(parset.get<double>("boundary.friction.initial_state")); + VectorType u_initial(finestSize); u_initial = 0.0; VectorType v_initial(finestSize); @@ -341,10 +345,6 @@ int main(int argc, char *argv[]) { // FIXME: This only happens to be correct VectorType a_initial(finestSize); a_initial = 0.0; - - SingletonVectorType alpha_initial(finestSize); - alpha_initial = - std::log(parset.get<double>("boundary.friction.initial_state")); // }}} // Set up TNNMG solver