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 maxstep=1.000000e+08, 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 using NONE norm type for convergence test 73 PC Object: (fieldsplit_x_velocity_) 3 MPI processes 74 type: bjacobi 75 number of blocks = 3 76 Local solver is the same for all blocks, as in the following KSP and PC objects on rank 0: 77 KSP Object: (fieldsplit_x_velocity_sub_) 1 MPI processes 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 processes 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 processes 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_) 1 MPI processes 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 is the same for all blocks, as in the following KSP and PC objects on rank 0: 122 KSP Object: (fieldsplit_y_velocity_sub_) 1 MPI processes 123 type: preonly 124 maximum iterations=10000, initial guess is zero 125 tolerances: relative=1e-05, absolute=1e-50, divergence=10000. 126 left preconditioning 127 using NONE norm type for convergence test 128 PC Object: (fieldsplit_y_velocity_sub_) 1 MPI processes 129 type: ilu 130 out-of-place factorization 131 0 levels of fill 132 tolerance for zero pivot 2.22045e-14 133 matrix ordering: natural 134 factor fill ratio given 1., needed 1. 135 Factored matrix follows: 136 Mat Object: 1 MPI processes 137 type: seqaij 138 rows=21, cols=21 139 package used to perform factorization: petsc 140 total: nonzeros=85, allocated nonzeros=85 141 not using I-node routines 142 linear system matrix = precond matrix: 143 Mat Object: (fieldsplit_y_velocity_) 1 MPI processes 144 type: seqaij 145 rows=21, cols=21 146 total: nonzeros=85, allocated nonzeros=85 147 total number of mallocs used during MatSetValues calls=0 148 not using I-node routines 149 linear system matrix = precond matrix: 150 Mat Object: (fieldsplit_y_velocity_) 3 MPI processes 151 type: mpiaij 152 rows=49, cols=49 153 total: nonzeros=217, allocated nonzeros=217 154 total number of mallocs used during MatSetValues calls=0 155 not using I-node (on process 0) routines 156 Split number 2 Defined by IS 157 KSP Object: (fieldsplit_Omega_) 3 MPI processes 158 type: preonly 159 maximum iterations=10000, initial guess is zero 160 tolerances: relative=1e-05, absolute=1e-50, divergence=10000. 161 left preconditioning 162 using NONE norm type for convergence test 163 PC Object: (fieldsplit_Omega_) 3 MPI processes 164 type: bjacobi 165 number of blocks = 3 166 Local solver is the same for all blocks, as in the following KSP and PC objects on rank 0: 167 KSP Object: (fieldsplit_Omega_sub_) 1 MPI processes 168 type: preonly 169 maximum iterations=10000, initial guess is zero 170 tolerances: relative=1e-05, absolute=1e-50, divergence=10000. 171 left preconditioning 172 using NONE norm type for convergence test 173 PC Object: (fieldsplit_Omega_sub_) 1 MPI processes 174 type: ilu 175 out-of-place factorization 176 0 levels of fill 177 tolerance for zero pivot 2.22045e-14 178 matrix ordering: natural 179 factor fill ratio given 1., needed 1. 180 Factored matrix follows: 181 Mat Object: 1 MPI processes 182 type: seqaij 183 rows=21, cols=21 184 package used to perform factorization: petsc 185 total: nonzeros=85, allocated nonzeros=85 186 not using I-node routines 187 linear system matrix = precond matrix: 188 Mat Object: (fieldsplit_Omega_) 1 MPI processes 189 type: seqaij 190 rows=21, cols=21 191 total: nonzeros=85, allocated nonzeros=85 192 total number of mallocs used during MatSetValues calls=0 193 not using I-node routines 194 linear system matrix = precond matrix: 195 Mat Object: (fieldsplit_Omega_) 3 MPI processes 196 type: mpiaij 197 rows=49, cols=49 198 total: nonzeros=217, allocated nonzeros=217 199 total number of mallocs used during MatSetValues calls=0 200 not using I-node (on process 0) routines 201 Split number 3 Defined by IS 202 KSP Object: (fieldsplit_temperature_) 3 MPI processes 203 type: preonly 204 maximum iterations=10000, initial guess is zero 205 tolerances: relative=1e-05, absolute=1e-50, divergence=10000. 206 left preconditioning 207 using NONE norm type for convergence test 208 PC Object: (fieldsplit_temperature_) 3 MPI processes 209 type: bjacobi 210 number of blocks = 3 211 Local solver is the same for all blocks, as in the following KSP and PC objects on rank 0: 212 KSP Object: (fieldsplit_temperature_sub_) 1 MPI processes 213 type: preonly 214 maximum iterations=10000, initial guess is zero 215 tolerances: relative=1e-05, absolute=1e-50, divergence=10000. 216 left preconditioning 217 using NONE norm type for convergence test 218 PC Object: (fieldsplit_temperature_sub_) 1 MPI processes 219 type: ilu 220 out-of-place factorization 221 0 levels of fill 222 tolerance for zero pivot 2.22045e-14 223 matrix ordering: natural 224 factor fill ratio given 1., needed 1. 225 Factored matrix follows: 226 Mat Object: 1 MPI processes 227 type: seqaij 228 rows=21, cols=21 229 package used to perform factorization: petsc 230 total: nonzeros=85, allocated nonzeros=85 231 not using I-node routines 232 linear system matrix = precond matrix: 233 Mat Object: (fieldsplit_temperature_) 1 MPI processes 234 type: seqaij 235 rows=21, cols=21 236 total: nonzeros=85, allocated nonzeros=85 237 total number of mallocs used during MatSetValues calls=0 238 not using I-node routines 239 linear system matrix = precond matrix: 240 Mat Object: (fieldsplit_temperature_) 3 MPI processes 241 type: mpiaij 242 rows=49, cols=49 243 total: nonzeros=217, allocated nonzeros=217 244 total number of mallocs used during MatSetValues calls=0 245 not using I-node (on process 0) routines 246 linear system matrix followed by preconditioner matrix: 247 Mat Object: 3 MPI processes 248 type: mffd 249 rows=196, cols=196 250 Matrix-free approximation: 251 err=1.49012e-08 (relative error in function evaluation) 252 Using wp compute h routine 253 Does not compute normU 254 Mat Object: 3 MPI processes 255 type: mpiaij 256 rows=196, cols=196, bs=4 257 total: nonzeros=3472, allocated nonzeros=3472 258 total number of mallocs used during MatSetValues calls=0 259Number of SNES iterations = 2 260