<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Wed, May 29, 2013 at 1:37 PM, Jed Brown <span dir="ltr"><<a href="mailto:jedbrown@mcs.anl.gov" target="_blank">jedbrown@mcs.anl.gov</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div class="im">Matthew Emmett <<a href="mailto:memmett@gmail.com">memmett@gmail.com</a>> writes:<br>> Ideally I think it would be nice if the user could just add "-ts_type sdc"<br>

> and have the code figure out which SDC variant (explicit, implicit, or<br>
> imex) to use.  The trouble I am having is that a users RHSFunction can be<br>
> recast into implicit form by TSComputeIFunction if an implicit time-stepper<br>
> is used... which means that I can't just check whether RHSFunction has been<br>
> registered to switch between explicit only and implicit only.<br>
<br>
</div>TSGetRHSFunction() should return NULL if the user did not set an<br>
RHSFunction, however, many stiff systems are described using RHSFunction<br>
as a matter of convenience, and because in some parameter regimes, they<br>
may be non-stiff.<br></blockquote><div><br></div><div style>OK, understood.</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div class="im">
> Is there a preferred way of doing this, or should I add another option<br>
> to select which SDC variant to use?<br>
<br>
</div>I would use an option for this.  The user likely wants to try both for<br>
moderately stiff systems.<br>
<br>
If you would like comments on your code and general approach, please<br>
push it in a fork on bitbucket or github.<br></blockquote><div><br></div><div style>Will do -- I will tidy up a few things and make the history a bit more sane first.</div><div style><br></div><div style>Thanks for your comments,</div>
<div style>Matt </div></div><br></div></div>