[petsc-dev] preallocation checks

[Stefano Zampini]

I was able to reproduce the error using
src/ksp/ksp/examples/tutorials/ex2.c using the patch attached.
Please see the error log obtained by executing with -n 10 -m 10 by command
line.
If you uncomment MatGetOwnershipRange, the code will work fine

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

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



-- 
Stefano
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20120128/e718b604/attachment.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: errlog_n10m10
Type: application/octet-stream
Size: 2945 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20120128/e718b604/attachment.obj>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: patch_ex2
Type: application/octet-stream
Size: 1775 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20120128/e718b604/attachment-0001.obj>