Actual source code: ex2f.F90
1: ! Synchronized printing: Fortran Example
2: #include <petsc/finclude/petscsys.h>
3: program main
4: use petscsys
6: implicit none
7: PetscErrorCode :: ierr
8: PetscMPIInt :: rank, size
9: character(len=PETSC_MAX_PATH_LEN) :: outputString
11: ! Every PETSc program should begin with the PetscInitialize() routine.
13: PetscCallA(PetscInitialize(ierr))
15: ! The following MPI calls return the number of processes
16: ! being used and the rank of this process in the group
18: PetscCallMPIA(MPI_Comm_size(MPI_COMM_WORLD, size, ierr))
19: PetscCallMPIA(MPI_Comm_rank(MPI_COMM_WORLD, rank, ierr))
21: ! Here we would like to print only one message that represents
22: ! all the processes in the group
23: write (outputString, *) 'No of Processors = ', size, ', rank = ', rank, '\n'
24: PetscCallA(PetscPrintf(PETSC_COMM_WORLD, outputString, ierr))
26: write (outputString, *) rank, 'Synchronized Hello World\n'
27: PetscCallA(PetscSynchronizedPrintf(PETSC_COMM_WORLD, outputString, ierr))
29: write (outputString, *) rank, 'Synchronized Hello World - Part II\n'
30: PetscCallA(PetscSynchronizedPrintf(PETSC_COMM_WORLD, outputString, ierr))
31: PetscCallA(PetscSynchronizedFlush(PETSC_COMM_WORLD, PETSC_STDOUT, ierr))
33: ! Here a barrier is used to separate the two program states.
34: PetscCallMPIA(MPI_Barrier(PETSC_COMM_WORLD, ierr))
36: write (outputString, *) rank, 'Jumbled Hello World\n'
37: PetscCallA(PetscPrintf(PETSC_COMM_SELF, outputString, ierr))
39: PetscCallA(PetscFinalize(ierr))
40: end program main
42: !/*TEST
43: !
44: ! test:
45: !
46: !TEST*/