[petsc-dev] TSSolve web page JED and EMIL READ THIS!

[Matthew Knepley]

On Fri, Feb 6, 2015 at 9:35 PM, Jed Brown <jed at jedbrown.org> wrote:

> Barry Smith <bsmith at mcs.anl.gov> writes:
> >    I am not proposing removing the solution from TSSolve() argument
> >    list, far from it. I am proposing removing the TSSetSolution(),
> >    keeping the solution in TSSolve() as the standard approach, but
> >    also allowing a callback for when the user wants the DM passed to
> >    TS to be in charge of creating all the vectors.
>
> Oh, that's fine if/when you have a use case for it.
>

I needed SNESSetSolution() so I put it in. Here is the use case.

I have a TS (magma dynamics), for which I would like to test the solver. I
have an explicit solution when I turn some terms off, and I test all these.
Sometimes I turn off time, so I TSGetSNES() and SNESSetSolution(), and
then call SNESComputeFunction(snes, exactSol, res). If I do not
SetSolution()
first, things are not setup correctly and it bombs.

   Matt

-- 
What most experimenters take for granted before they begin their
experiments is infinitely more interesting than any results to which their
experiments lead.
-- Norbert Wiener
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20150207/6d257520/attachment.html>