1lid velocity = 0.000106281, prandtl # = 1., grashof # = 1. 2 0 SNES Function norm 1.036007954337e-02 3 0 KSP Residual norm 4.888657920112e-01 4 1 KSP Residual norm 3.141627616286e-01 5 2 KSP Residual norm 2.358812378525e-01 6 3 KSP Residual norm 1.259045437573e-01 7 4 KSP Residual norm 7.408575133258e-02 8 5 KSP Residual norm 4.159658325629e-02 9 6 KSP Residual norm 2.181861584091e-02 10 7 KSP Residual norm 1.076244407256e-02 11 8 KSP Residual norm 5.242924144023e-03 12 9 KSP Residual norm 2.588726461570e-03 13 10 KSP Residual norm 1.400700093556e-03 14 11 KSP Residual norm 7.759158101290e-04 15 12 KSP Residual norm 4.655744671687e-04 16 13 KSP Residual norm 3.244240900196e-04 17 14 KSP Residual norm 2.327144448481e-04 18 15 KSP Residual norm 1.500951563875e-04 19 16 KSP Residual norm 8.327955132177e-05 20 17 KSP Residual norm 4.672468799190e-05 21 18 KSP Residual norm 2.889321996959e-05 22 19 KSP Residual norm 1.903514112532e-05 23 20 KSP Residual norm 1.298123631775e-05 24 21 KSP Residual norm 8.775705379722e-06 25 22 KSP Residual norm 5.330175482523e-06 26 23 KSP Residual norm 2.869775929854e-06 27 1 SNES Function norm 1.370527153606e-06 28 0 KSP Residual norm 1.177060567864e-05 29 1 KSP Residual norm 6.928128579778e-06 30 2 KSP Residual norm 4.304601930959e-06 31 3 KSP Residual norm 2.454875323346e-06 32 4 KSP Residual norm 1.512790232898e-06 33 5 KSP Residual norm 1.001915525185e-06 34 6 KSP Residual norm 7.387623107528e-07 35 7 KSP Residual norm 4.953711915641e-07 36 8 KSP Residual norm 3.144821641519e-07 37 9 KSP Residual norm 2.134819078071e-07 38 10 KSP Residual norm 1.518551249827e-07 39 11 KSP Residual norm 1.014185303569e-07 40 12 KSP Residual norm 5.608673008996e-08 41 13 KSP Residual norm 3.030519067648e-08 42 14 KSP Residual norm 1.786947998191e-08 43 15 KSP Residual norm 1.177816744649e-08 44 16 KSP Residual norm 8.997335866987e-09 45 17 KSP Residual norm 6.356106556680e-09 46 18 KSP Residual norm 4.147037399694e-09 47 19 KSP Residual norm 2.711154831563e-09 48 20 KSP Residual norm 1.750865415567e-09 49 21 KSP Residual norm 1.127681613823e-09 50 22 KSP Residual norm 7.783034009647e-10 51 23 KSP Residual norm 5.047971816405e-10 52 24 KSP Residual norm 3.137090395920e-10 53 25 KSP Residual norm 2.182629803057e-10 54 26 KSP Residual norm 1.265836293994e-10 55 27 KSP Residual norm 7.756916170515e-11 56 2 SNES Function norm 2.339734806142e-11 57SNES Object: 4 MPI processes 58 type: newtonls 59 maximum iterations=50, maximum function evaluations=10000 60 tolerances: relative=1e-08, absolute=1e-50, solution=1e-08 61 total number of linear solver iterations=50 62 total number of function evaluations=3 63 norm schedule ALWAYS 64 Jacobian is built using colored finite differences on a DMDA 65 SNESLineSearch Object: 4 MPI processes 66 type: bt 67 interpolation: cubic 68 alpha=1.000000e-04 69 maxstep=1.000000e+08, minlambda=1.000000e-12 70 tolerances: relative=1.000000e-08, absolute=1.000000e-15, lambda=1.000000e-08 71 maximum iterations=40 72 KSP Object: 4 MPI processes 73 type: gmres 74 restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement 75 happy breakdown tolerance 1e-30 76 maximum iterations=10000, initial guess is zero 77 tolerances: relative=1e-05, absolute=1e-50, divergence=10000. 78 left preconditioning 79 using PRECONDITIONED norm type for convergence test 80 PC Object: 4 MPI processes 81 type: ml 82 type is MULTIPLICATIVE, levels=5 cycles=v 83 Cycles per PCApply=1 84 Using externally compute Galerkin coarse grid matrices 85 Coarse grid solver -- level 0 ------------------------------- 86 KSP Object: (mg_coarse_) 4 MPI processes 87 type: preonly 88 maximum iterations=10000, initial guess is zero 89 tolerances: relative=1e-05, absolute=1e-50, divergence=10000. 90 left preconditioning 91 using NONE norm type for convergence test 92 PC Object: (mg_coarse_) 4 MPI processes 93 type: redundant 94 First (color=0) of 4 PCs follows 95 KSP Object: (mg_coarse_redundant_) 1 MPI process 96 type: preonly 97 maximum iterations=10000, initial guess is zero 98 tolerances: relative=1e-05, absolute=1e-50, divergence=10000. 99 left preconditioning 100 using NONE norm type for convergence test 101 PC Object: (mg_coarse_redundant_) 1 MPI process 102 type: lu 103 out-of-place factorization 104 tolerance for zero pivot 2.22045e-14 105 using diagonal shift on blocks to prevent zero pivot [INBLOCKS] 106 matrix ordering: nd 107 factor fill ratio given 5., needed 1. 108 Factored matrix follows: 109 Mat Object: (mg_coarse_redundant_) 1 MPI process 110 type: seqaij 111 rows=16, cols=16, bs=4 112 package used to perform factorization: petsc 113 total: nonzeros=256, allocated nonzeros=256 114 using I-node routines: found 4 nodes, limit used is 5 115 linear system matrix = precond matrix: 116 Mat Object: 1 MPI process 117 type: seqaij 118 rows=16, cols=16, bs=4 119 total: nonzeros=256, allocated nonzeros=256 120 total number of mallocs used during MatSetValues calls=0 121 using I-node routines: found 4 nodes, limit used is 5 122 linear system matrix = precond matrix: 123 Mat Object: 4 MPI processes 124 type: mpiaij 125 rows=16, cols=16, bs=4 126 total: nonzeros=256, allocated nonzeros=256 127 total number of mallocs used during MatSetValues calls=0 128 using I-node (on process 0) routines: found 1 nodes, limit used is 5 129 Down solver (pre-smoother) on level 1 ------------------------------- 130 KSP Object: (mg_levels_1_) 4 MPI processes 131 type: richardson 132 damping factor=1. 133 maximum iterations=2, nonzero initial guess 134 tolerances: relative=1e-05, absolute=1e-50, divergence=10000. 135 left preconditioning 136 using NONE norm type for convergence test 137 PC Object: (mg_levels_1_) 4 MPI processes 138 type: sor 139 type = local_symmetric, iterations = 1, local iterations = 1, omega = 1. 140 linear system matrix = precond matrix: 141 Mat Object: 4 MPI processes 142 type: mpiaij 143 rows=88, cols=88, bs=4 144 total: nonzeros=4960, allocated nonzeros=4960 145 total number of mallocs used during MatSetValues calls=0 146 using I-node (on process 0) routines: found 4 nodes, limit used is 5 147 Up solver (post-smoother) same as down solver (pre-smoother) 148 Down solver (pre-smoother) on level 2 ------------------------------- 149 KSP Object: (mg_levels_2_) 4 MPI processes 150 type: richardson 151 damping factor=1. 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=764, cols=764, bs=4 163 total: nonzeros=35672, allocated nonzeros=35672 164 total number of mallocs used during MatSetValues calls=0 165 using I-node (on process 0) routines: found 50 nodes, limit used is 5 166 Up solver (post-smoother) same as down solver (pre-smoother) 167 Down solver (pre-smoother) on level 3 ------------------------------- 168 KSP Object: (mg_levels_3_) 4 MPI processes 169 type: richardson 170 damping factor=1. 171 maximum iterations=2, nonzero initial guess 172 tolerances: relative=1e-05, absolute=1e-50, divergence=10000. 173 left preconditioning 174 using NONE norm type for convergence test 175 PC Object: (mg_levels_3_) 4 MPI processes 176 type: sor 177 type = local_symmetric, iterations = 1, local iterations = 1, omega = 1. 178 linear system matrix = precond matrix: 179 Mat Object: 4 MPI processes 180 type: mpiaij 181 rows=7156, cols=7156, bs=4 182 total: nonzeros=259236, allocated nonzeros=259236 183 total number of mallocs used during MatSetValues calls=0 184 not using I-node (on process 0) routines 185 Up solver (post-smoother) same as down solver (pre-smoother) 186 Down solver (pre-smoother) on level 4 ------------------------------- 187 KSP Object: (mg_levels_4_) 4 MPI processes 188 type: richardson 189 damping factor=1. 190 maximum iterations=2, nonzero initial guess 191 tolerances: relative=1e-05, absolute=1e-50, divergence=10000. 192 left preconditioning 193 using NONE norm type for convergence test 194 PC Object: (mg_levels_4_) 4 MPI processes 195 type: sor 196 type = local_symmetric, iterations = 1, local iterations = 1, omega = 1. 197 linear system matrix = precond matrix: 198 Mat Object: 4 MPI processes 199 type: mpiaij 200 rows=37636, cols=37636, bs=4 201 total: nonzeros=746512, allocated nonzeros=746512 202 total number of mallocs used during MatSetValues calls=0 203 Up solver (post-smoother) same as down solver (pre-smoother) 204 linear system matrix = precond matrix: 205 Mat Object: 4 MPI processes 206 type: mpiaij 207 rows=37636, cols=37636, bs=4 208 total: nonzeros=746512, allocated nonzeros=746512 209 total number of mallocs used during MatSetValues calls=0 210Number of SNES iterations = 2 211