From 82edb3bab07d4f1cc47d26dcbe4c9faac6413aa4 Mon Sep 17 00:00:00 2001
From: Elias Pipping <elias.pipping@fu-berlin.de>
Date: Wed, 14 Mar 2012 14:38:45 +0100
Subject: [PATCH] Add a means to compute the steady state
---
src/one-body-sample.cc | 4 ++++
1 file changed, 4 insertions(+)
diff --git a/src/one-body-sample.cc b/src/one-body-sample.cc
index a2b5b42b..2067be0e 100644
--- a/src/one-body-sample.cc
+++ b/src/one-body-sample.cc
@@ -269,6 +269,10 @@ int main(int argc, char *argv[]) {
double const L = parset.get<double>("boundary.friction.ruina.L");
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]);
}
}
--
GitLab