<div class="gmail_quote">On Wed, Oct 26, 2011 at 23:40, Paul Sicsic <span dir="ltr"><<a href="mailto:paul.sicsic@polytechnique.edu">paul.sicsic@polytechnique.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
I have attached the output of valgrind for a run with no options</blockquote><div><br></div><div>You can do this run in a debugger with the option -fp_trap. It will break at the first place where a denormal or NaN is computed.</div>
<div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;"> and for a run with the options you have recommended earlier. </blockquote></div><br><div>The X you passed in is not a valid Mat (we could see that from the earlier call). More specifically, according to valgrind, it is uninitialized memory. The shift is still infinite coming into this function, which I don't understand.</div>
<div><br></div><div>From the singularity at 0, it looks like perhaps the initial time step size is being set to zero. Next time you get a stack trace, can you on ts->time_step?</div>