# -*- mode:conf -*-
gravity         = 9.81     # [m/s^2]

[body0]
length          = 0.5      # [m]
height          = 0.5     # [m]
depth           = 0.12     # [m]
bulkModulus     = 1.5e5    # [Pa] #2190
poissonRatio    = 0.11     # [1]  #0.11
[body0.elastic]
distFromDiag    = 0.2
density         = 1300      # [kg/m^3] #750
shearViscosity  = 0     # [Pas]
bulkViscosity   = 0     # [Pas]
[body0.viscoelastic]
density         = 1300     # [kg/m^3]
shearViscosity  = 1e4     # [Pas]
bulkViscosity   = 1e4     # [Pas]

[body1]
length          = 0.5     # [m]
height          = 0.5     # [m]
depth           = 0.12     # [m]
bulkModulus     = 1.5e5    # [Pa]
poissonRatio    = 0.11     # [1]
[body1.elastic]
distFromDiag    = 0.2
density         = 1300      # [kg/m^3]
shearViscosity  = 0     # [Pas]
bulkViscosity   = 0     # [Pas]
[body1.viscoelastic]
density         = 1300     # [kg/m^3]
shearViscosity  = 1e4     # [Pas]
bulkViscosity   = 1e4     # [Pas]

[boundary.friction]
C               = 6       # [Pa]
mu0             = 0.48      # [ ]
V0              = 1e-3     # [m/s]
L               = 1e-6  # [m]
initialAlpha    = 0        # [ ]
stateModel      = AgeingLaw
frictionModel   = Truncated #Regularised
[boundary.friction.weakening]
a               = 0.054    # [ ]
b               = 0.074    # [ ]
[boundary.friction.strengthening]
a               = 0.054    # [ ]
b               = 0.074    # [ ]

[boundary.neumann]
sigmaN          = 200.0      # [Pa]

[boundary.dirichlet]
finalVelocity   = 1e-4     # [m/s]

[io]
data.write      = true
printProgress   = true
restarts.first  = 0
restarts.spacing= 1
restarts.write  = true #true
vtk.write       = true

[problem]
finalTime       = 100     # [s] #1000
bodyCount       = 2

[initialTime]
timeStep = 0
relativeTime = 0.0
relativeTau = 2e-2 # 1e-6

[timeSteps]
scheme = newmark
timeSteps = 5

[u0.solver]
maximumIterations = 100
verbosity         = full

[a0.solver]
maximumIterations = 100
verbosity         = full

[v.solver]
maximumIterations = 100
verbosity         = quiet

[v.fpi]
maximumIterations = 10000
lambda            = 0.5

[solver.tnnmg.preconditioner]
mode         = additive
patchDepth   = 1
maximumIterations = 2
verbosity         = quiet
[solver.tnnmg.preconditioner.patchsolver]
maximumIterations = 100
verbosity         = quiet
[solver.tnnmg.preconditioner.basesolver]
maximumIterations = 10000
verbosity         = quiet

[solver.tnnmg.main]
pre   = 1
multi = 5 # number of multigrid steps
post  = 0