<div class="gmail_quote">On Fri, Dec 16, 2011 at 13:04, Konstantinos Kontzialis <span dir="ltr"><<a href="mailto:ckontzialis@lycos.com">ckontzialis@lycos.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
I use ts_type arkimex, and I want at the end of every RK stage to apply a limiter function. How can I set TS so that to apply the function at the end of a RK stage?</blockquote></div><br><div>If you are doing something at every stage, it's best to write it into your residual/right-hand-side. This is how limiters work in TVD finite volume and discontinuous Galerkin methods.</div>
<div><br></div><div>Accessing at each stage is bad because:</div><div><br></div><div>1. The meaning of what is solved in each stage varies between methods.</div><div>2. It may spoil the order of accuracy and/or other stability properties of the method.</div>
<div>3. The stages may not all be increasing or located inside the time interval.</div><div><br></div><div>If you must, you can hook into SNES directly (e.g. by defining your own "line search" that applies the filter), but I don't recommend it.</div>