1lid velocity = 0.000106281, prandtl # = 1., grashof # = 1. 2 0 SNES Function norm 1.036007954337e-02 3 0 KSP Residual norm 2.388589583549e+00 4 1 KSP Residual norm 5.715829806981e-01 5 2 KSP Residual norm 4.623679005936e-02 6 3 KSP Residual norm 1.143381177646e-02 7 4 KSP Residual norm 2.015139840224e-03 8 5 KSP Residual norm 4.356196119798e-04 9 6 KSP Residual norm 4.240953066710e-05 10 7 KSP Residual norm 8.848315297175e-06 11 1 SNES Function norm 9.854304971115e-06 12 0 KSP Residual norm 3.868049496775e-05 13 1 KSP Residual norm 7.693574326868e-06 14 2 KSP Residual norm 1.059429116239e-06 15 3 KSP Residual norm 4.004524784804e-07 16 4 KSP Residual norm 1.050186948327e-07 17 5 KSP Residual norm 5.073180513583e-08 18 6 KSP Residual norm 2.510513776297e-08 19 7 KSP Residual norm 1.211886495400e-08 20 8 KSP Residual norm 1.911963112131e-09 21 9 KSP Residual norm 3.005260864225e-10 22 2 SNES Function norm 3.117674497824e-10 23 0 KSP Residual norm 3.005042584730e-10 24 1 KSP Residual norm 1.120673922713e-10 25 2 KSP Residual norm 3.288439453292e-11 26 3 KSP Residual norm 5.822504321413e-12 27 4 KSP Residual norm 2.486684466178e-12 28 5 KSP Residual norm 1.198858055503e-12 29 6 KSP Residual norm 6.255669709502e-13 30 7 KSP Residual norm 1.544647758005e-13 31 8 KSP Residual norm 4.592122224907e-14 32 9 KSP Residual norm 4.984149547392e-15 33 10 KSP Residual norm 8.905129652955e-16 34 3 SNES Function norm 1.045594761851e-14 35SNES Object: 4 MPI processes 36 type: newtonls 37 maximum iterations=50, maximum function evaluations=10000 38 tolerances: relative=1e-08, absolute=1e-50, solution=1e-08 39 total number of linear solver iterations=26 40 total number of function evaluations=4 41 norm schedule ALWAYS 42 Jacobian is built using colored finite differences on a DMDA 43 SNESLineSearch Object: 4 MPI processes 44 type: bt 45 interpolation: cubic 46 alpha=1.000000e-04 47 maxstep=1.000000e+08, minlambda=1.000000e-12 48 tolerances: relative=1.000000e-08, absolute=1.000000e-15, lambda=1.000000e-08 49 maximum iterations=40 50 KSP Object: 4 MPI processes 51 type: gmres 52 restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement 53 happy breakdown tolerance 1e-30 54 maximum iterations=10000, initial guess is zero 55 tolerances: relative=1e-05, absolute=1e-50, divergence=10000. 56 left preconditioning 57 using PRECONDITIONED norm type for convergence test 58 PC Object: 4 MPI processes 59 type: mg 60 type is MULTIPLICATIVE, levels=6 cycles=v 61 Cycles per PCApply=1 62 Not using Galerkin computed coarse grid matrices 63 Coarse grid solver -- level 0 ------------------------------- 64 KSP Object: (mg_coarse_) 4 MPI processes 65 type: preonly 66 maximum iterations=10000, initial guess is zero 67 tolerances: relative=1e-05, absolute=1e-50, divergence=10000. 68 left preconditioning 69 using NONE norm type for convergence test 70 PC Object: (mg_coarse_) 4 MPI processes 71 type: redundant 72 First (color=0) of 4 PCs follows 73 KSP Object: (mg_coarse_redundant_) 1 MPI process 74 type: preonly 75 maximum iterations=10000, initial guess is zero 76 tolerances: relative=1e-05, absolute=1e-50, divergence=10000. 77 left preconditioning 78 using NONE norm type for convergence test 79 PC Object: (mg_coarse_redundant_) 1 MPI process 80 type: lu 81 out-of-place factorization 82 tolerance for zero pivot 2.22045e-14 83 using diagonal shift on blocks to prevent zero pivot [INBLOCKS] 84 matrix ordering: nd 85 factor fill ratio given 5., needed 1.875 86 Factored matrix follows: 87 Mat Object: (mg_coarse_redundant_) 1 MPI process 88 type: seqaij 89 rows=64, cols=64, bs=4 90 package used to perform factorization: petsc 91 total: nonzeros=1920, allocated nonzeros=1920 92 using I-node routines: found 16 nodes, limit used is 5 93 linear system matrix = precond matrix: 94 Mat Object: 1 MPI process 95 type: seqaij 96 rows=64, cols=64, bs=4 97 total: nonzeros=1024, allocated nonzeros=1024 98 total number of mallocs used during MatSetValues calls=0 99 using I-node routines: found 16 nodes, limit used is 5 100 linear system matrix = precond matrix: 101 Mat Object: 4 MPI processes 102 type: mpiaij 103 rows=64, cols=64, bs=4 104 total: nonzeros=1024, allocated nonzeros=1024 105 total number of mallocs used during MatSetValues calls=0 106 Down solver (pre-smoother) on level 1 ------------------------------- 107 KSP Object: (mg_levels_1_) 4 MPI processes 108 type: chebyshev 109 Chebyshev polynomial of first kind 110 eigenvalue targets used: min 0.148269, max 1.63095 111 eigenvalues estimated via gmres: min 0.144902, max 1.48269 112 eigenvalues estimated using gmres with transform: [0. 0.1; 0. 1.1] 113 KSP Object: (mg_levels_1_esteig_) 4 MPI processes 114 type: gmres 115 restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement 116 happy breakdown tolerance 1e-30 117 maximum iterations=10, initial guess is zero 118 tolerances: relative=1e-12, absolute=1e-50, divergence=10000. 119 left preconditioning 120 using PRECONDITIONED norm type for convergence test 121 estimating eigenvalues using a noisy random number generated right-hand side 122 maximum iterations=2, nonzero initial guess 123 tolerances: relative=1e-05, absolute=1e-50, divergence=10000. 124 left preconditioning 125 using NONE norm type for convergence test 126 PC Object: (mg_levels_1_) 4 MPI processes 127 type: sor 128 type = local_symmetric, iterations = 1, local iterations = 1, omega = 1. 129 linear system matrix = precond matrix: 130 Mat Object: 4 MPI processes 131 type: mpiaij 132 rows=196, cols=196, bs=4 133 total: nonzeros=3472, allocated nonzeros=3472 134 total number of mallocs used during MatSetValues calls=0 135 Up solver (post-smoother) same as down solver (pre-smoother) 136 Down solver (pre-smoother) on level 2 ------------------------------- 137 KSP Object: (mg_levels_2_) 4 MPI processes 138 type: chebyshev 139 Chebyshev polynomial of first kind 140 eigenvalue targets used: min 0.149178, max 1.64096 141 eigenvalues estimated via gmres: min 0.0843938, max 1.49178 142 eigenvalues estimated using gmres with transform: [0. 0.1; 0. 1.1] 143 KSP Object: (mg_levels_2_esteig_) 4 MPI processes 144 type: gmres 145 restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement 146 happy breakdown tolerance 1e-30 147 maximum iterations=10, initial guess is zero 148 tolerances: relative=1e-12, absolute=1e-50, divergence=10000. 149 left preconditioning 150 using PRECONDITIONED norm type for convergence test 151 estimating eigenvalues using a noisy random number generated right-hand side 152 maximum iterations=2, nonzero initial guess 153 tolerances: relative=1e-05, absolute=1e-50, divergence=10000. 154 left preconditioning 155 using NONE norm type for convergence test 156 PC Object: (mg_levels_2_) 4 MPI processes 157 type: sor 158 type = local_symmetric, iterations = 1, local iterations = 1, omega = 1. 159 linear system matrix = precond matrix: 160 Mat Object: 4 MPI processes 161 type: mpiaij 162 rows=676, cols=676, bs=4 163 total: nonzeros=12688, allocated nonzeros=12688 164 total number of mallocs used during MatSetValues calls=0 165 Up solver (post-smoother) same as down solver (pre-smoother) 166 Down solver (pre-smoother) on level 3 ------------------------------- 167 KSP Object: (mg_levels_3_) 4 MPI processes 168 type: chebyshev 169 Chebyshev polynomial of first kind 170 eigenvalue targets used: min 0.146454, max 1.61099 171 eigenvalues estimated via gmres: min 0.0659153, max 1.46454 172 eigenvalues estimated using gmres with transform: [0. 0.1; 0. 1.1] 173 KSP Object: (mg_levels_3_esteig_) 4 MPI processes 174 type: gmres 175 restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement 176 happy breakdown tolerance 1e-30 177 maximum iterations=10, initial guess is zero 178 tolerances: relative=1e-12, absolute=1e-50, divergence=10000. 179 left preconditioning 180 using PRECONDITIONED norm type for convergence test 181 estimating eigenvalues using a noisy random number generated right-hand side 182 maximum iterations=2, nonzero initial guess 183 tolerances: relative=1e-05, absolute=1e-50, divergence=10000. 184 left preconditioning 185 using NONE norm type for convergence test 186 PC Object: (mg_levels_3_) 4 MPI processes 187 type: sor 188 type = local_symmetric, iterations = 1, local iterations = 1, omega = 1. 189 linear system matrix = precond matrix: 190 Mat Object: 4 MPI processes 191 type: mpiaij 192 rows=2500, cols=2500, bs=4 193 total: nonzeros=48400, allocated nonzeros=48400 194 total number of mallocs used during MatSetValues calls=0 195 Up solver (post-smoother) same as down solver (pre-smoother) 196 Down solver (pre-smoother) on level 4 ------------------------------- 197 KSP Object: (mg_levels_4_) 4 MPI processes 198 type: chebyshev 199 Chebyshev polynomial of first kind 200 eigenvalue targets used: min 0.141089, max 1.55197 201 eigenvalues estimated via gmres: min 0.044097, max 1.41089 202 eigenvalues estimated using gmres with transform: [0. 0.1; 0. 1.1] 203 KSP Object: (mg_levels_4_esteig_) 4 MPI processes 204 type: gmres 205 restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement 206 happy breakdown tolerance 1e-30 207 maximum iterations=10, initial guess is zero 208 tolerances: relative=1e-12, absolute=1e-50, divergence=10000. 209 left preconditioning 210 using PRECONDITIONED norm type for convergence test 211 estimating eigenvalues using a noisy random number generated right-hand side 212 maximum iterations=2, nonzero initial guess 213 tolerances: relative=1e-05, absolute=1e-50, divergence=10000. 214 left preconditioning 215 using NONE norm type for convergence test 216 PC Object: (mg_levels_4_) 4 MPI processes 217 type: sor 218 type = local_symmetric, iterations = 1, local iterations = 1, omega = 1. 219 linear system matrix = precond matrix: 220 Mat Object: 4 MPI processes 221 type: mpiaij 222 rows=9604, cols=9604, bs=4 223 total: nonzeros=188944, allocated nonzeros=188944 224 total number of mallocs used during MatSetValues calls=0 225 Up solver (post-smoother) same as down solver (pre-smoother) 226 Down solver (pre-smoother) on level 5 ------------------------------- 227 KSP Object: (mg_levels_5_) 4 MPI processes 228 type: chebyshev 229 Chebyshev polynomial of first kind 230 eigenvalue targets used: min 0.127956, max 1.40751 231 eigenvalues estimated via gmres: min 0.0380398, max 1.27956 232 eigenvalues estimated using gmres with transform: [0. 0.1; 0. 1.1] 233 KSP Object: (mg_levels_5_esteig_) 4 MPI processes 234 type: gmres 235 restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement 236 happy breakdown tolerance 1e-30 237 maximum iterations=10, initial guess is zero 238 tolerances: relative=1e-12, absolute=1e-50, divergence=10000. 239 left preconditioning 240 using PRECONDITIONED norm type for convergence test 241 estimating eigenvalues using a noisy random number generated right-hand side 242 maximum iterations=2, nonzero initial guess 243 tolerances: relative=1e-05, absolute=1e-50, divergence=10000. 244 left preconditioning 245 using NONE norm type for convergence test 246 PC Object: (mg_levels_5_) 4 MPI processes 247 type: sor 248 type = local_symmetric, iterations = 1, local iterations = 1, omega = 1. 249 linear system matrix = precond matrix: 250 Mat Object: 4 MPI processes 251 type: mpiaij 252 rows=37636, cols=37636, bs=4 253 total: nonzeros=746512, allocated nonzeros=746512 254 total number of mallocs used during MatSetValues calls=0 255 Up solver (post-smoother) same as down solver (pre-smoother) 256 linear system matrix = precond matrix: 257 Mat Object: 4 MPI processes 258 type: mpiaij 259 rows=37636, cols=37636, bs=4 260 total: nonzeros=746512, allocated nonzeros=746512 261 total number of mallocs used during MatSetValues calls=0 262Number of SNES iterations = 3 263