[petsc-dev] snes/examples/tutorials/ex27.c

[Matthew Knepley]

On Tue, Nov 15, 2011 at 6:06 PM, Jed Brown <jedbrown at mcs.anl.gov> wrote:

> On Tue, Nov 15, 2011 at 18:03, Matthew Knepley <knepley at gmail.com> wrote:
>
>> What are you on? It is F(x) = 0. This embedding will not work for all
>> problems, but it will work for many, and be easy. If you
>> have a more complicated system, by all means us TS, however I am pointing
>> out that there are useful regimes where this
>> kind of embedding would work.
>>
>
> Is the underlying transient problem xdot+F(x)=0 or xdot=F(x)?
>
> What happens if you enforce boundary conditions without eliminating them
> (like almost every PETSc example)?
>
> Just use the TS interface if time might be involved. Then you can do
> everything.
>

If you have a hammer, everything looks like a nail. I think you are
exaggerating the difficulty of doing this for simple cases.

   Matt


> Making assumptions about what the user meant when they were using a
> lower-level interface doesn't help. And just to get an "enclosing TS", they
> had to write some setup code. At that point, it's no more difficult to
> implement the TS callbacks than to implement the SNES callbacks.
>



-- 
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/20111115/6c172f5b/attachment.html>