<div>I take over from Blaise. </div><br><div class="gmail_quote">2011/10/26 Jed Brown <span dir="ltr"><<a href="mailto:jedbrown@mcs.anl.gov" target="_blank">jedbrown@mcs.anl.gov</a>></span><br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div class="gmail_quote"><div>On Wed, Oct 26, 2011 at 22:21, Blaise Bourdin <span dir="ltr"><<a href="mailto:bourdin@lsu.edu" target="_blank">bourdin@lsu.edu</a>></span> wrote:<br>


<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div>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:</div>



<div><br></div><div style="font-family:arial, sans-serif;font-size:13px;background-color:rgb(255, 255, 255)">(gdb) </div><div style="font-family:arial, sans-serif;font-size:13px;background-color:rgb(255, 255, 255)">713       ierr = VecAXPBYPCZ(ark->Ydot,-ark->shift,ark->shift,0,ark->Z,X);CHKERRQ(ierr); /* Ydot = shift*(X-Z) */</div>



</blockquote><div><br></div></div><div>Can you show the value of ark->shift and <span class="Apple-style-span" style="font-family: arial, sans-serif; font-size: 14px; background-color: rgb(255, 255, 255); ">arkimex.c:713</span>, </div>

</div></blockquote><div>It is inf  </div><div><div>(gdb) p ark->shift</div><div>$3 = inf</div></div><div><br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div class="gmail_quote">
<div>VecView(X,0)? Is the system singular?</div></div>
</blockquote></div><br>VecView(X,0)? gives 1 and 0 <div><br></div><div>VecAXPBYPCZ (z=0x1582630, alpha=-inf, beta=inf, gamma=0, x=0x15892b0, y=0x155c230)<br clear="all"><div><br></div><br>
</div>