[petsc-dev] preallocation checks

[Stefano Zampini]

Hi Jed,

you can use any code which calls the following sequence

ierr = MatCreate(PETSC_COMM_SELF,&mat);CHKERRQ(ierr);
ierr = MatSetSize(mat,n,m,n,m);CHKERRQ(ierr);
ierr = MatSetType(mat,MATSEQDENSE);CHKERRQ(ierr);
-> Then a loop to set values
ierr = MatAssemblyBegin(mat);CHKERRQ(ierr);
ierr = MatAssemblyEnd(mat);CHKERRQ(ierr);

If you don't insert a "dummy" call to
MatSeqDenseSetPreallocation(mat,PETSC_NULL) before doing the insertion
loop, it will give an error.

Stefano

2012/1/28 Jed Brown <jedbrown at mcs.anl.gov>

> On Sat, Jan 28, 2012 at 05:40, Stefano Zampini <stefano.zampini at gmail.com>wrote:
>
>> I'm getting errors when not preallocating for MATSEQDENSE matrices. Is it
>> what you want with dense matrices too?
>
>
> What can I run to see the error?
>



-- 
Stefano
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20120128/4740f6b7/attachment.html>