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 total number of mallocs used during MatSetValues calls=0 106c4762a1bSJed Brown using I-node routines: found 16 nodes, limit used is 5 107c4762a1bSJed Brown linear system matrix = precond matrix: 108c4762a1bSJed Brown Mat Object: 1 MPI processes 109c4762a1bSJed Brown type: seqaij 110c4762a1bSJed Brown rows=64, cols=64, bs=2 111c4762a1bSJed Brown total: nonzeros=2304, allocated nonzeros=2304 112c4762a1bSJed Brown total number of mallocs used during MatSetValues calls=0 113c4762a1bSJed Brown using I-node routines: found 16 nodes, limit used is 5 114c4762a1bSJed Brown linear system matrix = precond matrix: 115c4762a1bSJed Brown Mat Object: 6 MPI processes 116c4762a1bSJed Brown type: mpiaij 117c4762a1bSJed Brown rows=64, cols=64, bs=2 118c4762a1bSJed Brown total: nonzeros=2304, allocated nonzeros=2304 119c4762a1bSJed Brown total number of mallocs used during MatSetValues calls=0 120c4762a1bSJed Brown Down solver (pre-smoother) on level 1 ------------------------------- 121c4762a1bSJed Brown KSP Object: (mg_levels_1_) 6 MPI processes 122c4762a1bSJed Brown type: gmres 123c4762a1bSJed Brown restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement 124c4762a1bSJed Brown happy breakdown tolerance 1e-30 125c4762a1bSJed Brown maximum iterations=1, nonzero initial guess 126c4762a1bSJed Brown tolerances: relative=1e-05, absolute=1e-50, divergence=10000. 127c4762a1bSJed Brown left preconditioning 128c4762a1bSJed Brown using NONE norm type for convergence test 129c4762a1bSJed Brown PC Object: (mg_levels_1_) 6 MPI processes 130c4762a1bSJed Brown type: bjacobi 131c4762a1bSJed Brown number of blocks = 6 132*8d76b567SPatrick Sanan Local solver is the same for all blocks, as in the following KSP and PC objects on rank 0: 133c4762a1bSJed Brown KSP Object: (mg_levels_1_sub_) 1 MPI processes 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 138c4762a1bSJed Brown using NONE norm type for convergence test 139c4762a1bSJed Brown PC Object: (mg_levels_1_sub_) 1 MPI processes 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 145c4762a1bSJed Brown Factored matrix follows: 146c4762a1bSJed Brown Mat Object: 1 MPI processes 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 151c4762a1bSJed Brown total number of mallocs used during MatSetValues calls=0 152c4762a1bSJed Brown block size is 2 153c4762a1bSJed Brown linear system matrix = precond matrix: 154c4762a1bSJed Brown Mat Object: 1 MPI processes 155c4762a1bSJed Brown type: seqsbaij 156c4762a1bSJed Brown rows=48, cols=48, bs=2 157c4762a1bSJed Brown total: nonzeros=608, allocated nonzeros=608 158c4762a1bSJed Brown total number of mallocs used during MatSetValues calls=0 159c4762a1bSJed Brown block size is 2 160c4762a1bSJed Brown linear system matrix = precond matrix: 161c4762a1bSJed Brown Mat Object: 6 MPI processes 162c4762a1bSJed Brown type: mpisbaij 163c4762a1bSJed Brown rows=256, cols=256, bs=2 164c4762a1bSJed Brown total: nonzeros=4864, allocated nonzeros=4864 165c4762a1bSJed Brown total number of mallocs used during MatSetValues calls=0 166c4762a1bSJed Brown Up solver (post-smoother) same as down solver (pre-smoother) 167c4762a1bSJed Brown Down solver (pre-smoother) on level 2 ------------------------------- 168c4762a1bSJed Brown KSP Object: (mg_levels_2_) 6 MPI processes 169c4762a1bSJed Brown type: gmres 170c4762a1bSJed Brown restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement 171c4762a1bSJed Brown happy breakdown tolerance 1e-30 172c4762a1bSJed Brown maximum iterations=1, nonzero initial guess 173c4762a1bSJed Brown tolerances: relative=1e-05, absolute=1e-50, divergence=10000. 174c4762a1bSJed Brown left preconditioning 175c4762a1bSJed Brown using NONE norm type for convergence test 176c4762a1bSJed Brown PC Object: (mg_levels_2_) 6 MPI processes 177c4762a1bSJed Brown type: bjacobi 178c4762a1bSJed Brown number of blocks = 6 179*8d76b567SPatrick Sanan Local solver is the same for all blocks, as in the following KSP and PC objects on rank 0: 180c4762a1bSJed Brown KSP Object: (mg_levels_2_sub_) 1 MPI processes 181c4762a1bSJed Brown type: preonly 182c4762a1bSJed Brown maximum iterations=10000, initial guess is zero 183c4762a1bSJed Brown tolerances: relative=1e-05, absolute=1e-50, divergence=10000. 184c4762a1bSJed Brown left preconditioning 185c4762a1bSJed Brown using NONE norm type for convergence test 186c4762a1bSJed Brown PC Object: (mg_levels_2_sub_) 1 MPI processes 187c4762a1bSJed Brown type: icc 188c4762a1bSJed Brown out-of-place factorization 189c4762a1bSJed Brown 0 levels of fill 190c4762a1bSJed Brown tolerance for zero pivot 2.22045e-14 191c4762a1bSJed Brown using Manteuffel shift [POSITIVE_DEFINITE] 192c4762a1bSJed Brown matrix ordering: natural 193c4762a1bSJed Brown factor fill ratio given 1., needed 1.00133 194c4762a1bSJed Brown Factored matrix follows: 195c4762a1bSJed Brown Mat Object: 1 MPI processes 196c4762a1bSJed Brown type: seqsbaij 197c4762a1bSJed Brown rows=192, cols=192, bs=2 198c4762a1bSJed Brown package used to perform factorization: petsc 199c4762a1bSJed Brown total: nonzeros=3012, allocated nonzeros=3012 200c4762a1bSJed Brown total number of mallocs used during MatSetValues calls=0 201c4762a1bSJed Brown block size is 2 202c4762a1bSJed Brown linear system matrix = precond matrix: 203c4762a1bSJed Brown Mat Object: 1 MPI processes 204c4762a1bSJed Brown type: seqsbaij 205c4762a1bSJed Brown rows=192, cols=192, bs=2 206c4762a1bSJed Brown total: nonzeros=3008, allocated nonzeros=3008 207c4762a1bSJed Brown total number of mallocs used during MatSetValues calls=0 208c4762a1bSJed Brown block size is 2 209c4762a1bSJed Brown linear system matrix = precond matrix: 210c4762a1bSJed Brown Mat Object: 6 MPI processes 211c4762a1bSJed Brown type: mpisbaij 212c4762a1bSJed Brown rows=1024, cols=1024, bs=2 213c4762a1bSJed Brown total: nonzeros=19456, allocated nonzeros=19456 214c4762a1bSJed Brown total number of mallocs used during MatSetValues calls=0 215c4762a1bSJed Brown Up solver (post-smoother) same as down solver (pre-smoother) 216c4762a1bSJed Brown Down solver (pre-smoother) on level 3 ------------------------------- 217c4762a1bSJed Brown KSP Object: (mg_levels_3_) 6 MPI processes 218c4762a1bSJed Brown type: gmres 219c4762a1bSJed Brown restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement 220c4762a1bSJed Brown happy breakdown tolerance 1e-30 221c4762a1bSJed Brown maximum iterations=1, nonzero initial guess 222c4762a1bSJed Brown tolerances: relative=1e-05, absolute=1e-50, divergence=10000. 223c4762a1bSJed Brown left preconditioning 224c4762a1bSJed Brown using NONE norm type for convergence test 225c4762a1bSJed Brown PC Object: (mg_levels_3_) 6 MPI processes 226c4762a1bSJed Brown type: bjacobi 227c4762a1bSJed Brown number of blocks = 6 228*8d76b567SPatrick Sanan Local solver is the same for all blocks, as in the following KSP and PC objects on rank 0: 229c4762a1bSJed Brown KSP Object: (mg_levels_3_sub_) 1 MPI processes 230c4762a1bSJed Brown type: preonly 231c4762a1bSJed Brown maximum iterations=10000, initial guess is zero 232c4762a1bSJed Brown tolerances: relative=1e-05, absolute=1e-50, divergence=10000. 233c4762a1bSJed Brown left preconditioning 234c4762a1bSJed Brown using NONE norm type for convergence test 235c4762a1bSJed Brown PC Object: (mg_levels_3_sub_) 1 MPI processes 236c4762a1bSJed Brown type: icc 237c4762a1bSJed Brown out-of-place factorization 238c4762a1bSJed Brown 0 levels of fill 239c4762a1bSJed Brown tolerance for zero pivot 2.22045e-14 240c4762a1bSJed Brown using Manteuffel shift [POSITIVE_DEFINITE] 241c4762a1bSJed Brown matrix ordering: natural 242c4762a1bSJed Brown factor fill ratio given 1., needed 1.00054 243c4762a1bSJed Brown Factored matrix follows: 244c4762a1bSJed Brown Mat Object: 1 MPI processes 245c4762a1bSJed Brown type: seqsbaij 246c4762a1bSJed Brown rows=384, cols=384, bs=2 247c4762a1bSJed Brown package used to perform factorization: petsc 248c4762a1bSJed Brown total: nonzeros=7428, allocated nonzeros=7428 249c4762a1bSJed Brown total number of mallocs used during MatSetValues calls=0 250c4762a1bSJed Brown block size is 2 251c4762a1bSJed Brown linear system matrix = precond matrix: 252c4762a1bSJed Brown Mat Object: 1 MPI processes 253c4762a1bSJed Brown type: seqsbaij 254c4762a1bSJed Brown rows=384, cols=384, bs=2 255c4762a1bSJed Brown total: nonzeros=7424, allocated nonzeros=7424 256c4762a1bSJed Brown total number of mallocs used during MatSetValues calls=0 257c4762a1bSJed Brown block size is 2 258c4762a1bSJed Brown linear system matrix = precond matrix: 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 Brown Up solver (post-smoother) same as down solver (pre-smoother) 265c4762a1bSJed Brown linear system matrix = precond matrix: 266c4762a1bSJed Brown Mat Object: 6 MPI processes 267c4762a1bSJed Brown type: mpisbaij 268c4762a1bSJed Brown rows=2048, cols=2048, bs=2 269c4762a1bSJed Brown total: nonzeros=48128, allocated nonzeros=48128 270c4762a1bSJed Brown total number of mallocs used during MatSetValues calls=0 271c4762a1bSJed BrownSolution statistics after solve: Full 272c4762a1bSJed BrownCONVERGED_FNORM_RELATIVE: Number of SNES iterations = 5, total linear iterations = 22 273c4762a1bSJed Brown|X|_2 564267. -187.04 <= u <= 23814.8 -979.938 <= v <= 5720.23 3.24766 <= c <= 23904.5 274c4762a1bSJed BrownSurface statistics: u in [-6.945709e+01, 2.381387e+04] mean 1.520632e+04 275c4762a1bSJed BrownGlobal eta range 2.50736e+09 to 9.2273e+12 converged range 2.50736e+09 to 3.8249e+10 276c4762a1bSJed BrownGlobal beta2 range 0. to 4.97763e+09 converged range 0. to 4.97763e+09 277