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 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 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 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 91c4762a1bSJed Brown using NONE norm type for convergence test 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 99c4762a1bSJed Brown Factored matrix follows: 100*26cc229bSBarry 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 106c4762a1bSJed Brown linear system matrix = precond matrix: 107*26cc229bSBarry 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 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 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 138c4762a1bSJed Brown using NONE norm type for convergence test 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 145c4762a1bSJed Brown Factored matrix follows: 146*26cc229bSBarry 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 151c4762a1bSJed Brown block size is 2 152c4762a1bSJed Brown linear system matrix = precond matrix: 1538cc725e6SPierre Jolivet Mat Object: (mg_levels_1_sub_) 1 MPI process 154c4762a1bSJed Brown type: seqsbaij 155c4762a1bSJed Brown rows=48, cols=48, bs=2 156c4762a1bSJed Brown total: nonzeros=608, allocated nonzeros=608 157c4762a1bSJed Brown total number of mallocs used during MatSetValues calls=0 158c4762a1bSJed Brown block size is 2 159c4762a1bSJed Brown linear system matrix = precond matrix: 160c4762a1bSJed Brown Mat Object: 6 MPI processes 161c4762a1bSJed Brown type: mpisbaij 162c4762a1bSJed Brown rows=256, cols=256, bs=2 163c4762a1bSJed Brown total: nonzeros=4864, allocated nonzeros=4864 164c4762a1bSJed Brown total number of mallocs used during MatSetValues calls=0 165c4762a1bSJed Brown Up solver (post-smoother) same as down solver (pre-smoother) 166c4762a1bSJed Brown Down solver (pre-smoother) on level 2 ------------------------------- 167c4762a1bSJed Brown KSP Object: (mg_levels_2_) 6 MPI processes 168c4762a1bSJed Brown type: gmres 169c4762a1bSJed Brown restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement 170c4762a1bSJed Brown happy breakdown tolerance 1e-30 171c4762a1bSJed Brown maximum iterations=1, 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_2_) 6 MPI processes 176c4762a1bSJed Brown type: bjacobi 177c4762a1bSJed Brown number of blocks = 6 178020d6619SPierre Jolivet Local solver information for first block is in the following KSP and PC objects on rank 0: 179020d6619SPierre Jolivet Use -mg_levels_2_ksp_view ::ascii_info_detail to display information for all blocks 1808cc725e6SPierre Jolivet KSP Object: (mg_levels_2_sub_) 1 MPI process 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 1868cc725e6SPierre Jolivet PC Object: (mg_levels_2_sub_) 1 MPI process 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: 195*26cc229bSBarry Smith Mat Object: (mg_levels_2_sub_) 1 MPI process 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 block size is 2 201c4762a1bSJed Brown linear system matrix = precond matrix: 2028cc725e6SPierre Jolivet Mat Object: (mg_levels_2_sub_) 1 MPI process 203c4762a1bSJed Brown type: seqsbaij 204c4762a1bSJed Brown rows=192, cols=192, bs=2 205c4762a1bSJed Brown total: nonzeros=3008, allocated nonzeros=3008 206c4762a1bSJed Brown total number of mallocs used during MatSetValues calls=0 207c4762a1bSJed Brown block size is 2 208c4762a1bSJed Brown linear system matrix = precond matrix: 209c4762a1bSJed Brown Mat Object: 6 MPI processes 210c4762a1bSJed Brown type: mpisbaij 211c4762a1bSJed Brown rows=1024, cols=1024, bs=2 212c4762a1bSJed Brown total: nonzeros=19456, allocated nonzeros=19456 213c4762a1bSJed Brown total number of mallocs used during MatSetValues calls=0 214c4762a1bSJed Brown Up solver (post-smoother) same as down solver (pre-smoother) 215c4762a1bSJed Brown Down solver (pre-smoother) on level 3 ------------------------------- 216c4762a1bSJed Brown KSP Object: (mg_levels_3_) 6 MPI processes 217c4762a1bSJed Brown type: gmres 218c4762a1bSJed Brown restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement 219c4762a1bSJed Brown happy breakdown tolerance 1e-30 220c4762a1bSJed Brown maximum iterations=1, nonzero initial guess 221c4762a1bSJed Brown tolerances: relative=1e-05, absolute=1e-50, divergence=10000. 222c4762a1bSJed Brown left preconditioning 223c4762a1bSJed Brown using NONE norm type for convergence test 224c4762a1bSJed Brown PC Object: (mg_levels_3_) 6 MPI processes 225c4762a1bSJed Brown type: bjacobi 226c4762a1bSJed Brown number of blocks = 6 227020d6619SPierre Jolivet Local solver information for first block is in the following KSP and PC objects on rank 0: 228020d6619SPierre Jolivet Use -mg_levels_3_ksp_view ::ascii_info_detail to display information for all blocks 2298cc725e6SPierre Jolivet KSP Object: (mg_levels_3_sub_) 1 MPI process 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 2358cc725e6SPierre Jolivet PC Object: (mg_levels_3_sub_) 1 MPI process 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: 244*26cc229bSBarry Smith Mat Object: (mg_levels_3_sub_) 1 MPI process 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 block size is 2 250c4762a1bSJed Brown linear system matrix = precond matrix: 2518cc725e6SPierre Jolivet Mat Object: (mg_levels_3_sub_) 1 MPI process 252c4762a1bSJed Brown type: seqsbaij 253c4762a1bSJed Brown rows=384, cols=384, bs=2 254c4762a1bSJed Brown total: nonzeros=7424, allocated nonzeros=7424 255c4762a1bSJed Brown total number of mallocs used during MatSetValues calls=0 256c4762a1bSJed Brown block size is 2 257c4762a1bSJed Brown linear system matrix = precond matrix: 258c4762a1bSJed Brown Mat Object: 6 MPI processes 259c4762a1bSJed Brown type: mpisbaij 260c4762a1bSJed Brown rows=2048, cols=2048, bs=2 261c4762a1bSJed Brown total: nonzeros=48128, allocated nonzeros=48128 262c4762a1bSJed Brown total number of mallocs used during MatSetValues calls=0 263c4762a1bSJed Brown Up solver (post-smoother) same as down solver (pre-smoother) 264c4762a1bSJed Brown linear system matrix = precond matrix: 265c4762a1bSJed Brown Mat Object: 6 MPI processes 266c4762a1bSJed Brown type: mpisbaij 267c4762a1bSJed Brown rows=2048, cols=2048, bs=2 268c4762a1bSJed Brown total: nonzeros=48128, allocated nonzeros=48128 269c4762a1bSJed Brown total number of mallocs used during MatSetValues calls=0 270c4762a1bSJed BrownSolution statistics after solve: Full 271c4762a1bSJed BrownCONVERGED_FNORM_RELATIVE: Number of SNES iterations = 5, total linear iterations = 22 272c4762a1bSJed Brown|X|_2 564267. -187.04 <= u <= 23814.8 -979.938 <= v <= 5720.23 3.24766 <= c <= 23904.5 273c4762a1bSJed BrownSurface statistics: u in [-6.945709e+01, 2.381387e+04] mean 1.520632e+04 274c4762a1bSJed BrownGlobal eta range 2.50736e+09 to 9.2273e+12 converged range 2.50736e+09 to 3.8249e+10 275c4762a1bSJed BrownGlobal beta2 range 0. to 4.97763e+09 converged range 0. to 4.97763e+09 276