<div class="gmail_quote">On Wed, Oct 26, 2011 at 22:44, 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;">
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div class="gmail_quote"><div>Can you show the value of ark->shift and <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></blockquote></div><br><div>Something must have computed a step size of 0 (or we're doing the wrong thing with a method that has explicit stages; we are changing some things with that). Of the recent methods that we have been playing with, I would recommend starting with "-ts_type rosw -ts_rosw_type ra3pw" or "-ts_type rosw -ts_rosw_type ra34pw2", both of which have nice stability properties, good error estimates, and are more permissive than the ARKs about how you write your equations.</div>