[petsc-users] dmplex with block size error

[Matthew Knepley]

Run with valgrind and send the log.

  Thanks,

     Matt

On Thu, Aug 18, 2016 at 3:59 AM, Rongliang Chen <rongliang.chan at gmail.com>
wrote:

> 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/d
> ocumentation/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
>



-- 
What most experimenters take for granted before they begin their
experiments is infinitely more interesting than any results to which their
experiments lead.
-- Norbert Wiener
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20160818/b4d863b5/attachment.html>