1*c4762a1bSJed 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. 2*c4762a1bSJed Brown Linear solve converged due to CONVERGED_RTOL iterations 1 3*c4762a1bSJed Brown Linear solve converged due to CONVERGED_RTOL iterations 1 4*c4762a1bSJed Brown Linear solve converged due to CONVERGED_RTOL iterations 1 5*c4762a1bSJed Brown Linear solve converged due to CONVERGED_RTOL iterations 1 6*c4762a1bSJed Brown Linear solve converged due to CONVERGED_RTOL iterations 1 7*c4762a1bSJed Brown Linear solve converged due to CONVERGED_RTOL iterations 1 8*c4762a1bSJed Brown Linear solve converged due to CONVERGED_RTOL iterations 1 9*c4762a1bSJed Brown Linear solve converged due to CONVERGED_RTOL iterations 1 10*c4762a1bSJed Brown Linear solve converged due to CONVERGED_RTOL iterations 1 11*c4762a1bSJed Brown Nonlinear solve converged due to CONVERGED_FNORM_RELATIVE iterations 9 12*c4762a1bSJed 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. 13*c4762a1bSJed 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. 14*c4762a1bSJed Brown Linear solve converged due to CONVERGED_RTOL iterations 4 15*c4762a1bSJed Brown Linear solve converged due to CONVERGED_RTOL iterations 5 16*c4762a1bSJed Brown Linear solve converged due to CONVERGED_RTOL iterations 5 17*c4762a1bSJed Brown Linear solve converged due to CONVERGED_RTOL iterations 6 18*c4762a1bSJed Brown Linear solve converged due to CONVERGED_RTOL iterations 6 19*c4762a1bSJed Brown Linear solve converged due to CONVERGED_RTOL iterations 5 20*c4762a1bSJed Brown Linear solve converged due to CONVERGED_RTOL iterations 6 21*c4762a1bSJed Brown Linear solve converged due to CONVERGED_RTOL iterations 7 22*c4762a1bSJed Brown Nonlinear solve converged due to CONVERGED_FNORM_RELATIVE iterations 8 23*c4762a1bSJed 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. 24*c4762a1bSJed 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. 25*c4762a1bSJed 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. 26*c4762a1bSJed Brown Linear solve converged due to CONVERGED_RTOL iterations 5 27*c4762a1bSJed Brown Linear solve converged due to CONVERGED_RTOL iterations 6 28*c4762a1bSJed Brown Linear solve converged due to CONVERGED_RTOL iterations 7 29*c4762a1bSJed Brown Linear solve converged due to CONVERGED_RTOL iterations 6 30*c4762a1bSJed Brown Linear solve converged due to CONVERGED_RTOL iterations 2 31*c4762a1bSJed Brown Linear solve converged due to CONVERGED_RTOL iterations 2 32*c4762a1bSJed Brown Linear solve converged due to CONVERGED_RTOL iterations 3 33*c4762a1bSJed Brown Linear solve converged due to CONVERGED_RTOL iterations 5 34*c4762a1bSJed Brown Linear solve converged due to CONVERGED_RTOL iterations 6 35*c4762a1bSJed Brown Linear solve converged due to CONVERGED_RTOL iterations 7 36*c4762a1bSJed Brown Nonlinear solve converged due to CONVERGED_FNORM_RELATIVE iterations 10 37*c4762a1bSJed 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 38*c4762a1bSJed 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. 39*c4762a1bSJed 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. 40*c4762a1bSJed 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. 41*c4762a1bSJed Brown Linear solve converged due to CONVERGED_RTOL iterations 2 42*c4762a1bSJed Brown Linear solve converged due to CONVERGED_RTOL iterations 2 43*c4762a1bSJed Brown Linear solve converged due to CONVERGED_RTOL iterations 4 44*c4762a1bSJed Brown Linear solve converged due to CONVERGED_RTOL iterations 7 45*c4762a1bSJed Brown Linear solve converged due to CONVERGED_RTOL iterations 7 46*c4762a1bSJed BrownNonlinear solve converged due to CONVERGED_FNORM_RELATIVE iterations 5 47*c4762a1bSJed BrownSNES Object: 6 MPI processes 48*c4762a1bSJed Brown type: newtonls 49*c4762a1bSJed Brown maximum iterations=50, maximum function evaluations=10000 50*c4762a1bSJed Brown tolerances: relative=1e-08, absolute=1e-50, solution=1e-12 51*c4762a1bSJed Brown total number of linear solver iterations=22 52*c4762a1bSJed Brown total number of function evaluations=6 53*c4762a1bSJed Brown norm schedule ALWAYS 54*c4762a1bSJed Brown total number of grid sequence refinements=3 55*c4762a1bSJed Brown Jacobian is built using a DMDA local Jacobian 56*c4762a1bSJed Brown SNESLineSearch Object: 6 MPI processes 57*c4762a1bSJed Brown type: bt 58*c4762a1bSJed Brown interpolation: cubic 59*c4762a1bSJed Brown alpha=1.000000e-04 60*c4762a1bSJed Brown maxstep=1.000000e+08, minlambda=1.000000e-12 61*c4762a1bSJed Brown tolerances: relative=1.000000e-08, absolute=1.000000e-15, lambda=1.000000e-08 62*c4762a1bSJed Brown maximum iterations=40 63*c4762a1bSJed Brown KSP Object: 6 MPI processes 64*c4762a1bSJed Brown type: fgmres 65*c4762a1bSJed Brown restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement 66*c4762a1bSJed Brown happy breakdown tolerance 1e-30 67*c4762a1bSJed Brown maximum iterations=10000, initial guess is zero 68*c4762a1bSJed Brown tolerances: relative=0.01, absolute=-1., divergence=10000. 69*c4762a1bSJed Brown right preconditioning 70*c4762a1bSJed Brown using UNPRECONDITIONED norm type for convergence test 71*c4762a1bSJed Brown PC Object: 6 MPI processes 72*c4762a1bSJed Brown type: mg 73*c4762a1bSJed Brown type is MULTIPLICATIVE, levels=4 cycles=v 74*c4762a1bSJed Brown Cycles per PCApply=1 75*c4762a1bSJed Brown Not using Galerkin computed coarse grid matrices 76*c4762a1bSJed Brown Coarse grid solver -- level ------------------------------- 77*c4762a1bSJed Brown KSP Object: (mg_coarse_) 6 MPI processes 78*c4762a1bSJed Brown type: preonly 79*c4762a1bSJed Brown maximum iterations=10000, initial guess is zero 80*c4762a1bSJed Brown tolerances: relative=1e-05, absolute=1e-50, divergence=10000. 81*c4762a1bSJed Brown left preconditioning 82*c4762a1bSJed Brown using NONE norm type for convergence test 83*c4762a1bSJed Brown PC Object: (mg_coarse_) 6 MPI processes 84*c4762a1bSJed Brown type: redundant 85*c4762a1bSJed Brown First (color=0) of 6 PCs follows 86*c4762a1bSJed Brown KSP Object: (mg_coarse_redundant_) 1 MPI processes 87*c4762a1bSJed Brown type: preonly 88*c4762a1bSJed Brown maximum iterations=10000, initial guess is zero 89*c4762a1bSJed Brown tolerances: relative=1e-05, absolute=1e-50, divergence=10000. 90*c4762a1bSJed Brown left preconditioning 91*c4762a1bSJed Brown using NONE norm type for convergence test 92*c4762a1bSJed Brown PC Object: (mg_coarse_redundant_) 1 MPI processes 93*c4762a1bSJed Brown type: lu 94*c4762a1bSJed Brown out-of-place factorization 95*c4762a1bSJed Brown tolerance for zero pivot 2.22045e-14 96*c4762a1bSJed Brown using diagonal shift on blocks to prevent zero pivot [INBLOCKS] 97*c4762a1bSJed Brown matrix ordering: nd 98*c4762a1bSJed Brown factor fill ratio given 5., needed 1.41667 99*c4762a1bSJed Brown Factored matrix follows: 100*c4762a1bSJed Brown Mat Object: 1 MPI processes 101*c4762a1bSJed Brown type: seqaij 102*c4762a1bSJed Brown rows=64, cols=64, bs=2 103*c4762a1bSJed Brown package used to perform factorization: petsc 104*c4762a1bSJed Brown total: nonzeros=3264, allocated nonzeros=3264 105*c4762a1bSJed Brown total number of mallocs used during MatSetValues calls=0 106*c4762a1bSJed Brown using I-node routines: found 16 nodes, limit used is 5 107*c4762a1bSJed Brown linear system matrix = precond matrix: 108*c4762a1bSJed Brown Mat Object: 1 MPI processes 109*c4762a1bSJed Brown type: seqaij 110*c4762a1bSJed Brown rows=64, cols=64, bs=2 111*c4762a1bSJed Brown total: nonzeros=2304, allocated nonzeros=2304 112*c4762a1bSJed Brown total number of mallocs used during MatSetValues calls=0 113*c4762a1bSJed Brown using I-node routines: found 16 nodes, limit used is 5 114*c4762a1bSJed Brown linear system matrix = precond matrix: 115*c4762a1bSJed Brown Mat Object: 6 MPI processes 116*c4762a1bSJed Brown type: mpiaij 117*c4762a1bSJed Brown rows=64, cols=64, bs=2 118*c4762a1bSJed Brown total: nonzeros=2304, allocated nonzeros=2304 119*c4762a1bSJed Brown total number of mallocs used during MatSetValues calls=0 120*c4762a1bSJed Brown Down solver (pre-smoother) on level 1 ------------------------------- 121*c4762a1bSJed Brown KSP Object: (mg_levels_1_) 6 MPI processes 122*c4762a1bSJed Brown type: gmres 123*c4762a1bSJed Brown restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement 124*c4762a1bSJed Brown happy breakdown tolerance 1e-30 125*c4762a1bSJed Brown maximum iterations=1, nonzero initial guess 126*c4762a1bSJed Brown tolerances: relative=1e-05, absolute=1e-50, divergence=10000. 127*c4762a1bSJed Brown left preconditioning 128*c4762a1bSJed Brown using NONE norm type for convergence test 129*c4762a1bSJed Brown PC Object: (mg_levels_1_) 6 MPI processes 130*c4762a1bSJed Brown type: bjacobi 131*c4762a1bSJed Brown number of blocks = 6 132*c4762a1bSJed Brown Local solve is same for all blocks, in the following KSP and PC objects: 133*c4762a1bSJed Brown KSP Object: (mg_levels_1_sub_) 1 MPI processes 134*c4762a1bSJed Brown type: preonly 135*c4762a1bSJed Brown maximum iterations=10000, initial guess is zero 136*c4762a1bSJed Brown tolerances: relative=1e-05, absolute=1e-50, divergence=10000. 137*c4762a1bSJed Brown left preconditioning 138*c4762a1bSJed Brown using NONE norm type for convergence test 139*c4762a1bSJed Brown PC Object: (mg_levels_1_sub_) 1 MPI processes 140*c4762a1bSJed Brown type: cholesky 141*c4762a1bSJed Brown out-of-place factorization 142*c4762a1bSJed Brown tolerance for zero pivot 2.22045e-14 143*c4762a1bSJed Brown matrix ordering: natural 144*c4762a1bSJed Brown factor fill ratio given 5., needed 1.32237 145*c4762a1bSJed Brown Factored matrix follows: 146*c4762a1bSJed Brown Mat Object: 1 MPI processes 147*c4762a1bSJed Brown type: seqsbaij 148*c4762a1bSJed Brown rows=48, cols=48, bs=2 149*c4762a1bSJed Brown package used to perform factorization: petsc 150*c4762a1bSJed Brown total: nonzeros=804, allocated nonzeros=804 151*c4762a1bSJed Brown total number of mallocs used during MatSetValues calls=0 152*c4762a1bSJed Brown block size is 2 153*c4762a1bSJed Brown linear system matrix = precond matrix: 154*c4762a1bSJed Brown Mat Object: 1 MPI processes 155*c4762a1bSJed Brown type: seqsbaij 156*c4762a1bSJed Brown rows=48, cols=48, bs=2 157*c4762a1bSJed Brown total: nonzeros=608, allocated nonzeros=608 158*c4762a1bSJed Brown total number of mallocs used during MatSetValues calls=0 159*c4762a1bSJed Brown block size is 2 160*c4762a1bSJed Brown linear system matrix = precond matrix: 161*c4762a1bSJed Brown Mat Object: 6 MPI processes 162*c4762a1bSJed Brown type: mpisbaij 163*c4762a1bSJed Brown rows=256, cols=256, bs=2 164*c4762a1bSJed Brown total: nonzeros=4864, allocated nonzeros=4864 165*c4762a1bSJed Brown total number of mallocs used during MatSetValues calls=0 166*c4762a1bSJed Brown Up solver (post-smoother) same as down solver (pre-smoother) 167*c4762a1bSJed Brown Down solver (pre-smoother) on level 2 ------------------------------- 168*c4762a1bSJed Brown KSP Object: (mg_levels_2_) 6 MPI processes 169*c4762a1bSJed Brown type: gmres 170*c4762a1bSJed Brown restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement 171*c4762a1bSJed Brown happy breakdown tolerance 1e-30 172*c4762a1bSJed Brown maximum iterations=1, nonzero initial guess 173*c4762a1bSJed Brown tolerances: relative=1e-05, absolute=1e-50, divergence=10000. 174*c4762a1bSJed Brown left preconditioning 175*c4762a1bSJed Brown using NONE norm type for convergence test 176*c4762a1bSJed Brown PC Object: (mg_levels_2_) 6 MPI processes 177*c4762a1bSJed Brown type: bjacobi 178*c4762a1bSJed Brown number of blocks = 6 179*c4762a1bSJed Brown Local solve is same for all blocks, in the following KSP and PC objects: 180*c4762a1bSJed Brown KSP Object: (mg_levels_2_sub_) 1 MPI processes 181*c4762a1bSJed Brown type: preonly 182*c4762a1bSJed Brown maximum iterations=10000, initial guess is zero 183*c4762a1bSJed Brown tolerances: relative=1e-05, absolute=1e-50, divergence=10000. 184*c4762a1bSJed Brown left preconditioning 185*c4762a1bSJed Brown using NONE norm type for convergence test 186*c4762a1bSJed Brown PC Object: (mg_levels_2_sub_) 1 MPI processes 187*c4762a1bSJed Brown type: icc 188*c4762a1bSJed Brown out-of-place factorization 189*c4762a1bSJed Brown 0 levels of fill 190*c4762a1bSJed Brown tolerance for zero pivot 2.22045e-14 191*c4762a1bSJed Brown using Manteuffel shift [POSITIVE_DEFINITE] 192*c4762a1bSJed Brown matrix ordering: natural 193*c4762a1bSJed Brown factor fill ratio given 1., needed 1.00133 194*c4762a1bSJed Brown Factored matrix follows: 195*c4762a1bSJed Brown Mat Object: 1 MPI processes 196*c4762a1bSJed Brown type: seqsbaij 197*c4762a1bSJed Brown rows=192, cols=192, bs=2 198*c4762a1bSJed Brown package used to perform factorization: petsc 199*c4762a1bSJed Brown total: nonzeros=3012, allocated nonzeros=3012 200*c4762a1bSJed Brown total number of mallocs used during MatSetValues calls=0 201*c4762a1bSJed Brown block size is 2 202*c4762a1bSJed Brown linear system matrix = precond matrix: 203*c4762a1bSJed Brown Mat Object: 1 MPI processes 204*c4762a1bSJed Brown type: seqsbaij 205*c4762a1bSJed Brown rows=192, cols=192, bs=2 206*c4762a1bSJed Brown total: nonzeros=3008, allocated nonzeros=3008 207*c4762a1bSJed Brown total number of mallocs used during MatSetValues calls=0 208*c4762a1bSJed Brown block size is 2 209*c4762a1bSJed Brown linear system matrix = precond matrix: 210*c4762a1bSJed Brown Mat Object: 6 MPI processes 211*c4762a1bSJed Brown type: mpisbaij 212*c4762a1bSJed Brown rows=1024, cols=1024, bs=2 213*c4762a1bSJed Brown total: nonzeros=19456, allocated nonzeros=19456 214*c4762a1bSJed Brown total number of mallocs used during MatSetValues calls=0 215*c4762a1bSJed Brown Up solver (post-smoother) same as down solver (pre-smoother) 216*c4762a1bSJed Brown Down solver (pre-smoother) on level 3 ------------------------------- 217*c4762a1bSJed Brown KSP Object: (mg_levels_3_) 6 MPI processes 218*c4762a1bSJed Brown type: gmres 219*c4762a1bSJed Brown restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement 220*c4762a1bSJed Brown happy breakdown tolerance 1e-30 221*c4762a1bSJed Brown maximum iterations=1, nonzero initial guess 222*c4762a1bSJed Brown tolerances: relative=1e-05, absolute=1e-50, divergence=10000. 223*c4762a1bSJed Brown left preconditioning 224*c4762a1bSJed Brown using NONE norm type for convergence test 225*c4762a1bSJed Brown PC Object: (mg_levels_3_) 6 MPI processes 226*c4762a1bSJed Brown type: bjacobi 227*c4762a1bSJed Brown number of blocks = 6 228*c4762a1bSJed Brown Local solve is same for all blocks, in the following KSP and PC objects: 229*c4762a1bSJed Brown KSP Object: (mg_levels_3_sub_) 1 MPI processes 230*c4762a1bSJed Brown type: preonly 231*c4762a1bSJed Brown maximum iterations=10000, initial guess is zero 232*c4762a1bSJed Brown tolerances: relative=1e-05, absolute=1e-50, divergence=10000. 233*c4762a1bSJed Brown left preconditioning 234*c4762a1bSJed Brown using NONE norm type for convergence test 235*c4762a1bSJed Brown PC Object: (mg_levels_3_sub_) 1 MPI processes 236*c4762a1bSJed Brown type: icc 237*c4762a1bSJed Brown out-of-place factorization 238*c4762a1bSJed Brown 0 levels of fill 239*c4762a1bSJed Brown tolerance for zero pivot 2.22045e-14 240*c4762a1bSJed Brown using Manteuffel shift [POSITIVE_DEFINITE] 241*c4762a1bSJed Brown matrix ordering: natural 242*c4762a1bSJed Brown factor fill ratio given 1., needed 1.00054 243*c4762a1bSJed Brown Factored matrix follows: 244*c4762a1bSJed Brown Mat Object: 1 MPI processes 245*c4762a1bSJed Brown type: seqsbaij 246*c4762a1bSJed Brown rows=384, cols=384, bs=2 247*c4762a1bSJed Brown package used to perform factorization: petsc 248*c4762a1bSJed Brown total: nonzeros=7428, allocated nonzeros=7428 249*c4762a1bSJed Brown total number of mallocs used during MatSetValues calls=0 250*c4762a1bSJed Brown block size is 2 251*c4762a1bSJed Brown linear system matrix = precond matrix: 252*c4762a1bSJed Brown Mat Object: 1 MPI processes 253*c4762a1bSJed Brown type: seqsbaij 254*c4762a1bSJed Brown rows=384, cols=384, bs=2 255*c4762a1bSJed Brown total: nonzeros=7424, allocated nonzeros=7424 256*c4762a1bSJed Brown total number of mallocs used during MatSetValues calls=0 257*c4762a1bSJed Brown block size is 2 258*c4762a1bSJed Brown linear system matrix = precond matrix: 259*c4762a1bSJed Brown Mat Object: 6 MPI processes 260*c4762a1bSJed Brown type: mpisbaij 261*c4762a1bSJed Brown rows=2048, cols=2048, bs=2 262*c4762a1bSJed Brown total: nonzeros=48128, allocated nonzeros=48128 263*c4762a1bSJed Brown total number of mallocs used during MatSetValues calls=0 264*c4762a1bSJed Brown Up solver (post-smoother) same as down solver (pre-smoother) 265*c4762a1bSJed Brown linear system matrix = precond matrix: 266*c4762a1bSJed Brown Mat Object: 6 MPI processes 267*c4762a1bSJed Brown type: mpisbaij 268*c4762a1bSJed Brown rows=2048, cols=2048, bs=2 269*c4762a1bSJed Brown total: nonzeros=48128, allocated nonzeros=48128 270*c4762a1bSJed Brown total number of mallocs used during MatSetValues calls=0 271*c4762a1bSJed BrownSolution statistics after solve: Full 272*c4762a1bSJed BrownCONVERGED_FNORM_RELATIVE: Number of SNES iterations = 5, total linear iterations = 22 273*c4762a1bSJed Brown|X|_2 564267. -187.04 <= u <= 23814.8 -979.938 <= v <= 5720.23 3.24766 <= c <= 23904.5 274*c4762a1bSJed BrownSurface statistics: u in [-6.945709e+01, 2.381387e+04] mean 1.520632e+04 275*c4762a1bSJed BrownGlobal eta range 2.50736e+09 to 9.2273e+12 converged range 2.50736e+09 to 3.8249e+10 276*c4762a1bSJed BrownGlobal beta2 range 0. to 4.97763e+09 converged range 0. to 4.97763e+09 277