[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