1---------------------------- 2Simulation parameters: 3 dim : 2 4 r : 0.1 5 eps : 0.001 6 alpha: 0.75 7 gamma: 0.75 8 D : 0.01 9 IC : 3 10 snum : 1 11 x0 : (0.25,0.25) 12 scals: (500.,0.,0.,1.) 13---------------------------- 14 ---------- Testing Jacobian ------------- 15 Run with -snes_test_jacobian_view and optionally -snes_test_jacobian <threshold> to show difference 16 of hand-coded and finite difference Jacobian entries greater than <threshold>. 17 Testing hand-coded Jacobian, if (for double precision runs) ||J - Jfd||_F/||J||_F is 18 O(1.e-8), the hand-coded Jacobian is probably correct. 19 ||J - Jfd||_F/||J||_F = 7.63695e-14, ||J - Jfd||_F = 7.42506e-13 20 0 TS: time 0., energy 0.796955, intp -3.07049e-21, ell 0. 21 ---------- Testing Jacobian ------------- 22 ||J - Jfd||_F/||J||_F = 0.999818, ||J - Jfd||_F = 7. 23 step 0 stage 0 nonlinear its 0, linear its 0 24 ---------- Testing Jacobian ------------- 25 ||J - Jfd||_F/||J||_F = 1.79e-08, ||J - Jfd||_F = 1.92973e-07 26 step 0 stage 1 nonlinear its 1, linear its 1 27 ---------- Testing Jacobian ------------- 28 ||J - Jfd||_F/||J||_F = 1.82482e-08, ||J - Jfd||_F = 1.88397e-07 29 step 0 stage 2 nonlinear its 1, linear its 1 30 1 TS: time 0.1, energy 0.728119, intp -9.52912e-22, ell 0. 31