Soft body; new tolerance

src/one-body-sample.parset

@@ -17,7 +17,7 @@ enable_timer = false
 refinements = 4
 
 [body]
-E = 1e8
+E = 2e7
 # "Most steels and rigid polymers when used within their design limits
 # (before yield) exhibit values of about 0.3"
 # -- http://en.wikipedia.org/wiki/Poisson's_ratio
@@ -25,12 +25,12 @@ nu = 0.3 # Needs to be strictly below 0.5; the closer we get, the
          # tougher the problem. But we do have convergence!
 
 [solver]
-tolerance = 1e-10
+tolerance = 1e-11
 
 [solver.tnnmg]
 maxiterations = 1000000
-fixed_point_iterations = 1000
-fixed_point_tolerance = 1e-12
+fixed_point_iterations = 1000000
+fixed_point_tolerance = 1e-11
 
 [solver.tnnmg.linear]
 maxiterations = 1000000
@@ -60,7 +60,7 @@ acceptFactor = 1.0
 requiredResidual = 1e-12
 
 [boundary.friction]
-normalstress = 5000000
+normalstress = 10000000
 # "Most dry materials in combination have friction coefficient values
 #  between 0.3 and 0.6"
 # -- http://en.wikipedia.org/wiki/Friction#Coefficient_of_friction
@@ -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
+L = 1e-3
 
 # 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

@@ -14,7 +14,7 @@ class neumannCondition:
 ## FIXME: If we change time discretisation, we currently change the problem
 class dirichletCondition:
     def __call__(self, x):
-        return 3e-3 / 16
+        return 2e-3
         # return 0
         fst = 3e-4
         snd = 5e-4