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>
