  Residual norms for stokes_ solve.
  0 KSP Residual norm 0.00337927
  1 KSP Residual norm 0.00100033
  2 KSP Residual norm 0.00198036
  3 KSP Residual norm 0.00138965
  4 KSP Residual norm 0.00117457
  5 KSP Residual norm 0.000844377
  6 KSP Residual norm 0.000464969
  7 KSP Residual norm 0.000344846
  8 KSP Residual norm 0.000173329
  9 KSP Residual norm 0.000138069
 10 KSP Residual norm 0.000181073
 11 KSP Residual norm 0.000197309
 12 KSP Residual norm 0.000122466
 13 KSP Residual norm 8.0686e-05
 14 KSP Residual norm 0.000129132
 15 KSP Residual norm 0.000192349
 16 KSP Residual norm 0.000142142
 17 KSP Residual norm 6.46078e-05
 18 KSP Residual norm 4.28934e-05
 19 KSP Residual norm 3.97016e-05
 20 KSP Residual norm 2.41842e-05
 21 KSP Residual norm 1.28061e-05
 22 KSP Residual norm 8.62207e-06
 23 KSP Residual norm 4.79149e-06
 24 KSP Residual norm 4.08285e-06
 25 KSP Residual norm 3.2439e-06
 26 KSP Residual norm 2.47285e-06
 27 KSP Residual norm 1.70736e-06
 28 KSP Residual norm 1.22027e-06
 29 KSP Residual norm 6.35357e-07
 30 KSP Residual norm 3.67964e-07
 31 KSP Residual norm 2.37451e-07
 32 KSP Residual norm 1.71305e-07
 33 KSP Residual norm 1.10396e-07
 34 KSP Residual norm 6.10948e-08
 35 KSP Residual norm 3.47249e-08
 36 KSP Residual norm 2.72656e-08
KSP Object: (stokes_) 4 MPI processes
  type: pipegcr
    max previous directions = 60
    preallocated 5 directions
    Using Notay's truncation strategy
    w unrolling = TRUE 
    restarts performed = 0 
  maximum iterations=10000, initial guess is zero
  tolerances: relative=1e-05, absolute=1e-50, divergence=10000.
  left preconditioning
  using NATURAL 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 0 -------------------------------
    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
      not checking for convergence
    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 process
          type: preonly
          maximum iterations=10000, initial guess is zero
          tolerances: relative=1e-05, absolute=1e-50, divergence=10000.
          left preconditioning
          not checking for convergence
        PC Object: (stokes_mg_coarse_redundant_) 1 MPI process
          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:
                Mat Object: (stokes_mg_coarse_redundant_) 1 MPI process
                  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, which is also used to construct the preconditioner:
          Mat Object: 1 MPI process
            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, which is also used to construct the preconditioner:
      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
        Chebyshev polynomial of first kind
        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 a noisy random number generated right-hand side
      maximum iterations=1, nonzero initial guess
      tolerances: relative=1e-05, absolute=1e-50, divergence=10000.
      left preconditioning
      not checking for convergence
    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
        not checking for convergence
      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, which is also used to construct the preconditioner:
        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
        not checking for convergence
      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, which is also used to construct the preconditioner:
        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, which is also used to construct the preconditioner:
      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
        Chebyshev polynomial of first kind
        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 a noisy random number generated right-hand side
      maximum iterations=1, nonzero initial guess
      tolerances: relative=1e-05, absolute=1e-50, divergence=10000.
      left preconditioning
      not checking for convergence
    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
        not checking for convergence
      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, which is also used to construct the preconditioner:
        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
        not checking for convergence
      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, which is also used to construct the preconditioner:
        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, which is also used to construct the preconditioner:
      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
        Chebyshev polynomial of first kind
        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 a noisy random number generated right-hand side
      maximum iterations=1, nonzero initial guess
      tolerances: relative=1e-05, absolute=1e-50, divergence=10000.
      left preconditioning
      not checking for convergence
    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
        not checking for convergence
      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, which is also used to construct the preconditioner:
        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
        not checking for convergence
      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, which is also used to construct the preconditioner:
        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, which is also used to construct the preconditioner:
      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 the matrix used to construct the preconditioner:
  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
