atol=1e-50, rtol=1e-08, stol=1e-08, maxit=50, maxf=10000
  0 SNES Function norm 5.41468
  ---------- Testing Jacobian -------------
  Testing hand-coded Jacobian, if (for double precision runs) ||J - Jfd||_F/||J||_F is
    O(1.e-8), the hand-coded Jacobian is probably correct.
  ||J - Jfd||_F/||J||_F = 0.481271, ||J - Jfd||_F = 32.
  Hand-coded Jacobian ----------
Mat Object: 1 MPI process
  type: seqaij
row 0: (0, 1.)  (1, 16.) 
row 1: (1, -31.)  (2, 16.) 
row 2: (1, 16.)  (2, -31.)  (3, 16.) 
row 3: (2, 16.)  (3, -31.) 
row 4: (3, 16.)  (4, 1.) 
  Finite difference Jacobian ----------
Mat Object: 1 MPI process
  type: seqaij
row 0: (0, 1.) 
row 1: (0, 16.)  (1, -31.)  (2, 16.) 
row 2: (1, 16.)  (2, -31.)  (3, 16.) 
row 3: (2, 16.)  (3, -31.)  (4, 16.) 
row 4: (4, 1.) 
  Hand-coded minus finite-difference Jacobian with tolerance 1e-05 ----------
Mat Object: 1 MPI process
  type: seqaij
row 0: (1, 16.) 
row 1: (0, -16.) 
row 2:
row 3: (4, -16.) 
row 4: (3, 16.) 
    0 KSP Residual norm 6.945929944146e+00
    1 KSP Residual norm 2.792520156365e-15
Vec Object: 1 MPI process
  type: seq
0.5
-1.25024
-2.76563
-4.42798
-0.5
Vec Object: 1 MPI process
  type: seq
-3.39308
0.243317
0.393287
0.345825
-6.0332
Mat Object: 1 MPI process
  type: seqaij
row 0: (0, 1.)  (1, 16.) 
row 1: (1, -31.)  (2, 16.) 
row 2: (1, 16.)  (2, -31.)  (3, 16.) 
row 3: (2, 16.)  (3, -31.) 
row 4: (3, 16.)  (4, 1.) 
Mat Object: 1 MPI process
  type: seqaij
row 0: (0, 1.) 
row 1: (1, 1.) 
row 2: (2, 1.) 
row 3: (3, 1.) 
row 4: (4, 1.) 
  1 SNES Function norm 111.09
SNES Object: 1 MPI process
  type: ksptransposeonly
  maximum iterations=50, maximum function evaluations=10000
  tolerances: relative=1e-08, absolute=1e-50, solution=1e-08
  total number of linear solver iterations=1
  total number of function evaluations=4
  norm schedule ALWAYS
  KSP Object: 1 MPI process
    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 process
    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:
          Mat Object: 1 MPI process
            type: seqaij
            rows=5, cols=5, bs=5
            package used to perform factorization: petsc
            total: nonzeros=11, allocated nonzeros=11
              not using I-node routines
    linear system matrix, which is also used to construct the preconditioner:
    Mat Object: 1 MPI process
      type: seqaij
      rows=5, cols=5, bs=5 variable blocks set
      total: nonzeros=11, allocated nonzeros=15
      total number of mallocs used during MatSetValues calls=0
        not using I-node routines
number of SNES iterations = 1

Norm of error 6.77514, Iterations 1