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