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

Matthew Knepley knepley at gmail.com
Tue Nov 15 17:34:31 CST 2011


On Tue, Nov 15, 2011 at 3:44 PM, Barry Smith <bsmith at mcs.anl.gov> wrote:

>
> On Nov 15, 2011, at 3:40 PM, Jed Brown wrote:
>
> > On Tue, Nov 15, 2011 at 15:36, Barry Smith <bsmith at mcs.anl.gov> wrote:
> > I think it is best if the TS example be cleaned up to do everything the
> SNES example did and not have both around. People will stumble upon the
> SNES one and use it as a template for their work when they should use the
> TS one.
> >
> > Okay, I'll document the options to reproduce the paper's results.
> >
> >
> >   I would like the TS pseudo time-stepping stuff to actually be a SNES
> solver. That is SNESType SNESPSEUDOTS or something and not be part of TS.
> This is controversial hence nothing has been done to make it this way. Some
> people may argue it is the wrong approach.
> >
> > Well, the way in which the transient term is added is really quite
> important. Especially since Ptc often uses some sort of local time
> stepping, I don't see a good way to produce that without creating the TS
> interface. Since the interface is TS, I think it belongs there.
>
>   How about a SNES level interface that backs into the TS for the
> implementation? The thing is "if one is solving a nonlinear problem they
> should use the nonlinear solver interface" not a completely different
> interface that exists for a different purpose.


I am not sure I like this since it creates a circular dependency in the
libraries.

   Matt


>
>   Barry
>
>
>


-- 
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/2edb4a21/attachment.html>


More information about the petsc-dev mailing list