1204aa523SStefano Zampini---------------------------- 2204aa523SStefano ZampiniSimulation parameters: 3*66edf50cSStefano Zampini dim : 2 4204aa523SStefano Zampini r : 0.1 5204aa523SStefano Zampini eps : 0.001 6204aa523SStefano Zampini alpha: 0.75 7204aa523SStefano Zampini gamma: 0.75 8204aa523SStefano Zampini D : 0.01 9*66edf50cSStefano Zampini IC : 3 10*66edf50cSStefano Zampini snum : 1 11204aa523SStefano Zampini x0 : (0.25,0.25) 12*66edf50cSStefano Zampini scals: (500.,0.,0.,1.) 13204aa523SStefano Zampini---------------------------- 14204aa523SStefano Zampini ---------- Testing Jacobian ------------- 15204aa523SStefano Zampini Run with -snes_test_jacobian_view and optionally -snes_test_jacobian <threshold> to show difference 16204aa523SStefano Zampini of hand-coded and finite difference Jacobian entries greater than <threshold>. 17204aa523SStefano Zampini Testing hand-coded Jacobian, if (for double precision runs) ||J - Jfd||_F/||J||_F is 18204aa523SStefano 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. 21204aa523SStefano Zampini ---------- Testing Jacobian ------------- 22204aa523SStefano Zampini ||J - Jfd||_F/||J||_F = 0.999818, ||J - Jfd||_F = 7. 23204aa523SStefano Zampini step 0 stage 0 nonlinear its 0, linear its 0 24204aa523SStefano Zampini ---------- Testing Jacobian ------------- 25*66edf50cSStefano Zampini ||J - Jfd||_F/||J||_F = 1.78926e-08, ||J - Jfd||_F = 1.92892e-07 26204aa523SStefano Zampini step 0 stage 1 nonlinear its 1, linear its 1 27204aa523SStefano Zampini ---------- Testing Jacobian ------------- 28*66edf50cSStefano Zampini ||J - Jfd||_F/||J||_F = 1.82468e-08, ||J - Jfd||_F = 1.89051e-07 29204aa523SStefano Zampini step 0 stage 2 nonlinear its 1, linear its 1 30*66edf50cSStefano Zampini 1 TS: time 0.1, energy 0.733515, intp 4.65868e-21, ell 0. 31