## Outer solver on the full block system
-stokes_ksp_type gcr
-stokes_pc_type mg
-levels 3
-stokes_mg_levels_ksp_type fgmres
-stokes_mg_levels_ksp_max_it 4
-stokes_mg_levels_pc_type fieldsplit 

  ## Set three blocks for fieldsplit - u,v,w,p
  -stokes_mg_levels_pc_fieldsplit_block_size 4 
#  -stokes_mg_levels_pc_fieldsplit_type SYMMETRIC_MULTIPLICATIVE
  -stokes_mg_levels_pc_fieldsplit_type ADDITIVE
    ## group u,v,w together into field_0
    -stokes_mg_levels_pc_fieldsplit_0_fields 0,1,2
    ## place p in field 1
    -stokes_mg_levels_pc_fieldsplit_1_fields 3

      ## --------------------------------------------------
      ## configure the A_UU^{-1} block to be a iterative solve
      -stokes_mg_levels_fieldsplit_0_ksp_type fgmres
#      -stokes_mg_levels_fieldsplit_0_ksp_max_it 5
      -stokes_mg_levels_fieldsplit_0_ksp_rtol 1.0e-2
      -stokes_mg_levels_fieldsplit_0_pc_type ilu 

      ## -------------------------------------------------------------
      ## configure the A_pp^{-1} block to just be the diagonal of A_pp
      -stokes_mg_levels_fieldsplit_1_ksp_type fgmres
      -stokes_mg_levels_fieldsplit_1_ksp_max_it 5
      -stokes_mg_levels_fieldsplit_1_pc_type jacobi 
    ## end fieldsplit block configuration
  ## end fieldsplit configuration
## end solver configuration
##
-stokes_ksp_monitor_blocks
-mx 16
-model 3

