1lid velocity = 0.000106281, prandtl # = 1., grashof # = 1. 2 0 SNES Function norm 1.036007954337e-02 3 0 KSP Residual norm 4.386611885521e-01 4 1 KSP Residual norm 3.127838547648e-01 5 2 KSP Residual norm 2.334238799976e-01 6 3 KSP Residual norm 1.318350769046e-01 7 4 KSP Residual norm 8.143004995277e-02 8 5 KSP Residual norm 4.813457707703e-02 9 6 KSP Residual norm 2.733719484710e-02 10 7 KSP Residual norm 1.462151166610e-02 11 8 KSP Residual norm 8.331275159714e-03 12 9 KSP Residual norm 5.275884606772e-03 13 10 KSP Residual norm 3.768690746652e-03 14 11 KSP Residual norm 2.699365250849e-03 15 12 KSP Residual norm 1.824681696989e-03 16 13 KSP Residual norm 1.179126381375e-03 17 14 KSP Residual norm 7.107306905422e-04 18 15 KSP Residual norm 4.155820345074e-04 19 16 KSP Residual norm 2.408438645714e-04 20 17 KSP Residual norm 1.476855736514e-04 21 18 KSP Residual norm 9.907578845784e-05 22 19 KSP Residual norm 6.797328938557e-05 23 20 KSP Residual norm 4.675944545290e-05 24 21 KSP Residual norm 3.352583819988e-05 25 22 KSP Residual norm 2.441653843583e-05 26 23 KSP Residual norm 1.768579774951e-05 27 24 KSP Residual norm 1.411236998154e-05 28 25 KSP Residual norm 1.062751951795e-05 29 26 KSP Residual norm 7.402652132828e-06 30 27 KSP Residual norm 5.376101942565e-06 31 28 KSP Residual norm 3.936750059583e-06 32 1 SNES Function norm 1.539897001872e-06 33 0 KSP Residual norm 1.219799734452e-05 34 1 KSP Residual norm 9.351061917402e-06 35 2 KSP Residual norm 5.158749230924e-06 36 3 KSP Residual norm 3.800583255559e-06 37 4 KSP Residual norm 3.358178466662e-06 38 5 KSP Residual norm 2.704711328549e-06 39 6 KSP Residual norm 1.965512422061e-06 40 7 KSP Residual norm 1.383231224585e-06 41 8 KSP Residual norm 1.016838426400e-06 42 9 KSP Residual norm 7.817076853999e-07 43 10 KSP Residual norm 5.714787985538e-07 44 11 KSP Residual norm 4.070517298179e-07 45 12 KSP Residual norm 3.199366855723e-07 46 13 KSP Residual norm 2.498338700906e-07 47 14 KSP Residual norm 1.774654138694e-07 48 15 KSP Residual norm 1.220894443119e-07 49 16 KSP Residual norm 8.191467113032e-08 50 17 KSP Residual norm 5.471425416521e-08 51 18 KSP Residual norm 3.800399846106e-08 52 19 KSP Residual norm 2.624332673801e-08 53 20 KSP Residual norm 1.817665454829e-08 54 21 KSP Residual norm 1.193934944948e-08 55 22 KSP Residual norm 7.058845517318e-09 56 23 KSP Residual norm 4.801132366324e-09 57 24 KSP Residual norm 3.136888134931e-09 58 25 KSP Residual norm 1.908878704398e-09 59 26 KSP Residual norm 1.179297203589e-09 60 27 KSP Residual norm 7.450436712605e-10 61 28 KSP Residual norm 4.727201847673e-10 62 29 KSP Residual norm 3.071103410289e-10 63 30 KSP Residual norm 1.954558918512e-10 64 31 KSP Residual norm 1.476035656576e-10 65 32 KSP Residual norm 9.544641423593e-11 66 2 SNES Function norm 2.155536219955e-11 67SNES Object: 4 MPI processes 68 type: newtonls 69 maximum iterations=50, maximum function evaluations=10000 70 tolerances: relative=1e-08, absolute=1e-50, solution=1e-08 71 total number of linear solver iterations=60 72 total number of function evaluations=3 73 norm schedule ALWAYS 74 Jacobian is built using colored finite differences on a DMDA 75 SNESLineSearch Object: 4 MPI processes 76 type: bt 77 interpolation: cubic 78 alpha=1.000000e-04 79 maxstep=1.000000e+08, minlambda=1.000000e-12 80 tolerances: relative=1.000000e-08, absolute=1.000000e-15, lambda=1.000000e-08 81 maximum iterations=40 82 KSP Object: 4 MPI processes 83 type: gmres 84 restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement 85 happy breakdown tolerance 1e-30 86 maximum iterations=10000, initial guess is zero 87 tolerances: relative=1e-05, absolute=1e-50, divergence=10000. 88 left preconditioning 89 using PRECONDITIONED norm type for convergence test 90 PC Object: 4 MPI processes 91 type: ml 92 type is MULTIPLICATIVE, levels=5 cycles=v 93 Cycles per PCApply=1 94 Using externally compute Galerkin coarse grid matrices 95 Coarse grid solver -- level ------------------------------- 96 KSP Object: (mg_coarse_) 4 MPI processes 97 type: preonly 98 maximum iterations=10000, initial guess is zero 99 tolerances: relative=1e-05, absolute=1e-50, divergence=10000. 100 left preconditioning 101 using NONE norm type for convergence test 102 PC Object: (mg_coarse_) 4 MPI processes 103 type: redundant 104 First (color=0) of 4 PCs follows 105 KSP Object: (mg_coarse_redundant_) 1 MPI processes 106 type: preonly 107 maximum iterations=10000, initial guess is zero 108 tolerances: relative=1e-05, absolute=1e-50, divergence=10000. 109 left preconditioning 110 using NONE norm type for convergence test 111 PC Object: (mg_coarse_redundant_) 1 MPI processes 112 type: lu 113 out-of-place factorization 114 tolerance for zero pivot 2.22045e-14 115 using diagonal shift on blocks to prevent zero pivot [INBLOCKS] 116 matrix ordering: nd 117 factor fill ratio given 5., needed 1. 118 Factored matrix follows: 119 Mat Object: 1 MPI processes 120 type: seqaij 121 rows=16, cols=16, bs=4 122 package used to perform factorization: petsc 123 total: nonzeros=256, allocated nonzeros=256 124 total number of mallocs used during MatSetValues calls=0 125 using I-node routines: found 4 nodes, limit used is 5 126 linear system matrix = precond matrix: 127 Mat Object: 1 MPI processes 128 type: seqaij 129 rows=16, cols=16, bs=4 130 total: nonzeros=256, allocated nonzeros=256 131 total number of mallocs used during MatSetValues calls=0 132 using I-node routines: found 4 nodes, limit used is 5 133 linear system matrix = precond matrix: 134 Mat Object: 4 MPI processes 135 type: mpiaij 136 rows=16, cols=16, bs=4 137 total: nonzeros=256, allocated nonzeros=256 138 total number of mallocs used during MatSetValues calls=0 139 using I-node (on process 0) routines: found 1 nodes, limit used is 5 140 Down solver (pre-smoother) on level 1 ------------------------------- 141 KSP Object: (mg_levels_1_) 4 MPI processes 142 type: richardson 143 damping factor=1. 144 maximum iterations=2, nonzero initial guess 145 tolerances: relative=1e-05, absolute=1e-50, divergence=10000. 146 left preconditioning 147 using NONE norm type for convergence test 148 PC Object: (mg_levels_1_) 4 MPI processes 149 type: sor 150 type = local_symmetric, iterations = 1, local iterations = 1, omega = 1. 151 linear system matrix = precond matrix: 152 Mat Object: 4 MPI processes 153 type: mpiaij 154 rows=88, cols=88, bs=4 155 total: nonzeros=4896, allocated nonzeros=4896 156 total number of mallocs used during MatSetValues calls=0 157 using I-node (on process 0) routines: found 5 nodes, limit used is 5 158 Up solver (post-smoother) same as down solver (pre-smoother) 159 Down solver (pre-smoother) on level 2 ------------------------------- 160 KSP Object: (mg_levels_2_) 4 MPI processes 161 type: richardson 162 damping factor=1. 163 maximum iterations=2, nonzero initial guess 164 tolerances: relative=1e-05, absolute=1e-50, divergence=10000. 165 left preconditioning 166 using NONE norm type for convergence test 167 PC Object: (mg_levels_2_) 4 MPI processes 168 type: sor 169 type = local_symmetric, iterations = 1, local iterations = 1, omega = 1. 170 linear system matrix = precond matrix: 171 Mat Object: 4 MPI processes 172 type: mpiaij 173 rows=772, cols=772, bs=4 174 total: nonzeros=35364, allocated nonzeros=35364 175 total number of mallocs used during MatSetValues calls=0 176 using I-node (on process 0) routines: found 56 nodes, limit used is 5 177 Up solver (post-smoother) same as down solver (pre-smoother) 178 Down solver (pre-smoother) on level 3 ------------------------------- 179 KSP Object: (mg_levels_3_) 4 MPI processes 180 type: richardson 181 damping factor=1. 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=7156, cols=7156, bs=4 193 total: nonzeros=258241, allocated nonzeros=258241 194 total number of mallocs used during MatSetValues calls=0 195 not using I-node (on process 0) routines 196 Up solver (post-smoother) same as down solver (pre-smoother) 197 Down solver (pre-smoother) on level 4 ------------------------------- 198 KSP Object: (mg_levels_4_) 4 MPI processes 199 type: richardson 200 damping factor=1. 201 maximum iterations=2, nonzero initial guess 202 tolerances: relative=1e-05, absolute=1e-50, divergence=10000. 203 left preconditioning 204 using NONE norm type for convergence test 205 PC Object: (mg_levels_4_) 4 MPI processes 206 type: sor 207 type = local_symmetric, iterations = 1, local iterations = 1, omega = 1. 208 linear system matrix = precond matrix: 209 Mat Object: 4 MPI processes 210 type: mpiaij 211 rows=37636, cols=37636, bs=4 212 total: nonzeros=746512, allocated nonzeros=746512 213 total number of mallocs used during MatSetValues calls=0 214 Up solver (post-smoother) same as down solver (pre-smoother) 215 linear system matrix = precond matrix: 216 Mat Object: 4 MPI processes 217 type: mpiaij 218 rows=37636, cols=37636, bs=4 219 total: nonzeros=746512, allocated nonzeros=746512 220 total number of mallocs used during MatSetValues calls=0 221Number of SNES iterations = 2 222