New problem

src/one-body-sample.parset

@@ -13,7 +13,7 @@ enable = false
 iterations = 50
 
 [grid]
-refinements = 7
+refinements = 4
 
 [body]
 E = 1e6
@@ -41,7 +41,7 @@ maxiterations = 1000000
 [solver.tnnmg]
 use = true
 maxiterations = 1000000
-fixed_point_iterations = 5
+fixed_point_iterations = 50
 
 [solver.tnnmg.linear]
 maxiterations = 1000000
@@ -80,12 +80,13 @@ normalstress = 0.1 # laursen depends a lot more on this
 #    constitutive properties for earthquake prediction
 #    http://earthquake.usgs.gov/research/physics/lab/prediction.pdf
 mu = 0.5
-eta = 1
+eta = 1e-8
 model = Ruina
 
 [boundary.friction.state]
 evolve = true
-initial = 1e-4
+# log(alpha(0))
+initial = 0
 
 [boundary.friction.ruina]
 # "For rocks, typical values of A and B range from 0.005 to 0.015"
@@ -93,7 +94,7 @@ initial = 1e-4
 # Continuum mathematical modeling of slip weakening in geological systems
 a = 0.015
 b = 0.005
-L = 1e-10
+L = 1e-9
 
 # The second-order a and b terms in [..] are very important in
 # determining whether sliding is stable or unstable and (in appli-