1*c4762a1bSJed Brownatol=1e-50, rtol=1e-08, stol=1e-08, maxit=50, maxf=10000 2*c4762a1bSJed Browniter = 0,SNES Function norm 5.41468 3*c4762a1bSJed Brown ---------- Testing Jacobian ------------- 4*c4762a1bSJed Brown Testing hand-coded Jacobian, if (for double precision runs) ||J - Jfd||_F/||J||_F is 5*c4762a1bSJed Brown O(1.e-8), the hand-coded Jacobian is probably correct. 6*c4762a1bSJed Brown ||J - Jfd||_F/||J||_F = 1.40605e-09, ||J - Jfd||_F = 9.34893e-08 7*c4762a1bSJed Brown Hand-coded Jacobian ---------- 8*c4762a1bSJed BrownMat Object: 2 MPI processes 9*c4762a1bSJed Brown type: mpiaij 10*c4762a1bSJed Brownrow 0: (0, 1.) 11*c4762a1bSJed Brownrow 1: (0, 16.) (1, -31.) (2, 16.) 12*c4762a1bSJed Brownrow 2: (1, 16.) (2, -31.) (3, 16.) 13*c4762a1bSJed Brownrow 3: (2, 16.) (3, -31.) (4, 16.) 14*c4762a1bSJed Brownrow 4: (4, 1.) 15*c4762a1bSJed Brown Finite difference Jacobian ---------- 16*c4762a1bSJed BrownMat Object: 2 MPI processes 17*c4762a1bSJed Brown type: mpiaij 18*c4762a1bSJed Brownrow 0: (0, 1.) 19*c4762a1bSJed Brownrow 1: (0, 16.) (1, -31.) (2, 16.) 20*c4762a1bSJed Brownrow 2: (1, 16.) (2, -31.) (3, 16.) 21*c4762a1bSJed Brownrow 3: (2, 16.) (3, -31.) (4, 16.) 22*c4762a1bSJed Brownrow 4: (4, 1.) 23*c4762a1bSJed Brown Hand-coded minus finite-difference Jacobian with tolerance 1e-05 ---------- 24*c4762a1bSJed BrownMat Object: 2 MPI processes 25*c4762a1bSJed Brown type: mpiaij 26*c4762a1bSJed Brownrow 0: 27*c4762a1bSJed Brownrow 1: 28*c4762a1bSJed Brownrow 2: 29*c4762a1bSJed Brownrow 3: 30*c4762a1bSJed Brownrow 4: 31*c4762a1bSJed Browniter = 1,SNES Function norm 0.295258 32*c4762a1bSJed Brown ---------- Testing Jacobian ------------- 33*c4762a1bSJed Brown ||J - Jfd||_F/||J||_F = 1.30859e-09, ||J - Jfd||_F = 8.82153e-08 34*c4762a1bSJed Brown Hand-coded Jacobian ---------- 35*c4762a1bSJed BrownMat Object: 2 MPI processes 36*c4762a1bSJed Brown type: mpiaij 37*c4762a1bSJed Brownrow 0: (0, 1.) 38*c4762a1bSJed Brownrow 1: (0, 16.) (1, -32.0028) (2, 16.) 39*c4762a1bSJed Brownrow 2: (1, 16.) (2, -31.7866) (3, 16.) 40*c4762a1bSJed Brownrow 3: (2, 16.) (3, -31.1755) (4, 16.) 41*c4762a1bSJed Brownrow 4: (4, 1.) 42*c4762a1bSJed Brown Finite difference Jacobian ---------- 43*c4762a1bSJed BrownMat Object: 2 MPI processes 44*c4762a1bSJed Brown type: mpiaij 45*c4762a1bSJed Brownrow 0: (0, 1.) 46*c4762a1bSJed Brownrow 1: (0, 16.) (1, -32.0028) (2, 16.) 47*c4762a1bSJed Brownrow 2: (1, 16.) (2, -31.7866) (3, 16.) 48*c4762a1bSJed Brownrow 3: (2, 16.) (3, -31.1755) (4, 16.) 49*c4762a1bSJed Brownrow 4: (4, 1.) 50*c4762a1bSJed Brown Hand-coded minus finite-difference Jacobian with tolerance 1e-05 ---------- 51*c4762a1bSJed BrownMat Object: 2 MPI processes 52*c4762a1bSJed Brown type: mpiaij 53*c4762a1bSJed Brownrow 0: 54*c4762a1bSJed Brownrow 1: 55*c4762a1bSJed Brownrow 2: 56*c4762a1bSJed Brownrow 3: 57*c4762a1bSJed Brownrow 4: 58*c4762a1bSJed Browniter = 2,SNES Function norm 0.000450229 59*c4762a1bSJed Brown ---------- Testing Jacobian ------------- 60*c4762a1bSJed Brown ||J - Jfd||_F/||J||_F = 9.03654e-10, ||J - Jfd||_F = 6.08796e-08 61*c4762a1bSJed Brown Hand-coded Jacobian ---------- 62*c4762a1bSJed BrownMat Object: 2 MPI processes 63*c4762a1bSJed Brown type: mpiaij 64*c4762a1bSJed Brownrow 0: (0, 1.) 65*c4762a1bSJed Brownrow 1: (0, 16.) (1, -31.9688) (2, 16.) 66*c4762a1bSJed Brownrow 2: (1, 16.) (2, -31.7501) (3, 16.) 67*c4762a1bSJed Brownrow 3: (2, 16.) (3, -31.1563) (4, 16.) 68*c4762a1bSJed Brownrow 4: (4, 1.) 69*c4762a1bSJed Brown Finite difference Jacobian ---------- 70*c4762a1bSJed BrownMat Object: 2 MPI processes 71*c4762a1bSJed Brown type: mpiaij 72*c4762a1bSJed Brownrow 0: (0, 1.) 73*c4762a1bSJed Brownrow 1: (0, 16.) (1, -31.9688) (2, 16.) 74*c4762a1bSJed Brownrow 2: (1, 16.) (2, -31.7501) (3, 16.) 75*c4762a1bSJed Brownrow 3: (2, 16.) (3, -31.1563) (4, 16.) 76*c4762a1bSJed Brownrow 4: (4, 1.) 77*c4762a1bSJed Brown Hand-coded minus finite-difference Jacobian with tolerance 1e-05 ---------- 78*c4762a1bSJed BrownMat Object: 2 MPI processes 79*c4762a1bSJed Brown type: mpiaij 80*c4762a1bSJed Brownrow 0: 81*c4762a1bSJed Brownrow 1: 82*c4762a1bSJed Brownrow 2: 83*c4762a1bSJed Brownrow 3: 84*c4762a1bSJed Brownrow 4: 85*c4762a1bSJed Browniter = 3,SNES Function norm 1.38967e-09 86*c4762a1bSJed BrownNumber of SNES iterations = 3 87*c4762a1bSJed BrownNorm of error 1.49751e-10 Iterations 3 88