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 VecView(ark->Z,0), VecView(X,0)? Is the system singular? -------------- next part -------------- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20111026/d3afc6ec/attachment.html>
