Works fine for me: $ petscmpiexec -n 2 ./ex1f Vec Object: 2 MPI processes type: mpi Process [0] 0 0 0 0 0 Process [1] 0 0 0 0 0 ~/Src/petsc/test-dir (next=) arch-debug $ cat ./ex1f.F90 program test
implicit none #include <petsc-finclude/petscsys.h> #include <petsc-finclude/petscvec.h> MPI_Comm :: comm PetscErrorCode :: ierr Vec :: v PetscInt, parameter :: dim = 10 call PetscInitialize(PETSC_NULL_CHARACTER, ierr); CHKERRQ(ierr) comm = PETSC_COMM_WORLD call VecCreateMPI(comm, PETSC_DECIDE, dim, v, ierr); CHKERRQ(ierr) call VecView(v, PETSC_VIEWER_STDOUT_WORLD, ierr); CHKERRQ(ierr) call VecDestroy(v, ierr); CHKERRQ(ierr) call PetscFinalize(ierr); CHKERRQ(ierr) end program test ~/Src/petsc/test-dir (next=) arch-debug I cannot explain why it doesn't work for you. What about one process? Can you run in the debugger and put a break point in VecView() Barry > On Mar 16, 2015, at 8:00 PM, Adrian Croucher <[email protected]> > wrote: > > hi, > > The attached minimal C test program to create a Vec and view it produces the > output I expect: > > acro018@des108:~/software/PETSc/tests/dmda$ mpirun -np 2 testc > Vec Object: 2 MPI processes > type: mpi > Process [0] > 0 > 0 > 0 > 0 > 0 > Process [1] > 0 > 0 > 0 > 0 > 0 > > But the attached Fortran program, which I thought should do the same, > produces no output at all. Am I missing something? I'm running 'next' branch. > > Cheers, Adrian > > -- > Dr Adrian Croucher > Senior Research Fellow > Department of Engineering Science > University of Auckland, New Zealand > email: [email protected] > tel: +64 (0)9 923 84611 > > <test.F90><testc.c>
