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>