1c4762a1bSJed Brownatol=1e-50, rtol=1e-08, stol=1e-08, maxit=50, maxf=10000 2c4762a1bSJed Brown 0 SNES Function norm 5.41468 3c4762a1bSJed Brown ---------- Testing Jacobian ------------- 4c4762a1bSJed Brown Testing hand-coded Jacobian, if (for double precision runs) ||J - Jfd||_F/||J||_F is 5c4762a1bSJed Brown O(1.e-8), the hand-coded Jacobian is probably correct. 6c4762a1bSJed Brown ||J - Jfd||_F/||J||_F = 0.481271, ||J - Jfd||_F = 32. 7c4762a1bSJed Brown Hand-coded Jacobian ---------- 88cc725e6SPierre JolivetMat Object: 1 MPI process 9c4762a1bSJed Brown type: seqaij 10c4762a1bSJed Brownrow 0: (0, 1.) (1, 16.) 11c4762a1bSJed Brownrow 1: (1, -31.) (2, 16.) 12c4762a1bSJed Brownrow 2: (1, 16.) (2, -31.) (3, 16.) 13c4762a1bSJed Brownrow 3: (2, 16.) (3, -31.) 14c4762a1bSJed Brownrow 4: (3, 16.) (4, 1.) 15c4762a1bSJed Brown Finite difference Jacobian ---------- 168cc725e6SPierre JolivetMat Object: 1 MPI process 17c4762a1bSJed Brown type: seqaij 18c4762a1bSJed Brownrow 0: (0, 1.) 19c4762a1bSJed Brownrow 1: (0, 16.) (1, -31.) (2, 16.) 20c4762a1bSJed Brownrow 2: (1, 16.) (2, -31.) (3, 16.) 21c4762a1bSJed Brownrow 3: (2, 16.) (3, -31.) (4, 16.) 22c4762a1bSJed Brownrow 4: (4, 1.) 23c4762a1bSJed Brown Hand-coded minus finite-difference Jacobian with tolerance 1e-05 ---------- 248cc725e6SPierre JolivetMat Object: 1 MPI process 25c4762a1bSJed Brown type: seqaij 26c4762a1bSJed Brownrow 0: (1, 16.) 27c4762a1bSJed Brownrow 1: (0, -16.) 28c4762a1bSJed Brownrow 2: 29c4762a1bSJed Brownrow 3: (4, -16.) 30c4762a1bSJed Brownrow 4: (3, 16.) 31c4762a1bSJed Brown 0 KSP Residual norm 6.945929944146e+00 32c4762a1bSJed Brown 1 KSP Residual norm 2.792520156365e-15 338cc725e6SPierre JolivetVec Object: 1 MPI process 34c4762a1bSJed Brown type: seq 35c4762a1bSJed Brown0.5 36c4762a1bSJed Brown-1.25024 37c4762a1bSJed Brown-2.76563 38c4762a1bSJed Brown-4.42798 39c4762a1bSJed Brown-0.5 408cc725e6SPierre JolivetVec Object: 1 MPI process 41c4762a1bSJed Brown type: seq 42c4762a1bSJed Brown-3.39308 43c4762a1bSJed Brown0.243317 44c4762a1bSJed Brown0.393287 45c4762a1bSJed Brown0.345825 46c4762a1bSJed Brown-6.0332 478cc725e6SPierre JolivetMat Object: 1 MPI process 48c4762a1bSJed Brown type: seqaij 49c4762a1bSJed Brownrow 0: (0, 1.) (1, 16.) 50c4762a1bSJed Brownrow 1: (1, -31.) (2, 16.) 51c4762a1bSJed Brownrow 2: (1, 16.) (2, -31.) (3, 16.) 52c4762a1bSJed Brownrow 3: (2, 16.) (3, -31.) 53c4762a1bSJed Brownrow 4: (3, 16.) (4, 1.) 548cc725e6SPierre JolivetMat Object: 1 MPI process 55c4762a1bSJed Brown type: seqaij 56c4762a1bSJed Brownrow 0: (0, 1.) 57c4762a1bSJed Brownrow 1: (1, 1.) 58c4762a1bSJed Brownrow 2: (2, 1.) 59c4762a1bSJed Brownrow 3: (3, 1.) 60c4762a1bSJed Brownrow 4: (4, 1.) 61c4762a1bSJed Brown 1 SNES Function norm 111.09 628cc725e6SPierre JolivetSNES Object: 1 MPI process 63c4762a1bSJed Brown type: ksptransposeonly 64c4762a1bSJed Brown maximum iterations=50, maximum function evaluations=10000 65c4762a1bSJed Brown tolerances: relative=1e-08, absolute=1e-50, solution=1e-08 66c4762a1bSJed Brown total number of linear solver iterations=1 670df40c35SBarry Smith total number of function evaluations=4 68c4762a1bSJed Brown norm schedule ALWAYS 698cc725e6SPierre Jolivet KSP Object: 1 MPI process 70c4762a1bSJed Brown type: gmres 71f971d498SPierre Jolivet restart=30, using classical (unmodified) Gram-Schmidt orthogonalization with no iterative refinement 72*143f2514SPierre Jolivet happy breakdown tolerance=1e-30 73c4762a1bSJed Brown maximum iterations=10000, initial guess is zero 74c4762a1bSJed Brown tolerances: relative=1e-05, absolute=1e-50, divergence=10000. 75c4762a1bSJed Brown left preconditioning 76c4762a1bSJed Brown using PRECONDITIONED norm type for convergence test 778cc725e6SPierre Jolivet PC Object: 1 MPI process 78c4762a1bSJed Brown type: ilu 79c4762a1bSJed Brown out-of-place factorization 80c4762a1bSJed Brown 0 levels of fill 81c4762a1bSJed Brown tolerance for zero pivot 2.22045e-14 82c4762a1bSJed Brown matrix ordering: natural 83c4762a1bSJed Brown factor fill ratio given 1., needed 1. 84ecf3d421SBarry Smith Factored matrix: 858cc725e6SPierre Jolivet Mat Object: 1 MPI process 86c4762a1bSJed Brown type: seqaij 87d7ff4d9fSStefano Zampini rows=5, cols=5, bs=5 88c4762a1bSJed Brown package used to perform factorization: petsc 89c4762a1bSJed Brown total: nonzeros=11, allocated nonzeros=11 90c4762a1bSJed Brown not using I-node routines 91ecf3d421SBarry Smith linear system matrix, which is also used to construct the preconditioner: 928cc725e6SPierre Jolivet Mat Object: 1 MPI process 93c4762a1bSJed Brown type: seqaij 94d7ff4d9fSStefano Zampini rows=5, cols=5, bs=5 variable blocks set 95c4762a1bSJed Brown total: nonzeros=11, allocated nonzeros=15 96c4762a1bSJed Brown total number of mallocs used during MatSetValues calls=0 97c4762a1bSJed Brown not using I-node routines 98c4762a1bSJed Brownnumber of SNES iterations = 1 99c4762a1bSJed Brown 100c4762a1bSJed BrownNorm of error 6.77514, Iterations 1 101