I take over from Blaise. 2011/10/26 Jed Brown <jedbrown at mcs.anl.gov>
> On Wed, Oct 26, 2011 at 22:21, Blaise Bourdin <bourdin at lsu.edu> wrote: > >> At this point we have overcome the original issue but are still having >> difficulties. We followed the lines of ex22 to setup our TS, with the most >> notable difference being that our DM is a DMMesh and not a DMDA. For some >> reason, the second time the Ifunction is called, the Udot argument is a Vec >> of NaN. We checked the TS in the debugger and dan't catch anything obviously >> wrong. We followed the flow of the code in the debugger, and the NaN get >> generated here: >> >> (gdb) >> 713 ierr = >> VecAXPBYPCZ(ark->Ydot,-ark->shift,ark->shift,0,ark->Z,X);CHKERRQ(ierr); /* >> Ydot = shift*(X-Z) */ >> > > Can you show the value of ark->shift and arkimex.c:713, > It is inf (gdb) p ark->shift $3 = inf VecView(X,0)? Is the system singular? > VecView(X,0)? gives 1 and 0 VecAXPBYPCZ (z=0x1582630, alpha=-inf, beta=inf, gamma=0, x=0x15892b0, y=0x155c230) -------------- next part -------------- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20111026/03d64fb7/attachment.html>
