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

Matthew Knepley knepley at gmail.com
Sat Feb 7 00:53:33 CST 2015

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
first, things are not setup correctly and it bombs.


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>

More information about the petsc-dev mailing list