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