<div dir="ltr">Ok, I made progress. But:<br><br><div> 1 - whatever I do, I have very slightly negative values, and therefore all my steps get rejected (values like 1e-16)</div><div> 2 - As I expected, SNES is only used with implicit methods. So if I use explicit Runge-Kutta, then there is no solution vector stored by the SNES object.</div><div><br></div><div>Reading the code for the Runge-Kutta solver, it seems that TSPostStage is where I can retrieve the current state, and TSAdaptCheckStage where I can reject it. But is this something I can rely on?</div><div><br></div><div>Thanks,</div><div><br></div><div>Pierre</div><br><div class="gmail_quote">On Wed Feb 11 2015 at 02:45:26 Jed Brown <<a href="mailto:jed@jedbrown.org" target="_blank">jed@jedbrown.org</a>> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Pierre Barbier de Reuille <<a href="mailto:pierre.barbierdereuille@gmail.com" target="_blank">pierre.barbierdereuille@<u></u>gmail<u></u>.com</a>> writes:<br>
<br>
> Ok, it seems if I set the domain error from the rhs function, it will<br>
> indeed fail and backtrack. I hope that is what was intended?<br>
<br>
Yes, the PostStep (or PostStage) callbacks are not intended for this.<br>
</blockquote></div></div>