[petsc-users] PCFieldSplit with MatNest

[Manuel Colera Rico]

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:111
     #1 0x7f1ddbef1213 in PetscMallocAlign 
(/opt/PETSc_library/petsc-3.10.4/mcr_20190313/lib/libpetsc.so.3.10+0x150213)

Thank you,

Manuel

---

On 3/12/19 7:08 PM, Zhang, Junchao wrote:
> 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 
> <petsc-users at mcs.anl.gov <mailto:petsc-users at mcs.anl.gov>> wrote:
>
>     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 the PC, trying to follow the procedure described here:
>     https://www.mcs.anl.gov/petsc/petsc-current/src/snes/examples/tutorials/ex70.c.html.
>
>     However, when I run the code with Leak Sanitizer, I get the
>     following error:
>
>     =================================================================
>     ==54927==ERROR: AddressSanitizer: attempting free on address which
>     was
>     not malloc()-ed: 0x627000051ab8 in thread T0
>          #0 0x7fbd95c08f30 in __interceptor_free
>     ../../../../gcc-8.1.0/libsanitizer/asan/asan_malloc_linux.cc:66
>          #1 0x7fbd92b99dcd in PetscFreeAlign
>     (/opt/PETSc_library/petsc/manuel_OpenBLAS_petsc/lib/libpetsc.so.3.8+0x146dcd)
>          #2 0x7fbd92ce0178 in VecRestoreArray_Nest
>     (/opt/PETSc_library/petsc/manuel_OpenBLAS_petsc/lib/libpetsc.so.3.8+0x28d178)
>          #3 0x7fbd92cd627d in VecRestoreArrayRead
>     (/opt/PETSc_library/petsc/manuel_OpenBLAS_petsc/lib/libpetsc.so.3.8+0x28327d)
>          #4 0x7fbd92d1189e in VecScatterBegin_SSToSS
>     (/opt/PETSc_library/petsc/manuel_OpenBLAS_petsc/lib/libpetsc.so.3.8+0x2be89e)
>          #5 0x7fbd92d1a414 in VecScatterBegin
>     (/opt/PETSc_library/petsc/manuel_OpenBLAS_petsc/lib/libpetsc.so.3.8+0x2c7414)
>          #6 0x7fbd934a999c in PCApply_FieldSplit
>     (/opt/PETSc_library/petsc/manuel_OpenBLAS_petsc/lib/libpetsc.so.3.8+0xa5699c)
>          #7 0x7fbd93369071 in PCApply
>     (/opt/PETSc_library/petsc/manuel_OpenBLAS_petsc/lib/libpetsc.so.3.8+0x916071)
>          #8 0x7fbd934efe77 in KSPInitialResidual
>     (/opt/PETSc_library/petsc/manuel_OpenBLAS_petsc/lib/libpetsc.so.3.8+0xa9ce77)
>          #9 0x7fbd9350272c in KSPSolve_GMRES
>     (/opt/PETSc_library/petsc/manuel_OpenBLAS_petsc/lib/libpetsc.so.3.8+0xaaf72c)
>          #10 0x7fbd934e3c01 in KSPSolve
>     (/opt/PETSc_library/petsc/manuel_OpenBLAS_petsc/lib/libpetsc.so.3.8+0xa90c01)
>
>     Disabling Leak Sanitizer also outputs an "invalid pointer" error.
>
>     Did I forget something when writing the code?
>
>     Thank you,
>
>     Manuel
>
>     ---
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20190313/680c60eb/attachment.html>