Re: [petsc-users] PCFieldSplit with MatNest

2019-03-13 Thread Manuel Colera Rico via petsc-users
After adding that line the problem gets fixed. Regards, Manuel --- On 3/13/19 3:13 PM, Zhang, Junchao wrote: Manuel,   Could you try to add this line      sbaij->free_imax_ilen = PETSC_TRUE;  after line 2431 in /opt/PETSc_library/petsc-3.10.4/src/mat/impls/sbaij/seq/sbaij.c  PS: Matt, this

Re: [petsc-users] PCFieldSplit with MatNest

2019-03-13 Thread Zhang, Junchao via petsc-users
Manuel, Could you try to add this line sbaij->free_imax_ilen = PETSC_TRUE; after line 2431 in /opt/PETSc_library/petsc-3.10.4/src/mat/impls/sbaij/seq/sbaij.c PS: Matt, this bug looks unrelated to my VecRestoreArrayRead_Nest fix. --Junchao Zhang On Wed, Mar 13, 2019 at 9:05 AM Matthew

Re: [petsc-users] PCFieldSplit with MatNest

2019-03-13 Thread Lawrence Mitchell via petsc-users
> On 13 Mar 2019, at 14:04, Matthew Knepley via petsc-users > wrote: > > On Wed, Mar 13, 2019 at 9:44 AM Manuel Colera Rico via petsc-users > wrote: > Yes: > > [ 0]8416 bytes MatCreateSeqSBAIJWithArrays() line 2431 in > /opt/PETSc_library/petsc-3.10.4/src/mat/impls/sbaij/seq/sbaij.c > [ 0

Re: [petsc-users] PCFieldSplit with MatNest

2019-03-13 Thread Manuel Colera Rico via petsc-users
The warning that Leak Sanitizer gives me is not what I wrote two messages before (I apologize). It is: Direct leak of 25920 byte(s) in 4 object(s) allocated from:     #0 0x7fa97e35aa38 in __interceptor_memalign ../../../../gcc-8.1.0/libsanitizer/asan/asan_malloc_linux.cc:111     #1 0x7fa97cc812

Re: [petsc-users] PCFieldSplit with MatNest

2019-03-13 Thread Manuel Colera Rico via petsc-users
Yes: [ 0]8416 bytes MatCreateSeqSBAIJWithArrays() line 2431 in /opt/PETSc_library/petsc-3.10.4/src/mat/impls/sbaij/seq/sbaij.c [ 0]8416 bytes MatCreateSeqSBAIJWithArrays() line 2431 in /opt/PETSc_library/petsc-3.10.4/src/mat/impls/sbaij/seq/sbaij.c [ 0]4544 bytes MatCreateSeqSBAIJWithArrays() l

Re: [petsc-users] PCFieldSplit with MatNest

2019-03-13 Thread Jed Brown via petsc-users
Is there any output if you run with -malloc_dump? Manuel Colera Rico via petsc-users writes: > Hi, Junchao, > > I have installed the newest version of PETSc and it works fine. I just > get the following memory leak warning: > > Direct leak of 28608 byte(s) in 12 object(s) allocated from: >    

Re: [petsc-users] PCFieldSplit with MatNest

2019-03-13 Thread Manuel Colera Rico via petsc-users
Hi, Junchao, I have installed the newest version of PETSc and it works fine. I just get the following memory leak warning: Direct leak of 28608 byte(s) in 12 object(s) allocated from:     #0 0x7f1ddd5caa38 in __interceptor_memalign ../../../../gcc-8.1.0/libsanitizer/asan/asan_malloc_linux.cc:

Re: [petsc-users] PCFieldSplit with MatNest

2019-03-12 Thread Zhang, Junchao via petsc-users
Hi, Manuel, I recently fixed a problem in VecRestoreArrayRead. Basically, I added VecRestoreArrayRead_Nest. Could you try the master branch of PETSc to see if it fixes your problem? Thanks. --Junchao Zhang On Mon, Mar 11, 2019 at 6:56 AM Manuel Colera Rico via petsc-users mailto:petsc-use

[petsc-users] PCFieldSplit with MatNest

2019-03-11 Thread Manuel Colera Rico via petsc-users
Hello, I need to solve a 2*2 block linear system. The matrices A_00, A_01, A_10, A_11 are constructed separately via MatCreateSeqAIJWithArrays and MatCreateSeqSBAIJWithArrays. Then, I construct the full system matrix with MatCreateNest, and use MatNestGetISs and PCFieldSplitSetIS to set up th