Re: [petsc-users] IS from MatNestGetISs

2019-06-15 Thread Smith, Barry F. via petsc-users
From the manual page: "The reference count is not increased on the returned ISs", so you should not delete them Note that the rule in PETSc is items obtained with a call that contains the work Create must destroy the object, if it contains Get then it should not destroy the object.

[petsc-users] IS from MatNestGetISs

2019-06-15 Thread tangqi via petsc-users
Dear all, Do I need to delete index_set from MatNestGetISs? I get my index_set from MatNestGetISs(P, index_set, NULL); //P is a nest matrix of 3x3 When I delete index_set later for (int i=0; i<3; i++) { ISDestroy(_set[i]); } It gives me error [0]PETSC ERROR:

Re: [petsc-users] Local sizes when creating PETSC MatIS

2019-06-15 Thread Jed Brown via petsc-users
José Lorenzo via petsc-users writes: > I'm using PETSC 3.10 with 64 bits indices. > > When I run valgrind I get the following message at the end of the report, > which I don't know how to interpret: > > ==399672== Invalid read of size 8 > ==399672==at 0x5627D05: MatZeroEntries_SeqAIJ