#pc-type=fieldsplit #gasm,lu,fieldsplit,ml,ilu
ksp-type=fgmres
pc-type=fieldsplit

#####################################################
# gasm
pc-gasm-overlap=1
sub-pc-type=lu #ilu
#sub-pc-factor-levels=3

#####################################################
# ilu
pc-factor-levels=3
#####################################################
# ml
pc-mg-levels=2
pc-mg-type=kaskade #multiplicative, additive, full, kaskade

mg-coarse.pc-type=lu

mg-levels1.pc-type=gasm
mg-levels1.sub-pc-type=lu #ilu
mg-levels1.pc-gasm-overlap=0
mg-levels1.ksp-maxit=50
mg-levels1.ksp-rtol=1e-5
#mg-levels1.ksp-monitor=true

######################################################
# fieldsplit
fieldsplit-type=schur #additive, multiplicative, schur
fieldsplit-schur-fact-type=full #diag, lower, upper, full

fieldsplit-0.pc-type=bjacobi #gasm,lu
fieldsplit-0.sub-pc-type=lu
fieldsplit-0.pc-gasm-overlap=1
fieldsplit-0.ksp-type=preonly #gmres #preonly
fieldsplit-0.reuse-prec=true

# D^-1 L^-1
fieldsplit-schur-inner-solver.use-outer-solver=true
fieldsplit-schur-inner-solver.pc-type=jacobi
fieldsplit-schur-inner-solver.ksp-type=preonly

# U^-1
fieldsplit-schur-upper-solver.use-outer-solver=true
fieldsplit-schur-upper-solver.pc-type=jacobi
fieldsplit-schur-upper-solver.sub-pc-type=lu
fieldsplit-schur-upper-solver.ksp-type=preonly

# Schur complement
fieldsplit-schur-precondition=user
fieldsplit-1.pc-type=ssor
fieldsplit-1.sub-pc-type=lu
fieldsplit-1.ksp-type=fgmres
fieldsplit-1.ksp-maxit=15
fieldsplit-1.rtol=1e-3

#fieldsplit-1.pc-type=lsc  #lsc,jacobi
#fieldsplit-1.ksp-type=gmres #preonly, gmres,minres
#fieldsplit-1.lsc.pc-type=lu # ml,lu,gasm

fieldsplit-1.lsc.ksp-type=gmres
fieldsplit-1.lsc.pc-type=gamg
fieldsplit-1.lsc.ksp-rtol=1e-10
fieldsplit-1.lsc.ksp-monitor=false
fieldsplit-1.lsc.reuse-prec=true

functions.g=(y+0.4)*(2.4-y):x:y

[gmsh]
#filename=$top_srcdir/data/gmsh/feelpp/feel2d-ns.geo
filename=$top_srcdir/data/gmsh/feelpp/feel3d-ns.geo
hsize=0.2

[bdf]
time-step=0.1
time-final=1

[exporter]
geometry=static
format=ensightgold