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