[petsc-users] MATIS and DMDA in petsc4py
Loic Gouarin
loic.gouarin at math.u-psud.fr
Thu Apr 26 03:05:31 CDT 2018
The problem is solved with the new 3.9.0 release on conda-forge.
Thanks Lisandro !!
Loic
Le 23/04/2018 à 13:25, Loic Gouarin a écrit :
> Thanks Lawrence for your reply.
>
> Le 23/04/2018 à 12:36, Lawrence Mitchell a écrit :
>>> On 23 Apr 2018, at 11:32, Loic Gouarin<loic.gouarin at math.u-psud.fr> wrote:
>>>
>>> Hi,
>>>
>>> I try to use MATIS from petsc4py but I am not able to use setValuesLocal (I have a segmentation fault). And I don't understand why.
>>>
>>> Could you tell me if I did a mistake or if it's a bug ?
>> This code works for me with:
>>
>> PETSc: 7bfdb6492ce
>> petsc4py: 9137bf78e
>>
>> Perhaps you have a mismatching petsc4py/PETSc combo?
> The problem probably comes from the conda package. I reinstalled a new
> environment and tried again but I still have this issue.
>
> Did you compile petsc and petsc4py ?
>
> The versions that I use are those given on conda-forge : petsc-3.8.4
> and petsc4py-3.8.1.
>
> I did the following thing
>
> conda create -n petsc_tmp
> source activate petsc_tmp
> conda install petsc4py -c conda-forge
>
> and execute my script.
>
> Loic
>> Lawrence
>>
>>> My example
>>> =======
>>> from petsc4py import PETSc
>>> import numpy as np
>>>
>>> n = 10
>>> da = PETSc.DMDA().create([n, n], dof=1, stencil_width=1)
>>> da.setMatType(PETSc.Mat.Type.IS)
>>>
>>> A = da.createMatrix()
>>> #ltog = da.getLGMap()
>>> #A.setLGMap(ltog, ltog)
>>>
>>> elem = da.getElements()
>>> melem = np.ones((4, 4))
>>>
>>> for e in elem:
>>> A.setValuesLocal(e, e, melem, PETSc.InsertMode.ADD_VALUES)
>>>
>>> Thanks,
>>> Loic
>>>
>>> --
>>> Tel: 01 69 15 60 14
>>> http://www.math.u-psud.fr/~gouarin
>>> https://github.com/gouarin
>>>
>
> --
> Tel: 01 69 15 60 14
> http://www.math.u-psud.fr/~gouarin
> https://github.com/gouarin
--
Tel: 01 69 15 60 14
http://www.math.u-psud.fr/~gouarin
https://github.com/gouarin
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20180426/a21a0db0/attachment.html>
More information about the petsc-users
mailing list