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.0554997 2 KSP Residual norm 0.0103331 3 KSP Residual norm 0.00349351 4 KSP Residual norm 0.00162771 5 KSP Residual norm 0.000420559 6 KSP Residual norm 0.000131423 7 KSP Residual norm 4.9738e-05 8 KSP Residual norm 1.91653e-05 9 KSP Residual norm 5.05572e-06 10 KSP Residual norm 1.51319e-06 11 KSP Residual norm 5.43836e-07 1 SNES Function norm 2.50119e-05 0 KSP Residual norm 2.50119e-05 1 KSP Residual norm 5.42803e-06 2 KSP Residual norm 1.19801e-06 3 KSP Residual norm 7.19919e-07 4 KSP Residual norm 3.23248e-07 5 KSP Residual norm 8.65301e-08 6 KSP Residual norm 3.40266e-08 7 KSP Residual norm 1.01576e-08 8 KSP Residual norm 1.77734e-09 9 KSP Residual norm 3.403e-10 10 KSP Residual norm 1.094e-10 2 SNES Function norm 1.093e-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=21 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: 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: 1 MPI process 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_) 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_) 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_) 1 MPI process type: seqaij rows=49, cols=49 package used to perform factorization: petsc total: nonzeros=217, allocated nonzeros=217 not using I-node routines linear system matrix, which is also used to construct the preconditioner: Mat Object: (fieldsplit_x_velocity_) 1 MPI process type: seqaij rows=49, cols=49 total: nonzeros=217, allocated nonzeros=217 total number of mallocs used during MatSetValues calls=0 not using I-node routines Split number 1 Defined by IS KSP Object: (fieldsplit_y_velocity_) 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_) 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_) 1 MPI process type: seqaij rows=49, cols=49 package used to perform factorization: petsc total: nonzeros=217, allocated nonzeros=217 not using I-node routines linear system matrix, which is also used to construct the preconditioner: Mat Object: (fieldsplit_y_velocity_) 1 MPI process type: seqaij rows=49, cols=49 total: nonzeros=217, allocated nonzeros=217 total number of mallocs used during MatSetValues calls=0 not using I-node routines Split number 2 Defined by IS KSP Object: (fieldsplit_Omega_) 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_) 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_) 1 MPI process type: seqaij rows=49, cols=49 package used to perform factorization: petsc total: nonzeros=217, allocated nonzeros=217 not using I-node routines linear system matrix, which is also used to construct the preconditioner: Mat Object: (fieldsplit_Omega_) 1 MPI process type: seqaij rows=49, cols=49 total: nonzeros=217, allocated nonzeros=217 total number of mallocs used during MatSetValues calls=0 not using I-node routines Split number 3 Defined by IS KSP Object: (fieldsplit_temperature_) 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_) 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_) 1 MPI process type: seqaij rows=49, cols=49 package used to perform factorization: petsc total: nonzeros=217, allocated nonzeros=217 not using I-node routines linear system matrix, which is also used to construct the preconditioner: Mat Object: (fieldsplit_temperature_) 1 MPI process type: seqaij rows=49, cols=49 total: nonzeros=217, allocated nonzeros=217 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: 1 MPI process type: seqaij rows=196, cols=196, bs=4 total: nonzeros=3472, allocated nonzeros=3472 total number of mallocs used during MatSetValues calls=0 using I-node routines: found 49 nodes, limit used is 5 Number of SNES iterations = 2