  Residual norms for stokes_ solve.
  0 KSP Residual norm 0.00454925 
  1 KSP Residual norm 0.00448491 
  2 KSP Residual norm 0.00422213 
  3 KSP Residual norm 0.00344387 
  4 KSP Residual norm 0.00274066 
  5 KSP Residual norm 0.00224447 
  6 KSP Residual norm 0.00190167 
  7 KSP Residual norm 0.00142513 
  8 KSP Residual norm 0.00104045 
  9 KSP Residual norm 0.000919474 
 10 KSP Residual norm 0.000843643 
 11 KSP Residual norm 0.000693814 
 12 KSP Residual norm 0.00051672 
 13 KSP Residual norm 0.000437398 
 14 KSP Residual norm 0.000338089 
 15 KSP Residual norm 0.000228793 
 16 KSP Residual norm 0.000126806 
 17 KSP Residual norm 8.46533e-05 
 18 KSP Residual norm 7.04128e-05 
 19 KSP Residual norm 5.58247e-05 
 20 KSP Residual norm 4.64416e-05 
 21 KSP Residual norm 3.71667e-05 
 22 KSP Residual norm 2.76327e-05 
 23 KSP Residual norm 1.86617e-05 
 24 KSP Residual norm 1.05427e-05 
 25 KSP Residual norm 6.7899e-06 
 26 KSP Residual norm 4.08236e-06 
 27 KSP Residual norm 2.65597e-06 
 28 KSP Residual norm 1.95963e-06 
 29 KSP Residual norm 1.28355e-06 
 30 KSP Residual norm 7.71982e-07 
 31 KSP Residual norm 5.11651e-07 
 32 KSP Residual norm 2.94772e-07 
 33 KSP Residual norm 1.56689e-07 
 34 KSP Residual norm 8.65988e-08 
 35 KSP Residual norm 4.03615e-08 
 36 KSP Residual norm 2.52198e-08 
 37 KSP Residual norm 1.32618e-08 
 38 KSP Residual norm 7.19101e-09 
 39 KSP Residual norm 4.00515e-09 
 40 KSP Residual norm 1.84714e-09 
 41 KSP Residual norm 9.348e-10 
 42 KSP Residual norm 4.520e-10 
 43 KSP Residual norm 2.092e-10 
 44 KSP Residual norm 1.045e-10 
 45 KSP Residual norm 5.122e-11 
 46 KSP Residual norm 2.221e-11 
KSP Object:(stokes_) 4 MPI processes
  type: gcr
    GCR: restart = 60 
    GCR: restarts performed = 1 
  maximum iterations=10000, initial guess is zero
  tolerances:  relative=1e-08, 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
          matrix ordering: nd
          factor fill ratio given 5, needed 3.27988
            Factored matrix follows:
              Matrix Object:               1 MPI processes
                type: seqaij
                rows=867, cols=867
                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:
        Matrix Object:         1 MPI processes
          type: seqaij
          rows=867, cols=867
          total: nonzeros=21609, allocated nonzeros=23409
          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:
      Matrix Object:       4 MPI processes
        type: mpiaij
        rows=867, cols=867
        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.306745, max = 1.6871
        Chebyshev: estimated using:  [0 0.2; 0 1.1]
        KSP Object:        (stokes_mg_levels_1_est_)         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-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
                SOR: type = local_symmetric, iterations = 1, local iterations = 1, omega = 1
              linear system matrix = precond matrix:
              Matrix Object:               4 MPI processes
                type: mpiaij
                rows=2178, cols=2178
                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:
              Matrix Object:               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:
          Matrix Object:           4 MPI processes
            type: mpiaij
            rows=3267, cols=3267
            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
      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:
          Matrix Object:           4 MPI processes
            type: mpiaij
            rows=2178, cols=2178
            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:
          Matrix Object:           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:
      Matrix Object:       4 MPI processes
        type: mpiaij
        rows=3267, cols=3267
        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.291361, max = 1.60249
        Chebyshev: estimated using:  [0 0.2; 0 1.1]
        KSP Object:        (stokes_mg_levels_2_est_)         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-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
                SOR: type = local_symmetric, iterations = 1, local iterations = 1, omega = 1
              linear system matrix = precond matrix:
              Matrix Object:               4 MPI processes
                type: mpiaij
                rows=8450, cols=8450
                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:
              Matrix Object:               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:
          Matrix Object:           4 MPI processes
            type: mpiaij
            rows=12675, cols=12675
            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
      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:
          Matrix Object:           4 MPI processes
            type: mpiaij
            rows=8450, cols=8450
            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:
          Matrix Object:           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:
      Matrix Object:       4 MPI processes
        type: mpiaij
        rows=12675, cols=12675
        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.356873, max = 1.9628
        Chebyshev: estimated using:  [0 0.2; 0 1.1]
        KSP Object:        (stokes_mg_levels_3_est_)         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-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
                SOR: type = local_symmetric, iterations = 1, local iterations = 1, omega = 1
              linear system matrix = precond matrix:
              Matrix Object:               4 MPI processes
                type: mpiaij
                rows=33282, cols=33282
                total: nonzeros=587256, allocated nonzeros=587256
                total number of mallocs used during MatSetValues calls =0
                  using I-node (on process 0) routines: found 4354 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:
              Matrix Object:               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:
          Matrix Object:           4 MPI processes
            type: mpiaij
            rows=49923, cols=49923, bs=3
            total: nonzeros=1325301, allocated nonzeros=1334025
            total number of mallocs used during MatSetValues calls =0
          Matrix Object:           4 MPI processes
            type: mpiaij
            rows=49923, cols=49923, bs=3
            total: nonzeros=1325301, allocated nonzeros=1334025
            total number of mallocs used during MatSetValues calls =0
      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:
          Matrix Object:           4 MPI processes
            type: mpiaij
            rows=33282, cols=33282
            total: nonzeros=587256, allocated nonzeros=587256
            total number of mallocs used during MatSetValues calls =0
              using I-node (on process 0) routines: found 4354 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:
          Matrix Object:           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:
      Matrix Object:       4 MPI processes
        type: mpiaij
        rows=49923, cols=49923, bs=3
        total: nonzeros=1325301, allocated nonzeros=1334025
        total number of mallocs used during MatSetValues calls =0
      Matrix Object:       4 MPI processes
        type: mpiaij
        rows=49923, cols=49923, bs=3
        total: nonzeros=1325301, 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:
  Matrix Object:   4 MPI processes
    type: mpiaij
    rows=49923, cols=49923, bs=3
    total: nonzeros=1325301, allocated nonzeros=1334025
    total number of mallocs used during MatSetValues calls =0
  Matrix Object:   4 MPI processes
    type: mpiaij
    rows=49923, cols=49923, bs=3
    total: nonzeros=1325301, allocated nonzeros=1334025
    total number of mallocs used during MatSetValues calls =0
