<div dir="ltr">Run with valgrind and send the log.<div><br></div><div>  Thanks,</div><div><br></div><div>     Matt</div></div><div class="gmail_extra"><br><div class="gmail_quote">On Thu, Aug 18, 2016 at 3:59 AM, Rongliang Chen <span dir="ltr"><<a href="mailto:rongliang.chan@gmail.com" target="_blank">rongliang.chan@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Dear All,<br>
<br>
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.<br>
<br>
------------------------------<wbr>----------------------<br>
[1]PETSC ERROR: PetscMallocValidate: error detected at VecAXPY_Seq() line 89 in /home/rlchen/soft/petsc-3.6.3/<wbr>src/vec/vec/impls/seq/bvec1.c<br>
[1]PETSC ERROR: Memory at address 0x1332571 is corrupted<br>
[1]PETSC ERROR: Probably write past beginning or end of array<br>
[1]PETSC ERROR: Last intact block allocated in PetscStrallocpy() line 188 in /home/rlchen/soft/petsc-3.6.3/<wbr>src/sys/utils/str.c<br>
[1]PETSC ERROR: --------------------- Error Message ------------------------------<wbr>------------------------------<wbr>--<br>
[1]PETSC ERROR: Memory corruption: <a href="http://www.mcs.anl.gov/petsc/documentation/installation.html#valgrind" rel="noreferrer" target="_blank">http://www.mcs.anl.gov/petsc/d<wbr>ocumentation/installation.html<wbr>#valgrind</a><br>
[1]PETSC ERROR:<br>
[1]PETSC ERROR: See <a href="http://www.mcs.anl.gov/petsc/documentation/faq.html" rel="noreferrer" target="_blank">http://www.mcs.anl.gov/petsc/d<wbr>ocumentation/faq.html</a> for trouble shooting.<br>
[1]PETSC ERROR: Petsc Release Version 3.6.3, Dec, 03, 2015<br>
[1]PETSC ERROR: ./fsi on a 64bit-debug named rlchen by rlchen Thu Aug 18 16:42:34 2016<br>
[1]PETSC ERROR: Configure options --download-blacs --download-scalapack --download-metis --download-parmetis --download-exodusii --download-netcdf --download-hdf5 --with-mpi-dir=/home/rlchen/so<wbr>ft/Program/mpich2-shared --with-debugging=1 --download-fblaslapack --with-64-bit-indices<br>
[1]PETSC ERROR: #1 PetscMallocValidate() line 136 in /home/rlchen/soft/petsc-3.6.3/<wbr>src/sys/memory/mtr.c<br>
[1]PETSC ERROR: #2 VecAXPY_Seq() line 89 in /home/rlchen/soft/petsc-3.6.3/<wbr>src/vec/vec/impls/seq/bvec1.c<br>
[1]PETSC ERROR: #3 VecAXPY() line 640 in /home/rlchen/soft/petsc-3.6.3/<wbr>src/vec/vec/interface/rvector.<wbr>c<br>
------------------------------<wbr>---------------------<br>
<br>
Best regards,<br>
Rongliang<br>
</blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature" data-smartmail="gmail_signature">What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>-- Norbert Wiener</div>
</div>