lid velocity = 0.0625, prandtl # = 1., grashof # = 1.
  0 SNES Function norm 2.391552133017e-01
    0 KSP Residual norm 2.897657436847e-01
    1 KSP Residual norm 1.707173075176e-02
    2 KSP Residual norm 6.405218762144e-03
    3 KSP Residual norm 1.314118438538e-03
    4 KSP Residual norm 2.530537469982e-04
    5 KSP Residual norm 8.384858565615e-05
    6 KSP Residual norm 1.667581529554e-05
    7 KSP Residual norm 3.047472053480e-06
    8 KSP Residual norm 3.176319008581e-07
  1 SNES Function norm 6.833625209978e-05
    0 KSP Residual norm 2.171349867503e-05
    1 KSP Residual norm 1.876727209396e-06
    2 KSP Residual norm 5.049137701209e-07
    3 KSP Residual norm 1.270156346290e-07
    4 KSP Residual norm 2.524572539168e-08
    5 KSP Residual norm 3.953358266575e-09
    6 KSP Residual norm 4.707984640669e-10
    7 KSP Residual norm 9.073835182022e-11
  2 SNES Function norm 3.064910055804e-10
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=15
  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: sor
      type = local_symmetric, iterations = 1, local iterations = 1, omega = 1.
    linear system matrix = precond matrix:
    Mat Object: 1 MPI process
      type: seqsell
      rows=64, cols=64, bs=4
      total: nonzeros=1152, allocated nonzeros=1152
      total number of mallocs used during MatSetValues calls=0
Number of SNES iterations = 2