1c4762a1bSJed Brownlid velocity = 0.000106281, prandtl # = 1., grashof # = 1. 2c4762a1bSJed Brown 0 SNES Function norm 1.036007954337e-02 32d157150SStefano Zampini 0 KSP Residual norm 4.888657920112e-01 42d157150SStefano Zampini 1 KSP Residual norm 3.141627616286e-01 52d157150SStefano Zampini 2 KSP Residual norm 2.358812378525e-01 62d157150SStefano Zampini 3 KSP Residual norm 1.259045437573e-01 72d157150SStefano Zampini 4 KSP Residual norm 7.408575133258e-02 82d157150SStefano Zampini 5 KSP Residual norm 4.159658325629e-02 92d157150SStefano Zampini 6 KSP Residual norm 2.181861584091e-02 102d157150SStefano Zampini 7 KSP Residual norm 1.076244407256e-02 112d157150SStefano Zampini 8 KSP Residual norm 5.242924144023e-03 122d157150SStefano Zampini 9 KSP Residual norm 2.588726461570e-03 132d157150SStefano Zampini 10 KSP Residual norm 1.400700093556e-03 142d157150SStefano Zampini 11 KSP Residual norm 7.759158101290e-04 152d157150SStefano Zampini 12 KSP Residual norm 4.655744671687e-04 162d157150SStefano Zampini 13 KSP Residual norm 3.244240900196e-04 172d157150SStefano Zampini 14 KSP Residual norm 2.327144448481e-04 182d157150SStefano Zampini 15 KSP Residual norm 1.500951563875e-04 192d157150SStefano Zampini 16 KSP Residual norm 8.327955132177e-05 202d157150SStefano Zampini 17 KSP Residual norm 4.672468799190e-05 212d157150SStefano Zampini 18 KSP Residual norm 2.889321996959e-05 222d157150SStefano Zampini 19 KSP Residual norm 1.903514112532e-05 232d157150SStefano Zampini 20 KSP Residual norm 1.298123631775e-05 242d157150SStefano Zampini 21 KSP Residual norm 8.775705379722e-06 252d157150SStefano Zampini 22 KSP Residual norm 5.330175482523e-06 262d157150SStefano Zampini 23 KSP Residual norm 2.869775929854e-06 272d157150SStefano Zampini 1 SNES Function norm 1.370527153606e-06 282d157150SStefano Zampini 0 KSP Residual norm 1.177060567864e-05 292d157150SStefano Zampini 1 KSP Residual norm 6.928128579778e-06 302d157150SStefano Zampini 2 KSP Residual norm 4.304601930959e-06 312d157150SStefano Zampini 3 KSP Residual norm 2.454875323346e-06 322d157150SStefano Zampini 4 KSP Residual norm 1.512790232898e-06 332d157150SStefano Zampini 5 KSP Residual norm 1.001915525185e-06 342d157150SStefano Zampini 6 KSP Residual norm 7.387623107528e-07 352d157150SStefano Zampini 7 KSP Residual norm 4.953711915641e-07 362d157150SStefano Zampini 8 KSP Residual norm 3.144821641519e-07 372d157150SStefano Zampini 9 KSP Residual norm 2.134819078071e-07 382d157150SStefano Zampini 10 KSP Residual norm 1.518551249827e-07 392d157150SStefano Zampini 11 KSP Residual norm 1.014185303569e-07 402d157150SStefano Zampini 12 KSP Residual norm 5.608673008996e-08 412d157150SStefano Zampini 13 KSP Residual norm 3.030519067648e-08 422d157150SStefano Zampini 14 KSP Residual norm 1.786947998191e-08 432d157150SStefano Zampini 15 KSP Residual norm 1.177816744649e-08 442d157150SStefano Zampini 16 KSP Residual norm 8.997335866987e-09 452d157150SStefano Zampini 17 KSP Residual norm 6.356106556680e-09 462d157150SStefano Zampini 18 KSP Residual norm 4.147037399694e-09 472d157150SStefano Zampini 19 KSP Residual norm 2.711154831563e-09 482d157150SStefano Zampini 20 KSP Residual norm 1.750865415567e-09 492d157150SStefano Zampini 21 KSP Residual norm 1.127681613823e-09 502d157150SStefano Zampini 22 KSP Residual norm 7.783034009647e-10 512d157150SStefano Zampini 23 KSP Residual norm 5.047971816405e-10 522d157150SStefano Zampini 24 KSP Residual norm 3.137090395920e-10 532d157150SStefano Zampini 25 KSP Residual norm 2.182629803057e-10 542d157150SStefano Zampini 26 KSP Residual norm 1.265836293994e-10 552d157150SStefano Zampini 27 KSP Residual norm 7.756916170515e-11 562d157150SStefano 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 612d157150SStefano 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 69a99ef635SJonas Heinzmann maxlambda=1.000000e+00, 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 74f971d498SPierre Jolivet restart=30, using classical (unmodified) Gram-Schmidt orthogonalization with no iterative refinement 75*143f2514SPierre Jolivet 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 918405740aSBarry Smith not checking for convergence 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 1008405740aSBarry Smith not checking for convergence 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. 108ecf3d421SBarry Smith Factored matrix: 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 115ecf3d421SBarry Smith linear system matrix, which is also used to construct the preconditioner: 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 122ecf3d421SBarry Smith linear system matrix, which is also used to construct the preconditioner: 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 1368405740aSBarry Smith not checking for convergence 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. 140ecf3d421SBarry Smith linear system matrix, which is also used to construct the preconditioner: 141c4762a1bSJed Brown Mat Object: 4 MPI processes 142c4762a1bSJed Brown type: mpiaij 143c4762a1bSJed Brown rows=88, cols=88, bs=4 1442d157150SStefano Zampini total: nonzeros=4960, allocated nonzeros=4960 145c4762a1bSJed Brown total number of mallocs used during MatSetValues calls=0 1462d157150SStefano 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 1558405740aSBarry Smith not checking for convergence 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. 159ecf3d421SBarry Smith linear system matrix, which is also used to construct the preconditioner: 160c4762a1bSJed Brown Mat Object: 4 MPI processes 161c4762a1bSJed Brown type: mpiaij 1622d157150SStefano Zampini rows=764, cols=764, bs=4 1632d157150SStefano Zampini total: nonzeros=35672, allocated nonzeros=35672 164c4762a1bSJed Brown total number of mallocs used during MatSetValues calls=0 1652d157150SStefano 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 1748405740aSBarry Smith not checking for convergence 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. 178ecf3d421SBarry Smith linear system matrix, which is also used to construct the preconditioner: 179c4762a1bSJed Brown Mat Object: 4 MPI processes 180c4762a1bSJed Brown type: mpiaij 181c4762a1bSJed Brown rows=7156, cols=7156, bs=4 1822d157150SStefano 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 1938405740aSBarry Smith not checking for convergence 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. 197ecf3d421SBarry Smith linear system matrix, which is also used to construct the preconditioner: 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) 204ecf3d421SBarry Smith linear system matrix, which is also used to construct the preconditioner: 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