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.48231e-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 = 0.999818, ||J - Jfd||_F = 7. 23*811d88abSStefano Zampini ---------- Testing Jacobian ------------- 24*811d88abSStefano Zampini ||J - Jfd||_F/||J||_F = 3.45224e-07, ||J - Jfd||_F = 6.16181e-06 25*811d88abSStefano Zampini ---------- Testing Jacobian ------------- 26*811d88abSStefano Zampini ||J - Jfd||_F/||J||_F = 3.49161e-07, ||J - Jfd||_F = 6.08119e-06 27*811d88abSStefano Zampini 1 TS: time 0.1, energy 1.24896, intp 2.91168e-21, ell 0.291667 28*811d88abSStefano Zampini nonlinear its 1, linear its 1 29