geometry = examples/beam.geo
mesh = examples/beam.vol

define coefficient e
1,

define coefficient nu
0.2,

define coefficient rho
1

define coefficient penalty
1e6, 0, 0, 0, 0, 0

define coefficient coef_source
1, 0, 0


define fespace v -dim=3
define fespace vp -dim=6 -l2

define gridfunction u -fespace=v
define gridfunction p -fespace=vp

define bilinearform a -fespace=v 
elasticity e nu
sysrobin penalty

define bilinearform m -fespace=v 
sysmass rho

define preconditioner c -type=multigrid -bilinearform=a


define linearform f -fespace=v
force forcex forcey forcez

numproc evp np1 -bilinearforma=a -bilinearformm=m  -gridfunction=u -preconditioner=c
# numproc calcflux np2 -bilinearform=a -solution=u -flux=p
# numproc zzerrorestimator np3 -bilinearform=a -linearform=f  -flux=p