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 43SNESLineSearch 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 ------------------------------- 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 processes 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 processes 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: 1 MPI processes 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 total number of mallocs used during MatSetValues calls=0 93 using I-node routines: found 16 nodes, limit used is 5 94 linear system matrix = precond matrix: 95 Mat Object: 1 MPI processes 96 type: seqaij 97 rows=64, cols=64, bs=4 98 total: nonzeros=1024, allocated nonzeros=1024 99 total number of mallocs used during MatSetValues calls=0 100 using I-node routines: found 16 nodes, limit used is 5 101 linear system matrix = precond matrix: 102 Mat Object: 4 MPI processes 103 type: mpiaij 104 rows=64, cols=64, bs=4 105 total: nonzeros=1024, allocated nonzeros=1024 106 total number of mallocs used during MatSetValues calls=0 107 Down solver (pre-smoother) on level 1 ------------------------------- 108 KSP Object: (mg_levels_1_) 4 MPI processes 109 type: chebyshev 110 eigenvalue estimates used: min = 0.148269, max = 1.63095 111 eigenvalues estimate via gmres min 0.144902, max 1.48269 112 eigenvalues estimated using gmres with translations [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 noisy 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 eigenvalue estimates used: min = 0.149178, max = 1.64096 140 eigenvalues estimate via gmres min 0.0843938, max 1.49178 141 eigenvalues estimated using gmres with translations [0. 0.1; 0. 1.1] 142 KSP Object: (mg_levels_2_esteig_) 4 MPI processes 143 type: gmres 144 restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement 145 happy breakdown tolerance 1e-30 146 maximum iterations=10, initial guess is zero 147 tolerances: relative=1e-12, absolute=1e-50, divergence=10000. 148 left preconditioning 149 using PRECONDITIONED norm type for convergence test 150 estimating eigenvalues using noisy right hand side 151 maximum iterations=2, nonzero initial guess 152 tolerances: relative=1e-05, absolute=1e-50, divergence=10000. 153 left preconditioning 154 using NONE norm type for convergence test 155 PC Object: (mg_levels_2_) 4 MPI processes 156 type: sor 157 type = local_symmetric, iterations = 1, local iterations = 1, omega = 1. 158 linear system matrix = precond matrix: 159 Mat Object: 4 MPI processes 160 type: mpiaij 161 rows=676, cols=676, bs=4 162 total: nonzeros=12688, allocated nonzeros=12688 163 total number of mallocs used during MatSetValues calls=0 164 Up solver (post-smoother) same as down solver (pre-smoother) 165 Down solver (pre-smoother) on level 3 ------------------------------- 166 KSP Object: (mg_levels_3_) 4 MPI processes 167 type: chebyshev 168 eigenvalue estimates used: min = 0.146454, max = 1.61099 169 eigenvalues estimate via gmres min 0.0659153, max 1.46454 170 eigenvalues estimated using gmres with translations [0. 0.1; 0. 1.1] 171 KSP Object: (mg_levels_3_esteig_) 4 MPI processes 172 type: gmres 173 restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement 174 happy breakdown tolerance 1e-30 175 maximum iterations=10, initial guess is zero 176 tolerances: relative=1e-12, absolute=1e-50, divergence=10000. 177 left preconditioning 178 using PRECONDITIONED norm type for convergence test 179 estimating eigenvalues using noisy right hand side 180 maximum iterations=2, nonzero initial guess 181 tolerances: relative=1e-05, absolute=1e-50, divergence=10000. 182 left preconditioning 183 using NONE norm type for convergence test 184 PC Object: (mg_levels_3_) 4 MPI processes 185 type: sor 186 type = local_symmetric, iterations = 1, local iterations = 1, omega = 1. 187 linear system matrix = precond matrix: 188 Mat Object: 4 MPI processes 189 type: mpiaij 190 rows=2500, cols=2500, bs=4 191 total: nonzeros=48400, allocated nonzeros=48400 192 total number of mallocs used during MatSetValues calls=0 193 Up solver (post-smoother) same as down solver (pre-smoother) 194 Down solver (pre-smoother) on level 4 ------------------------------- 195 KSP Object: (mg_levels_4_) 4 MPI processes 196 type: chebyshev 197 eigenvalue estimates used: min = 0.141089, max = 1.55197 198 eigenvalues estimate via gmres min 0.044097, max 1.41089 199 eigenvalues estimated using gmres with translations [0. 0.1; 0. 1.1] 200 KSP Object: (mg_levels_4_esteig_) 4 MPI processes 201 type: gmres 202 restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement 203 happy breakdown tolerance 1e-30 204 maximum iterations=10, initial guess is zero 205 tolerances: relative=1e-12, absolute=1e-50, divergence=10000. 206 left preconditioning 207 using PRECONDITIONED norm type for convergence test 208 estimating eigenvalues using noisy right hand side 209 maximum iterations=2, nonzero initial guess 210 tolerances: relative=1e-05, absolute=1e-50, divergence=10000. 211 left preconditioning 212 using NONE norm type for convergence test 213 PC Object: (mg_levels_4_) 4 MPI processes 214 type: sor 215 type = local_symmetric, iterations = 1, local iterations = 1, omega = 1. 216 linear system matrix = precond matrix: 217 Mat Object: 4 MPI processes 218 type: mpiaij 219 rows=9604, cols=9604, bs=4 220 total: nonzeros=188944, allocated nonzeros=188944 221 total number of mallocs used during MatSetValues calls=0 222 Up solver (post-smoother) same as down solver (pre-smoother) 223 Down solver (pre-smoother) on level 5 ------------------------------- 224 KSP Object: (mg_levels_5_) 4 MPI processes 225 type: chebyshev 226 eigenvalue estimates used: min = 0.127956, max = 1.40751 227 eigenvalues estimate via gmres min 0.0380398, max 1.27956 228 eigenvalues estimated using gmres with translations [0. 0.1; 0. 1.1] 229 KSP Object: (mg_levels_5_esteig_) 4 MPI processes 230 type: gmres 231 restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement 232 happy breakdown tolerance 1e-30 233 maximum iterations=10, initial guess is zero 234 tolerances: relative=1e-12, absolute=1e-50, divergence=10000. 235 left preconditioning 236 using PRECONDITIONED norm type for convergence test 237 estimating eigenvalues using noisy right hand side 238 maximum iterations=2, nonzero initial guess 239 tolerances: relative=1e-05, absolute=1e-50, divergence=10000. 240 left preconditioning 241 using NONE norm type for convergence test 242 PC Object: (mg_levels_5_) 4 MPI processes 243 type: sor 244 type = local_symmetric, iterations = 1, local iterations = 1, omega = 1. 245 linear system matrix = precond matrix: 246 Mat Object: 4 MPI processes 247 type: mpiaij 248 rows=37636, cols=37636, bs=4 249 total: nonzeros=746512, allocated nonzeros=746512 250 total number of mallocs used during MatSetValues calls=0 251 Up solver (post-smoother) same as down solver (pre-smoother) 252 linear system matrix = precond matrix: 253 Mat Object: 4 MPI processes 254 type: mpiaij 255 rows=37636, cols=37636, bs=4 256 total: nonzeros=746512, allocated nonzeros=746512 257 total number of mallocs used during MatSetValues calls=0 258Number of SNES iterations = 3 259