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