On Wed, Jul 25, 2012 at 4:35 PM, Kirk, Benjamin (JSC-EG311) <span dir="ltr"><<a href="mailto:benjamin.kirk-1@nasa.gov" target="_blank">benjamin.kirk-1@nasa.gov</a>></span> wrote:<br><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Hello -<br>
<br>
I've been using PETSc 3.1 for quite a while now and have been hesitant to<br>
upgrade because of some new behavior I found in 3.2.  Let me explain...<br>
<br>
In petsc-3.1, if the KSP encountered a NaN it would return it to the<br>
application code.  We actually liked this feature because it gives us an<br>
opportunity to catch the NaN and attempt recovery, in our case by decreasing<br>
the time step and trying again.<br>
<br>
It seems in petsc-3.2, however, that PETSc itself aborts internally, so we<br>
are unable to recover from the situation.<br>
<br>
Is there any way to get the old behavior back?<br></blockquote><div><br></div><div>1) How exactly could the KSP generate a NaN if it was not injected in A or b?</div><div><br></div><div>2) You can always check the return value of KSPSolve() and do what you did last time.</div>
<div><br></div><div>   Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Thanks,<br>
<br>
-Ben<br>
<br>
</blockquote></div><br><br clear="all"><div><br></div>-- <br>What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>
-- Norbert Wiener<br>