Skip to content
Snippets Groups Projects
Commit c6fb60c0 authored by Elias Pipping's avatar Elias Pipping Committed by Elias Pipping
Browse files

[Cleanup] Move declaration/initialisation of alpha

parent 4f346fa3
No related branches found
No related tags found
No related merge requests found
......@@ -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
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment