lid velocity = 0.0204082, prandtl # = 1., grashof # = 1.
  0 SNES Function norm 0.146194
    0 KSP Residual norm 0.146194
    1 KSP Residual norm 0.0781702
    2 KSP Residual norm 0.0289896
    3 KSP Residual norm 0.00874626
    4 KSP Residual norm 0.00360523
    5 KSP Residual norm 0.00187927
    6 KSP Residual norm 0.000997932
    7 KSP Residual norm 0.000325883
    8 KSP Residual norm 0.000122254
    9 KSP Residual norm 6.21287e-05
   10 KSP Residual norm 4.04759e-05
   11 KSP Residual norm 2.1448e-05
   12 KSP Residual norm 9.92327e-06
   13 KSP Residual norm 4.95927e-06
   14 KSP Residual norm 2.41652e-06
   15 KSP Residual norm 1.39592e-06
  1 SNES Function norm 2.5489e-05
    0 KSP Residual norm 2.5489e-05
    1 KSP Residual norm 1.07362e-05
    2 KSP Residual norm 3.78935e-06
    3 KSP Residual norm 1.43754e-06
    4 KSP Residual norm 6.54959e-07
    5 KSP Residual norm 5.12221e-07
    6 KSP Residual norm 3.42561e-07
    7 KSP Residual norm 1.45139e-07
    8 KSP Residual norm 7.40916e-08
    9 KSP Residual norm 3.79621e-08
   10 KSP Residual norm 1.50283e-08
   11 KSP Residual norm 5.03342e-09
   12 KSP Residual norm 2.41715e-09
   13 KSP Residual norm 1.38641e-09
   14 KSP Residual norm 7.047e-10
   15 KSP Residual norm 3.897e-10
   16 KSP Residual norm 2.021e-10
  2 SNES Function norm 2.020e-10
SNES Object: 3 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=31
  total number of function evaluations=36
  norm schedule ALWAYS
  Jacobian is applied matrix-free with differencing
  Preconditioning Jacobian is built using colored finite differences on a DMDA
  SNESLineSearch Object: 3 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: 3 MPI processes
    type: fgmres
      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.
    right preconditioning
    using UNPRECONDITIONED norm type for convergence test
  PC Object: 3 MPI processes
    type: fieldsplit
      FieldSplit with MULTIPLICATIVE composition: total splits = 4
      Solver info for each split is in the following KSP objects:
    Split number 0 Defined by IS
    KSP Object: (fieldsplit_x_velocity_) 3 MPI processes
      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: (fieldsplit_x_velocity_) 3 MPI processes
      type: bjacobi
        number of blocks = 3
        Local solver information for first block is in the following KSP and PC objects on rank 0:
        Use -fieldsplit_x_velocity_ksp_view ::ascii_info_detail to display information for all blocks
        KSP Object: (fieldsplit_x_velocity_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: (fieldsplit_x_velocity_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: (fieldsplit_x_velocity_sub_) 1 MPI process
                  type: seqaij
                  rows=21, cols=21
                  package used to perform factorization: petsc
                  total: nonzeros=85, allocated nonzeros=85
                    not using I-node routines
          linear system matrix, which is also used to construct the preconditioner:
          Mat Object: (fieldsplit_x_velocity_sub_) 1 MPI process
            type: seqaij
            rows=21, cols=21
            total: nonzeros=85, allocated nonzeros=85
            total number of mallocs used during MatSetValues calls=0
              not using I-node routines
      linear system matrix, which is also used to construct the preconditioner:
      Mat Object: (fieldsplit_x_velocity_) 3 MPI processes
        type: mpiaij
        rows=49, cols=49
        total: nonzeros=217, allocated nonzeros=217
        total number of mallocs used during MatSetValues calls=0
          not using I-node (on process 0) routines
    Split number 1 Defined by IS
    KSP Object: (fieldsplit_y_velocity_) 3 MPI processes
      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: (fieldsplit_y_velocity_) 3 MPI processes
      type: bjacobi
        number of blocks = 3
        Local solver information for first block is in the following KSP and PC objects on rank 0:
        Use -fieldsplit_y_velocity_ksp_view ::ascii_info_detail to display information for all blocks
        KSP Object: (fieldsplit_y_velocity_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: (fieldsplit_y_velocity_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: (fieldsplit_y_velocity_sub_) 1 MPI process
                  type: seqaij
                  rows=21, cols=21
                  package used to perform factorization: petsc
                  total: nonzeros=85, allocated nonzeros=85
                    not using I-node routines
          linear system matrix, which is also used to construct the preconditioner:
          Mat Object: (fieldsplit_y_velocity_sub_) 1 MPI process
            type: seqaij
            rows=21, cols=21
            total: nonzeros=85, allocated nonzeros=85
            total number of mallocs used during MatSetValues calls=0
              not using I-node routines
      linear system matrix, which is also used to construct the preconditioner:
      Mat Object: (fieldsplit_y_velocity_) 3 MPI processes
        type: mpiaij
        rows=49, cols=49
        total: nonzeros=217, allocated nonzeros=217
        total number of mallocs used during MatSetValues calls=0
          not using I-node (on process 0) routines
    Split number 2 Defined by IS
    KSP Object: (fieldsplit_Omega_) 3 MPI processes
      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: (fieldsplit_Omega_) 3 MPI processes
      type: bjacobi
        number of blocks = 3
        Local solver information for first block is in the following KSP and PC objects on rank 0:
        Use -fieldsplit_Omega_ksp_view ::ascii_info_detail to display information for all blocks
        KSP Object: (fieldsplit_Omega_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: (fieldsplit_Omega_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: (fieldsplit_Omega_sub_) 1 MPI process
                  type: seqaij
                  rows=21, cols=21
                  package used to perform factorization: petsc
                  total: nonzeros=85, allocated nonzeros=85
                    not using I-node routines
          linear system matrix, which is also used to construct the preconditioner:
          Mat Object: (fieldsplit_Omega_sub_) 1 MPI process
            type: seqaij
            rows=21, cols=21
            total: nonzeros=85, allocated nonzeros=85
            total number of mallocs used during MatSetValues calls=0
              not using I-node routines
      linear system matrix, which is also used to construct the preconditioner:
      Mat Object: (fieldsplit_Omega_) 3 MPI processes
        type: mpiaij
        rows=49, cols=49
        total: nonzeros=217, allocated nonzeros=217
        total number of mallocs used during MatSetValues calls=0
          not using I-node (on process 0) routines
    Split number 3 Defined by IS
    KSP Object: (fieldsplit_temperature_) 3 MPI processes
      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: (fieldsplit_temperature_) 3 MPI processes
      type: bjacobi
        number of blocks = 3
        Local solver information for first block is in the following KSP and PC objects on rank 0:
        Use -fieldsplit_temperature_ksp_view ::ascii_info_detail to display information for all blocks
        KSP Object: (fieldsplit_temperature_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: (fieldsplit_temperature_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: (fieldsplit_temperature_sub_) 1 MPI process
                  type: seqaij
                  rows=21, cols=21
                  package used to perform factorization: petsc
                  total: nonzeros=85, allocated nonzeros=85
                    not using I-node routines
          linear system matrix, which is also used to construct the preconditioner:
          Mat Object: (fieldsplit_temperature_sub_) 1 MPI process
            type: seqaij
            rows=21, cols=21
            total: nonzeros=85, allocated nonzeros=85
            total number of mallocs used during MatSetValues calls=0
              not using I-node routines
      linear system matrix, which is also used to construct the preconditioner:
      Mat Object: (fieldsplit_temperature_) 3 MPI processes
        type: mpiaij
        rows=49, cols=49
        total: nonzeros=217, allocated nonzeros=217
        total number of mallocs used during MatSetValues calls=0
          not using I-node (on process 0) routines
    linear system matrix, followed by the matrix used to construct the preconditioner:
    Mat Object: 3 MPI processes
      type: mffd
      rows=196, cols=196
        Matrix-free approximation:
          err=1.49012e-08 (relative error in function evaluation)
          Using wp compute h routine
              Does not compute normU
    Mat Object: 3 MPI processes
      type: mpiaij
      rows=196, cols=196, bs=4
      total: nonzeros=3472, allocated nonzeros=3472
      total number of mallocs used during MatSetValues calls=0
Number of SNES iterations = 2
