lid velocity = 0.00591716, prandtl # = 1., grashof # = 1. 0 KSP Residual norm 0.22318 1 KSP Residual norm 0.0399771 2 KSP Residual norm 0.0263924 3 KSP Residual norm 0.0415116 4 KSP Residual norm 0.00334741 5 KSP Residual norm 0.00188077 6 KSP Residual norm 0.000591979 7 KSP Residual norm 0.000847763 8 KSP Residual norm 9.09501e-05 9 KSP Residual norm 4.07915e-05 10 KSP Residual norm 1.45997e-05 11 KSP Residual norm 2.32689e-05 12 KSP Residual norm 3.48513e-06 13 KSP Residual norm 4.00787e-06 14 KSP Residual norm 2.49517e-07 0 KSP Residual norm 6.24163e-06 1 KSP Residual norm 1.89562e-06 2 KSP Residual norm 9.1461e-07 3 KSP Residual norm 4.2732e-07 4 KSP Residual norm 4.10552e-07 5 KSP Residual norm 7.36292e-08 6 KSP Residual norm 4.22272e-08 7 KSP Residual norm 2.60372e-08 8 KSP Residual norm 1.40502e-08 9 KSP Residual norm 8.34563e-09 10 KSP Residual norm 9.87196e-09 11 KSP Residual norm 1.52953e-09 12 KSP Residual norm 2.726e-10 13 KSP Residual norm 2.133e-10 14 KSP Residual norm 9.213e-11 15 KSP Residual norm 7.141e-11 16 KSP Residual norm 3.286e-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 maxstep=1.000000e+08, 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: ibcgs 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 using NONE norm type for convergence test 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 follows: 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 = precond matrix: 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 = precond matrix: 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