[petsc-dev] preallocation

[Barry Smith]

On Sep 20, 2011, at 2:03 PM, Matthew Knepley wrote:

> On Tue, Sep 20, 2011 at 2:01 PM, Jed Brown <jedbrown at mcs.anl.gov> wrote:
> 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.
> 
> Yes, this is fine with me.
> 
>    Matt

   Ok, and make it support SBAIJ also?

   Barry

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