[petsc-dev] preallocation checks

[Jed Brown]

On Sat, Jan 28, 2012 at 10:11, Stefano Zampini <stefano.zampini at gmail.com>wrote:

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

I couldn't reproduce this with a modified
src/ksp/ksp/examples/tutorials/ex2.c

MatSetUpPreallocation() should be sufficient in other cases.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20120128/3bee4073/attachment.html>