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