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