---------------------------- 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 to show difference of hand-coded and finite difference Jacobian entries greater than . 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 = 7.63695e-14, ||J - Jfd||_F = 7.42506e-13 0 TS: time 0., energy 0.796955, intp -3.07049e-21, 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 = 1.78926e-08, ||J - Jfd||_F = 1.92892e-07 step 0 stage 1 nonlinear its 1, linear its 1 ---------- Testing Jacobian ------------- ||J - Jfd||_F/||J||_F = 1.82468e-08, ||J - Jfd||_F = 1.89051e-07 step 0 stage 2 nonlinear its 1, linear its 1 1 TS: time 0.1, energy 0.733515, intp 4.65868e-21, ell 0.