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 0 ------------------------------- 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 process 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 process 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: (mg_coarse_redundant_) 1 MPI process 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 using I-node routines: found 4 nodes, limit used is 5 125 linear system matrix = precond matrix: 126 Mat Object: 1 MPI process 127 type: seqaij 128 rows=16, cols=16, bs=4 129 total: nonzeros=256, allocated nonzeros=256 130 total number of mallocs used during MatSetValues calls=0 131 using I-node routines: found 4 nodes, limit used is 5 132 linear system matrix = precond matrix: 133 Mat Object: 4 MPI processes 134 type: mpiaij 135 rows=16, cols=16, bs=4 136 total: nonzeros=256, allocated nonzeros=256 137 total number of mallocs used during MatSetValues calls=0 138 using I-node (on process 0) routines: found 1 nodes, limit used is 5 139 Down solver (pre-smoother) on level 1 ------------------------------- 140 KSP Object: (mg_levels_1_) 4 MPI processes 141 type: richardson 142 damping factor=1. 143 maximum iterations=2, nonzero initial guess 144 tolerances: relative=1e-05, absolute=1e-50, divergence=10000. 145 left preconditioning 146 using NONE norm type for convergence test 147 PC Object: (mg_levels_1_) 4 MPI processes 148 type: sor 149 type = local_symmetric, iterations = 1, local iterations = 1, omega = 1. 150 linear system matrix = precond matrix: 151 Mat Object: 4 MPI processes 152 type: mpiaij 153 rows=88, cols=88, bs=4 154 total: nonzeros=4896, allocated nonzeros=4896 155 total number of mallocs used during MatSetValues calls=0 156 using I-node (on process 0) routines: found 5 nodes, limit used is 5 157 Up solver (post-smoother) same as down solver (pre-smoother) 158 Down solver (pre-smoother) on level 2 ------------------------------- 159 KSP Object: (mg_levels_2_) 4 MPI processes 160 type: richardson 161 damping factor=1. 162 maximum iterations=2, nonzero initial guess 163 tolerances: relative=1e-05, absolute=1e-50, divergence=10000. 164 left preconditioning 165 using NONE norm type for convergence test 166 PC Object: (mg_levels_2_) 4 MPI processes 167 type: sor 168 type = local_symmetric, iterations = 1, local iterations = 1, omega = 1. 169 linear system matrix = precond matrix: 170 Mat Object: 4 MPI processes 171 type: mpiaij 172 rows=772, cols=772, bs=4 173 total: nonzeros=35364, allocated nonzeros=35364 174 total number of mallocs used during MatSetValues calls=0 175 using I-node (on process 0) routines: found 56 nodes, limit used is 5 176 Up solver (post-smoother) same as down solver (pre-smoother) 177 Down solver (pre-smoother) on level 3 ------------------------------- 178 KSP Object: (mg_levels_3_) 4 MPI processes 179 type: richardson 180 damping factor=1. 181 maximum iterations=2, nonzero initial guess 182 tolerances: relative=1e-05, absolute=1e-50, divergence=10000. 183 left preconditioning 184 using NONE norm type for convergence test 185 PC Object: (mg_levels_3_) 4 MPI processes 186 type: sor 187 type = local_symmetric, iterations = 1, local iterations = 1, omega = 1. 188 linear system matrix = precond matrix: 189 Mat Object: 4 MPI processes 190 type: mpiaij 191 rows=7156, cols=7156, bs=4 192 total: nonzeros=258241, allocated nonzeros=258241 193 total number of mallocs used during MatSetValues calls=0 194 not using I-node (on process 0) routines 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: richardson 199 damping factor=1. 200 maximum iterations=2, nonzero initial guess 201 tolerances: relative=1e-05, absolute=1e-50, divergence=10000. 202 left preconditioning 203 using NONE norm type for convergence test 204 PC Object: (mg_levels_4_) 4 MPI processes 205 type: sor 206 type = local_symmetric, iterations = 1, local iterations = 1, omega = 1. 207 linear system matrix = precond matrix: 208 Mat Object: 4 MPI processes 209 type: mpiaij 210 rows=37636, cols=37636, bs=4 211 total: nonzeros=746512, allocated nonzeros=746512 212 total number of mallocs used during MatSetValues calls=0 213 Up solver (post-smoother) same as down solver (pre-smoother) 214 linear system matrix = precond matrix: 215 Mat Object: 4 MPI processes 216 type: mpiaij 217 rows=37636, cols=37636, bs=4 218 total: nonzeros=746512, allocated nonzeros=746512 219 total number of mallocs used during MatSetValues calls=0 220Number of SNES iterations = 2 221