1*66edf50cSStefano Zampini---------------------------- 2*66edf50cSStefano ZampiniSimulation parameters: 3*66edf50cSStefano Zampini dim : 3 4*66edf50cSStefano Zampini r : 0.1 5*66edf50cSStefano Zampini eps : 0.001 6*66edf50cSStefano Zampini alpha: 0.75 7*66edf50cSStefano Zampini gamma: 0.75 8*66edf50cSStefano Zampini D : 0.01 9*66edf50cSStefano Zampini IC : 3 10*66edf50cSStefano Zampini snum : 1 11*66edf50cSStefano Zampini x0 : (0.25,0.25,0.25) 12*66edf50cSStefano Zampini scals: (500.,0.,0.,1.) 13*66edf50cSStefano Zampini---------------------------- 14*66edf50cSStefano Zampini ---------- Testing Jacobian ------------- 15*66edf50cSStefano Zampini Run with -snes_test_jacobian_view and optionally -snes_test_jacobian <threshold> to show difference 16*66edf50cSStefano Zampini of hand-coded and finite difference Jacobian entries greater than <threshold>. 17*66edf50cSStefano Zampini Testing hand-coded Jacobian, if (for double precision runs) ||J - Jfd||_F/||J||_F is 18*66edf50cSStefano Zampini O(1.e-8), the hand-coded Jacobian is probably correct. 19*66edf50cSStefano Zampini ||J - Jfd||_F/||J||_F = 5.46812e-15, ||J - Jfd||_F = 1.24968e-14 20*66edf50cSStefano Zampini 0 TS: time 0., energy 0.94281, intp 0., ell 0. 21*66edf50cSStefano Zampini ---------- Testing Jacobian ------------- 22*66edf50cSStefano Zampini ||J - Jfd||_F/||J||_F = 4.58971e-09, ||J - Jfd||_F = 1.97484e-08 23*66edf50cSStefano Zampini step 0 stage 0 nonlinear its 1, linear its 1 24*66edf50cSStefano Zampini 1 TS: time 0.1, energy 0.914191, intp 1.32349e-23, ell 0. 25