Residual norms for sys_fieldsplit_1_ solve. 0 KSP preconditioned resid norm 5.474119726817e+00 true resid norm 6.271263985588e-02 ||r(i)||/||b|| 1.000000000000e+00 1 KSP preconditioned resid norm 3.558305192777e+00 true resid norm 1.058973705920e-01 ||r(i)||/||b|| 1.688612867124e+00 2 KSP preconditioned resid norm 3.298100273162e+00 true resid norm 1.156783348037e-01 ||r(i)||/||b|| 1.844577665197e+00 3 KSP preconditioned resid norm 3.204506681182e+00 true resid norm 1.224321236388e-01 ||r(i)||/||b|| 1.952271885224e+00 4 KSP preconditioned resid norm 3.200617351529e+00 true resid norm 1.226723763792e-01 ||r(i)||/||b|| 1.956102895064e+00 5 KSP preconditioned resid norm 3.193899612077e+00 true resid norm 1.230768731373e-01 ||r(i)||/||b|| 1.962552898748e+00 6 KSP preconditioned resid norm 3.192284751425e+00 true resid norm 1.230987039998e-01 ||r(i)||/||b|| 1.962901008197e+00 7 KSP preconditioned resid norm 3.182046417643e+00 true resid norm 1.251415509775e-01 ||r(i)||/||b|| 1.995475732884e+00 8 KSP preconditioned resid norm 3.175321474925e+00 true resid norm 1.253310625779e-01 ||r(i)||/||b|| 1.998497637254e+00 9 KSP preconditioned resid norm 3.168533786026e+00 true resid norm 1.264431488479e-01 ||r(i)||/||b|| 2.016230685529e+00 10 KSP preconditioned resid norm 3.168532845686e+00 true resid norm 1.264402013564e-01 ||r(i)||/||b|| 2.016183685569e+00 11 KSP preconditioned resid norm 3.165215670685e+00 true resid norm 1.274314672561e-01 ||r(i)||/||b|| 2.031990162572e+00 12 KSP preconditioned resid norm 3.147302100681e+00 true resid norm 1.334750838443e-01 ||r(i)||/||b|| 2.128360154365e+00 13 KSP preconditioned resid norm 3.147025696521e+00 true resid norm 1.335181380413e-01 ||r(i)||/||b|| 2.129046685773e+00 14 KSP preconditioned resid norm 3.144679052588e+00 true resid norm 1.342410792822e-01 ||r(i)||/||b|| 2.140574525179e+00 15 KSP preconditioned resid norm 3.143606723134e+00 true resid norm 1.338179140593e-01 ||r(i)||/||b|| 2.133826838844e+00 16 KSP preconditioned resid norm 3.141865092040e+00 true resid norm 1.344693829147e-01 ||r(i)||/||b|| 2.144214997546e+00 17 KSP preconditioned resid norm 3.141055662337e+00 true resid norm 1.349063443413e-01 ||r(i)||/||b|| 2.151182674678e+00 18 KSP preconditioned resid norm 3.140962810074e+00 true resid norm 1.346244714740e-01 ||r(i)||/||b|| 2.146688000751e+00 19 KSP preconditioned resid norm 3.140772607271e+00 true resid norm 1.352717996624e-01 ||r(i)||/||b|| 2.157010133416e+00 20 KSP preconditioned resid norm 3.139130384488e+00 true resid norm 1.332379377142e-01 ||r(i)||/||b|| 2.124578681751e+00 21 KSP preconditioned resid norm 3.139059160909e+00 true resid norm 1.346109015609e-01 ||r(i)||/||b|| 2.146471618325e+00 22 KSP preconditioned resid norm 3.133748617091e+00 true resid norm 4.191669300537e-01 ||r(i)||/||b|| 6.683930560362e+00 23 KSP preconditioned resid norm 3.133142611227e+00 true resid norm 1.346944344540e-01 ||r(i)||/||b|| 2.147803612854e+00 24 KSP preconditioned resid norm 3.132493649959e+00 true resid norm 1.771555031587e+00 ||r(i)||/||b|| 2.824877147029e+01 25 KSP preconditioned resid norm 3.130056283425e+00 true resid norm 1.741340963848e+01 ||r(i)||/||b|| 2.776698553673e+02 26 KSP preconditioned resid norm 3.128377028918e+00 true resid norm 2.525433101112e+01 ||r(i)||/||b|| 4.026992177201e+02 27 KSP preconditioned resid norm 3.128181250784e+00 true resid norm 8.709449961570e+01 ||r(i)||/||b|| 1.388787010336e+03 28 KSP preconditioned resid norm 3.125904782721e+00 true resid norm 9.807951100693e+02 ||r(i)||/||b|| 1.563951242243e+04 29 KSP preconditioned resid norm 3.125382631664e+00 true resid norm 2.585740254134e+03 ||r(i)||/||b|| 4.123156448327e+04 30 KSP preconditioned resid norm 4.740936896249e+05 true resid norm 1.566479652926e+04 ||r(i)||/||b|| 2.497869100274e+05 KSP Object:(sys_) 4 MPI processes type: fgmres GMRES: restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement GMRES: happy breakdown tolerance 1e-30 maximum iterations=1, initial guess is zero tolerances: relative=0.01, absolute=1e-50, divergence=10000. right preconditioning using UNPRECONDITIONED norm type for convergence test PC Object:(sys_) 4 MPI processes type: fieldsplit FieldSplit with Schur preconditioner, factorization UPPER Preconditioner for the Schur complement formed from Sp, an assembled approximation to S, which uses (lumped, if requested) A00's diagonal's inverse Split info: Split number 0 Defined by IS Split number 1 Defined by IS KSP solver for A00 block KSP Object: (sys_fieldsplit_0_) 4 MPI processes type: gmres GMRES: restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement GMRES: happy breakdown tolerance 1e-30 maximum iterations=10000, initial guess is zero tolerances: relative=1e-06, absolute=1e-50, divergence=10000. right preconditioning using UNPRECONDITIONED norm type for convergence test PC Object: (sys_fieldsplit_0_) 4 MPI processes type: bjacobi block Jacobi: number of blocks = 4 Local solve is same for all blocks, in the following KSP and PC objects: KSP Object: (sys_fieldsplit_0_sub_) 1 MPI processes type: preonly maximum iterations=10000, initial guess is zero tolerances: relative=1e-05, absolute=1e-50, divergence=10000. left preconditioning using NONE norm type for convergence test PC Object: (sys_fieldsplit_0_sub_) 1 MPI processes type: ilu ILU: out-of-place factorization 0 levels of fill tolerance for zero pivot 2.22045e-14 matrix ordering: natural factor fill ratio given 1., needed 1. Factored matrix follows: Mat Object: 1 MPI processes type: seqaij rows=192, cols=192 package used to perform factorization: petsc total: nonzeros=864, allocated nonzeros=864 total number of mallocs used during MatSetValues calls =0 not using I-node routines linear system matrix = precond matrix: Mat Object: (sys_fieldsplit_0_) 1 MPI processes type: seqaij rows=192, cols=192 total: nonzeros=864, allocated nonzeros=864 total number of mallocs used during MatSetValues calls =0 not using I-node routines linear system matrix = precond matrix: Mat Object: (sys_fieldsplit_0_) 4 MPI processes type: mpiaij rows=768, cols=768 total: nonzeros=3648, allocated nonzeros=3648 total number of mallocs used during MatSetValues calls =0 not using I-node (on process 0) routines KSP solver for upper A00 in upper triangular factor KSP Object: (sys_fieldsplit_1_upper_) 4 MPI processes type: preonly maximum iterations=10000, initial guess is zero tolerances: relative=1e-05, absolute=1e-50, divergence=10000. left preconditioning using DEFAULT norm type for convergence test PC Object: (sys_fieldsplit_1_upper_) 4 MPI processes type: jacobi PC has not been set up so information may be incomplete linear system matrix = precond matrix: Mat Object: (sys_fieldsplit_0_) 4 MPI processes type: mpiaij rows=768, cols=768 total: nonzeros=3648, allocated nonzeros=3648 total number of mallocs used during MatSetValues calls =0 not using I-node (on process 0) routines KSP solver for S = A11 - A10 inv(A00) A01 KSP Object: (sys_fieldsplit_1_) 4 MPI processes type: gmres GMRES: restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement GMRES: happy breakdown tolerance 1e-30 maximum iterations=10000, initial guess is zero tolerances: relative=1e-06, absolute=1e-50, divergence=10000. left preconditioning using PRECONDITIONED norm type for convergence test PC Object: (sys_fieldsplit_1_) 4 MPI processes type: bjacobi block Jacobi: number of blocks = 4 Local solve is same for all blocks, in the following KSP and PC objects: KSP Object: (sys_fieldsplit_1_sub_) 1 MPI processes type: preonly maximum iterations=10000, initial guess is zero tolerances: relative=1e-05, absolute=1e-50, divergence=10000. left preconditioning using NONE norm type for convergence test PC Object: (sys_fieldsplit_1_sub_) 1 MPI processes type: ilu ILU: out-of-place factorization 0 levels of fill tolerance for zero pivot 2.22045e-14 matrix ordering: natural factor fill ratio given 1., needed 1. Factored matrix follows: Mat Object: 1 MPI processes type: seqaij rows=64, cols=64 package used to perform factorization: petsc total: nonzeros=676, allocated nonzeros=676 total number of mallocs used during MatSetValues calls =0 not using I-node routines linear system matrix = precond matrix: Mat Object: 1 MPI processes type: seqaij rows=64, cols=64 total: nonzeros=676, allocated nonzeros=676 total number of mallocs used during MatSetValues calls =0 not using I-node routines linear system matrix followed by preconditioner matrix: Mat Object: (sys_fieldsplit_1_) 4 MPI processes type: schurcomplement rows=256, cols=256 Schur complement A11 - A10 inv(A00) A01 A11 Mat Object: (sys_fieldsplit_1_) 4 MPI processes type: mpiaij rows=256, cols=256 total: nonzeros=1216, allocated nonzeros=1216 total number of mallocs used during MatSetValues calls =0 not using I-node (on process 0) routines A10 Mat Object: (div_) 4 MPI processes type: mpiaij rows=256, cols=768 total: nonzeros=3648, allocated nonzeros=3648 total number of mallocs used during MatSetValues calls =0 not using I-node (on process 0) routines KSP of A00 KSP Object: (sys_fieldsplit_1_inner_) 4 MPI processes type: gmres GMRES: restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement GMRES: happy breakdown tolerance 1e-30 maximum iterations=10000, initial guess is zero tolerances: relative=1e-06, absolute=1e-50, divergence=10000. right preconditioning using UNPRECONDITIONED norm type for convergence test PC Object: (sys_fieldsplit_1_inner_) 4 MPI processes type: bjacobi block Jacobi: number of blocks = 4 Local solve is same for all blocks, in the following KSP and PC objects: KSP Object: (sys_fieldsplit_1_inner_sub_) 1 MPI processes type: preonly maximum iterations=10000, initial guess is zero tolerances: relative=1e-05, absolute=1e-50, divergence=10000. left preconditioning using NONE norm type for convergence test PC Object: (sys_fieldsplit_1_inner_sub_) 1 MPI processes type: ilu ILU: out-of-place factorization 0 levels of fill tolerance for zero pivot 2.22045e-14 matrix ordering: natural factor fill ratio given 1., needed 1. Factored matrix follows: Mat Object: 1 MPI processes type: seqaij rows=192, cols=192 package used to perform factorization: petsc total: nonzeros=864, allocated nonzeros=864 total number of mallocs used during MatSetValues calls =0 not using I-node routines linear system matrix = precond matrix: Mat Object: (sys_fieldsplit_0_) 1 MPI processes type: seqaij rows=192, cols=192 total: nonzeros=864, allocated nonzeros=864 total number of mallocs used during MatSetValues calls =0 not using I-node routines linear system matrix = precond matrix: Mat Object: (sys_fieldsplit_0_) 4 MPI processes type: mpiaij rows=768, cols=768 total: nonzeros=3648, allocated nonzeros=3648 total number of mallocs used during MatSetValues calls =0 not using I-node (on process 0) routines A01 Mat Object: (grad_) 4 MPI processes type: mpiaij rows=768, cols=256 total: nonzeros=3648, allocated nonzeros=3648 total number of mallocs used during MatSetValues calls =0 not using I-node (on process 0) routines Mat Object: 4 MPI processes type: mpiaij rows=256, cols=256 total: nonzeros=3012, allocated nonzeros=3012 total number of mallocs used during MatSetValues calls =0 not using I-node (on process 0) routines linear system matrix = precond matrix: Mat Object: 4 MPI processes type: nest rows=1024, cols=1024 has attached null space Matrix object: type=nest, rows=2, cols=2 MatNest structure: (0,0) : prefix="sys_fieldsplit_0_", type=mpiaij, rows=768, cols=768 (0,1) : prefix="grad_", type=mpiaij, rows=768, cols=256 (1,0) : prefix="div_", type=mpiaij, rows=256, cols=768 (1,1) : prefix="sys_fieldsplit_1_", type=mpiaij, rows=256, cols=256