lid velocity = 0.0625, prandtl # = 1., grashof # = 1. 0 KSP Residual norm 0.244212 1 KSP Residual norm 0.0263811 2 KSP Residual norm 0.001012 3 KSP Residual norm 0.000219671 4 KSP Residual norm 8.29849e-06 5 KSP Residual norm 1.63718e-06 0 KSP Residual norm 2.29703e-05 1 KSP Residual norm 2.80698e-06 2 KSP Residual norm 6.02504e-07 3 KSP Residual norm 1.59775e-07 4 KSP Residual norm 6.11071e-09 5 KSP Residual norm 7.777e-10 6 KSP Residual norm 8.330e-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=11 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 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: 2 MPI processes 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: 2 MPI processes type: parms global preconditioner: Restricted Additive Schwarz local preconditioner: ARMS non-zero elements/original non-zero entries: 0.75 Tolerance for local solve: 0.01 Number of levels: 10 Using nonsymmetric permutation Block size: 250 Tolerance for independent sets: 0.4 Inner Krylov dimension: 0 Maximum number of inner iterations: 0 amount of fill-in for ilut, iluk and arms: 20 amount of fill-in for schur: 20 amount of fill-in for ILUT L and U: 20 drop tolerance for L, U, L^{-1}F and EU^{-1}: 1e-05 drop tolerance for schur complement at each level: 0.001 drop tolerance for ILUT in last level schur complement: 0.001 linear system matrix, which is also used to construct the preconditioner: Mat Object: 2 MPI processes type: mpiaij 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