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 60c4762a1bSJed Brown maxstep=1.000000e+08, 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 65c4762a1bSJed Brown restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement 66c4762a1bSJed Brown 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 76c4762a1bSJed Brown Coarse grid solver -- level ------------------------------- 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 82c4762a1bSJed Brown using NONE norm type for convergence test 83c4762a1bSJed Brown PC Object: (mg_coarse_) 6 MPI processes 84c4762a1bSJed Brown type: redundant 85c4762a1bSJed Brown First (color=0) of 6 PCs follows 86c4762a1bSJed Brown KSP Object: (mg_coarse_redundant_) 1 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_redundant_) 1 MPI processes 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 99c4762a1bSJed Brown Factored matrix follows: 100c4762a1bSJed Brown Mat Object: 1 MPI processes 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 106c4762a1bSJed Brown linear system matrix = precond matrix: 107c4762a1bSJed Brown Mat Object: 1 MPI processes 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 113c4762a1bSJed Brown linear system matrix = precond matrix: 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 122c4762a1bSJed Brown restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement 123c4762a1bSJed Brown 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 127c4762a1bSJed Brown using NONE norm type for convergence test 128c4762a1bSJed Brown PC Object: (mg_levels_1_) 6 MPI processes 129c4762a1bSJed Brown type: bjacobi 130c4762a1bSJed Brown number of blocks = 6 1318d76b567SPatrick Sanan Local solver is the same for all blocks, as in the following KSP and PC objects on rank 0: 132c4762a1bSJed Brown KSP Object: (mg_levels_1_sub_) 1 MPI processes 133c4762a1bSJed Brown type: preonly 134c4762a1bSJed Brown maximum iterations=10000, initial guess is zero 135c4762a1bSJed Brown tolerances: relative=1e-05, absolute=1e-50, divergence=10000. 136c4762a1bSJed Brown left preconditioning 137c4762a1bSJed Brown using NONE norm type for convergence test 138c4762a1bSJed Brown PC Object: (mg_levels_1_sub_) 1 MPI processes 139c4762a1bSJed Brown type: cholesky 140c4762a1bSJed Brown out-of-place factorization 141c4762a1bSJed Brown tolerance for zero pivot 2.22045e-14 142c4762a1bSJed Brown matrix ordering: natural 143c4762a1bSJed Brown factor fill ratio given 5., needed 1.32237 144c4762a1bSJed Brown Factored matrix follows: 145c4762a1bSJed Brown Mat Object: 1 MPI processes 146c4762a1bSJed Brown type: seqsbaij 147c4762a1bSJed Brown rows=48, cols=48, bs=2 148c4762a1bSJed Brown package used to perform factorization: petsc 149c4762a1bSJed Brown total: nonzeros=804, allocated nonzeros=804 150c4762a1bSJed Brown block size is 2 151c4762a1bSJed Brown linear system matrix = precond matrix: 152*793a3527SPierre Jolivet Mat Object: (mg_levels_1_sub_) 1 MPI processes 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 157c4762a1bSJed Brown block size is 2 158c4762a1bSJed Brown linear system matrix = precond matrix: 159c4762a1bSJed Brown Mat Object: 6 MPI processes 160c4762a1bSJed Brown type: mpisbaij 161c4762a1bSJed Brown rows=256, cols=256, bs=2 162c4762a1bSJed Brown total: nonzeros=4864, allocated nonzeros=4864 163c4762a1bSJed Brown total number of mallocs used during MatSetValues calls=0 164c4762a1bSJed Brown Up solver (post-smoother) same as down solver (pre-smoother) 165c4762a1bSJed Brown Down solver (pre-smoother) on level 2 ------------------------------- 166c4762a1bSJed Brown KSP Object: (mg_levels_2_) 6 MPI processes 167c4762a1bSJed Brown type: gmres 168c4762a1bSJed Brown restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement 169c4762a1bSJed Brown happy breakdown tolerance 1e-30 170c4762a1bSJed Brown maximum iterations=1, nonzero initial guess 171c4762a1bSJed Brown tolerances: relative=1e-05, absolute=1e-50, divergence=10000. 172c4762a1bSJed Brown left preconditioning 173c4762a1bSJed Brown using NONE norm type for convergence test 174c4762a1bSJed Brown PC Object: (mg_levels_2_) 6 MPI processes 175c4762a1bSJed Brown type: bjacobi 176c4762a1bSJed Brown number of blocks = 6 1778d76b567SPatrick Sanan Local solver is the same for all blocks, as in the following KSP and PC objects on rank 0: 178c4762a1bSJed Brown KSP Object: (mg_levels_2_sub_) 1 MPI processes 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 183c4762a1bSJed Brown using NONE norm type for convergence test 184c4762a1bSJed Brown PC Object: (mg_levels_2_sub_) 1 MPI processes 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 192c4762a1bSJed Brown Factored matrix follows: 193c4762a1bSJed Brown Mat Object: 1 MPI processes 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 198c4762a1bSJed Brown block size is 2 199c4762a1bSJed Brown linear system matrix = precond matrix: 200*793a3527SPierre Jolivet Mat Object: (mg_levels_2_sub_) 1 MPI processes 201c4762a1bSJed Brown type: seqsbaij 202c4762a1bSJed Brown rows=192, cols=192, bs=2 203c4762a1bSJed Brown total: nonzeros=3008, allocated nonzeros=3008 204c4762a1bSJed Brown total number of mallocs used during MatSetValues calls=0 205c4762a1bSJed Brown block size is 2 206c4762a1bSJed Brown linear system matrix = precond matrix: 207c4762a1bSJed Brown Mat Object: 6 MPI processes 208c4762a1bSJed Brown type: mpisbaij 209c4762a1bSJed Brown rows=1024, cols=1024, bs=2 210c4762a1bSJed Brown total: nonzeros=19456, allocated nonzeros=19456 211c4762a1bSJed Brown total number of mallocs used during MatSetValues calls=0 212c4762a1bSJed Brown Up solver (post-smoother) same as down solver (pre-smoother) 213c4762a1bSJed Brown Down solver (pre-smoother) on level 3 ------------------------------- 214c4762a1bSJed Brown KSP Object: (mg_levels_3_) 6 MPI processes 215c4762a1bSJed Brown type: gmres 216c4762a1bSJed Brown restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement 217c4762a1bSJed Brown happy breakdown tolerance 1e-30 218c4762a1bSJed Brown maximum iterations=1, nonzero initial guess 219c4762a1bSJed Brown tolerances: relative=1e-05, absolute=1e-50, divergence=10000. 220c4762a1bSJed Brown left preconditioning 221c4762a1bSJed Brown using NONE norm type for convergence test 222c4762a1bSJed Brown PC Object: (mg_levels_3_) 6 MPI processes 223c4762a1bSJed Brown type: bjacobi 224c4762a1bSJed Brown number of blocks = 6 2258d76b567SPatrick Sanan Local solver is the same for all blocks, as in the following KSP and PC objects on rank 0: 226c4762a1bSJed Brown KSP Object: (mg_levels_3_sub_) 1 MPI processes 227c4762a1bSJed Brown type: preonly 228c4762a1bSJed Brown maximum iterations=10000, initial guess is zero 229c4762a1bSJed Brown tolerances: relative=1e-05, absolute=1e-50, divergence=10000. 230c4762a1bSJed Brown left preconditioning 231c4762a1bSJed Brown using NONE norm type for convergence test 232c4762a1bSJed Brown PC Object: (mg_levels_3_sub_) 1 MPI processes 233c4762a1bSJed Brown type: icc 234c4762a1bSJed Brown out-of-place factorization 235c4762a1bSJed Brown 0 levels of fill 236c4762a1bSJed Brown tolerance for zero pivot 2.22045e-14 237c4762a1bSJed Brown using Manteuffel shift [POSITIVE_DEFINITE] 238c4762a1bSJed Brown matrix ordering: natural 239c4762a1bSJed Brown factor fill ratio given 1., needed 1.00054 240c4762a1bSJed Brown Factored matrix follows: 241c4762a1bSJed Brown Mat Object: 1 MPI processes 242c4762a1bSJed Brown type: seqsbaij 243c4762a1bSJed Brown rows=384, cols=384, bs=2 244c4762a1bSJed Brown package used to perform factorization: petsc 245c4762a1bSJed Brown total: nonzeros=7428, allocated nonzeros=7428 246c4762a1bSJed Brown block size is 2 247c4762a1bSJed Brown linear system matrix = precond matrix: 248*793a3527SPierre Jolivet Mat Object: (mg_levels_3_sub_) 1 MPI processes 249c4762a1bSJed Brown type: seqsbaij 250c4762a1bSJed Brown rows=384, cols=384, bs=2 251c4762a1bSJed Brown total: nonzeros=7424, allocated nonzeros=7424 252c4762a1bSJed Brown total number of mallocs used during MatSetValues calls=0 253c4762a1bSJed Brown block size is 2 254c4762a1bSJed Brown linear system matrix = precond matrix: 255c4762a1bSJed Brown Mat Object: 6 MPI processes 256c4762a1bSJed Brown type: mpisbaij 257c4762a1bSJed Brown rows=2048, cols=2048, bs=2 258c4762a1bSJed Brown total: nonzeros=48128, allocated nonzeros=48128 259c4762a1bSJed Brown total number of mallocs used during MatSetValues calls=0 260c4762a1bSJed Brown Up solver (post-smoother) same as down solver (pre-smoother) 261c4762a1bSJed Brown linear system matrix = precond matrix: 262c4762a1bSJed Brown Mat Object: 6 MPI processes 263c4762a1bSJed Brown type: mpisbaij 264c4762a1bSJed Brown rows=2048, cols=2048, bs=2 265c4762a1bSJed Brown total: nonzeros=48128, allocated nonzeros=48128 266c4762a1bSJed Brown total number of mallocs used during MatSetValues calls=0 267c4762a1bSJed BrownSolution statistics after solve: Full 268c4762a1bSJed BrownCONVERGED_FNORM_RELATIVE: Number of SNES iterations = 5, total linear iterations = 22 269c4762a1bSJed Brown|X|_2 564267. -187.04 <= u <= 23814.8 -979.938 <= v <= 5720.23 3.24766 <= c <= 23904.5 270c4762a1bSJed BrownSurface statistics: u in [-6.945709e+01, 2.381387e+04] mean 1.520632e+04 271c4762a1bSJed BrownGlobal eta range 2.50736e+09 to 9.2273e+12 converged range 2.50736e+09 to 3.8249e+10 272c4762a1bSJed BrownGlobal beta2 range 0. to 4.97763e+09 converged range 0. to 4.97763e+09 273