Re: [petsc-dev] Using PCFieldSplitSetIS

18 Mar 2011


      Barry Smith wrote:
...
Did you get a chance to run as I requested with "Please run with -ksp_view and send the output."
Yes, the problem with the second iteration occurs only when I make use 
of bcgs for the outer iteration (which is set in my code as the standard 
solver). When using richardson or gmres, I get the expected output. Here 
is the output for -ksp_type bcgs and -ksp_view

        Residual norms for fieldsplit_boundary_ solve.
    0 KSP Residual norm 1.238748554482e+02
    1 KSP Residual norm 1.712708691574e+01
    2 KSP Residual norm 5.537715643639e+00
    3 KSP Residual norm 2.767431394864e+00
    4 KSP Residual norm 2.308410861046e+00
    5 KSP Residual norm 1.225474327957e+00
    6 KSP Residual norm 4.705888493777e-01
    7 KSP Residual norm 4.476449896581e-01
    8 KSP Residual norm 1.859190943730e-01
    9 KSP Residual norm 9.038347953900e-02
   10 KSP Residual norm 7.545891428899e-02
   11 KSP Residual norm 4.046165237615e-02
   12 KSP Residual norm 1.499755811468e-02
   13 KSP Residual norm 1.100288579779e-02
   14 KSP Residual norm 4.530437449505e-03
   15 KSP Residual norm 2.511112676997e-03
   16 KSP Residual norm 1.425956836799e-03
   17 KSP Residual norm 1.008189124008e-03
   18 KSP Residual norm 4.759971096889e-04
   19 KSP Residual norm 2.418445744584e-04
   20 KSP Residual norm 1.345530122184e-04
   21 KSP Residual norm 6.618130426631e-05
   22 KSP Residual norm 3.502075214549e-05
   23 KSP Residual norm 2.093145547362e-05
   24 KSP Residual norm 1.426936785120e-05
   25 KSP Residual norm 6.716866486925e-06
   26 KSP Residual norm 2.938650130672e-06
   27 KSP Residual norm 1.194333200607e-06
   28 KSP Residual norm 8.117544753927e-07
   29 KSP Residual norm 4.436526804864e-07
   30 KSP Residual norm 2.619889310859e-07
   31 KSP Residual norm 2.608365614137e-07
   32 KSP Residual norm 2.502840136682e-07
   33 KSP Residual norm 2.108068726014e-07
   34 KSP Residual norm 1.190606661068e-07
   35 KSP Residual norm 8.587022449664e-08
   36 KSP Residual norm 2.657505256794e-08
   37 KSP Residual norm 1.091795047578e-08
   38 KSP Residual norm 5.622627482150e-09
  Linear solve converged due to CONVERGED_ATOL iterations 38
  0 KSP preconditioned resid norm 1.039480141214e+02 true resid norm 
2.318397033201e+03 ||r(i)||/||b|| 1.000000000000e+00
    Residual norms for fieldsplit_boundary_ solve.
    0 KSP Residual norm 1.238748554488e+02
    1 KSP Residual norm 1.712708691920e+01
    2 KSP Residual norm 5.537715644318e+00
    3 KSP Residual norm 2.767431394929e+00
    4 KSP Residual norm 2.308410861005e+00
    5 KSP Residual norm 1.225474328165e+00
    6 KSP Residual norm 4.705888492572e-01
    7 KSP Residual norm 4.476449894826e-01
    8 KSP Residual norm 1.859190943717e-01
    9 KSP Residual norm 9.038347954447e-02
   10 KSP Residual norm 7.545891428241e-02
   11 KSP Residual norm 4.046165236351e-02
   12 KSP Residual norm 1.499755811079e-02
   13 KSP Residual norm 1.100288580023e-02
   14 KSP Residual norm 4.530437448754e-03
   15 KSP Residual norm 2.511112675547e-03
   16 KSP Residual norm 1.425956832972e-03
   17 KSP Residual norm 1.008189121211e-03
   18 KSP Residual norm 4.759971083063e-04
   19 KSP Residual norm 2.418445729000e-04
   20 KSP Residual norm 1.345530095639e-04
   21 KSP Residual norm 6.618130177611e-05
   22 KSP Residual norm 3.502074795980e-05
   23 KSP Residual norm 2.093145116748e-05
   24 KSP Residual norm 1.426936277614e-05
   25 KSP Residual norm 6.716863232176e-06
   26 KSP Residual norm 2.938647588875e-06
   27 KSP Residual norm 1.194330143586e-06
   28 KSP Residual norm 8.117422448346e-07
   29 KSP Residual norm 4.436425441849e-07
   30 KSP Residual norm 2.619889078643e-07
   31 KSP Residual norm 2.608368186135e-07
   32 KSP Residual norm 2.502849238080e-07
   33 KSP Residual norm 2.108123763958e-07
   34 KSP Residual norm 1.190912989234e-07
   35 KSP Residual norm 8.589705060827e-08
   36 KSP Residual norm 2.657654667253e-08
   37 KSP Residual norm 1.091743530704e-08
   38 KSP Residual norm 5.623018769513e-09
  Linear solve converged due to CONVERGED_ATOL iterations 38
    Residual norms for fieldsplit_boundary_ solve.
    0 KSP Residual norm 5.628142396422e-09
  Linear solve converged due to CONVERGED_ATOL iterations 0
  1 KSP preconditioned resid norm 2.834168769772e-08 true resid norm 
5.737810813629e-09 ||r(i)||/||b|| 2.474904311669e-12
    Residual norms for fieldsplit_boundary_ solve.
    0 KSP Residual norm 5.626350752150e-09
  Linear solve converged due to CONVERGED_ATOL iterations 0
    Residual norms for fieldsplit_boundary_ solve.
    0 KSP Residual norm 5.628142396422e-09
  Linear solve converged due to CONVERGED_ATOL iterations 0
  2 KSP preconditioned resid norm 2.834287716574e-08 true resid norm 
1.921317696012e-06 ||r(i)||/||b|| 8.287267747921e-10
    Residual norms for fieldsplit_boundary_ solve.
    0 KSP Residual norm 6.360593934279e-07
    1 KSP Residual norm 6.207657513778e-07
    2 KSP Residual norm 4.595521952584e-07
    3 KSP Residual norm 2.353791066544e-07
    4 KSP Residual norm 1.659055771510e-07
    5 KSP Residual norm 1.057381480368e-07
    6 KSP Residual norm 2.733159212163e-08
    7 KSP Residual norm 1.685966550943e-08
    8 KSP Residual norm 5.250226264235e-09
  Linear solve converged due to CONVERGED_ATOL iterations 8
    Residual norms for fieldsplit_boundary_ solve.
    0 KSP Residual norm 5.314638497757e-10
  Linear solve converged due to CONVERGED_ATOL iterations 0
  3 KSP preconditioned resid norm 2.418607617684e-09 true resid norm 
1.943729846844e-06 ||r(i)||/||b|| 8.383938639538e-10
KSP Object:
  type: bcgs
  maximum iterations=10000, initial guess is zero
  tolerances:  relative=0, absolute=1e-08, divergence=10000
  left preconditioning
  using PRECONDITIONED norm type for convergence test
PC Object:
  type: fieldsplit
    FieldSplit with Schur preconditioner, blocksize = -1, factorization FULL
    Split info:
    Split number 0 Defined by IS
    Split number 1 Defined by IS
    KSP solver for A00 block
      KSP Object:      (fieldsplit_interior_)     
        type: preonly
        maximum iterations=10000, initial guess is zero
        tolerances:  relative=1e-05, absolute=1e-50, divergence=10000
        left preconditioning
        using PRECONDITIONED norm type for convergence test
      PC Object:      (fieldsplit_interior_)     
        type: bjacobi
          block Jacobi: number of blocks = 2
          Local solve is same for all blocks, in the following KSP and 
PC objects:
        KSP Object:        (fieldsplit_interior_sub_)       
          type: preonly
          maximum iterations=10000, initial guess is zero
          tolerances:  relative=1e-05, absolute=1e-50, divergence=10000
          left preconditioning
          using PRECONDITIONED norm type for convergence test
        PC Object:        (fieldsplit_interior_sub_)       
          type: lu
            LU: out-of-place factorization
            tolerance for zero pivot 1e-12
            matrix ordering: nd
            factor fill ratio given 0, needed 0
              Factored matrix follows:
                Matrix Object:               
                  type: seqaij
                  rows=12336, cols=12336
                  package used to perform factorization: umfpack
                  total: nonzeros=0, allocated nonzeros=0
                  total number of mallocs used during MatSetValues calls =0
                    not using I-node routines
                    UMFPACK run parameters:
                      Control[UMFPACK_PRL]: 1
                      Control[UMFPACK_STRATEGY]: 0
                      Control[UMFPACK_DENSE_COL]: 0.2
                      Control[UMFPACK_DENSE_ROW]: 0.2
                      Control[UMFPACK_AMD_DENSE]: 10
                      Control[UMFPACK_BLOCK_SIZE]: 32
                      Control[UMFPACK_2BY2_TOLERANCE]: 0.01
                      Control[UMFPACK_FIXQ]: 0
                      Control[UMFPACK_AGGRESSIVE]: 1
                      Control[UMFPACK_PIVOT_TOLERANCE]: 0.1
                      Control[UMFPACK_SYM_PIVOT_TOLERANCE]: 0.001
                      Control[UMFPACK_SCALE]: 1
                      Control[UMFPACK_ALLOC_INIT]: 0.7
                      Control[UMFPACK_DROPTOL]: 0
                      Control[UMFPACK_IRSTEP]: 0
                      UMFPACK default matrix ordering is used (not the 
PETSc matrix ordering)
          linear system matrix = precond matrix:
          Matrix Object:         
            type: seqaij
            rows=12336, cols=12336
            total: nonzeros=173608, allocated nonzeros=173608
            total number of mallocs used during MatSetValues calls =0
              not using I-node routines
        linear system matrix = precond matrix:
        Matrix Object:       
          type: mpiaij
          rows=24672, cols=24672
          total: nonzeros=347216, allocated nonzeros=347216
          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:      (fieldsplit_boundary_)     
        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-50, absolute=1e-08, divergence=10000
        left preconditioning
        using PRECONDITIONED norm type for convergence test
      PC Object:      (fieldsplit_boundary_)     
        type: bjacobi
          block Jacobi: number of blocks = 2
          Local solve is same for all blocks, in the following KSP and 
PC objects:
        KSP Object:        (fieldsplit_boundary_sub_)       
          type: preonly
          maximum iterations=10000, initial guess is zero
          tolerances:  relative=1e-05, absolute=1e-50, divergence=10000
          left preconditioning
          using PRECONDITIONED norm type for convergence test
        PC Object:        (fieldsplit_boundary_sub_)       
          type: ilu
            ILU: out-of-place factorization
            0 levels of fill
            tolerance for zero pivot 1e-12
            using diagonal shift to prevent zero pivot
            matrix ordering: natural
            factor fill ratio given 1, needed 1
              Factored matrix follows:
                Matrix Object:               
                  type: seqaij
                  rows=0, cols=0
                  package used to perform factorization: petsc
                  total: nonzeros=0, allocated nonzeros=0
                  total number of mallocs used during MatSetValues calls =0
                    not using I-node routines
          linear system matrix = precond matrix:
          Matrix Object:         
            type: seqaij
            rows=0, cols=0
            total: nonzeros=0, allocated nonzeros=0
            total number of mallocs used during MatSetValues calls =0
              not using I-node routines
        linear system matrix followed by preconditioner matrix:
        Matrix Object:       
          type: schurcomplement
          rows=291, cols=291
            Schur complement A11 - A10 inv(A00) A01
            A11
              Matrix Object:             
                type: mpiaij
                rows=291, cols=291
                total: nonzeros=1831, allocated nonzeros=1831
                total number of mallocs used during MatSetValues calls =0
                  not using I-node (on process 0) routines
            A10
              Matrix Object:             
                type: mpiaij
                rows=291, cols=24672
                total: nonzeros=2568, allocated nonzeros=2568
                total number of mallocs used during MatSetValues calls =0
                  not using I-node (on process 0) routines
            KSP of A00
              KSP Object:              (fieldsplit_interior_)             
                type: preonly
                maximum iterations=10000, initial guess is zero
                tolerances:  relative=1e-05, absolute=1e-50, 
divergence=10000
                left preconditioning
                using PRECONDITIONED norm type for convergence test
              PC Object:              (fieldsplit_interior_)             
                type: bjacobi
                  block Jacobi: number of blocks = 2
                  Local solve is same for all blocks, in the following 
KSP and PC objects:
                KSP Object:                
(fieldsplit_interior_sub_)               
                  type: preonly
                  maximum iterations=10000, initial guess is zero
                  tolerances:  relative=1e-05, absolute=1e-50, 
divergence=10000
                  left preconditioning
                  using PRECONDITIONED norm type for convergence test
                PC Object:                
(fieldsplit_interior_sub_)               
                  type: lu
                    LU: out-of-place factorization
                    tolerance for zero pivot 1e-12
                    matrix ordering: nd
                    factor fill ratio given 0, needed 0
                      Factored matrix follows:
                        Matrix Object:                       
                          type: seqaij
                          rows=12336, cols=12336
                          package used to perform factorization: umfpack
                          total: nonzeros=0, allocated nonzeros=0
                          total number of mallocs used during 
MatSetValues calls =0
                            not using I-node routines
                            UMFPACK run parameters:
                              Control[UMFPACK_PRL]: 1
                              Control[UMFPACK_STRATEGY]: 0
                              Control[UMFPACK_DENSE_COL]: 0.2
                              Control[UMFPACK_DENSE_ROW]: 0.2
                              Control[UMFPACK_AMD_DENSE]: 10
                              Control[UMFPACK_BLOCK_SIZE]: 32
                              Control[UMFPACK_2BY2_TOLERANCE]: 0.01
                              Control[UMFPACK_FIXQ]: 0
                              Control[UMFPACK_AGGRESSIVE]: 1
                              Control[UMFPACK_PIVOT_TOLERANCE]: 0.1
                              Control[UMFPACK_SYM_PIVOT_TOLERANCE]: 0.001
                              Control[UMFPACK_SCALE]: 1
                              Control[UMFPACK_ALLOC_INIT]: 0.7
                              Control[UMFPACK_DROPTOL]: 0
                              Control[UMFPACK_IRSTEP]: 0
                              UMFPACK default matrix ordering is used 
(not the PETSc matrix ordering)
                  linear system matrix = precond matrix:
                  Matrix Object:                 
                    type: seqaij
                    rows=12336, cols=12336
                    total: nonzeros=173608, allocated nonzeros=173608
                    total number of mallocs used during MatSetValues 
calls =0
                      not using I-node routines
                linear system matrix = precond matrix:
                Matrix Object:               
                  type: mpiaij
                  rows=24672, cols=24672
                  total: nonzeros=347216, allocated nonzeros=347216
                  total number of mallocs used during MatSetValues calls =0
                    not using I-node (on process 0) routines
            A01
              Matrix Object:             
                type: mpiaij
                rows=24672, cols=291
                total: nonzeros=2568, allocated nonzeros=2568
                total number of mallocs used during MatSetValues calls =0
                  using I-node (on process 0) routines: found 2468 
nodes, limit used is 5
        Matrix Object:       
          type: mpiaij
          rows=291, cols=291
          total: nonzeros=1831, allocated nonzeros=1831
          total number of mallocs used during MatSetValues calls =0
            not using I-node (on process 0) routines
  linear system matrix = precond matrix:
  Matrix Object: 
    type: mpiaij
    rows=24963, cols=24963
    total: nonzeros=354183, allocated nonzeros=356014
    total number of mallocs used during MatSetValues calls =0
      not using I-node (on process 0) routines


Thomas
...
We are eager to determine the problem.
Barry
On Mar 16, 2011, at 11:03 AM, Barry Smith wrote:
...
On Mar 16, 2011, at 8:57 AM, Thomas Witkowski wrote:
-pc_fieldsplit_type schur -fieldsplit_interior_ksp_type preonly -fieldsplit_interior_pc_type bjacobi -fieldsplit_interior_sub_pc_type lu -fieldsplit_boundary_ksp_monitor -ksp_monitor_true_residual
...
Residual norms for fieldsplit_boundary_ solve.
 0 KSP Residual norm 1.790059331071e-04
.....
18 KSP Residual norm 1.593454221996e-09
0 KSP preconditioned resid norm 5.745043818120e-01 true resid norm 1.584249437360e-01 ||r(i)||/||b|| 1.000000000000e+00
This is the norms on the outer iteration variables
...
Residual norms for fieldsplit_boundary_ solve.
 0 KSP Residual norm 1.790059331000e-04
This is the norms on the Schur complement iteration
...
....
18 KSP Residual norm 1.593449602302e-09
It is unclear to me why it starts up another solve here, starting basically where it left off. I think it should just kick out to print the norm of the residual of the outer iteration.
...
Residual norms for fieldsplit_boundary_ solve.
 0 KSP Residual norm 1.593449603155e-09
...
18 KSP Residual norm 1.335744434252e-14
1 KSP preconditioned resid norm 2.608069232839e-13 true resid norm 1.449869994655e-13 ||r(i)||/||b|| 9.151778504469e-13
Are you running with petsc-dev?
Please run with -ksp_view and send the output.
When I run tutorials/multiphysics/mp.c with the options ./mp -couple -snes_monitor_short -pc_type fieldsplit -ksp_monitor -pc_fieldsplit_type schur -fieldsplit_ksp_monitor -fieldsplit_1_ksp_type fgmres -fieldsplit_0_ksp_type gmres   -snes_view I do not see the "double" iteration business.
Barry
...
So, are there now two or one outer iterations?
Thomas
Matthew Knepley wrote:
...
On Wed, Mar 16, 2011 at 8:27 AM, Thomas Witkowski <[email protected] <mailto:[email protected]>> wrote:
Jed Brown wrote:
On Wed, Mar 16, 2011 at 07:37, Thomas Witkowski
      <[email protected]
      <mailto:[email protected]>
      <mailto:[email protected]
      <mailto:[email protected]>>> wrote:
Thanks for explanations! It works fine in my code. But I
      have two
         questions about it, maybe you can help me with them:
         - To the first, is the LU factorization on block A_00 done
      only once?
Yes, one way to confirm this is to run with -log_summary and
      check MatLUFactor{Sym,Num}.
                  - I have run the code with -fieldsplit_1_ksp_monitor to get
      some
         more information about the internal solves. I expected to get
         information about one iterative solver (for solving the Schur
         complement system), but I got three, all of them need around 20
         iterations for my example. Could you explain to me what is
         actually solved there?
I'm going to guess that you still have an outer KSP that (in
      the global norm, rather than the partitioned norm used inside
      of splits) has a tighter tolerance, therefore it takes a few
      outer iterations. If you use loose inner tolerances then the
      preconditioner becomes nonlinear and you'll need to use FGMRES
      for the outer. When in doubt, run with -ksp_view and show us
      the results if you don't understand. Additionally, monitoring
      inner solves separately can be useful, e.g.
      -fieldsplit_0_ksp_converged_reason -fieldsplit_1_ksp_monitor
      -ksp_monitor_true_residual.
Okay, I run my code with the options "-pc_fieldsplit_type schur
  -fieldsplit_interior_ksp_type preonly -fieldsplit_interior_pc_type
  bjacobi -fieldsplit_interior_sub_pc_type lu
  -fieldsplit_boundary_ksp_monitor -ksp_monitor_true_residual
  -fieldsplit_interior_ksp_converged_reason". The splits are named
  "interior" and "boundary". The ksp output is as follows:
Linear solve converged due to CONVERGED_ITS iterations 1
    Residual norms for fieldsplit_boundary_ solve.
    0 KSP Residual norm 1.790059331071e-04
    Linear solve converged due to CONVERGED_ITS iterations 1
    1 KSP Residual norm 1.237356212928e-04
    Linear solve converged due to CONVERGED_ITS iterations 1
    2 KSP Residual norm 7.952220245101e-05
    Linear solve converged due to CONVERGED_ITS iterations 1
    3 KSP Residual norm 3.502858665398e-05
    Linear solve converged due to CONVERGED_ITS iterations 1
    4 KSP Residual norm 1.601085190129e-05
    Linear solve converged due to CONVERGED_ITS iterations 1
    5 KSP Residual norm 8.491161243057e-06
    Linear solve converged due to CONVERGED_ITS iterations 1
    6 KSP Residual norm 4.778227395989e-06
    Linear solve converged due to CONVERGED_ITS iterations 1
    7 KSP Residual norm 2.571919352633e-06
    Linear solve converged due to CONVERGED_ITS iterations 1
    8 KSP Residual norm 1.353979940120e-06
    Linear solve converged due to CONVERGED_ITS iterations 1
    9 KSP Residual norm 7.408377337496e-07
    Linear solve converged due to CONVERGED_ITS iterations 1
   10 KSP Residual norm 3.810694345834e-07
    Linear solve converged due to CONVERGED_ITS iterations 1
   11 KSP Residual norm 1.955889542056e-07
    Linear solve converged due to CONVERGED_ITS iterations 1
   12 KSP Residual norm 9.938906113920e-08
    Linear solve converged due to CONVERGED_ITS iterations 1
   13 KSP Residual norm 4.893445675773e-08
    Linear solve converged due to CONVERGED_ITS iterations 1
   14 KSP Residual norm 2.395128458062e-08
    Linear solve converged due to CONVERGED_ITS iterations 1
   15 KSP Residual norm 1.2464425896 <tel:1.2464425896>57e-08
    Linear solve converged due to CONVERGED_ITS iterations 1
   16 KSP Residual norm 6.316454504022e-09
    Linear solve converged due to CONVERGED_ITS iterations 1
   17 KSP Residual norm 3.088668836241e-09
    Linear solve converged due to CONVERGED_ITS iterations 1
   18 KSP Residual norm 1.593454221996e-09
    Linear solve converged due to CONVERGED_ITS iterations 1
  [0]  Petsc-Iteration 0: 0.574504
   0 KSP preconditioned resid norm 5.745043818120e-01 true resid
  norm 1.584249437360e-01 ||r(i)||/||b|| 1.000000000000e+00
    Linear solve converged due to CONVERGED_ITS iterations 1
    Residual norms for fieldsplit_boundary_ solve.
    0 KSP Residual norm 1.790059331000e-04
    Linear solve converged due to CONVERGED_ITS iterations 1
    1 KSP Residual norm 1.237356237744e-04
    Linear solve converged due to CONVERGED_ITS iterations 1
    2 KSP Residual norm 7.952220078415e-05
    Linear solve converged due to CONVERGED_ITS iterations 1
    3 KSP Residual norm 3.502859285561e-05
    Linear solve converged due to CONVERGED_ITS iterations 1
    4 KSP Residual norm 1.601086485810e-05
    Linear solve converged due to CONVERGED_ITS iterations 1
    5 KSP Residual norm 8.491169384185e-06
    Linear solve converged due to CONVERGED_ITS iterations 1
    6 KSP Residual norm 4.778229157872e-06
    Linear solve converged due to CONVERGED_ITS iterations 1
    7 KSP Residual norm 2.571917112191e-06
    Linear solve converged due to CONVERGED_ITS iterations 1
    8 KSP Residual norm 1.353977280183e-06
    Linear solve converged due to CONVERGED_ITS iterations 1
    9 KSP Residual norm 7.408352317360e-07
    Linear solve converged due to CONVERGED_ITS iterations 1
   10 KSP Residual norm 3.810710583336e-07
    Linear solve converged due to CONVERGED_ITS iterations 1
   11 KSP Residual norm 1.955898169394e-07
    Linear solve converged due to CONVERGED_ITS iterations 1
   12 KSP Residual norm 9.938892788465e-08
    Linear solve converged due to CONVERGED_ITS iterations 1
   13 KSP Residual norm 4.893428546840e-08
    Linear solve converged due to CONVERGED_ITS iterations 1
   14 KSP Residual norm 2.395120857135e-08
    Linear solve converged due to CONVERGED_ITS iterations 1
   15 KSP Residual norm 1.2464367927 <tel:1.2464367927>84e-08
    Linear solve converged due to CONVERGED_ITS iterations 1
   16 KSP Residual norm 6.316409779118e-09
    Linear solve converged due to CONVERGED_ITS iterations 1
   17 KSP Residual norm 3.088649355374e-09
    Linear solve converged due to CONVERGED_ITS iterations 1
   18 KSP Residual norm 1.593449602302e-09
    Linear solve converged due to CONVERGED_ITS iterations 1
    Linear solve converged due to CONVERGED_ITS iterations 1
    Residual norms for fieldsplit_boundary_ solve.
    0 KSP Residual norm 1.593449603155e-09
    Linear solve converged due to CONVERGED_ITS iterations 1
    1 KSP Residual norm 1.011857676662e-09
    Linear solve converged due to CONVERGED_ITS iterations 1
    2 KSP Residual norm 5.660892669788e-10
    Linear solve converged due to CONVERGED_ITS iterations 1
    3 KSP Residual norm 3.046799343903e-10
    Linear solve converged due to CONVERGED_ITS iterations 1
    4 KSP Residual norm 1.6004713993 <tel:1.6004713993>29e-10
    Linear solve converged due to CONVERGED_ITS iterations 1
    5 KSP Residual norm 8.083216601016e-11
    Linear solve converged due to CONVERGED_ITS iterations 1
    6 KSP Residual norm 4.308671648345e-11
    Linear solve converged due to CONVERGED_ITS iterations 1
    7 KSP Residual norm 2.381959817139e-11
    Linear solve converged due to CONVERGED_ITS iterations 1
    8 KSP Residual norm 1.264839790114e-11
    Linear solve converged due to CONVERGED_ITS iterations 1
    9 KSP Residual norm 6.435098853282e-12
    Linear solve converged due to CONVERGED_ITS iterations 1
   10 KSP Residual norm 3.184769579797e-12
    Linear solve converged due to CONVERGED_ITS iterations 1
   11 KSP Residual norm 1.5022019430 <tel:1.5022019430>16e-12
    Linear solve converged due to CONVERGED_ITS iterations 1
   12 KSP Residual norm 7.099486208425e-13
    Linear solve converged due to CONVERGED_ITS iterations 1
   13 KSP Residual norm 3.443222176388e-13
    Linear solve converged due to CONVERGED_ITS iterations 1
   14 KSP Residual norm 1.6478434989 <tel:1.6478434989>12e-13
    Linear solve converged due to CONVERGED_ITS iterations 1
   15 KSP Residual norm 8.208508173926e-14
    Linear solve converged due to CONVERGED_ITS iterations 1
   16 KSP Residual norm 4.376604355287e-14
    Linear solve converged due to CONVERGED_ITS iterations 1
   17 KSP Residual norm 2.596069152755e-14
    Linear solve converged due to CONVERGED_ITS iterations 1
   18 KSP Residual norm 1.335744434252e-14
    Linear solve converged due to CONVERGED_ITS iterations 1
   1 KSP preconditioned resid norm 2.608069232839e-13 true resid
  norm 1.449869994655e-13 ||r(i)||/||b|| 9.151778504469e-13
I have no idea how to interpret this output! Could you help me
  with it?
Get rid of -fieldsplit_interior_ksp_converged_reason. It is confusing things. Then you will see
the outer iteration and the inner boundary iteration which is indented.
Matt
Thomas
-- 
What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.
-- Norbert Wiener

Thomas Witkowski

tags

participants (1)