From dc6efef44f66ceba2c465c6c64daf450d1dcca6e Mon Sep 17 00:00:00 2001
From: Elias Pipping <elias.pipping@fu-berlin.de>
Date: Wed, 29 Aug 2012 09:50:28 +0200
Subject: [PATCH] Comments
---
src/one-body-sample.parset | 2 +-
src/one-body-sample.py | 1 +
2 files changed, 2 insertions(+), 1 deletion(-)
diff --git a/src/one-body-sample.parset b/src/one-body-sample.parset
index b5c63002..695ac2ca 100644
--- a/src/one-body-sample.parset
+++ b/src/one-body-sample.parset
@@ -84,7 +84,7 @@ model = Dieterich # Ruina
a = 0.0075
b = 0.015
# If L is too large, reaching a steady state takes too long
-L = 1e-8 # 1e-3 for the Dirichlet case, 1e-5 - 1e-7 for the Neumann case
+L = 1e-8
# The second-order a and b terms in [..] are very important in
# determining whether sliding is stable or unstable and (in appli-
diff --git a/src/one-body-sample.py b/src/one-body-sample.py
index ed940280..8f4759d4 100644
--- a/src/one-body-sample.py
+++ b/src/one-body-sample.py
@@ -11,6 +11,7 @@ class neumannCondition:
else:
return trd * (x - 2.0/3) + (fst + snd) * 1.0/3
+## FIXME: If we change time discretisation, we currently change the problem
class dirichletCondition:
def __call__(self, x):
return 3e-3 / 16
--
GitLab