  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 = 1 
    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 -------------------------------
    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
                  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 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 estimates used:  min = 0.305926, max = 1.68259
        eigenvalues estimate via gmres min 0.189158, max 1.52963
        eigenvalues estimated using gmres with translations  [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 estimates used:  min = 0.291447, max = 1.60296
        eigenvalues estimate via gmres min 0.394933, max 1.45724
        eigenvalues estimated using gmres with translations  [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 estimates used:  min = 0.301145, max = 1.6563
        eigenvalues estimate via gmres min 0.485531, max 1.50573
        eigenvalues estimated using gmres with translations  [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
