[petsc-dev] preallocation

[Matthew Knepley]

On Tue, Sep 20, 2011 at 2:44 AM, Jed Brown <jedbrown at mcs.anl.gov> wrote:

> I'm so tired of calling a different preallocation variant for each matrix
> format.
>
> MatMPIBAIJSetPreallocation() has enough information to preallocate
> {Seq,MPI}{AIJ,BAIJ}. Should we make this one call enough?


I disagree with this. I think the correct reduction in complexity is for the
DM to call all relevant versions.

   Matt


> Downsides:
>
> If the user wants to allocate some special pattern within blocks, this
> would fill it all in, so they would have to check the types and call only
> the version that they really want.
>
> Old code calls all these functions which could mean a lot of reallocation,
> they can change the code to only call this one function.
>



-- 
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/20110920/6b39d4ce/attachment.html>