Segmentation Violation with getNestLocalISs

[Benoit NENNIG]

Hi,

I continue to work with nest matrix and I have now another problem with getNestLocalISs. I would like to use this function to get the local index set of the submatrices. It works well with absolute index set getNestISs but crash with Segmentation Violation for the local index.

C=PETSc.Mat().createNest([[A, A], [None, 2*A]])
ind=C.getNestISs() 
ind[0][0].getIndices() 

yields array([0, 1, 2, 3], dtype=int32) which is ok, but

indl=C.getNestLocalISs() # seems ok, I got a tuple of length 2
indl[0][0] # yield <petsc4py.PETSc.IS object at 0x7f64cd6517d0>
indl[0][0].getIndices()  # Caught signal number 11 SEGV: Segmentation Violation, probably memory access out of range

Thank you for your help.

Benoit

info : petsc4py, commit 3619809 maint branch, arch-linux2-c-optim, complex