<div class="gmail_quote">On Fri, Mar 18, 2011 at 12:07, Gianluca Meneghello <span dir="ltr"><<a href="mailto:gianmail@gmail.com">gianmail@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
<div id=":32f">esting hand-coded Jacobian, if the ratio is<br>
O(1.e-8), the hand-coded Jacobian is probably correct.<br>
Run with -snes_test_display to show difference<br>
of hand-coded and finite difference Jacobian.<br>
Norm of matrix ratio 1.08967e-07 difference 2.43658e-07<br>
Norm of matrix ratio 7.6194e-08 difference 3.4075e-07<br>
Norm of matrix ratio 7.6211e-08 difference 3.40826e-07<br></div></blockquote><div><br></div><div>Your Jacobian looks good so you don't need to use TEST.</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
<div id=":32f">
[0]PETSC ERROR: SNESSolve() line 2360 in src/snes/interface/snes.c<br>
[0]PETSC ERROR: main() line 144 in src/snes/examples/tutorials/ex1.c<br>
application called MPI_Abort(MPI_COMM_WORLD, 73) - process 0<br>
[unset]: aborting job:<br>
application called MPI_Abort(MPI_COMM_WORLD, 73) - process 0</div></blockquote></div><br><div>This is intentional, the end of SNESSolve_Test() has</div><div><br></div><div><div> /*</div><div> Return error code cause Jacobian not good</div>
<div> */</div><div> PetscFunctionReturn(PETSC_ERR_ARG_WRONGSTATE);</div></div><div><br></div><div>Perhaps it should be done differently, maybe even solve the Newton step and continue, but that's not the current behavior.</div>