  Residual norms for stokes_ solve.
  0 KSP Residual norm 0.00454925 
  1 KSP Residual norm 0.00445265 
  2 KSP Residual norm 0.00390246 
  3 KSP Residual norm 0.00285626 
  4 KSP Residual norm 0.00217438 
  5 KSP Residual norm 0.00137538 
  6 KSP Residual norm 0.00100592 
  7 KSP Residual norm 0.000805012 
  8 KSP Residual norm 0.000703663 
  9 KSP Residual norm 0.000670825 
 10 KSP Residual norm 0.000619861 
 11 KSP Residual norm 0.000548141 
 12 KSP Residual norm 0.000467674 
 13 KSP Residual norm 0.00045039 
 14 KSP Residual norm 0.000428086 
 15 KSP Residual norm 0.000331618 
 16 KSP Residual norm 0.00020597 
 17 KSP Residual norm 0.000115577 
 18 KSP Residual norm 8.8882e-05 
 19 KSP Residual norm 6.37702e-05 
 20 KSP Residual norm 3.59527e-05 
KSP Object: (stokes_) 4 MPI processes
  type: gcr
    restart = 60 
    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
    type is MULTIPLICATIVE, levels=4 cycles=v
      Cycles per PCApply=1
      Using Galerkin computed coarse grid matrices for pmat
  Coarse grid solver -- level -------------------------------
    KSP Object: (stokes_mg_coarse_) 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_coarse_) 4 MPI processes
      type: redundant
        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
            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
                    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 nonscalable MatPtAP() implementation
          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
        eigenvalue targets used: min 0.305926, max 1.68259
        eigenvalues estimated via gmres: min 0.189158, max 1.52963
        eigenvalues estimated using gmres with transform: [0. 0.2; 0. 1.1]
        KSP Object: (stokes_mg_levels_1_esteig_) 4 MPI processes
          type: gmres
            restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement
            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
        estimating eigenvalues using noisy right hand side
      maximum iterations=1, nonzero initial guess
      tolerances:  relative=1e-05, absolute=1e-50, divergence=10000.
      left preconditioning
      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
          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
          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 nonscalable MatPtAP() implementation
          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
        eigenvalue targets used: min 0.291447, max 1.60296
        eigenvalues estimated via gmres: min 0.394933, max 1.45724
        eigenvalues estimated using gmres with transform: [0. 0.2; 0. 1.1]
        KSP Object: (stokes_mg_levels_2_esteig_) 4 MPI processes
          type: gmres
            restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement
            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
        estimating eigenvalues using noisy right hand side
      maximum iterations=1, nonzero initial guess
      tolerances:  relative=1e-05, absolute=1e-50, divergence=10000.
      left preconditioning
      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
          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
          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 nonscalable MatPtAP() implementation
          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
        eigenvalue targets used: min 0.301145, max 1.6563
        eigenvalues estimated via gmres: min 0.485531, max 1.50573
        eigenvalues estimated using gmres with transform: [0. 0.2; 0. 1.1]
        KSP Object: (stokes_mg_levels_3_esteig_) 4 MPI processes
          type: gmres
            restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement
            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
        estimating eigenvalues using noisy right hand side
      maximum iterations=1, nonzero initial guess
      tolerances:  relative=1e-05, absolute=1e-50, divergence=10000.
      left preconditioning
      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
          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
          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 = precond 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
  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
