[petsc-users] Problem with Mat.setpreallocationNNZ in petsc4py

[Xin Zhao]

yeap...sorry...
L=DA.createMat()

then
L.setPreallocationNNZ((3,2))

On Sat, Jul 7, 2012 at 6:01 PM, Matthew Knepley <knepley at gmail.com> wrote:

> On Sat, Jul 7, 2012 at 10:58 AM, Xin Zhao <sean.null at gmail.com> wrote:
>
>> Dear all,
>>
>> I generate a matrix L by
>> DA =PETSc.DA().create(...some...)
>> L = DA.create()
>>
>
> Is this createMatrix()? The matrix returned from a DA is already
> preallocated.
>
>    Matt
>
>
>> Then I want to preallocate memory for L
>> L.setPreallocationNNZ((3,2))
>> This works when for mpiexec -np 1
>> but it gives the error message below when mpiexec -np 4
>>
>> [3] MatAnyAIJSetPreallocation() line 311 in
>> petsc4py-1.2/src/include/custom.h
>> [3] Operation done in wrong order
>> [3] matrix is already preallocated
>>
>> How to solve this?
>>
>> Thanks in advance.
>>
>> Cheers,
>> Xin
>>
>
>
>
> --
> 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-users/attachments/20120707/167997b5/attachment.html>