On Sun, Nov 27, 2011 at 1:21 PM, Mike Sussman <span dir="ltr"><<a href="mailto:mmsussman@gmail.com">mmsussman@gmail.com</a>></span> wrote:<br><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
Hello,<br>
<br>
Recently I ran a very simple problem with a known solution using the DAE<br>
form in TS. When it got the wrong answer, I looked carefully at its<br>
progress and discovered that the nonlinear (snes) solves were failing<br>
because of my incorrect Jacobian. BUT these failures were not reported<br>
by TS, and there was no indication that the solution was suspect.<br>
<br>
I tried setting -ts_max_snes_failures 1, but the calculation still<br>
proceeded without halting. Am I missing some setting?<br>
<br>
I am using 3.2p5 from Fortran.<br>
</blockquote></div><br>1) Are you calling TSStep() or TSSolve()?<div><br></div><div>2) Do you check all return codes from PETSc calls?</div><div><br></div><div> Matt<br clear="all"><div><br></div>-- <br>What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>
-- Norbert Wiener<br>
</div>