<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Fri, Feb 6, 2015 at 9:35 PM, Jed Brown <span dir="ltr"><<a href="mailto:jed@jedbrown.org" target="_blank">jed@jedbrown.org</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><span class="">Barry Smith <<a href="mailto:bsmith@mcs.anl.gov">bsmith@mcs.anl.gov</a>> writes:<br>
>    I am not proposing removing the solution from TSSolve() argument<br>
>    list, far from it. I am proposing removing the TSSetSolution(),<br>
>    keeping the solution in TSSolve() as the standard approach, but<br>
>    also allowing a callback for when the user wants the DM passed to<br>
>    TS to be in charge of creating all the vectors.<br>
<br>
</span>Oh, that's fine if/when you have a use case for it.<br>
</blockquote></div><br>I needed SNESSetSolution() so I put it in. Here is the use case.</div><div class="gmail_extra"><br></div><div class="gmail_extra">I have a TS (magma dynamics), for which I would like to test the solver. I</div><div class="gmail_extra">have an explicit solution when I turn some terms off, and I test all these.</div><div class="gmail_extra">Sometimes I turn off time, so I TSGetSNES() and SNESSetSolution(), and</div><div class="gmail_extra">then call SNESComputeFunction(snes, exactSol, res). If I do not SetSolution()</div><div class="gmail_extra">first, things are not setup correctly and it bombs.</div><div class="gmail_extra"><br></div><div class="gmail_extra">   Matt<br clear="all"><div><br></div>-- <br><div class="gmail_signature">What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>-- Norbert Wiener</div>
</div></div>