1811d88abSStefano Zampini---------------------------- 2811d88abSStefano ZampiniSimulation parameters: 3811d88abSStefano Zampini r : 0.1 4811d88abSStefano Zampini eps : 0.001 5811d88abSStefano Zampini alpha: 0.75 6811d88abSStefano Zampini gamma: 0.75 7811d88abSStefano Zampini D : 0.01 8811d88abSStefano Zampini c : 5. 9811d88abSStefano Zampini IC : 0 10811d88abSStefano Zampini S : 0 11811d88abSStefano Zampini x0 : (0.25,0.25) 12811d88abSStefano Zampini---------------------------- 13811d88abSStefano Zampini ---------- Testing Jacobian ------------- 14811d88abSStefano Zampini Run with -snes_test_jacobian_view and optionally -snes_test_jacobian <threshold> to show difference 15811d88abSStefano Zampini of hand-coded and finite difference Jacobian entries greater than <threshold>. 16811d88abSStefano Zampini Testing hand-coded Jacobian, if (for double precision runs) ||J - Jfd||_F/||J||_F is 17811d88abSStefano Zampini O(1.e-8), the hand-coded Jacobian is probably correct. 18811d88abSStefano Zampini ||J - Jfd||_F/||J||_F = 5.45199e-14, ||J - Jfd||_F = 9.40197e-13 19*204aa523SStefano Zampini 0 TS: time 0., energy 1.29684, intp -8.47033e-22, ell 0. 20811d88abSStefano Zampini ---------- Testing Jacobian ------------- 21*204aa523SStefano Zampini ||J - Jfd||_F/||J||_F = 3.56169e-07, ||J - Jfd||_F = 6.16169e-06 22*204aa523SStefano Zampini step 0 stage 0 nonlinear its 1, linear its 1 23*204aa523SStefano Zampini 1 TS: time 0.1, energy 1.24867, intp 4.23516e-22, ell 0. 24