<div class="gmail_quote">On Wed, Aug 17, 2011 at 07:02, Richard Katz <span dir="ltr"><<a href="mailto:richard.katz@earth.ox.ac.uk">richard.katz@earth.ox.ac.uk</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
<div id=":94">Is there an example showing the use of these controllers?<br></div></blockquote><div><br></div><div>Simplest is to run with -ts_type alpha -ts_alpha_adapt. Run with -help |grep ts_alpha_adapt for the parameters.</div>
<div><br></div><div>src/ts/examples/tutorials/ex18.c uses this in an isogeometric method for Navier-Stokes-Korteweg</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
<div id=":94">
<div class="im"><br>
> I think having these interfaces is really undesirable and that we should unify it, but the information provided by the error estimates for TSGL are quite different from most methods (OTOH, they seem to me noisy which limits their utility, but there are other possible reasons for that, including "starting methods").<br>

<br>
</div>Is it possible for the user to push a custom controller?  I.e. a call-back that gets passed the SNES converged reason, the last SNES residual, the current time-step size, the solution vector, etc.<br></div></blockquote>
<div><br></div><div>That is how it was designed, except that the controller would probably pull out the SNES themselves and query the reason themselves (because it's hard to predict every piece of information that they might want to use).</div>
<div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;"><div id=":94"><div class="im">
<br>
</div>So now that I have TSBEULER working for my problem, should I try to upgrade to TSALPHA ?</div></blockquote></div><br><div>The source you sent is a DAE, which TSALPHA (not being L-stable) will not be robust for. You should still be able to "try it" by running with -ts_type alpha (or -ts_type theta which is closely related; I just tried alpha and TSInterpolate_Alpha is missing, will fix).</div>