Also suggest using latest version 3.7.
Satish
________________________________
From: Ilyas YILMAZ<mailto:[email protected]>
Sent: 4/30/2016 10:04 AM
To: PETSc users list<mailto:[email protected]>
Subject: [petsc-users] a crash due to memory issue
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.
