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