On Tue, Jul 5, 2011 at 09:04, Stefan Kolb <skolb at rocketmail.com> wrote:
> I tried to set the local size of the matrix by hand with the command > MatSetSizes(A,50, 100,N,M) but then i get an error from PETSc. How can i > tell PETSc to store the first 50 rows on process 0 and the next 50 rows on > process 1? What you did is correct, but MatGetOwnershipRangeColumn() returns the ownership range of a vector that the matrix can be applied to. See the this man page for details on how MPIAIJ matrices are stored: http://www.mcs.anl.gov/petsc/petsc-2/snapshots/petsc-dev/docs/manualpages/Mat/MatMPIAIJSetPreallocation.html -------------- next part -------------- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20110705/1f389a6b/attachment.htm>
