lid velocity = 0.0625, prandtl # = 1., grashof # = 1.
  0 SNES Function norm 2.391552133017e-01
    0 KSP Residual norm 2.515333131411e-01
    1 KSP Residual norm 1.396390175249e-02
    2 KSP Residual norm 3.903932692033e-03
    3 KSP Residual norm 7.056288033751e-04
    4 KSP Residual norm 4.669631283998e-05
    5 KSP Residual norm 5.322915676287e-06
    6 KSP Residual norm 9.940062885412e-07
  1 SNES Function norm 6.839858507068e-05
    0 KSP Residual norm 2.243458836024e-05
    1 KSP Residual norm 1.469665797093e-06
    2 KSP Residual norm 2.107807911905e-07
    3 KSP Residual norm 2.945467263227e-08
    4 KSP Residual norm 5.695843492021e-09
    5 KSP Residual norm 4.882870253526e-10
    6 KSP Residual norm 4.483853760249e-11
  2 SNES Function norm 8.557160032460e-11
SNES Object: 1 MPI process
  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=12
  total number of function evaluations=3
  norm schedule ALWAYS
  Jacobian is built using colored finite differences on a DMDA
  SNESLineSearch Object: 1 MPI process
    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: 1 MPI process
    type: gmres
      restart=30, using classical (unmodified) Gram-Schmidt orthogonalization with no iterative refinement
      happy breakdown tolerance=1e-30
    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: 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: 1 MPI process
            type: seqbaij
            rows=64, cols=64, bs=4
            package used to perform factorization: petsc
            total: nonzeros=1024, allocated nonzeros=1024
    linear system matrix, which is also used to construct the preconditioner:
    Mat Object: 1 MPI process
      type: seqbaij
      rows=64, cols=64, bs=4
      total: nonzeros=1024, allocated nonzeros=1024
      total number of mallocs used during MatSetValues calls=0
Number of SNES iterations = 2