1---------------------------- 2Simulation parameters: 3 r : 0.1 4 eps : 0.001 5 alpha: 0.75 6 gamma: 0.75 7 D : 0.01 8 c : 5. 9 IC : 0 10 S : 0 11 x0 : (0.25,0.25) 12---------------------------- 13 ---------- Testing Jacobian ------------- 14 Run with -snes_test_jacobian_view and optionally -snes_test_jacobian <threshold> to show difference 15 of hand-coded and finite difference Jacobian entries greater than <threshold>. 16 Testing hand-coded Jacobian, if (for double precision runs) ||J - Jfd||_F/||J||_F is 17 O(1.e-8), the hand-coded Jacobian is probably correct. 18 ||J - Jfd||_F/||J||_F = 5.45199e-14, ||J - Jfd||_F = 9.40197e-13 19 0 TS: time 0., energy 1.29684, intp -8.47033e-22, ell 0. 20 ---------- Testing Jacobian ------------- 21 ||J - Jfd||_F/||J||_F = 0.999818, ||J - Jfd||_F = 7. 22 ---------- Testing Jacobian ------------- 23 ||J - Jfd||_F/||J||_F = 0.998116, ||J - Jfd||_F = 4. 24 step 0 stage 0 nonlinear its 0, linear its 0 25 ---------- Testing Jacobian ------------- 26 ||J - Jfd||_F/||J||_F = 3.4522e-07, ||J - Jfd||_F = 6.16175e-06 27 ---------- Testing Jacobian ------------- 28 ||J - Jfd||_F/||J||_F = 1.94234e-07, ||J - Jfd||_F = 2.25123e-06 29 step 0 stage 1 nonlinear its 1, linear its 1 30 ---------- Testing Jacobian ------------- 31 ||J - Jfd||_F/||J||_F = 3.49164e-07, ||J - Jfd||_F = 6.08125e-06 32 ---------- Testing Jacobian ------------- 33 ||J - Jfd||_F/||J||_F = 1.94125e-07, ||J - Jfd||_F = 2.22328e-06 34 step 0 stage 2 nonlinear its 1, linear its 1 35 1 TS: time 0.1, energy 1.24896, intp -5.29396e-22, ell 0. 36