[petsc-users] Preallocation woes again

[Jed Brown]

On Fri, Jan 27, 2012 at 12:14, Matthew Knepley <knepley at gmail.com> wrote:

> Basically. You have to store indices while counting if you don't want to
> overcount. This is the part he glosses over.


No, he has structure implicitly through the mesh. The mesh overcounts some
points, but if you can easily determine how much it overcounts, then you
have an efficient way to compute a non-redundant count. It's actually
straightforward for low-order simplices, but not for more general bases.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20120127/5d42e851/attachment.htm>