lid velocity = 0.000106281, prandtl # = 1., grashof # = 1. 0 SNES Function norm 1.036007954337e-02 0 KSP Residual norm 4.386611885521e-01 1 KSP Residual norm 3.127838547648e-01 2 KSP Residual norm 2.334238799976e-01 3 KSP Residual norm 1.318350769046e-01 4 KSP Residual norm 8.143004995277e-02 5 KSP Residual norm 4.813457707703e-02 6 KSP Residual norm 2.733719484710e-02 7 KSP Residual norm 1.462151166610e-02 8 KSP Residual norm 8.331275159714e-03 9 KSP Residual norm 5.275884606772e-03 10 KSP Residual norm 3.768690746652e-03 11 KSP Residual norm 2.699365250849e-03 12 KSP Residual norm 1.824681696989e-03 13 KSP Residual norm 1.179126381375e-03 14 KSP Residual norm 7.107306905422e-04 15 KSP Residual norm 4.155820345074e-04 16 KSP Residual norm 2.408438645714e-04 17 KSP Residual norm 1.476855736514e-04 18 KSP Residual norm 9.907578845784e-05 19 KSP Residual norm 6.797328938557e-05 20 KSP Residual norm 4.675944545290e-05 21 KSP Residual norm 3.352583819988e-05 22 KSP Residual norm 2.441653843583e-05 23 KSP Residual norm 1.768579774951e-05 24 KSP Residual norm 1.411236998154e-05 25 KSP Residual norm 1.062751951795e-05 26 KSP Residual norm 7.402652132828e-06 27 KSP Residual norm 5.376101942565e-06 28 KSP Residual norm 3.936750059583e-06 1 SNES Function norm 1.539897001872e-06 0 KSP Residual norm 1.219799734452e-05 1 KSP Residual norm 9.351061917402e-06 2 KSP Residual norm 5.158749230924e-06 3 KSP Residual norm 3.800583255559e-06 4 KSP Residual norm 3.358178466662e-06 5 KSP Residual norm 2.704711328549e-06 6 KSP Residual norm 1.965512422061e-06 7 KSP Residual norm 1.383231224585e-06 8 KSP Residual norm 1.016838426400e-06 9 KSP Residual norm 7.817076853999e-07 10 KSP Residual norm 5.714787985538e-07 11 KSP Residual norm 4.070517298179e-07 12 KSP Residual norm 3.199366855723e-07 13 KSP Residual norm 2.498338700906e-07 14 KSP Residual norm 1.774654138694e-07 15 KSP Residual norm 1.220894443119e-07 16 KSP Residual norm 8.191467113032e-08 17 KSP Residual norm 5.471425416521e-08 18 KSP Residual norm 3.800399846106e-08 19 KSP Residual norm 2.624332673801e-08 20 KSP Residual norm 1.817665454829e-08 21 KSP Residual norm 1.193934944948e-08 22 KSP Residual norm 7.058845517318e-09 23 KSP Residual norm 4.801132366324e-09 24 KSP Residual norm 3.136888134931e-09 25 KSP Residual norm 1.908878704398e-09 26 KSP Residual norm 1.179297203589e-09 27 KSP Residual norm 7.450436712605e-10 28 KSP Residual norm 4.727201847673e-10 29 KSP Residual norm 3.071103410289e-10 30 KSP Residual norm 1.954558918512e-10 31 KSP Residual norm 1.476035656576e-10 32 KSP Residual norm 9.544641423593e-11 2 SNES Function norm 2.155536219955e-11 SNES Object: 4 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=60 total number of function evaluations=3 norm schedule ALWAYS Jacobian is built using colored finite differences on a DMDA SNESLineSearch Object: 4 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: 4 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: 4 MPI processes type: ml type is MULTIPLICATIVE, levels=5 cycles=v Cycles per PCApply=1 Using externally compute Galerkin coarse grid matrices Coarse grid solver -- level ------------------------------- KSP Object: (mg_coarse_) 4 MPI processes 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: (mg_coarse_) 4 MPI processes type: redundant First (color=0) of 4 PCs follows KSP Object: (mg_coarse_redundant_) 1 MPI processes 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: (mg_coarse_redundant_) 1 MPI processes type: lu out-of-place factorization tolerance for zero pivot 2.22045e-14 using diagonal shift on blocks to prevent zero pivot [INBLOCKS] matrix ordering: nd factor fill ratio given 5., needed 1. Factored matrix follows: Mat Object: 1 MPI processes type: seqaij rows=16, cols=16, bs=4 package used to perform factorization: petsc total: nonzeros=256, allocated nonzeros=256 using I-node routines: found 4 nodes, limit used is 5 linear system matrix = precond matrix: Mat Object: 1 MPI processes type: seqaij rows=16, cols=16, bs=4 total: nonzeros=256, allocated nonzeros=256 total number of mallocs used during MatSetValues calls=0 using I-node routines: found 4 nodes, limit used is 5 linear system matrix = precond matrix: Mat Object: 4 MPI processes type: mpiaij rows=16, cols=16, bs=4 total: nonzeros=256, allocated nonzeros=256 total number of mallocs used during MatSetValues calls=0 using I-node (on process 0) routines: found 1 nodes, limit used is 5 Down solver (pre-smoother) on level 1 ------------------------------- KSP Object: (mg_levels_1_) 4 MPI processes type: richardson damping factor=1. maximum iterations=2, nonzero initial guess tolerances: relative=1e-05, absolute=1e-50, divergence=10000. left preconditioning using NONE norm type for convergence test PC Object: (mg_levels_1_) 4 MPI processes type: sor type = local_symmetric, iterations = 1, local iterations = 1, omega = 1. linear system matrix = precond matrix: Mat Object: 4 MPI processes type: mpiaij rows=88, cols=88, bs=4 total: nonzeros=4896, allocated nonzeros=4896 total number of mallocs used during MatSetValues calls=0 using I-node (on process 0) routines: found 5 nodes, limit used is 5 Up solver (post-smoother) same as down solver (pre-smoother) Down solver (pre-smoother) on level 2 ------------------------------- KSP Object: (mg_levels_2_) 4 MPI processes type: richardson damping factor=1. maximum iterations=2, nonzero initial guess tolerances: relative=1e-05, absolute=1e-50, divergence=10000. left preconditioning using NONE norm type for convergence test PC Object: (mg_levels_2_) 4 MPI processes type: sor type = local_symmetric, iterations = 1, local iterations = 1, omega = 1. linear system matrix = precond matrix: Mat Object: 4 MPI processes type: mpiaij rows=772, cols=772, bs=4 total: nonzeros=35364, allocated nonzeros=35364 total number of mallocs used during MatSetValues calls=0 using I-node (on process 0) routines: found 56 nodes, limit used is 5 Up solver (post-smoother) same as down solver (pre-smoother) Down solver (pre-smoother) on level 3 ------------------------------- KSP Object: (mg_levels_3_) 4 MPI processes type: richardson damping factor=1. maximum iterations=2, nonzero initial guess tolerances: relative=1e-05, absolute=1e-50, divergence=10000. left preconditioning using NONE norm type for convergence test PC Object: (mg_levels_3_) 4 MPI processes type: sor type = local_symmetric, iterations = 1, local iterations = 1, omega = 1. linear system matrix = precond matrix: Mat Object: 4 MPI processes type: mpiaij rows=7156, cols=7156, bs=4 total: nonzeros=258241, allocated nonzeros=258241 total number of mallocs used during MatSetValues calls=0 not using I-node (on process 0) routines Up solver (post-smoother) same as down solver (pre-smoother) Down solver (pre-smoother) on level 4 ------------------------------- KSP Object: (mg_levels_4_) 4 MPI processes type: richardson damping factor=1. maximum iterations=2, nonzero initial guess tolerances: relative=1e-05, absolute=1e-50, divergence=10000. left preconditioning using NONE norm type for convergence test PC Object: (mg_levels_4_) 4 MPI processes type: sor type = local_symmetric, iterations = 1, local iterations = 1, omega = 1. linear system matrix = precond matrix: Mat Object: 4 MPI processes type: mpiaij rows=37636, cols=37636, bs=4 total: nonzeros=746512, allocated nonzeros=746512 total number of mallocs used during MatSetValues calls=0 Up solver (post-smoother) same as down solver (pre-smoother) linear system matrix = precond matrix: Mat Object: 4 MPI processes type: mpiaij rows=37636, cols=37636, bs=4 total: nonzeros=746512, allocated nonzeros=746512 total number of mallocs used during MatSetValues calls=0 Number of SNES iterations = 2