[petsc-users] vecview problem

[Klaij, Christiaan]

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