1lid velocity = 0.000106281, prandtl # = 1., grashof # = 1. 2 0 SNES Function norm 1.036007954337e-02 3 0 KSP Residual norm 4.888657920255e-01 4 1 KSP Residual norm 3.141627615859e-01 5 2 KSP Residual norm 2.358812378409e-01 6 3 KSP Residual norm 1.259045437422e-01 7 4 KSP Residual norm 7.408575133770e-02 8 5 KSP Residual norm 4.159658326097e-02 9 6 KSP Residual norm 2.181861584671e-02 10 7 KSP Residual norm 1.076244408006e-02 11 8 KSP Residual norm 5.242924152904e-03 12 9 KSP Residual norm 2.588726470721e-03 13 10 KSP Residual norm 1.400700101742e-03 14 11 KSP Residual norm 7.759158144031e-04 15 12 KSP Residual norm 4.655744672970e-04 16 13 KSP Residual norm 3.244240883869e-04 17 14 KSP Residual norm 2.327144434323e-04 18 15 KSP Residual norm 1.500951557724e-04 19 16 KSP Residual norm 8.327955099929e-05 20 17 KSP Residual norm 4.672468779586e-05 21 18 KSP Residual norm 2.889321997942e-05 22 19 KSP Residual norm 1.903514126270e-05 23 20 KSP Residual norm 1.298123639853e-05 24 21 KSP Residual norm 8.775705332110e-06 25 22 KSP Residual norm 5.330175356427e-06 26 23 KSP Residual norm 2.869775815387e-06 27 1 SNES Function norm 1.370527134130e-06 28 0 KSP Residual norm 1.177060566542e-05 29 1 KSP Residual norm 6.928128565305e-06 30 2 KSP Residual norm 4.304601882686e-06 31 3 KSP Residual norm 2.454875297598e-06 32 4 KSP Residual norm 1.512790218693e-06 33 5 KSP Residual norm 1.001915521115e-06 34 6 KSP Residual norm 7.387623119646e-07 35 7 KSP Residual norm 4.953711923012e-07 36 8 KSP Residual norm 3.144821617223e-07 37 9 KSP Residual norm 2.134819050748e-07 38 10 KSP Residual norm 1.518551228740e-07 39 11 KSP Residual norm 1.014185286269e-07 40 12 KSP Residual norm 5.608672865319e-08 41 13 KSP Residual norm 3.030518934915e-08 42 14 KSP Residual norm 1.786947902120e-08 43 15 KSP Residual norm 1.177816691797e-08 44 16 KSP Residual norm 8.997335693730e-09 45 17 KSP Residual norm 6.356106653561e-09 46 18 KSP Residual norm 4.147037531401e-09 47 19 KSP Residual norm 2.711154844415e-09 48 20 KSP Residual norm 1.750865364654e-09 49 21 KSP Residual norm 1.127681551733e-09 50 22 KSP Residual norm 7.783033361682e-10 51 23 KSP Residual norm 5.047971405840e-10 52 24 KSP Residual norm 3.137090060579e-10 53 25 KSP Residual norm 2.182629537231e-10 54 26 KSP Residual norm 1.265836094717e-10 55 27 KSP Residual norm 7.756914910531e-11 56 2 SNES Function norm 2.339735156000e-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 ------------------------------- 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 processes 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 processes 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: 1 MPI processes 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 total number of mallocs used during MatSetValues calls=0 115 using I-node routines: found 4 nodes, limit used is 5 116 linear system matrix = precond matrix: 117 Mat Object: 1 MPI processes 118 type: seqaij 119 rows=16, cols=16, bs=4 120 total: nonzeros=256, allocated nonzeros=256 121 total number of mallocs used during MatSetValues calls=0 122 using I-node routines: found 4 nodes, limit used is 5 123 linear system matrix = precond matrix: 124 Mat Object: 4 MPI processes 125 type: mpiaij 126 rows=16, cols=16, bs=4 127 total: nonzeros=256, allocated nonzeros=256 128 total number of mallocs used during MatSetValues calls=0 129 using I-node (on process 0) routines: found 1 nodes, limit used is 5 130 Down solver (pre-smoother) on level 1 ------------------------------- 131 KSP Object: (mg_levels_1_) 4 MPI processes 132 type: richardson 133 damping factor=1. 134 maximum iterations=2, nonzero initial guess 135 tolerances: relative=1e-05, absolute=1e-50, divergence=10000. 136 left preconditioning 137 using NONE norm type for convergence test 138 PC Object: (mg_levels_1_) 4 MPI processes 139 type: sor 140 type = local_symmetric, iterations = 1, local iterations = 1, omega = 1. 141 linear system matrix = precond matrix: 142 Mat Object: 4 MPI processes 143 type: mpiaij 144 rows=88, cols=88, bs=4 145 total: nonzeros=4960, allocated nonzeros=4960 146 total number of mallocs used during MatSetValues calls=0 147 using I-node (on process 0) routines: found 4 nodes, limit used is 5 148 Up solver (post-smoother) same as down solver (pre-smoother) 149 Down solver (pre-smoother) on level 2 ------------------------------- 150 KSP Object: (mg_levels_2_) 4 MPI processes 151 type: richardson 152 damping factor=1. 153 maximum iterations=2, nonzero initial guess 154 tolerances: relative=1e-05, absolute=1e-50, divergence=10000. 155 left preconditioning 156 using NONE norm type for convergence test 157 PC Object: (mg_levels_2_) 4 MPI processes 158 type: sor 159 type = local_symmetric, iterations = 1, local iterations = 1, omega = 1. 160 linear system matrix = precond matrix: 161 Mat Object: 4 MPI processes 162 type: mpiaij 163 rows=764, cols=764, bs=4 164 total: nonzeros=35672, allocated nonzeros=35672 165 total number of mallocs used during MatSetValues calls=0 166 using I-node (on process 0) routines: found 50 nodes, limit used is 5 167 Up solver (post-smoother) same as down solver (pre-smoother) 168 Down solver (pre-smoother) on level 3 ------------------------------- 169 KSP Object: (mg_levels_3_) 4 MPI processes 170 type: richardson 171 damping factor=1. 172 maximum iterations=2, nonzero initial guess 173 tolerances: relative=1e-05, absolute=1e-50, divergence=10000. 174 left preconditioning 175 using NONE norm type for convergence test 176 PC Object: (mg_levels_3_) 4 MPI processes 177 type: sor 178 type = local_symmetric, iterations = 1, local iterations = 1, omega = 1. 179 linear system matrix = precond matrix: 180 Mat Object: 4 MPI processes 181 type: mpiaij 182 rows=7156, cols=7156, bs=4 183 total: nonzeros=259236, allocated nonzeros=259236 184 total number of mallocs used during MatSetValues calls=0 185 not using I-node (on process 0) routines 186 Up solver (post-smoother) same as down solver (pre-smoother) 187 Down solver (pre-smoother) on level 4 ------------------------------- 188 KSP Object: (mg_levels_4_) 4 MPI processes 189 type: richardson 190 damping factor=1. 191 maximum iterations=2, nonzero initial guess 192 tolerances: relative=1e-05, absolute=1e-50, divergence=10000. 193 left preconditioning 194 using NONE norm type for convergence test 195 PC Object: (mg_levels_4_) 4 MPI processes 196 type: sor 197 type = local_symmetric, iterations = 1, local iterations = 1, omega = 1. 198 linear system matrix = precond matrix: 199 Mat Object: 4 MPI processes 200 type: mpiaij 201 rows=37636, cols=37636, bs=4 202 total: nonzeros=746512, allocated nonzeros=746512 203 total number of mallocs used during MatSetValues calls=0 204 Up solver (post-smoother) same as down solver (pre-smoother) 205 linear system matrix = precond matrix: 206 Mat Object: 4 MPI processes 207 type: mpiaij 208 rows=37636, cols=37636, bs=4 209 total: nonzeros=746512, allocated nonzeros=746512 210 total number of mallocs used during MatSetValues calls=0 211Number of SNES iterations = 2 212