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 is the same for all blocks, as in the following KSP and PC objects on rank 0: 76 KSP Object: (fieldsplit_x_velocity_sub_) 1 MPI processes 77 type: preonly 78 maximum iterations=10000, initial guess is zero 79 tolerances: relative=1e-05, absolute=1e-50, divergence=10000. 80 left preconditioning 81 using NONE norm type for convergence test 82 PC Object: (fieldsplit_x_velocity_sub_) 1 MPI processes 83 type: ilu 84 out-of-place factorization 85 0 levels of fill 86 tolerance for zero pivot 2.22045e-14 87 matrix ordering: natural 88 factor fill ratio given 1., needed 1. 89 Factored matrix follows: 90 Mat Object: 1 MPI processes 91 type: seqaij 92 rows=21, cols=21 93 package used to perform factorization: petsc 94 total: nonzeros=85, allocated nonzeros=85 95 not using I-node routines 96 linear system matrix = precond matrix: 97 Mat Object: (fieldsplit_x_velocity_) 1 MPI processes 98 type: seqaij 99 rows=21, cols=21 100 total: nonzeros=85, allocated nonzeros=85 101 total number of mallocs used during MatSetValues calls=0 102 not using I-node routines 103 linear system matrix = precond matrix: 104 Mat Object: (fieldsplit_x_velocity_) 3 MPI processes 105 type: mpiaij 106 rows=49, cols=49 107 total: nonzeros=217, allocated nonzeros=217 108 total number of mallocs used during MatSetValues calls=0 109 not using I-node (on process 0) routines 110 Split number 1 Defined by IS 111 KSP Object: (fieldsplit_y_velocity_) 3 MPI processes 112 type: preonly 113 maximum iterations=10000, initial guess is zero 114 tolerances: relative=1e-05, absolute=1e-50, divergence=10000. 115 left preconditioning 116 using NONE norm type for convergence test 117 PC Object: (fieldsplit_y_velocity_) 3 MPI processes 118 type: bjacobi 119 number of blocks = 3 120 Local solver is the same for all blocks, as in the following KSP and PC objects on rank 0: 121 KSP Object: (fieldsplit_y_velocity_sub_) 1 MPI processes 122 type: preonly 123 maximum iterations=10000, initial guess is zero 124 tolerances: relative=1e-05, absolute=1e-50, divergence=10000. 125 left preconditioning 126 using NONE norm type for convergence test 127 PC Object: (fieldsplit_y_velocity_sub_) 1 MPI processes 128 type: ilu 129 out-of-place factorization 130 0 levels of fill 131 tolerance for zero pivot 2.22045e-14 132 matrix ordering: natural 133 factor fill ratio given 1., needed 1. 134 Factored matrix follows: 135 Mat Object: 1 MPI processes 136 type: seqaij 137 rows=21, cols=21 138 package used to perform factorization: petsc 139 total: nonzeros=85, allocated nonzeros=85 140 not using I-node routines 141 linear system matrix = precond matrix: 142 Mat Object: (fieldsplit_y_velocity_) 1 MPI processes 143 type: seqaij 144 rows=21, cols=21 145 total: nonzeros=85, allocated nonzeros=85 146 total number of mallocs used during MatSetValues calls=0 147 not using I-node routines 148 linear system matrix = precond matrix: 149 Mat Object: (fieldsplit_y_velocity_) 3 MPI processes 150 type: mpiaij 151 rows=49, cols=49 152 total: nonzeros=217, allocated nonzeros=217 153 total number of mallocs used during MatSetValues calls=0 154 not using I-node (on process 0) routines 155 Split number 2 Defined by IS 156 KSP Object: (fieldsplit_Omega_) 3 MPI processes 157 type: preonly 158 maximum iterations=10000, initial guess is zero 159 tolerances: relative=1e-05, absolute=1e-50, divergence=10000. 160 left preconditioning 161 using NONE norm type for convergence test 162 PC Object: (fieldsplit_Omega_) 3 MPI processes 163 type: bjacobi 164 number of blocks = 3 165 Local solver is the same for all blocks, as in the following KSP and PC objects on rank 0: 166 KSP Object: (fieldsplit_Omega_sub_) 1 MPI processes 167 type: preonly 168 maximum iterations=10000, initial guess is zero 169 tolerances: relative=1e-05, absolute=1e-50, divergence=10000. 170 left preconditioning 171 using NONE norm type for convergence test 172 PC Object: (fieldsplit_Omega_sub_) 1 MPI processes 173 type: ilu 174 out-of-place factorization 175 0 levels of fill 176 tolerance for zero pivot 2.22045e-14 177 matrix ordering: natural 178 factor fill ratio given 1., needed 1. 179 Factored matrix follows: 180 Mat Object: 1 MPI processes 181 type: seqaij 182 rows=21, cols=21 183 package used to perform factorization: petsc 184 total: nonzeros=85, allocated nonzeros=85 185 not using I-node routines 186 linear system matrix = precond matrix: 187 Mat Object: (fieldsplit_Omega_) 1 MPI processes 188 type: seqaij 189 rows=21, cols=21 190 total: nonzeros=85, allocated nonzeros=85 191 total number of mallocs used during MatSetValues calls=0 192 not using I-node routines 193 linear system matrix = precond matrix: 194 Mat Object: (fieldsplit_Omega_) 3 MPI processes 195 type: mpiaij 196 rows=49, cols=49 197 total: nonzeros=217, allocated nonzeros=217 198 total number of mallocs used during MatSetValues calls=0 199 not using I-node (on process 0) routines 200 Split number 3 Defined by IS 201 KSP Object: (fieldsplit_temperature_) 3 MPI processes 202 type: preonly 203 maximum iterations=10000, initial guess is zero 204 tolerances: relative=1e-05, absolute=1e-50, divergence=10000. 205 left preconditioning 206 using NONE norm type for convergence test 207 PC Object: (fieldsplit_temperature_) 3 MPI processes 208 type: bjacobi 209 number of blocks = 3 210 Local solver is the same for all blocks, as in the following KSP and PC objects on rank 0: 211 KSP Object: (fieldsplit_temperature_sub_) 1 MPI processes 212 type: preonly 213 maximum iterations=10000, initial guess is zero 214 tolerances: relative=1e-05, absolute=1e-50, divergence=10000. 215 left preconditioning 216 using NONE norm type for convergence test 217 PC Object: (fieldsplit_temperature_sub_) 1 MPI processes 218 type: ilu 219 out-of-place factorization 220 0 levels of fill 221 tolerance for zero pivot 2.22045e-14 222 matrix ordering: natural 223 factor fill ratio given 1., needed 1. 224 Factored matrix follows: 225 Mat Object: 1 MPI processes 226 type: seqaij 227 rows=21, cols=21 228 package used to perform factorization: petsc 229 total: nonzeros=85, allocated nonzeros=85 230 not using I-node routines 231 linear system matrix = precond matrix: 232 Mat Object: (fieldsplit_temperature_) 1 MPI processes 233 type: seqaij 234 rows=21, cols=21 235 total: nonzeros=85, allocated nonzeros=85 236 total number of mallocs used during MatSetValues calls=0 237 not using I-node routines 238 linear system matrix = precond matrix: 239 Mat Object: (fieldsplit_temperature_) 3 MPI processes 240 type: mpiaij 241 rows=49, cols=49 242 total: nonzeros=217, allocated nonzeros=217 243 total number of mallocs used during MatSetValues calls=0 244 not using I-node (on process 0) routines 245 linear system matrix = precond matrix: 246 Mat Object: 3 MPI processes 247 type: mpiaij 248 rows=196, cols=196, bs=4 249 total: nonzeros=3472, allocated nonzeros=3472 250 total number of mallocs used during MatSetValues calls=0 251Number of SNES iterations = 2 252