[petsc-users] a crash due to memory issue

[Ilyas YILMAZ]

Hello,

The code segment I wrote based on "src/dm/da/examples/tutorials/ex2.c"
crashes when destroying things / freeing memory as given below.
I can't figure out what I'm missing? Any comments are welcome. (Petsc
3.1.p8, intel fortran compiler). The code and error are below.

        !...Create scatter from global DA parallel vector to local vector
that contains all entries
        call DAGlobalToNaturalAllCreate(da,tolocalall,ierr)
        call DAGlobalToNaturalAllCreate(da,tolocalall2,ierr)
        call DAGlobalToNaturalAllCreate(da,tolocalall3,ierr)

        call VecCreateSeq(PETSC_COMM_SELF,im*jm*km,CSRENOMlocalall,ierr)
        call
VecScatterBegin(tolocalall,CSRENOM,CSRENOMlocalall,INSERT_VALUES,SCATTER_FORWARD,ierr)

        call VecScatterEnd(
tolocalall,CSRENOM,CSRENOMlocalall,INSERT_VALUES,SCATTER_FORWARD,ierr)

        call VecCreateSeq(PETSC_COMM_SELF,im*jm*km,CSDENOMlocalall,ierr)
        call
VecScatterBegin(tolocalall2,CSDENOM,CSDENOMlocalall,INSERT_VALUES,SCATTER_FORWARD,ierr)

        call VecScatterEnd(
tolocalall2,CSDENOM,CSDENOMlocalall,INSERT_VALUES,SCATTER_FORWARD,ierr)

        call VecCreateSeq(PETSC_COMM_SELF,im*jm*km,CSlocalall,ierr)
        call
VecScatterBegin(tolocalall3,CS,CSlocalall,INSERT_VALUES,SCATTER_FORWARD,ierr)

        call VecScatterEnd(
tolocalall3,CS,CSlocalall,INSERT_VALUES,SCATTER_FORWARD,ierr)


        call
VecGetArray(CSRENOMlocalall,scaCSRENOMlocalall,idCSRENOMlocalall,ierr)
        call
VecGetArray(CSDENOMlocalall,scaCSDENOMlocalall,idCSDENOMlocalall,ierr)
        call VecGetArray(     CSlocalall,scaCSlocalall
,idCSlocalall     ,ierr)


                          SOME WORK HERE


        call VecRestoreArray(CSlocalall     ,scaCSlocalall
,idCSlocalall,ierr)
        call
VecRestoreArray(CSRENOMlocalall,scaCSRENOMlocalall,idCSRENOMlocalall,ierr)
        call
VecRestoreArray(CSDENOMlocalall,scaCSDENOMlocalall,idCSDENOMlocalall,ierr)


        !...scatter back to global vector
        call DANaturalAllToGlobalCreate(da,fromlocalall,ierr)
        call DANaturalAllToGlobalCreate(da,fromlocalall2,ierr)
        call DANaturalAllToGlobalCreate(da,fromlocalall3,ierr)

        call
VecScatterBegin(fromlocalall,CSRENOMlocalall,CSRENOM,INSERT_VALUES,SCATTER_FORWARD,ierr)

        call VecScatterEnd(
fromlocalall,CSRENOMlocalall,CSRENOM,INSERT_VALUES,SCATTER_FORWARD,ierr)

        call
VecScatterBegin(fromlocalall2,CSDENOMlocalall,CSDENOM,INSERT_VALUES,SCATTER_FORWARD,ierr)

        call VecScatterEnd(
fromlocalall2,CSDENOMlocalall,CSDENOM,INSERT_VALUES,SCATTER_FORWARD,ierr)

        call
VecScatterBegin(fromlocalall3,CSlocalall,CS,INSERT_VALUES,SCATTER_FORWARD,ierr)

        call VecScatterEnd(
fromlocalall3,CSlocalall,CS,INSERT_VALUES,SCATTER_FORWARD,ierr)

        !..free memory
       call VecScatterDestroy(fromlocalall,ierr)
       call VecScatterDestroy(fromlocalall2,ierr)
       call VecScatterDestroy(fromlocalall3,ierr)

      call VecScatterDestroy(tolocalall,ierr)
      call VecScatterDestroy(tolocalall2,ierr)
      call VecScatterDestroy(tolocalall3,ierr)

      call VecDestroy(CSDENOMlocalall,ierr)
      call VecDestroy(CSRENOMlocalall,ierr)
      call VecDestroy(CSlocalall,ierr)


.
.
.
.
.
.
 ...................... SGS model is called
 MUT computed
 MUT computed
 MUT computed
 MUT computed
[1]PETSC ERROR:
------------------------------------------------------------------------
[1]PETSC ERROR: [0]PETSC ERROR:
------------------------------------------------------------------------
[0]PETSC ERROR: Caught signal number 11 SEGV: Segmentation Violation,
probably memory access out of range
[0]PETSC ERROR: Try option -start_in_debugger or -on_error_attach_debugger
[0]PETSC ERROR: or see
http://www.mcs.anl.gov/petsc/petsc-as/documentation/troubleshooting.html#Signal[0]PETSC
ERROR: or try http://valgrind.org on GNU/linux and Apple Mac OS X to find
memory corruption errors
[0]PETSC ERROR: configure using --with-debugging=yes, recompile, link, and
run
[0]PETSC ERROR: to get more information on the crash.
.
.
.
.
.

Thank you

Ilyas.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20160430/d7cfb286/attachment.html>