lid velocity = 0.00591716, prandtl # = 1., grashof # = 1.
    0 KSP Residual norm 0.22318
    3 KSP Residual norm 0.0397558
    6 KSP Residual norm 0.000417387
    9 KSP Residual norm 2.72144e-05
   12 KSP Residual norm 2.60861e-06
   15 KSP Residual norm 1.0556e-07
    0 KSP Residual norm 6.21497e-06
    3 KSP Residual norm 4.01807e-07
    6 KSP Residual norm 3.02758e-08
    9 KSP Residual norm 5.54048e-09
   12 KSP Residual norm 1.540e-10
   15 KSP Residual norm 2.084e-11
SNES Object: 2 MPI processes
  type: newtonls
  maximum iterations=50, maximum function evaluations=10000
  tolerances: relative=1e-08, absolute=1e-50, solution=1e-08
  total number of linear solver iterations=30
  total number of function evaluations=3
  norm schedule ALWAYS
  Jacobian is built using colored finite differences on a DMDA
  SNESLineSearch Object: 2 MPI processes
    type: bt
      interpolation: cubic
      alpha=1.000000e-04
    maxlambda=1.000000e+00, minlambda=1.000000e-12
    tolerances: relative=1.000000e-08, absolute=1.000000e-15, lambda=1.000000e-08
    maximum iterations=40
  KSP Object: 2 MPI processes
    type: bcgsl
      Ell = 3
      Delta = 0.
    maximum iterations=10000, initial guess is zero
    tolerances: relative=1e-05, absolute=1e-50, divergence=10000.
    left preconditioning
    using PRECONDITIONED norm type for convergence test
  PC Object: 2 MPI processes
    type: bjacobi
      number of blocks = 2
      Local solver information for first block is in the following KSP and PC objects on rank 0:
      Use -ksp_view ::ascii_info_detail to display information for all blocks
      KSP Object: (sub_) 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: (sub_) 1 MPI process
        type: ilu
          out-of-place factorization
          0 levels of fill
          tolerance for zero pivot 2.22045e-14
          matrix ordering: natural
          factor fill ratio given 1., needed 1.
            Factored matrix:
              Mat Object: (sub_) 1 MPI process
                type: seqaij
                rows=364, cols=364, bs=4
                package used to perform factorization: petsc
                total: nonzeros=6640, allocated nonzeros=6640
                  using I-node routines: found 91 nodes, limit used is 5
        linear system matrix, which is also used to construct the preconditioner:
        Mat Object: (sub_) 1 MPI process
          type: seqaij
          rows=364, cols=364, bs=4
          total: nonzeros=6640, allocated nonzeros=6640
          total number of mallocs used during MatSetValues calls=0
            using I-node routines: found 91 nodes, limit used is 5
    linear system matrix, which is also used to construct the preconditioner:
    Mat Object: 2 MPI processes
      type: mpiaij
      rows=676, cols=676, bs=4
      total: nonzeros=12688, allocated nonzeros=12688
      total number of mallocs used during MatSetValues calls=0
Number of SNES iterations = 2
