[petsc-dev] Problem mixing sieve and TS

[Paul Sicsic]

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>