yjwu@yjwu-XPS-8910:~/petsc-3.10.1/src/snes/examples/tutorials$ mpiexec -n 1 ./ex216 -f wu-readtwogroups -snes_mf_operator -snes_view -snes_converged_reason -snes_monitor -ksp_converged_reason -ksp_monitor_true_residual
  0 SNES Function norm 8.235090086536e-02 
iter = 0, SNES Function norm 0.0823509
iter = 0, Keff ======= 1.
  Linear solve did not converge due to DIVERGED_PCSETUP_FAILED iterations 0
                 PCSETUP_FAILED due to FACTOR_NUMERIC_ZEROPIVOT 
Nonlinear solve did not converge due to DIVERGED_LINEAR_SOLVE iterations 0
SNES Object: 1 MPI processes
  type: newtonls
  maximum iterations=50, maximum function evaluations=10000
  tolerances: relative=1e-08, absolute=1e-50, solution=1e-08
  total number of linear solver iterations=0
  total number of function evaluations=1
  norm schedule ALWAYS
  SNESLineSearch Object: 1 MPI processes
    type: bt
      interpolation: cubic
      alpha=1.000000e-04
    maxstep=1.000000e+08, minlambda=1.000000e-12
    tolerances: relative=1.000000e-08, absolute=1.000000e-15, lambda=1.000000e-08
    maximum iterations=40
  KSP Object: 1 MPI processes
    type: gmres
      restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement
      happy breakdown tolerance 1e-30
    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: 1 MPI processes
    type: 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=961, cols=961
            package used to perform factorization: petsc
            total: nonzeros=4625, allocated nonzeros=4625
            total number of mallocs used during MatSetValues calls =0
              not using I-node routines
    linear system matrix followed by preconditioner matrix:
    Mat Object: 1 MPI processes
      type: mffd
      rows=961, cols=961
        Matrix-free approximation:
          err=1.49012e-08 (relative error in function evaluation)
          The compute h routine has not yet been set
    Mat Object: 1 MPI processes
      type: seqaij
      rows=961, cols=961
      total: nonzeros=4625, allocated nonzeros=4625
      total number of mallocs used during MatSetValues calls =0
        not using I-node routines