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