diff --git a/src/one-body-sample.parset b/src/one-body-sample.parset index b5c63002ea3394d5fa3b3b8a074f5fdc5dc6885e..695ac2ca305d70ff9fd1fc518db8f56c443d2595 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 ed9402805e755c99962fc7341a410b0f2285bc0f..8f4759d49f5b927db3b535e0058728724c4b0b4d 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