[petsc-dev] preallocation

Jed Brown jedbrown at mcs.anl.gov
Tue Sep 20 14:01:11 CDT 2011


On Tue, Sep 20, 2011 at 19:04, Barry Smith <bsmith at mcs.anl.gov> wrote:

> Why not write a wrapper preallocator that then calls all the known ones,
> for example
>
>    MatPreallocateJed(.....)
>    {
>        ierr = MatPreallocateSeqAIJ(....)
>        ierr = MatPreallocate(MPIAIJ(....)
>
>    }
>

And make this a public function, MatXAIJPreallocate(A,bs,dnz,dnnz,onz,onnz)
? That sounds fine to me.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20110920/b7f47232/attachment.html>


More information about the petsc-dev mailing list