[petsc-users] Preallocation woes again

[Matthew Knepley]

On Fri, Jan 27, 2012 at 1:12 PM, Jed Brown <jedbrown at mcs.anl.gov> wrote:

> 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.
>

I am doing completely general, in parallel.

   Matt

-- 
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-users/attachments/20120127/f614d725/attachment.htm>