1811d88abSStefano Zampini---------------------------- 2811d88abSStefano ZampiniSimulation parameters: 3*66edf50cSStefano Zampini dim : 2 4811d88abSStefano Zampini r : 0.1 5811d88abSStefano Zampini eps : 0.001 6811d88abSStefano Zampini alpha: 0.75 7811d88abSStefano Zampini gamma: 0.75 8811d88abSStefano Zampini D : 0.01 9*66edf50cSStefano Zampini IC : 3 10*66edf50cSStefano Zampini snum : 1 11811d88abSStefano Zampini x0 : (0.25,0.25) 12*66edf50cSStefano Zampini scals: (500.,0.,0.,1.) 13811d88abSStefano Zampini---------------------------- 14811d88abSStefano Zampini ---------- Testing Jacobian ------------- 15811d88abSStefano Zampini Run with -snes_test_jacobian_view and optionally -snes_test_jacobian <threshold> to show difference 16811d88abSStefano Zampini of hand-coded and finite difference Jacobian entries greater than <threshold>. 17811d88abSStefano Zampini Testing hand-coded Jacobian, if (for double precision runs) ||J - Jfd||_F/||J||_F is 18811d88abSStefano Zampini O(1.e-8), the hand-coded Jacobian is probably correct. 19*66edf50cSStefano Zampini ||J - Jfd||_F/||J||_F = 7.63695e-14, ||J - Jfd||_F = 7.42506e-13 20*66edf50cSStefano Zampini 0 TS: time 0., energy 0.796955, intp -3.07049e-21, ell 0. 21811d88abSStefano Zampini ---------- Testing Jacobian ------------- 22811d88abSStefano Zampini ||J - Jfd||_F/||J||_F = 0.999818, ||J - Jfd||_F = 7. 23811d88abSStefano Zampini ---------- Testing Jacobian ------------- 24811d88abSStefano Zampini ||J - Jfd||_F/||J||_F = 0.998116, ||J - Jfd||_F = 4. 25204aa523SStefano Zampini step 0 stage 0 nonlinear its 0, linear its 0 26811d88abSStefano Zampini ---------- Testing Jacobian ------------- 27*66edf50cSStefano Zampini ||J - Jfd||_F/||J||_F = 1.79202e-08, ||J - Jfd||_F = 1.9319e-07 28811d88abSStefano Zampini ---------- Testing Jacobian ------------- 29*66edf50cSStefano Zampini ||J - Jfd||_F/||J||_F = 7.71182e-09, ||J - Jfd||_F = 7.37102e-08 30204aa523SStefano Zampini step 0 stage 1 nonlinear its 1, linear its 1 31811d88abSStefano Zampini ---------- Testing Jacobian ------------- 32*66edf50cSStefano Zampini ||J - Jfd||_F/||J||_F = 1.82401e-08, ||J - Jfd||_F = 1.88314e-07 33811d88abSStefano Zampini ---------- Testing Jacobian ------------- 34*66edf50cSStefano Zampini ||J - Jfd||_F/||J||_F = 7.51701e-09, ||J - Jfd||_F = 7.11916e-08 35204aa523SStefano Zampini step 0 stage 2 nonlinear its 1, linear its 1 36*66edf50cSStefano Zampini 1 TS: time 0.1, energy 0.728179, intp 2.85874e-21, ell 0. 37