---------- 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 = 1.32937e-08, ||J - Jfd||_F = 9.15765e-07
    ---------- Testing Jacobian -------------
    ||J - Jfd||_F/||J||_F = 1.32668e-08, ||J - Jfd||_F = 9.13916e-07
    ---------- Testing Jacobian -------------
    ||J - Jfd||_F/||J||_F = 1.32465e-08, ||J - Jfd||_F = 9.12519e-07
    ---------- Testing Jacobian -------------
    ||J - Jfd||_F/||J||_F = 4.53613e-05, ||J - Jfd||_F = 0.0030338