  Residual norms for stokes_ solve.
  0 KSP Residual norm 0.00454925 
  1 KSP Residual norm 0.00448497 
  2 KSP Residual norm 0.00422229 
  3 KSP Residual norm 0.00344459 
  4 KSP Residual norm 0.00274557 
  5 KSP Residual norm 0.00224997 
  6 KSP Residual norm 0.00190515 
  7 KSP Residual norm 0.00142908 
  8 KSP Residual norm 0.00104305 
  9 KSP Residual norm 0.000922675 
 10 KSP Residual norm 0.000846683 
 11 KSP Residual norm 0.000694187 
 12 KSP Residual norm 0.000517029 
 13 KSP Residual norm 0.000437923 
 14 KSP Residual norm 0.000340909 
 15 KSP Residual norm 0.000231052 
 16 KSP Residual norm 0.000127722 
 17 KSP Residual norm 8.49558e-05 
 18 KSP Residual norm 6.9969e-05 
 19 KSP Residual norm 5.5701e-05 
 20 KSP Residual norm 4.66091e-05 
 21 KSP Residual norm 3.71995e-05 
KSP Object:(stokes_) 4 MPI processes
  type: gcr
    GCR: restart = 60 
    GCR: restarts performed = 1 
  maximum iterations=10000, initial guess is zero
  tolerances:  relative=0.01, absolute=1e-50, divergence=10000.
  right preconditioning
  using UNPRECONDITIONED norm type for convergence test
PC Object:(stokes_) 4 MPI processes
  type: mg
    MG: type is MULTIPLICATIVE, levels=4 cycles=v
      Cycles per PCApply=1
      Using Galerkin computed coarse grid matrices
  Coarse grid solver -- level -------------------------------
    KSP Object:    (stokes_mg_coarse_)     4 MPI processes
      type: preonly
      maximum iterations=1, initial guess is zero
      tolerances:  relative=1e-05, absolute=1e-50, divergence=10000.
      left preconditioning
      using NONE norm type for convergence test
    PC Object:    (stokes_mg_coarse_)     4 MPI processes
      type: redundant
        Redundant preconditioner: First (color=0) of 4 PCs follows
        KSP Object:        (stokes_mg_coarse_redundant_)         1 MPI processes
          type: preonly
          maximum iterations=10000, initial guess is zero
          tolerances:  relative=1e-05, absolute=1e-50, divergence=10000.
          left preconditioning
          using NONE norm type for convergence test
        PC Object:        (stokes_mg_coarse_redundant_)         1 MPI processes
          type: lu
            LU: out-of-place factorization
            tolerance for zero pivot 2.22045e-14
            using diagonal shift on blocks to prevent zero pivot [INBLOCKS]
            matrix ordering: nd
            factor fill ratio given 5., needed 3.27988
              Factored matrix follows:
                Mat Object:                 1 MPI processes
                  type: seqaij
                  rows=867, cols=867, bs=3
                  package used to perform factorization: petsc
                  total: nonzeros=70875, allocated nonzeros=70875
                  total number of mallocs used during MatSetValues calls =0
                    using I-node routines: found 287 nodes, limit used is 5
          linear system matrix = precond matrix:
          Mat Object:           1 MPI processes
            type: seqaij
            rows=867, cols=867, bs=3
            total: nonzeros=21609, allocated nonzeros=21609
            total number of mallocs used during MatSetValues calls =0
              using I-node routines: found 289 nodes, limit used is 5
      linear system matrix = precond matrix:
      Mat Object:       4 MPI processes
        type: mpiaij
        rows=867, cols=867, bs=3
        total: nonzeros=21609, allocated nonzeros=21609
        total number of mallocs used during MatSetValues calls =0
          using I-node (on process 0) routines: found 81 nodes, limit used is 5
  Down solver (pre-smoother) on level 1 -------------------------------
    KSP Object:    (stokes_mg_levels_1_)     4 MPI processes
      type: chebyshev
        Chebyshev: eigenvalue estimates:  min = 0.306009, max = 1.68305
        Chebyshev: eigenvalues estimated using gmres with translations  [0. 0.2; 0. 1.1]
        KSP Object:        (stokes_mg_levels_1_esteig_)         4 MPI processes
          type: gmres
            GMRES: restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement
            GMRES: happy breakdown tolerance 1e-30
          maximum iterations=10, initial guess is zero
          tolerances:  relative=1e-12, absolute=1e-50, divergence=10000.
          left preconditioning
          using PRECONDITIONED norm type for convergence test
      maximum iterations=1
      tolerances:  relative=1e-05, absolute=1e-50, divergence=10000.
      left preconditioning
      using nonzero initial guess
      using NONE norm type for convergence test
    PC Object:    (stokes_mg_levels_1_)     4 MPI processes
      type: fieldsplit
        FieldSplit with MULTIPLICATIVE composition: total splits = 2, blocksize = 3
        Solver info for each split is in the following KSP objects:
        Split number 0 Fields  0, 1
        KSP Object:        (stokes_mg_levels_1_fieldsplit_0_)         4 MPI processes
          type: preonly
          maximum iterations=10000, initial guess is zero
          tolerances:  relative=1e-05, absolute=1e-50, divergence=10000.
          left preconditioning
          using NONE norm type for convergence test
        PC Object:        (stokes_mg_levels_1_fieldsplit_0_)         4 MPI processes
          type: sor
            SOR: type = local_symmetric, iterations = 1, local iterations = 1, omega = 1.
          linear system matrix = precond matrix:
          Mat Object:          (stokes_mg_levels_1_fieldsplit_0_)           4 MPI processes
            type: mpiaij
            rows=2178, cols=2178, bs=2
            total: nonzeros=37636, allocated nonzeros=37636
            total number of mallocs used during MatSetValues calls =0
              using I-node (on process 0) routines: found 289 nodes, limit used is 5
        Split number 1 Fields  2
        KSP Object:        (stokes_mg_levels_1_fieldsplit_1_)         4 MPI processes
          type: preonly
          maximum iterations=10000, initial guess is zero
          tolerances:  relative=1e-05, absolute=1e-50, divergence=10000.
          left preconditioning
          using NONE norm type for convergence test
        PC Object:        (stokes_mg_levels_1_fieldsplit_1_)         4 MPI processes
          type: sor
            SOR: type = local_symmetric, iterations = 1, local iterations = 1, omega = 1.
          linear system matrix = precond matrix:
          Mat Object:          (stokes_mg_levels_1_fieldsplit_1_)           4 MPI processes
            type: mpiaij
            rows=1089, cols=1089
            total: nonzeros=9409, allocated nonzeros=9409
            total number of mallocs used during MatSetValues calls =0
              not using I-node (on process 0) routines
      linear system matrix = precond matrix:
      Mat Object:       4 MPI processes
        type: mpiaij
        rows=3267, cols=3267, bs=3
        total: nonzeros=84681, allocated nonzeros=84681
        total number of mallocs used during MatSetValues calls =0
          using I-node (on process 0) routines: found 289 nodes, limit used is 5
  Up solver (post-smoother) same as down solver (pre-smoother)
  Down solver (pre-smoother) on level 2 -------------------------------
    KSP Object:    (stokes_mg_levels_2_)     4 MPI processes
      type: chebyshev
        Chebyshev: eigenvalue estimates:  min = 0.291297, max = 1.60213
        Chebyshev: eigenvalues estimated using gmres with translations  [0. 0.2; 0. 1.1]
        KSP Object:        (stokes_mg_levels_2_esteig_)         4 MPI processes
          type: gmres
            GMRES: restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement
            GMRES: happy breakdown tolerance 1e-30
          maximum iterations=10, initial guess is zero
          tolerances:  relative=1e-12, absolute=1e-50, divergence=10000.
          left preconditioning
          using PRECONDITIONED norm type for convergence test
      maximum iterations=1
      tolerances:  relative=1e-05, absolute=1e-50, divergence=10000.
      left preconditioning
      using nonzero initial guess
      using NONE norm type for convergence test
    PC Object:    (stokes_mg_levels_2_)     4 MPI processes
      type: fieldsplit
        FieldSplit with MULTIPLICATIVE composition: total splits = 2, blocksize = 3
        Solver info for each split is in the following KSP objects:
        Split number 0 Fields  0, 1
        KSP Object:        (stokes_mg_levels_2_fieldsplit_0_)         4 MPI processes
          type: preonly
          maximum iterations=10000, initial guess is zero
          tolerances:  relative=1e-05, absolute=1e-50, divergence=10000.
          left preconditioning
          using NONE norm type for convergence test
        PC Object:        (stokes_mg_levels_2_fieldsplit_0_)         4 MPI processes
          type: sor
            SOR: type = local_symmetric, iterations = 1, local iterations = 1, omega = 1.
          linear system matrix = precond matrix:
          Mat Object:          (stokes_mg_levels_2_fieldsplit_0_)           4 MPI processes
            type: mpiaij
            rows=8450, cols=8450, bs=2
            total: nonzeros=148996, allocated nonzeros=148996
            total number of mallocs used during MatSetValues calls =0
              using I-node (on process 0) routines: found 1089 nodes, limit used is 5
        Split number 1 Fields  2
        KSP Object:        (stokes_mg_levels_2_fieldsplit_1_)         4 MPI processes
          type: preonly
          maximum iterations=10000, initial guess is zero
          tolerances:  relative=1e-05, absolute=1e-50, divergence=10000.
          left preconditioning
          using NONE norm type for convergence test
        PC Object:        (stokes_mg_levels_2_fieldsplit_1_)         4 MPI processes
          type: sor
            SOR: type = local_symmetric, iterations = 1, local iterations = 1, omega = 1.
          linear system matrix = precond matrix:
          Mat Object:          (stokes_mg_levels_2_fieldsplit_1_)           4 MPI processes
            type: mpiaij
            rows=4225, cols=4225
            total: nonzeros=37249, allocated nonzeros=37249
            total number of mallocs used during MatSetValues calls =0
              not using I-node (on process 0) routines
      linear system matrix = precond matrix:
      Mat Object:       4 MPI processes
        type: mpiaij
        rows=12675, cols=12675, bs=3
        total: nonzeros=335241, allocated nonzeros=335241
        total number of mallocs used during MatSetValues calls =0
          using I-node (on process 0) routines: found 1089 nodes, limit used is 5
  Up solver (post-smoother) same as down solver (pre-smoother)
  Down solver (pre-smoother) on level 3 -------------------------------
    KSP Object:    (stokes_mg_levels_3_)     4 MPI processes
      type: chebyshev
        Chebyshev: eigenvalue estimates:  min = 0.356152, max = 1.95884
        Chebyshev: eigenvalues estimated using gmres with translations  [0. 0.2; 0. 1.1]
        KSP Object:        (stokes_mg_levels_3_esteig_)         4 MPI processes
          type: gmres
            GMRES: restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement
            GMRES: happy breakdown tolerance 1e-30
          maximum iterations=10, initial guess is zero
          tolerances:  relative=1e-12, absolute=1e-50, divergence=10000.
          left preconditioning
          using PRECONDITIONED norm type for convergence test
      maximum iterations=1
      tolerances:  relative=1e-05, absolute=1e-50, divergence=10000.
      left preconditioning
      using nonzero initial guess
      using NONE norm type for convergence test
    PC Object:    (stokes_mg_levels_3_)     4 MPI processes
      type: fieldsplit
        FieldSplit with MULTIPLICATIVE composition: total splits = 2, blocksize = 3
        Solver info for each split is in the following KSP objects:
        Split number 0 Fields  0, 1
        KSP Object:        (stokes_mg_levels_3_fieldsplit_0_)         4 MPI processes
          type: preonly
          maximum iterations=10000, initial guess is zero
          tolerances:  relative=1e-05, absolute=1e-50, divergence=10000.
          left preconditioning
          using NONE norm type for convergence test
        PC Object:        (stokes_mg_levels_3_fieldsplit_0_)         4 MPI processes
          type: sor
            SOR: type = local_symmetric, iterations = 1, local iterations = 1, omega = 1.
          linear system matrix = precond matrix:
          Mat Object:          (stokes_mg_levels_3_fieldsplit_0_)           4 MPI processes
            type: mpiaij
            rows=33282, cols=33282, bs=2
            total: nonzeros=592900, allocated nonzeros=592900
            total number of mallocs used during MatSetValues calls =0
              using I-node (on process 0) routines: found 4225 nodes, limit used is 5
        Split number 1 Fields  2
        KSP Object:        (stokes_mg_levels_3_fieldsplit_1_)         4 MPI processes
          type: preonly
          maximum iterations=10000, initial guess is zero
          tolerances:  relative=1e-05, absolute=1e-50, divergence=10000.
          left preconditioning
          using NONE norm type for convergence test
        PC Object:        (stokes_mg_levels_3_fieldsplit_1_)         4 MPI processes
          type: sor
            SOR: type = local_symmetric, iterations = 1, local iterations = 1, omega = 1.
          linear system matrix = precond matrix:
          Mat Object:          (stokes_mg_levels_3_fieldsplit_1_)           4 MPI processes
            type: mpiaij
            rows=16641, cols=16641
            total: nonzeros=148225, allocated nonzeros=148225
            total number of mallocs used during MatSetValues calls =0
              not using I-node (on process 0) routines
      linear system matrix followed by preconditioner matrix:
      Mat Object:       4 MPI processes
        type: mpiaij
        rows=49923, cols=49923, bs=3
        total: nonzeros=1334025, allocated nonzeros=1334025
        total number of mallocs used during MatSetValues calls =0
      Mat Object:       4 MPI processes
        type: mpiaij
        rows=49923, cols=49923, bs=3
        total: nonzeros=1334025, allocated nonzeros=1334025
        total number of mallocs used during MatSetValues calls =0
  Up solver (post-smoother) same as down solver (pre-smoother)
  linear system matrix followed by preconditioner matrix:
  Mat Object:   4 MPI processes
    type: mpiaij
    rows=49923, cols=49923, bs=3
    total: nonzeros=1334025, allocated nonzeros=1334025
    total number of mallocs used during MatSetValues calls =0
  Mat Object:   4 MPI processes
    type: mpiaij
    rows=49923, cols=49923, bs=3
    total: nonzeros=1334025, allocated nonzeros=1334025
    total number of mallocs used during MatSetValues calls =0
