1c4762a1bSJed BrownLevel 0 domain size (m) 8e+04 x 8e+04 x 1e+03, num elements 4 x 4 x 2 (32), size (m) 20000. x 20000. x 1000. 2c4762a1bSJed Brown Linear solve converged due to CONVERGED_RTOL iterations 1 3c4762a1bSJed Brown Linear solve converged due to CONVERGED_RTOL iterations 1 4c4762a1bSJed Brown Linear solve converged due to CONVERGED_RTOL iterations 1 5c4762a1bSJed Brown Linear solve converged due to CONVERGED_RTOL iterations 1 6c4762a1bSJed Brown Linear solve converged due to CONVERGED_RTOL iterations 1 7c4762a1bSJed Brown Linear solve converged due to CONVERGED_RTOL iterations 1 8c4762a1bSJed Brown Linear solve converged due to CONVERGED_RTOL iterations 1 9c4762a1bSJed Brown Linear solve converged due to CONVERGED_RTOL iterations 1 10c4762a1bSJed Brown Linear solve converged due to CONVERGED_RTOL iterations 1 11c4762a1bSJed Brown Nonlinear solve converged due to CONVERGED_FNORM_RELATIVE iterations 9 12c4762a1bSJed BrownLevel 1 domain size (m) 8e+04 x 8e+04 x 1e+03, num elements 8 x 8 x 2 (128), size (m) 10000. x 10000. x 1000. 13c4762a1bSJed BrownLevel 0 domain size (m) 8e+04 x 8e+04 x 1e+03, num elements 4 x 4 x 2 (32), size (m) 20000. x 20000. x 1000. 14c4762a1bSJed Brown Linear solve converged due to CONVERGED_RTOL iterations 4 15c4762a1bSJed Brown Linear solve converged due to CONVERGED_RTOL iterations 5 16c4762a1bSJed Brown Linear solve converged due to CONVERGED_RTOL iterations 5 17c4762a1bSJed Brown Linear solve converged due to CONVERGED_RTOL iterations 6 18c4762a1bSJed Brown Linear solve converged due to CONVERGED_RTOL iterations 6 19c4762a1bSJed Brown Linear solve converged due to CONVERGED_RTOL iterations 5 20c4762a1bSJed Brown Linear solve converged due to CONVERGED_RTOL iterations 6 21c4762a1bSJed Brown Linear solve converged due to CONVERGED_RTOL iterations 7 22c4762a1bSJed Brown Nonlinear solve converged due to CONVERGED_FNORM_RELATIVE iterations 8 23c4762a1bSJed BrownLevel 2 domain size (m) 8e+04 x 8e+04 x 1e+03, num elements 16 x 16 x 2 (512), size (m) 5000. x 5000. x 1000. 24c4762a1bSJed BrownLevel 1 domain size (m) 8e+04 x 8e+04 x 1e+03, num elements 8 x 8 x 2 (128), size (m) 10000. x 10000. x 1000. 25c4762a1bSJed BrownLevel 0 domain size (m) 8e+04 x 8e+04 x 1e+03, num elements 4 x 4 x 2 (32), size (m) 20000. x 20000. x 1000. 26c4762a1bSJed Brown Linear solve converged due to CONVERGED_RTOL iterations 5 27c4762a1bSJed Brown Linear solve converged due to CONVERGED_RTOL iterations 6 28c4762a1bSJed Brown Linear solve converged due to CONVERGED_RTOL iterations 7 29c4762a1bSJed Brown Linear solve converged due to CONVERGED_RTOL iterations 6 30c4762a1bSJed Brown Linear solve converged due to CONVERGED_RTOL iterations 2 31c4762a1bSJed Brown Linear solve converged due to CONVERGED_RTOL iterations 2 32c4762a1bSJed Brown Linear solve converged due to CONVERGED_RTOL iterations 3 33c4762a1bSJed Brown Linear solve converged due to CONVERGED_RTOL iterations 5 34c4762a1bSJed Brown Linear solve converged due to CONVERGED_RTOL iterations 6 35c4762a1bSJed Brown Linear solve converged due to CONVERGED_RTOL iterations 7 36c4762a1bSJed Brown Nonlinear solve converged due to CONVERGED_FNORM_RELATIVE iterations 10 37c4762a1bSJed BrownLevel 3 domain size (m) 8e+04 x 8e+04 x 1e+03, num elements 16 x 16 x 4 (1024), size (m) 5000. x 5000. x 333.333 38c4762a1bSJed BrownLevel 2 domain size (m) 8e+04 x 8e+04 x 1e+03, num elements 16 x 16 x 2 (512), size (m) 5000. x 5000. x 1000. 39c4762a1bSJed BrownLevel 1 domain size (m) 8e+04 x 8e+04 x 1e+03, num elements 8 x 8 x 2 (128), size (m) 10000. x 10000. x 1000. 40c4762a1bSJed BrownLevel 0 domain size (m) 8e+04 x 8e+04 x 1e+03, num elements 4 x 4 x 2 (32), size (m) 20000. x 20000. x 1000. 41c4762a1bSJed Brown Linear solve converged due to CONVERGED_RTOL iterations 2 42c4762a1bSJed Brown Linear solve converged due to CONVERGED_RTOL iterations 2 43c4762a1bSJed Brown Linear solve converged due to CONVERGED_RTOL iterations 4 44c4762a1bSJed Brown Linear solve converged due to CONVERGED_RTOL iterations 7 45c4762a1bSJed Brown Linear solve converged due to CONVERGED_RTOL iterations 7 46c4762a1bSJed BrownNonlinear solve converged due to CONVERGED_FNORM_RELATIVE iterations 5 47c4762a1bSJed BrownSNES Object: 6 MPI processes 48c4762a1bSJed Brown type: newtonls 49c4762a1bSJed Brown maximum iterations=50, maximum function evaluations=10000 50c4762a1bSJed Brown tolerances: relative=1e-08, absolute=1e-50, solution=1e-12 51c4762a1bSJed Brown total number of linear solver iterations=22 52c4762a1bSJed Brown total number of function evaluations=6 53c4762a1bSJed Brown norm schedule ALWAYS 54c4762a1bSJed Brown total number of grid sequence refinements=3 55c4762a1bSJed Brown Jacobian is built using a DMDA local Jacobian 56c4762a1bSJed Brown SNESLineSearch Object: 6 MPI processes 57c4762a1bSJed Brown type: bt 58c4762a1bSJed Brown interpolation: cubic 59c4762a1bSJed Brown alpha=1.000000e-04 60a99ef635SJonas Heinzmann maxlambda=1.000000e+00, minlambda=1.000000e-12 61c4762a1bSJed Brown tolerances: relative=1.000000e-08, absolute=1.000000e-15, lambda=1.000000e-08 62c4762a1bSJed Brown maximum iterations=40 63c4762a1bSJed Brown KSP Object: 6 MPI processes 64c4762a1bSJed Brown type: fgmres 65f971d498SPierre Jolivet restart=30, using classical (unmodified) Gram-Schmidt orthogonalization with no iterative refinement 66*143f2514SPierre Jolivet happy breakdown tolerance=1e-30 67c4762a1bSJed Brown maximum iterations=10000, initial guess is zero 68c4762a1bSJed Brown tolerances: relative=0.01, absolute=-1., divergence=10000. 69c4762a1bSJed Brown right preconditioning 70c4762a1bSJed Brown using UNPRECONDITIONED norm type for convergence test 71c4762a1bSJed Brown PC Object: 6 MPI processes 72c4762a1bSJed Brown type: mg 73c4762a1bSJed Brown type is MULTIPLICATIVE, levels=4 cycles=v 74c4762a1bSJed Brown Cycles per PCApply=1 75c4762a1bSJed Brown Not using Galerkin computed coarse grid matrices 7663a3b9bcSJacob Faibussowitsch Coarse grid solver -- level 0 ------------------------------- 77c4762a1bSJed Brown KSP Object: (mg_coarse_) 6 MPI processes 78c4762a1bSJed Brown type: preonly 79c4762a1bSJed Brown maximum iterations=10000, initial guess is zero 80c4762a1bSJed Brown tolerances: relative=1e-05, absolute=1e-50, divergence=10000. 81c4762a1bSJed Brown left preconditioning 828405740aSBarry Smith not checking for convergence 83c4762a1bSJed Brown PC Object: (mg_coarse_) 6 MPI processes 84c4762a1bSJed Brown type: redundant 85c4762a1bSJed Brown First (color=0) of 6 PCs follows 868cc725e6SPierre Jolivet KSP Object: (mg_coarse_redundant_) 1 MPI process 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 928cc725e6SPierre Jolivet PC Object: (mg_coarse_redundant_) 1 MPI process 93c4762a1bSJed Brown type: lu 94c4762a1bSJed Brown out-of-place factorization 95c4762a1bSJed Brown tolerance for zero pivot 2.22045e-14 96c4762a1bSJed Brown using diagonal shift on blocks to prevent zero pivot [INBLOCKS] 97c4762a1bSJed Brown matrix ordering: nd 98c4762a1bSJed Brown factor fill ratio given 5., needed 1.41667 99ecf3d421SBarry Smith Factored matrix: 10026cc229bSBarry Smith Mat Object: (mg_coarse_redundant_) 1 MPI process 101c4762a1bSJed Brown type: seqaij 102c4762a1bSJed Brown rows=64, cols=64, bs=2 103c4762a1bSJed Brown package used to perform factorization: petsc 104c4762a1bSJed Brown total: nonzeros=3264, allocated nonzeros=3264 105c4762a1bSJed Brown using I-node routines: found 16 nodes, limit used is 5 106ecf3d421SBarry Smith linear system matrix, which is also used to construct the preconditioner: 10726cc229bSBarry Smith Mat Object: 1 MPI process 108c4762a1bSJed Brown type: seqaij 109c4762a1bSJed Brown rows=64, cols=64, bs=2 110c4762a1bSJed Brown total: nonzeros=2304, allocated nonzeros=2304 111c4762a1bSJed Brown total number of mallocs used during MatSetValues calls=0 112c4762a1bSJed Brown using I-node routines: found 16 nodes, limit used is 5 113ecf3d421SBarry Smith linear system matrix, which is also used to construct the preconditioner: 114c4762a1bSJed Brown Mat Object: 6 MPI processes 115c4762a1bSJed Brown type: mpiaij 116c4762a1bSJed Brown rows=64, cols=64, bs=2 117c4762a1bSJed Brown total: nonzeros=2304, allocated nonzeros=2304 118c4762a1bSJed Brown total number of mallocs used during MatSetValues calls=0 119c4762a1bSJed Brown Down solver (pre-smoother) on level 1 ------------------------------- 120c4762a1bSJed Brown KSP Object: (mg_levels_1_) 6 MPI processes 121c4762a1bSJed Brown type: gmres 122f971d498SPierre Jolivet restart=30, using classical (unmodified) Gram-Schmidt orthogonalization with no iterative refinement 123*143f2514SPierre Jolivet happy breakdown tolerance=1e-30 124c4762a1bSJed Brown maximum iterations=1, nonzero initial guess 125c4762a1bSJed Brown tolerances: relative=1e-05, absolute=1e-50, divergence=10000. 126c4762a1bSJed Brown left preconditioning 1278405740aSBarry Smith not checking for convergence 128c4762a1bSJed Brown PC Object: (mg_levels_1_) 6 MPI processes 129c4762a1bSJed Brown type: bjacobi 130c4762a1bSJed Brown number of blocks = 6 131020d6619SPierre Jolivet Local solver information for first block is in the following KSP and PC objects on rank 0: 132020d6619SPierre Jolivet Use -mg_levels_1_ksp_view ::ascii_info_detail to display information for all blocks 1338cc725e6SPierre Jolivet KSP Object: (mg_levels_1_sub_) 1 MPI process 134c4762a1bSJed Brown type: preonly 135c4762a1bSJed Brown maximum iterations=10000, initial guess is zero 136c4762a1bSJed Brown tolerances: relative=1e-05, absolute=1e-50, divergence=10000. 137c4762a1bSJed Brown left preconditioning 1388405740aSBarry Smith not checking for convergence 1398cc725e6SPierre Jolivet PC Object: (mg_levels_1_sub_) 1 MPI process 140c4762a1bSJed Brown type: cholesky 141c4762a1bSJed Brown out-of-place factorization 142c4762a1bSJed Brown tolerance for zero pivot 2.22045e-14 143c4762a1bSJed Brown matrix ordering: natural 144c4762a1bSJed Brown factor fill ratio given 5., needed 1.32237 145ecf3d421SBarry Smith Factored matrix: 14626cc229bSBarry Smith Mat Object: (mg_levels_1_sub_) 1 MPI process 147c4762a1bSJed Brown type: seqsbaij 148c4762a1bSJed Brown rows=48, cols=48, bs=2 149c4762a1bSJed Brown package used to perform factorization: petsc 150c4762a1bSJed Brown total: nonzeros=804, allocated nonzeros=804 151ecf3d421SBarry Smith linear system matrix, which is also used to construct the preconditioner: 1528cc725e6SPierre Jolivet Mat Object: (mg_levels_1_sub_) 1 MPI process 153c4762a1bSJed Brown type: seqsbaij 154c4762a1bSJed Brown rows=48, cols=48, bs=2 155c4762a1bSJed Brown total: nonzeros=608, allocated nonzeros=608 156c4762a1bSJed Brown total number of mallocs used during MatSetValues calls=0 157ecf3d421SBarry Smith linear system matrix, which is also used to construct the preconditioner: 158c4762a1bSJed Brown Mat Object: 6 MPI processes 159c4762a1bSJed Brown type: mpisbaij 160c4762a1bSJed Brown rows=256, cols=256, bs=2 161c4762a1bSJed Brown total: nonzeros=4864, allocated nonzeros=4864 162c4762a1bSJed Brown total number of mallocs used during MatSetValues calls=0 163c4762a1bSJed Brown Up solver (post-smoother) same as down solver (pre-smoother) 164c4762a1bSJed Brown Down solver (pre-smoother) on level 2 ------------------------------- 165c4762a1bSJed Brown KSP Object: (mg_levels_2_) 6 MPI processes 166c4762a1bSJed Brown type: gmres 167f971d498SPierre Jolivet restart=30, using classical (unmodified) Gram-Schmidt orthogonalization with no iterative refinement 168*143f2514SPierre Jolivet happy breakdown tolerance=1e-30 169c4762a1bSJed Brown maximum iterations=1, nonzero initial guess 170c4762a1bSJed Brown tolerances: relative=1e-05, absolute=1e-50, divergence=10000. 171c4762a1bSJed Brown left preconditioning 1728405740aSBarry Smith not checking for convergence 173c4762a1bSJed Brown PC Object: (mg_levels_2_) 6 MPI processes 174c4762a1bSJed Brown type: bjacobi 175c4762a1bSJed Brown number of blocks = 6 176020d6619SPierre Jolivet Local solver information for first block is in the following KSP and PC objects on rank 0: 177020d6619SPierre Jolivet Use -mg_levels_2_ksp_view ::ascii_info_detail to display information for all blocks 1788cc725e6SPierre Jolivet KSP Object: (mg_levels_2_sub_) 1 MPI process 179c4762a1bSJed Brown type: preonly 180c4762a1bSJed Brown maximum iterations=10000, initial guess is zero 181c4762a1bSJed Brown tolerances: relative=1e-05, absolute=1e-50, divergence=10000. 182c4762a1bSJed Brown left preconditioning 1838405740aSBarry Smith not checking for convergence 1848cc725e6SPierre Jolivet PC Object: (mg_levels_2_sub_) 1 MPI process 185c4762a1bSJed Brown type: icc 186c4762a1bSJed Brown out-of-place factorization 187c4762a1bSJed Brown 0 levels of fill 188c4762a1bSJed Brown tolerance for zero pivot 2.22045e-14 189c4762a1bSJed Brown using Manteuffel shift [POSITIVE_DEFINITE] 190c4762a1bSJed Brown matrix ordering: natural 191c4762a1bSJed Brown factor fill ratio given 1., needed 1.00133 192ecf3d421SBarry Smith Factored matrix: 19326cc229bSBarry Smith Mat Object: (mg_levels_2_sub_) 1 MPI process 194c4762a1bSJed Brown type: seqsbaij 195c4762a1bSJed Brown rows=192, cols=192, bs=2 196c4762a1bSJed Brown package used to perform factorization: petsc 197c4762a1bSJed Brown total: nonzeros=3012, allocated nonzeros=3012 198ecf3d421SBarry Smith linear system matrix, which is also used to construct the preconditioner: 1998cc725e6SPierre Jolivet Mat Object: (mg_levels_2_sub_) 1 MPI process 200c4762a1bSJed Brown type: seqsbaij 201c4762a1bSJed Brown rows=192, cols=192, bs=2 202c4762a1bSJed Brown total: nonzeros=3008, allocated nonzeros=3008 203c4762a1bSJed Brown total number of mallocs used during MatSetValues calls=0 204ecf3d421SBarry Smith linear system matrix, which is also used to construct the preconditioner: 205c4762a1bSJed Brown Mat Object: 6 MPI processes 206c4762a1bSJed Brown type: mpisbaij 207c4762a1bSJed Brown rows=1024, cols=1024, bs=2 208c4762a1bSJed Brown total: nonzeros=19456, allocated nonzeros=19456 209c4762a1bSJed Brown total number of mallocs used during MatSetValues calls=0 210c4762a1bSJed Brown Up solver (post-smoother) same as down solver (pre-smoother) 211c4762a1bSJed Brown Down solver (pre-smoother) on level 3 ------------------------------- 212c4762a1bSJed Brown KSP Object: (mg_levels_3_) 6 MPI processes 213c4762a1bSJed Brown type: gmres 214f971d498SPierre Jolivet restart=30, using classical (unmodified) Gram-Schmidt orthogonalization with no iterative refinement 215*143f2514SPierre Jolivet happy breakdown tolerance=1e-30 216c4762a1bSJed Brown maximum iterations=1, nonzero initial guess 217c4762a1bSJed Brown tolerances: relative=1e-05, absolute=1e-50, divergence=10000. 218c4762a1bSJed Brown left preconditioning 2198405740aSBarry Smith not checking for convergence 220c4762a1bSJed Brown PC Object: (mg_levels_3_) 6 MPI processes 221c4762a1bSJed Brown type: bjacobi 222c4762a1bSJed Brown number of blocks = 6 223020d6619SPierre Jolivet Local solver information for first block is in the following KSP and PC objects on rank 0: 224020d6619SPierre Jolivet Use -mg_levels_3_ksp_view ::ascii_info_detail to display information for all blocks 2258cc725e6SPierre Jolivet KSP Object: (mg_levels_3_sub_) 1 MPI process 226c4762a1bSJed Brown type: preonly 227c4762a1bSJed Brown maximum iterations=10000, initial guess is zero 228c4762a1bSJed Brown tolerances: relative=1e-05, absolute=1e-50, divergence=10000. 229c4762a1bSJed Brown left preconditioning 2308405740aSBarry Smith not checking for convergence 2318cc725e6SPierre Jolivet PC Object: (mg_levels_3_sub_) 1 MPI process 232c4762a1bSJed Brown type: icc 233c4762a1bSJed Brown out-of-place factorization 234c4762a1bSJed Brown 0 levels of fill 235c4762a1bSJed Brown tolerance for zero pivot 2.22045e-14 236c4762a1bSJed Brown using Manteuffel shift [POSITIVE_DEFINITE] 237c4762a1bSJed Brown matrix ordering: natural 238c4762a1bSJed Brown factor fill ratio given 1., needed 1.00054 239ecf3d421SBarry Smith Factored matrix: 24026cc229bSBarry Smith Mat Object: (mg_levels_3_sub_) 1 MPI process 241c4762a1bSJed Brown type: seqsbaij 242c4762a1bSJed Brown rows=384, cols=384, bs=2 243c4762a1bSJed Brown package used to perform factorization: petsc 244c4762a1bSJed Brown total: nonzeros=7428, allocated nonzeros=7428 245ecf3d421SBarry Smith linear system matrix, which is also used to construct the preconditioner: 2468cc725e6SPierre Jolivet Mat Object: (mg_levels_3_sub_) 1 MPI process 247c4762a1bSJed Brown type: seqsbaij 248c4762a1bSJed Brown rows=384, cols=384, bs=2 249c4762a1bSJed Brown total: nonzeros=7424, allocated nonzeros=7424 250c4762a1bSJed Brown total number of mallocs used during MatSetValues calls=0 251ecf3d421SBarry Smith linear system matrix, which is also used to construct the preconditioner: 252c4762a1bSJed Brown Mat Object: 6 MPI processes 253c4762a1bSJed Brown type: mpisbaij 254c4762a1bSJed Brown rows=2048, cols=2048, bs=2 255c4762a1bSJed Brown total: nonzeros=48128, allocated nonzeros=48128 256c4762a1bSJed Brown total number of mallocs used during MatSetValues calls=0 257c4762a1bSJed Brown Up solver (post-smoother) same as down solver (pre-smoother) 258ecf3d421SBarry Smith linear system matrix, which is also used to construct the preconditioner: 259c4762a1bSJed Brown Mat Object: 6 MPI processes 260c4762a1bSJed Brown type: mpisbaij 261c4762a1bSJed Brown rows=2048, cols=2048, bs=2 262c4762a1bSJed Brown total: nonzeros=48128, allocated nonzeros=48128 263c4762a1bSJed Brown total number of mallocs used during MatSetValues calls=0 264c4762a1bSJed BrownSolution statistics after solve: Full 265c4762a1bSJed BrownCONVERGED_FNORM_RELATIVE: Number of SNES iterations = 5, total linear iterations = 22 266c4762a1bSJed Brown|X|_2 564267. -187.04 <= u <= 23814.8 -979.938 <= v <= 5720.23 3.24766 <= c <= 23904.5 267c4762a1bSJed BrownSurface statistics: u in [-6.945709e+01, 2.381387e+04] mean 1.520632e+04 268c4762a1bSJed BrownGlobal eta range 2.50736e+09 to 9.2273e+12 converged range 2.50736e+09 to 3.8249e+10 269c4762a1bSJed BrownGlobal beta2 range 0. to 4.97763e+09 converged range 0. to 4.97763e+09 270