1*811d88abSStefano Zampini---------------------------- 2*811d88abSStefano ZampiniSimulation parameters: 3*811d88abSStefano Zampini r : 0.1 4*811d88abSStefano Zampini eps : 0.001 5*811d88abSStefano Zampini alpha: 0.75 6*811d88abSStefano Zampini gamma: 0.75 7*811d88abSStefano Zampini D : 0.01 8*811d88abSStefano Zampini c : 5. 9*811d88abSStefano Zampini IC : 0 10*811d88abSStefano Zampini S : 0 11*811d88abSStefano Zampini x0 : (0.25,0.25) 12*811d88abSStefano Zampini---------------------------- 13*811d88abSStefano Zampini ---------- Testing Jacobian ------------- 14*811d88abSStefano Zampini Run with -snes_test_jacobian_view and optionally -snes_test_jacobian <threshold> to show difference 15*811d88abSStefano Zampini of hand-coded and finite difference Jacobian entries greater than <threshold>. 16*811d88abSStefano Zampini Testing hand-coded Jacobian, if (for double precision runs) ||J - Jfd||_F/||J||_F is 17*811d88abSStefano Zampini O(1.e-8), the hand-coded Jacobian is probably correct. 18*811d88abSStefano Zampini ||J - Jfd||_F/||J||_F = 5.45199e-14, ||J - Jfd||_F = 9.40197e-13 19*811d88abSStefano Zampini 0 TS: time 0., energy 1.29684, intp 1.69407e-21, ell 0. 20*811d88abSStefano Zampini nonlinear its 0, linear its 0 21*811d88abSStefano Zampini ---------- Testing Jacobian ------------- 22*811d88abSStefano Zampini ||J - Jfd||_F/||J||_F = 3.56176e-07, ||J - Jfd||_F = 6.16181e-06 23*811d88abSStefano Zampini 1 TS: time 0.1, energy 1.24867, intp -2.11758e-22, ell 0.25 24*811d88abSStefano Zampini nonlinear its 1, linear its 1 25