[petsc-dev] SNES for 2x2 problems

Jed Brown jedbrown at mcs.anl.gov
Wed Jan 2 16:19:37 CST 2013


On Wed, Jan 2, 2013 at 4:16 PM, Matthew Knepley <knepley at gmail.com> wrote:

> On Wed, Jan 2, 2013 at 3:21 PM, Jed Brown <jedbrown at mcs.anl.gov> wrote:
>
>> Matt, you know how we always tell people not to use KSP (or even Lapack)
>> for 2x2 problems? Surely the same advice applies to SNES, but your point
>> location code is doing it.
>>
>
> Yep. The wrapping overhead is far less with some Newton iterations
>

Bogus


> , and we have no good alternative.
>

I posit that it's fewer lines of code (and much faster) to just write out
the Newton iteration on the spot. You don't need the full algorithmic
infrastructure for point location.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20130102/9572e815/attachment.html>


More information about the petsc-dev mailing list