I'm having a segmentation fault with vecview in fortran90
with this program. Am I doing something wrong?
$ cat bug.F90
module bug
use petscksp
implicit none
#include "finclude/petsckspdef.h"
PetscErrorCode, public :: ierr
Vec, private :: x
public bugSubroutine
contains
subroutine bugSubroutine()
call setVector()
call VecView(x,PETSC_VIEWER_DEFAULT,ierr); CHKERRQ(ierr)
end subroutine bugSubroutine
subroutine setVector()
call VecCreate(PETSC_COMM_WORLD,x,ierr); CHKERRQ(ierr)
call VecSetSizes(x,PETSC_DECIDE,5,ierr); CHKERRQ(ierr)
call VecSetType(x,VECMPI,ierr); CHKERRQ(ierr)
call VecView(x,PETSC_VIEWER_DEFAULT,ierr); CHKERRQ(ierr)
end subroutine setVector
end module bug
program testBug
use bug
use petscksp
implicit none
#include "finclude/petsckspdef.h"
call PetscInitialize(PETSC_NULL_CHARACTER,ierr)
call bugSubroutine();
call PetscFinalize(ierr)
end program testBug
$ mpiexec -n 1 ./bug
Vector Object: 1 MPI processes
type: mpi
Process [0]
0
0
0
0
0
[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/documentation/faq.html#valgrind[0]PETSC ERROR: or
try http://valgrind.org on GNU/linux and Apple Mac OS X to find memory
corruption errors
[0]PETSC ERROR: likely location of problem given in stack below
[0]PETSC ERROR: --------------------- Stack Frames
------------------------------------
[0]PETSC ERROR: Note: The EXACT line numbers in the stack are not available,
[0]PETSC ERROR: INSTEAD the line number of the start of the function
[0]PETSC ERROR: is given.
[0]PETSC ERROR: [0] VecView line 747
/home/CKlaij/ReFRESCO/Dev/trunk/Libs/build/petsc-3.3-p1/src/vec/vec/interface/vector.c
[0]PETSC ERROR: --------------------- Error Message
------------------------------------
[0]PETSC ERROR: Signal received!
[0]PETSC ERROR:
------------------------------------------------------------------------
[0]PETSC ERROR: Petsc Release Version 3.3.0, Patch 1, Fri Jun 15 09:30:49 CDT
2012
[0]PETSC ERROR: See docs/changes/index.html for recent updates.
[0]PETSC ERROR: See docs/faq.html for hints about trouble shooting.
[0]PETSC ERROR: See docs/index.html for manual pages.
[0]PETSC ERROR:
------------------------------------------------------------------------
[0]PETSC ERROR: ./bug on a linux_64b named lin0133 by cklaij Mon Jul 9
15:46:34 2012
[0]PETSC ERROR: Libraries linked from
/home/CKlaij/ReFRESCO/Dev/trunk/Libs/build/petsc-3.3-p1/linux_64bit_debug/lib
[0]PETSC ERROR: Configure run at Wed Jun 20 12:08:20 2012
[0]PETSC ERROR: Configure options
--with-mpi-dir=/opt/refresco/libraries_cklaij/openmpi-1.4.5
--with-clanguage=c++ --with-x=1 --with-debugging=1
--with-hypre-include=/opt/refresco/libraries_cklaij/hypre-2.7.0b/include
--with-hypre-lib=/opt/refresco/libraries_cklaij/hypre-2.7.0b/lib/libHYPRE.a
--with-ml-include=/opt/refresco/libraries_cklaij/ml-6.2/include
--with-ml-lib=/opt/refresco/libraries_cklaij/ml-6.2/lib/libml.a
--with-blas-lapack-dir=/opt/intel/mkl
[0]PETSC ERROR:
------------------------------------------------------------------------
[0]PETSC ERROR: User provided function() line 0 in unknown directory unknown
file
--------------------------------------------------------------------------
MPI_ABORT was invoked on rank 0 in communicator MPI_COMM_WORLD
with errorcode 59.
NOTE: invoking MPI_ABORT causes Open MPI to kill all MPI processes.
You may or may not see output from other processes, depending on
exactly when Open MPI kills them.
--------------------------------------------------------------------------
(gdb) run
Starting program: /home/CKlaij/Programming/PETSc/Laplace/bug
[Thread debugging using libthread_db enabled]
Vector Object: 1 MPI processes
type: mpi
Process [0]
0
0
0
0
0
Program received signal SIGSEGV, Segmentation fault.
0x000000000043f7ee in VecView (vec=0xd8b1f0, viewer=0x69706d00000000)
at
/home/CKlaij/ReFRESCO/Dev/trunk/Libs/build/petsc-3.3-p1/src/vec/vec/interface/vector.c:753
753 PetscValidHeaderSpecific(viewer,PETSC_VIEWER_CLASSID,2);
(gdb) bt
#0 0x000000000043f7ee in VecView (vec=0xd8b1f0, viewer=0x69706d00000000)
at
/home/CKlaij/ReFRESCO/Dev/trunk/Libs/build/petsc-3.3-p1/src/vec/vec/interface/vector.c:753
#1 0x00000000004319e8 in vecview_ (x=0xb56e20, vin=0x822788, ierr=0xb56e38)
at
/home/CKlaij/ReFRESCO/Dev/trunk/Libs/build/petsc-3.3-p1/src/vec/vec/interface/ftn-custom/zvectorf.c:56
#2 0x000000000042c936 in bug_mp_bugsubroutine_ ()
#3 0x000000000042cafb in testbug () at bug.F90:37
dr. ir. Christiaan Klaij
CFD Researcher
Research & Development
E mailto:C.Klaij at marin.nl
T +31 317 49 33 44
MARIN
2, Haagsteeg, P.O. Box 28, 6700 AA Wageningen, The Netherlands
T +31 317 49 39 11, F +31 317 49 32 45, I www.marin.nl
