Residual norms for sys_fieldsplit_1_ solve. 0 KSP preconditioned resid norm 2.651426524461e-01 true resid norm 1.947268006892e-02 ||r(i)||/||b|| 1.000000000000e+00 1 KSP preconditioned resid norm 8.503233836734e-02 true resid norm 5.617906297269e-03 ||r(i)||/||b|| 2.885019564531e-01 2 KSP preconditioned resid norm 3.238469701444e-02 true resid norm 9.800097635323e-04 ||r(i)||/||b|| 5.032742077946e-02 3 KSP preconditioned resid norm 1.669239032022e-02 true resid norm 6.260563099774e-04 ||r(i)||/||b|| 3.215049534844e-02 4 KSP preconditioned resid norm 9.076875144325e-03 true resid norm 3.153918045795e-04 ||r(i)||/||b|| 1.619663053381e-02 5 KSP preconditioned resid norm 4.653388201951e-03 true resid norm 1.779714224717e-04 ||r(i)||/||b|| 9.139544317565e-03 6 KSP preconditioned resid norm 1.825996916687e-03 true resid norm 8.495216428016e-05 ||r(i)||/||b|| 4.362633390960e-03 7 KSP preconditioned resid norm 5.466414574239e-04 true resid norm 3.537829111281e-05 ||r(i)||/||b|| 1.816816739534e-03 8 KSP preconditioned resid norm 1.683622226756e-04 true resid norm 1.019217974220e-05 ||r(i)||/||b|| 5.234091920641e-04 9 KSP preconditioned resid norm 6.251221637485e-05 true resid norm 2.953723588582e-06 ||r(i)||/||b|| 1.516855193085e-04 10 KSP preconditioned resid norm 2.077476920262e-05 true resid norm 1.304632634893e-06 ||r(i)||/||b|| 6.699810351096e-05 11 KSP preconditioned resid norm 9.124600530800e-06 true resid norm 4.218353839242e-07 ||r(i)||/||b|| 2.166293404047e-05 12 KSP preconditioned resid norm 3.166475857787e-06 true resid norm 1.789078563249e-07 ||r(i)||/||b|| 9.187633941076e-06 13 KSP preconditioned resid norm 1.308705119819e-06 true resid norm 7.703574807697e-08 ||r(i)||/||b|| 3.956093758246e-06 14 KSP preconditioned resid norm 7.249147427007e-07 true resid norm 2.657083328289e-08 ||r(i)||/||b|| 1.364518555682e-06 15 KSP preconditioned resid norm 4.284522202550e-07 true resid norm 2.050262990218e-08 ||r(i)||/||b|| 1.052892043089e-06 16 KSP preconditioned resid norm 1.914519076262e-07 true resid norm 1.265931988789e-08 ||r(i)||/||b|| 6.501067055530e-07 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 has attached null space 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 has attached null space 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 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