1c4762a1bSJed Brownlid velocity = 0.0204082, prandtl # = 1., grashof # = 1. 2c4762a1bSJed Brown 0 SNES Function norm 0.146194 3c4762a1bSJed Brown 0 KSP Residual norm 0.146194 4c4762a1bSJed Brown 1 KSP Residual norm 0.0781702 5c4762a1bSJed Brown 2 KSP Residual norm 0.0289896 6c4762a1bSJed Brown 3 KSP Residual norm 0.00874626 7c4762a1bSJed Brown 4 KSP Residual norm 0.00360523 8c4762a1bSJed Brown 5 KSP Residual norm 0.00187927 9c4762a1bSJed Brown 6 KSP Residual norm 0.000997932 10c4762a1bSJed Brown 7 KSP Residual norm 0.000325883 11c4762a1bSJed Brown 8 KSP Residual norm 0.000122254 12c4762a1bSJed Brown 9 KSP Residual norm 6.21287e-05 13c4762a1bSJed Brown 10 KSP Residual norm 4.04759e-05 14c4762a1bSJed Brown 11 KSP Residual norm 2.1448e-05 15c4762a1bSJed Brown 12 KSP Residual norm 9.92327e-06 16c4762a1bSJed Brown 13 KSP Residual norm 4.95927e-06 17c4762a1bSJed Brown 14 KSP Residual norm 2.41652e-06 18c4762a1bSJed Brown 15 KSP Residual norm 1.39592e-06 19c4762a1bSJed Brown 1 SNES Function norm 2.5489e-05 20c4762a1bSJed Brown 0 KSP Residual norm 2.5489e-05 21c4762a1bSJed Brown 1 KSP Residual norm 1.07362e-05 22c4762a1bSJed Brown 2 KSP Residual norm 3.78935e-06 23c4762a1bSJed Brown 3 KSP Residual norm 1.43754e-06 24c4762a1bSJed Brown 4 KSP Residual norm 6.54959e-07 25c4762a1bSJed Brown 5 KSP Residual norm 5.12221e-07 26c4762a1bSJed Brown 6 KSP Residual norm 3.42561e-07 27c4762a1bSJed Brown 7 KSP Residual norm 1.45139e-07 28c4762a1bSJed Brown 8 KSP Residual norm 7.40916e-08 29c4762a1bSJed Brown 9 KSP Residual norm 3.79621e-08 30c4762a1bSJed Brown 10 KSP Residual norm 1.50283e-08 31c4762a1bSJed Brown 11 KSP Residual norm 5.03342e-09 32c4762a1bSJed Brown 12 KSP Residual norm 2.41715e-09 33c4762a1bSJed Brown 13 KSP Residual norm 1.38641e-09 34c4762a1bSJed Brown 14 KSP Residual norm 7.047e-10 35c4762a1bSJed Brown 15 KSP Residual norm 3.897e-10 36c4762a1bSJed Brown 16 KSP Residual norm 2.021e-10 37c4762a1bSJed Brown 2 SNES Function norm 2.020e-10 38c4762a1bSJed BrownSNES Object: 3 MPI processes 39c4762a1bSJed Brown type: newtonls 40c4762a1bSJed Brown maximum iterations=50, maximum function evaluations=10000 41c4762a1bSJed Brown tolerances: relative=1e-08, absolute=1e-50, solution=1e-08 42c4762a1bSJed Brown total number of linear solver iterations=31 43c4762a1bSJed Brown total number of function evaluations=36 44c4762a1bSJed Brown norm schedule ALWAYS 45c4762a1bSJed Brown Jacobian is applied matrix-free with differencing 46c4762a1bSJed Brown Preconditioning Jacobian is built using colored finite differences on a DMDA 47c4762a1bSJed Brown SNESLineSearch Object: 3 MPI processes 48c4762a1bSJed Brown type: bt 49c4762a1bSJed Brown interpolation: cubic 50c4762a1bSJed Brown alpha=1.000000e-04 51c4762a1bSJed Brown maxstep=1.000000e+08, minlambda=1.000000e-12 52c4762a1bSJed Brown tolerances: relative=1.000000e-08, absolute=1.000000e-15, lambda=1.000000e-08 53c4762a1bSJed Brown maximum iterations=40 54c4762a1bSJed Brown KSP Object: 3 MPI processes 55c4762a1bSJed Brown type: fgmres 56c4762a1bSJed Brown restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement 57c4762a1bSJed Brown happy breakdown tolerance 1e-30 58c4762a1bSJed Brown maximum iterations=10000, initial guess is zero 59c4762a1bSJed Brown tolerances: relative=1e-05, absolute=1e-50, divergence=10000. 60c4762a1bSJed Brown right preconditioning 61c4762a1bSJed Brown using UNPRECONDITIONED norm type for convergence test 62c4762a1bSJed Brown PC Object: 3 MPI processes 63c4762a1bSJed Brown type: fieldsplit 64c4762a1bSJed Brown FieldSplit with MULTIPLICATIVE composition: total splits = 4 65c4762a1bSJed Brown Solver info for each split is in the following KSP objects: 66c4762a1bSJed Brown Split number 0 Defined by IS 67c4762a1bSJed Brown KSP Object: (fieldsplit_x_velocity_) 3 MPI processes 68c4762a1bSJed Brown type: preonly 69c4762a1bSJed Brown maximum iterations=10000, initial guess is zero 70c4762a1bSJed Brown tolerances: relative=1e-05, absolute=1e-50, divergence=10000. 71c4762a1bSJed Brown left preconditioning 72c4762a1bSJed Brown using NONE norm type for convergence test 73c4762a1bSJed Brown PC Object: (fieldsplit_x_velocity_) 3 MPI processes 74c4762a1bSJed Brown type: bjacobi 75c4762a1bSJed Brown number of blocks = 3 768d76b567SPatrick Sanan Local solver is the same for all blocks, as in the following KSP and PC objects on rank 0: 77c4762a1bSJed Brown KSP Object: (fieldsplit_x_velocity_sub_) 1 MPI processes 78c4762a1bSJed Brown type: preonly 79c4762a1bSJed Brown maximum iterations=10000, initial guess is zero 80c4762a1bSJed Brown tolerances: relative=1e-05, absolute=1e-50, divergence=10000. 81c4762a1bSJed Brown left preconditioning 82c4762a1bSJed Brown using NONE norm type for convergence test 83c4762a1bSJed Brown PC Object: (fieldsplit_x_velocity_sub_) 1 MPI processes 84c4762a1bSJed Brown type: ilu 85c4762a1bSJed Brown out-of-place factorization 86c4762a1bSJed Brown 0 levels of fill 87c4762a1bSJed Brown tolerance for zero pivot 2.22045e-14 88c4762a1bSJed Brown matrix ordering: natural 89c4762a1bSJed Brown factor fill ratio given 1., needed 1. 90c4762a1bSJed Brown Factored matrix follows: 91c4762a1bSJed Brown Mat Object: 1 MPI processes 92c4762a1bSJed Brown type: seqaij 93c4762a1bSJed Brown rows=21, cols=21 94c4762a1bSJed Brown package used to perform factorization: petsc 95c4762a1bSJed Brown total: nonzeros=85, allocated nonzeros=85 96c4762a1bSJed Brown not using I-node routines 97c4762a1bSJed Brown linear system matrix = precond matrix: 98*793a3527SPierre Jolivet Mat Object: (fieldsplit_x_velocity_sub_) 1 MPI processes 99c4762a1bSJed Brown type: seqaij 100c4762a1bSJed Brown rows=21, cols=21 101c4762a1bSJed Brown total: nonzeros=85, allocated nonzeros=85 102c4762a1bSJed Brown total number of mallocs used during MatSetValues calls=0 103c4762a1bSJed Brown not using I-node routines 104c4762a1bSJed Brown linear system matrix = precond matrix: 105c4762a1bSJed Brown Mat Object: (fieldsplit_x_velocity_) 3 MPI processes 106c4762a1bSJed Brown type: mpiaij 107c4762a1bSJed Brown rows=49, cols=49 108c4762a1bSJed Brown total: nonzeros=217, allocated nonzeros=217 109c4762a1bSJed Brown total number of mallocs used during MatSetValues calls=0 110c4762a1bSJed Brown not using I-node (on process 0) routines 111c4762a1bSJed Brown Split number 1 Defined by IS 112c4762a1bSJed Brown KSP Object: (fieldsplit_y_velocity_) 3 MPI processes 113c4762a1bSJed Brown type: preonly 114c4762a1bSJed Brown maximum iterations=10000, initial guess is zero 115c4762a1bSJed Brown tolerances: relative=1e-05, absolute=1e-50, divergence=10000. 116c4762a1bSJed Brown left preconditioning 117c4762a1bSJed Brown using NONE norm type for convergence test 118c4762a1bSJed Brown PC Object: (fieldsplit_y_velocity_) 3 MPI processes 119c4762a1bSJed Brown type: bjacobi 120c4762a1bSJed Brown number of blocks = 3 1218d76b567SPatrick Sanan Local solver is the same for all blocks, as in the following KSP and PC objects on rank 0: 122c4762a1bSJed Brown KSP Object: (fieldsplit_y_velocity_sub_) 1 MPI processes 123c4762a1bSJed Brown type: preonly 124c4762a1bSJed Brown maximum iterations=10000, initial guess is zero 125c4762a1bSJed Brown tolerances: relative=1e-05, absolute=1e-50, divergence=10000. 126c4762a1bSJed Brown left preconditioning 127c4762a1bSJed Brown using NONE norm type for convergence test 128c4762a1bSJed Brown PC Object: (fieldsplit_y_velocity_sub_) 1 MPI processes 129c4762a1bSJed Brown type: ilu 130c4762a1bSJed Brown out-of-place factorization 131c4762a1bSJed Brown 0 levels of fill 132c4762a1bSJed Brown tolerance for zero pivot 2.22045e-14 133c4762a1bSJed Brown matrix ordering: natural 134c4762a1bSJed Brown factor fill ratio given 1., needed 1. 135c4762a1bSJed Brown Factored matrix follows: 136c4762a1bSJed Brown Mat Object: 1 MPI processes 137c4762a1bSJed Brown type: seqaij 138c4762a1bSJed Brown rows=21, cols=21 139c4762a1bSJed Brown package used to perform factorization: petsc 140c4762a1bSJed Brown total: nonzeros=85, allocated nonzeros=85 141c4762a1bSJed Brown not using I-node routines 142c4762a1bSJed Brown linear system matrix = precond matrix: 143*793a3527SPierre Jolivet Mat Object: (fieldsplit_y_velocity_sub_) 1 MPI processes 144c4762a1bSJed Brown type: seqaij 145c4762a1bSJed Brown rows=21, cols=21 146c4762a1bSJed Brown total: nonzeros=85, allocated nonzeros=85 147c4762a1bSJed Brown total number of mallocs used during MatSetValues calls=0 148c4762a1bSJed Brown not using I-node routines 149c4762a1bSJed Brown linear system matrix = precond matrix: 150c4762a1bSJed Brown Mat Object: (fieldsplit_y_velocity_) 3 MPI processes 151c4762a1bSJed Brown type: mpiaij 152c4762a1bSJed Brown rows=49, cols=49 153c4762a1bSJed Brown total: nonzeros=217, allocated nonzeros=217 154c4762a1bSJed Brown total number of mallocs used during MatSetValues calls=0 155c4762a1bSJed Brown not using I-node (on process 0) routines 156c4762a1bSJed Brown Split number 2 Defined by IS 157c4762a1bSJed Brown KSP Object: (fieldsplit_Omega_) 3 MPI processes 158c4762a1bSJed Brown type: preonly 159c4762a1bSJed Brown maximum iterations=10000, initial guess is zero 160c4762a1bSJed Brown tolerances: relative=1e-05, absolute=1e-50, divergence=10000. 161c4762a1bSJed Brown left preconditioning 162c4762a1bSJed Brown using NONE norm type for convergence test 163c4762a1bSJed Brown PC Object: (fieldsplit_Omega_) 3 MPI processes 164c4762a1bSJed Brown type: bjacobi 165c4762a1bSJed Brown number of blocks = 3 1668d76b567SPatrick Sanan Local solver is the same for all blocks, as in the following KSP and PC objects on rank 0: 167c4762a1bSJed Brown KSP Object: (fieldsplit_Omega_sub_) 1 MPI processes 168c4762a1bSJed Brown type: preonly 169c4762a1bSJed Brown maximum iterations=10000, initial guess is zero 170c4762a1bSJed Brown tolerances: relative=1e-05, absolute=1e-50, divergence=10000. 171c4762a1bSJed Brown left preconditioning 172c4762a1bSJed Brown using NONE norm type for convergence test 173c4762a1bSJed Brown PC Object: (fieldsplit_Omega_sub_) 1 MPI processes 174c4762a1bSJed Brown type: ilu 175c4762a1bSJed Brown out-of-place factorization 176c4762a1bSJed Brown 0 levels of fill 177c4762a1bSJed Brown tolerance for zero pivot 2.22045e-14 178c4762a1bSJed Brown matrix ordering: natural 179c4762a1bSJed Brown factor fill ratio given 1., needed 1. 180c4762a1bSJed Brown Factored matrix follows: 181c4762a1bSJed Brown Mat Object: 1 MPI processes 182c4762a1bSJed Brown type: seqaij 183c4762a1bSJed Brown rows=21, cols=21 184c4762a1bSJed Brown package used to perform factorization: petsc 185c4762a1bSJed Brown total: nonzeros=85, allocated nonzeros=85 186c4762a1bSJed Brown not using I-node routines 187c4762a1bSJed Brown linear system matrix = precond matrix: 188*793a3527SPierre Jolivet Mat Object: (fieldsplit_Omega_sub_) 1 MPI processes 189c4762a1bSJed Brown type: seqaij 190c4762a1bSJed Brown rows=21, cols=21 191c4762a1bSJed Brown total: nonzeros=85, allocated nonzeros=85 192c4762a1bSJed Brown total number of mallocs used during MatSetValues calls=0 193c4762a1bSJed Brown not using I-node routines 194c4762a1bSJed Brown linear system matrix = precond matrix: 195c4762a1bSJed Brown Mat Object: (fieldsplit_Omega_) 3 MPI processes 196c4762a1bSJed Brown type: mpiaij 197c4762a1bSJed Brown rows=49, cols=49 198c4762a1bSJed Brown total: nonzeros=217, allocated nonzeros=217 199c4762a1bSJed Brown total number of mallocs used during MatSetValues calls=0 200c4762a1bSJed Brown not using I-node (on process 0) routines 201c4762a1bSJed Brown Split number 3 Defined by IS 202c4762a1bSJed Brown KSP Object: (fieldsplit_temperature_) 3 MPI processes 203c4762a1bSJed Brown type: preonly 204c4762a1bSJed Brown maximum iterations=10000, initial guess is zero 205c4762a1bSJed Brown tolerances: relative=1e-05, absolute=1e-50, divergence=10000. 206c4762a1bSJed Brown left preconditioning 207c4762a1bSJed Brown using NONE norm type for convergence test 208c4762a1bSJed Brown PC Object: (fieldsplit_temperature_) 3 MPI processes 209c4762a1bSJed Brown type: bjacobi 210c4762a1bSJed Brown number of blocks = 3 2118d76b567SPatrick Sanan Local solver is the same for all blocks, as in the following KSP and PC objects on rank 0: 212c4762a1bSJed Brown KSP Object: (fieldsplit_temperature_sub_) 1 MPI processes 213c4762a1bSJed Brown type: preonly 214c4762a1bSJed Brown maximum iterations=10000, initial guess is zero 215c4762a1bSJed Brown tolerances: relative=1e-05, absolute=1e-50, divergence=10000. 216c4762a1bSJed Brown left preconditioning 217c4762a1bSJed Brown using NONE norm type for convergence test 218c4762a1bSJed Brown PC Object: (fieldsplit_temperature_sub_) 1 MPI processes 219c4762a1bSJed Brown type: ilu 220c4762a1bSJed Brown out-of-place factorization 221c4762a1bSJed Brown 0 levels of fill 222c4762a1bSJed Brown tolerance for zero pivot 2.22045e-14 223c4762a1bSJed Brown matrix ordering: natural 224c4762a1bSJed Brown factor fill ratio given 1., needed 1. 225c4762a1bSJed Brown Factored matrix follows: 226c4762a1bSJed Brown Mat Object: 1 MPI processes 227c4762a1bSJed Brown type: seqaij 228c4762a1bSJed Brown rows=21, cols=21 229c4762a1bSJed Brown package used to perform factorization: petsc 230c4762a1bSJed Brown total: nonzeros=85, allocated nonzeros=85 231c4762a1bSJed Brown not using I-node routines 232c4762a1bSJed Brown linear system matrix = precond matrix: 233*793a3527SPierre Jolivet Mat Object: (fieldsplit_temperature_sub_) 1 MPI processes 234c4762a1bSJed Brown type: seqaij 235c4762a1bSJed Brown rows=21, cols=21 236c4762a1bSJed Brown total: nonzeros=85, allocated nonzeros=85 237c4762a1bSJed Brown total number of mallocs used during MatSetValues calls=0 238c4762a1bSJed Brown not using I-node routines 239c4762a1bSJed Brown linear system matrix = precond matrix: 240c4762a1bSJed Brown Mat Object: (fieldsplit_temperature_) 3 MPI processes 241c4762a1bSJed Brown type: mpiaij 242c4762a1bSJed Brown rows=49, cols=49 243c4762a1bSJed Brown total: nonzeros=217, allocated nonzeros=217 244c4762a1bSJed Brown total number of mallocs used during MatSetValues calls=0 245c4762a1bSJed Brown not using I-node (on process 0) routines 246c4762a1bSJed Brown linear system matrix followed by preconditioner matrix: 247c4762a1bSJed Brown Mat Object: 3 MPI processes 248c4762a1bSJed Brown type: mffd 249c4762a1bSJed Brown rows=196, cols=196 250c4762a1bSJed Brown Matrix-free approximation: 251c4762a1bSJed Brown err=1.49012e-08 (relative error in function evaluation) 252c4762a1bSJed Brown Using wp compute h routine 253c4762a1bSJed Brown Does not compute normU 254c4762a1bSJed Brown Mat Object: 3 MPI processes 255c4762a1bSJed Brown type: mpiaij 256c4762a1bSJed Brown rows=196, cols=196, bs=4 257c4762a1bSJed Brown total: nonzeros=3472, allocated nonzeros=3472 258c4762a1bSJed Brown total number of mallocs used during MatSetValues calls=0 259c4762a1bSJed BrownNumber of SNES iterations = 2 260