[petsc-users] dmplex with block size error

[Rongliang Chen]

Dear All,

I try to use the block matrix (BAIJ) for the dmplex data structure with 
the option "-dm_mat_type baij" (the block size is 7). The code works 
fine when np = 1 but the following error comes up when np>1. And the 
code also works fine for np>1 if I set the block size to be 1. Any 
suggestions are highly appreciated.

----------------------------------------------------
[1]PETSC ERROR: PetscMallocValidate: error detected at VecAXPY_Seq() 
line 89 in /home/rlchen/soft/petsc-3.6.3/src/vec/vec/impls/seq/bvec1.c
[1]PETSC ERROR: Memory at address 0x1332571 is corrupted
[1]PETSC ERROR: Probably write past beginning or end of array
[1]PETSC ERROR: Last intact block allocated in PetscStrallocpy() line 
188 in /home/rlchen/soft/petsc-3.6.3/src/sys/utils/str.c
[1]PETSC ERROR: --------------------- Error Message 
--------------------------------------------------------------
[1]PETSC ERROR: Memory corruption: 
http://www.mcs.anl.gov/petsc/documentation/installation.html#valgrind
[1]PETSC ERROR:
[1]PETSC ERROR: See http://www.mcs.anl.gov/petsc/documentation/faq.html 
for trouble shooting.
[1]PETSC ERROR: Petsc Release Version 3.6.3, Dec, 03, 2015
[1]PETSC ERROR: ./fsi on a 64bit-debug named rlchen by rlchen Thu Aug 18 
16:42:34 2016
[1]PETSC ERROR: Configure options --download-blacs --download-scalapack 
--download-metis --download-parmetis --download-exodusii 
--download-netcdf --download-hdf5 
--with-mpi-dir=/home/rlchen/soft/Program/mpich2-shared 
--with-debugging=1 --download-fblaslapack --with-64-bit-indices
[1]PETSC ERROR: #1 PetscMallocValidate() line 136 in 
/home/rlchen/soft/petsc-3.6.3/src/sys/memory/mtr.c
[1]PETSC ERROR: #2 VecAXPY_Seq() line 89 in 
/home/rlchen/soft/petsc-3.6.3/src/vec/vec/impls/seq/bvec1.c
[1]PETSC ERROR: #3 VecAXPY() line 640 in 
/home/rlchen/soft/petsc-3.6.3/src/vec/vec/interface/rvector.c
---------------------------------------------------

Best regards,
Rongliang