From cefcac789bd65c3b9ff7e92d9b3b00306153b940 Mon Sep 17 00:00:00 2001
From: Elias Pipping <elias.pipping@fu-berlin.de>
Date: Thu, 22 Nov 2012 16:32:38 +0100
Subject: [PATCH] Pass theta rather than alpha initially
---
src/one-body-sample.cc | 3 ++-
src/one-body-sample.parset | 2 +-
2 files changed, 3 insertions(+), 2 deletions(-)
diff --git a/src/one-body-sample.cc b/src/one-body-sample.cc
index 624e94fc..b5c4c0cd 100644
--- a/src/one-body-sample.cc
+++ b/src/one-body-sample.cc
@@ -316,7 +316,8 @@ int main(int argc, char *argv[]) {
udd_initial = 0.0;
SingletonVectorType alpha_initial(finestSize);
- alpha_initial = parset.get<double>("boundary.friction.initial_log_state");
+ alpha_initial =
+ std::log(parset.get<double>("boundary.friction.initial_state"));
// }}}
// Set up TNNMG solver
diff --git a/src/one-body-sample.parset b/src/one-body-sample.parset
index b5c30ef1..94d982ac 100644
--- a/src/one-body-sample.parset
+++ b/src/one-body-sample.parset
@@ -53,5 +53,5 @@ a = 0.010
b = 0.015
V0 = 1
L = 1e-5
-initial_log_state = -10
+initial_state = 4.54e-05 # exp(-10)
state_model = Dieterich
--
GitLab