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.0479011 2 KSP Residual norm 0.010778 3 KSP Residual norm 0.00272661 4 KSP Residual norm 0.000690594 5 KSP Residual norm 0.000210687 6 KSP Residual norm 0.000105322 7 KSP Residual norm 6.5975e-05 8 KSP Residual norm 2.80726e-05 9 KSP Residual norm 1.87879e-05 10 KSP Residual norm 8.24594e-06 11 KSP Residual norm 1.72955e-06 12 KSP Residual norm 5.17053e-07 1 SNES Function norm 2.5023e-05 0 KSP Residual norm 2.5023e-05 1 KSP Residual norm 5.03001e-06 2 KSP Residual norm 1.40061e-06 3 KSP Residual norm 5.68423e-07 4 KSP Residual norm 1.98909e-07 5 KSP Residual norm 9.6167e-08 6 KSP Residual norm 4.64061e-08 7 KSP Residual norm 9.70042e-09 8 KSP Residual norm 2.76341e-09 9 KSP Residual norm 9.531e-10 10 KSP Residual norm 4.491e-10 11 KSP Residual norm 2.978e-10 12 KSP Residual norm 1.484e-10 2 SNES Function norm 1.483e-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=24 total number of function evaluations=3 norm schedule ALWAYS Jacobian is built using colored finite differences on a DMDA SNESLineSearch Object: 3 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: 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 SYMMETRIC_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 using NONE norm type for convergence test 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 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: (fieldsplit_x_velocity_sub_) 1 MPI processes 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: 1 MPI processes 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 = precond matrix: Mat Object: (fieldsplit_x_velocity_sub_) 1 MPI processes 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 = precond matrix: 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 using NONE norm type for convergence test 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 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: (fieldsplit_y_velocity_sub_) 1 MPI processes 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: 1 MPI processes 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 = precond matrix: Mat Object: (fieldsplit_y_velocity_sub_) 1 MPI processes 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 = precond matrix: 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 using NONE norm type for convergence test 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 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: (fieldsplit_Omega_sub_) 1 MPI processes 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: 1 MPI processes 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 = precond matrix: Mat Object: (fieldsplit_Omega_sub_) 1 MPI processes 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 = precond matrix: 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 using NONE norm type for convergence test 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 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: (fieldsplit_temperature_sub_) 1 MPI processes 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: 1 MPI processes 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 = precond matrix: Mat Object: (fieldsplit_temperature_sub_) 1 MPI processes 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 = precond matrix: 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 = precond matrix: 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