[petsc-dev] SNES for 2x2 problems

[Jed Brown]

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.

Also, please remember to test that complex actually works before declaring
it fixed. I just cleaned up a bunch more stuff, including pain invoked by
SNES operating in PetscScalar for those 2x2 point location problems.

https://bitbucket.org/petsc/petsc-dev/commits/0213b4504a412b7304d6032d6047079603a8c2c8


To be vectorization-friendly, I would write point location to operate
directly on a batch of points, wrapped in the little Newton loop (not
SNES). Maybe this implementation of yours is a reference that you're
planning to replace?
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20130102/8107de22/attachment.html>