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

Add a means to compute the steady state

parent ad66c3a5
No related branches found
No related tags found
No related merge requests found
...@@ -269,6 +269,10 @@ int main(int argc, char *argv[]) { ...@@ -269,6 +269,10 @@ int main(int argc, char *argv[]) {
double const L = parset.get<double>("boundary.friction.ruina.L"); double const L = parset.get<double>("boundary.friction.ruina.L");
double const unorm = u4_diff[i].two_norm(); double const unorm = u4_diff[i].two_norm();
// the (logarithmic) steady state corresponding to the current
// velocity
// std::cout << std::log(L/unorm * h) << std::endl;
(*s4_new)[i] = compute_state_update(h, unorm, L, s4_old[i]); (*s4_new)[i] = compute_state_update(h, unorm, L, s4_old[i]);
} }
} }
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment