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