[petsc-dev] preallocation

[Matthew Knepley]

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

-- 
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/d076056e/attachment.html>