Comments

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-

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