[petsc-users] TS snes failure

[Matthew Knepley]

On Sun, Nov 27, 2011 at 1:21 PM, Mike Sussman <mmsussman at gmail.com> wrote:

> Hello,
>
> Recently I ran a very simple problem with a known solution using the DAE
> form in TS.  When it got the wrong answer, I looked carefully at its
> progress and discovered that the nonlinear (snes) solves were failing
> because of my incorrect Jacobian.  BUT these failures were not reported
> by TS, and there was no indication that the solution was suspect.
>
> I tried setting -ts_max_snes_failures 1, but the calculation still
> proceeded without halting.  Am I missing some setting?
>
> I am using 3.2p5 from Fortran.
>

1) Are you calling TSStep() or TSSolve()?

2) Do you check all  return codes from PETSc calls?

  Matt

-- 
What most experimenters take for granted before they begin their
experiments is infinitely more interesting than any results to which their
experiments lead.
-- Norbert Wiener
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20111127/ad7cd2f7/attachment.htm>