----------------------------
Simulation parameters:
  dim  : 2
  r    : 0.1
  eps  : 0.001
  alpha: 0.75
  gamma: 0.75
  D    : 0.01
  IC   : 3
  snum : 1
  x0   : (0.25,0.25)
  scals: (500.,0.,0.,1.)
----------------------------
  ---------- 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 = 2.27208e-12, ||J - Jfd||_F = 4.24161e-11
   0 TS: time 0., energy 0.745385, intp 3.04932e-20, ell 0.
    ---------- Testing Jacobian -------------
    ||J - Jfd||_F/||J||_F = 2.75773e-08, ||J - Jfd||_F = 5.15249e-07
    ---------- Testing Jacobian -------------
    ||J - Jfd||_F/||J||_F = 2.02598e-08, ||J - Jfd||_F = 1.98015e-07
    ---------- Testing Jacobian -------------
    ||J - Jfd||_F/||J||_F = 1.68913e-08, ||J - Jfd||_F = 8.00889e-08
         step 0 stage 0 nonlinear its 1, linear its 1
   1 TS: time 0.1, energy 0.670449, intp 4.40457e-20, ell 0.