[petsc-dev] TSSolve web page JED and EMIL READ THIS!
Barry Smith
bsmith at mcs.anl.gov
Fri Feb 6 21:23:38 CST 2015
> On Feb 6, 2015, at 9:16 PM, Jed Brown <jed at jedbrown.org> wrote:
>
> Barry Smith <bsmith at mcs.anl.gov> writes:
>
>>> On Feb 6, 2015, at 8:59 PM, Jed Brown <jed at jedbrown.org> wrote:
>>>
>>> Barry Smith <bsmith at mcs.anl.gov> writes:
>>>> So going back to my email, now with more clarity, do we get rid of
>>>> the XXXSetSolution/RightHandSide() paradigm; this would void
>>>> questions like Mark's about what the differences are between
>>>> them**. But do we lose any useful functionality?
>>>
>>> I think the main possible loss is if we want an automatic
>>> spatially-adaptive TSSolve. The way I see it, the caller must provide
>>> an initial condition (on some grid)
>>
>> They could instead provide a "set initial conditions" function that the TS solver calls instead in that case.
>
> It's ugly and complicated state management if we were to require that
> the user provide their initial conditions in a callback.
Nonsense. They could always put the initial conditions in an appropriate vector and the callback just be a vec copy, then you get back the effect of TSSet"Initial"Solution()
Maybe I'll just mark the XXXSetSolution() routines as depreciated :-)
More information about the petsc-dev
mailing list