----------------------------
Simulation parameters:
  r    : 0.1
  eps  : 0.001
  alpha: 0.75
  gamma: 0.75
  D    : 0.01
  c    : 5.
  IC   : 0
  S    : 0
  x0   : (0.25,0.25)
----------------------------
  ---------- Testing Jacobian -------------
  Run with -snes_test_jacobian_view and optionally -snes_test_jacobian <threshold> to show difference
    of hand-coded and finite difference Jacobian entries greater than <threshold>.
  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 = 5.45199e-14, ||J - Jfd||_F = 9.40197e-13
   0 TS: time 0., energy 1.29684, intp -8.47033e-22, ell 0.
      ---------- Testing Jacobian -------------
      ||J - Jfd||_F/||J||_F = 0.999818, ||J - Jfd||_F = 7.
         step 0 stage 0 nonlinear its 0, linear its 0
    ---------- Testing Jacobian -------------
    ||J - Jfd||_F/||J||_F = 3.45221e-07, ||J - Jfd||_F = 6.16175e-06
         step 0 stage 1 nonlinear its 1, linear its 1
    ---------- Testing Jacobian -------------
    ||J - Jfd||_F/||J||_F = 3.49164e-07, ||J - Jfd||_F = 6.08126e-06
         step 0 stage 2 nonlinear its 1, linear its 1
   1 TS: time 0.1, energy 1.24896, intp -1.05879e-21, ell 0.
