C******************************************************************************
|
C******************************************************************************
|
C FILE: omp_workshare2.f
|
C FILE: omp_workshare2.f
|
C DESCRIPTION:
|
C DESCRIPTION:
|
C OpenMP Example - Sections Work-sharing - Fortran Version
|
C OpenMP Example - Sections Work-sharing - Fortran Version
|
C In this example, the OpenMP SECTION directive is used to assign
|
C In this example, the OpenMP SECTION directive is used to assign
|
C different array operations to threads that execute a SECTION. Each
|
C different array operations to threads that execute a SECTION. Each
|
C thread receives its own copy of the result array to work with.
|
C thread receives its own copy of the result array to work with.
|
C AUTHOR: Blaise Barney 5/99
|
C AUTHOR: Blaise Barney 5/99
|
C LAST REVISED: 01/09/04
|
C LAST REVISED: 01/09/04
|
C******************************************************************************
|
C******************************************************************************
|
|
|
PROGRAM WORKSHARE2
|
PROGRAM WORKSHARE2
|
|
|
INTEGER N, I, NTHREADS, TID, OMP_GET_NUM_THREADS,
|
INTEGER N, I, NTHREADS, TID, OMP_GET_NUM_THREADS,
|
+ OMP_GET_THREAD_NUM
|
+ OMP_GET_THREAD_NUM
|
PARAMETER (N=50)
|
PARAMETER (N=50)
|
REAL A(N), B(N), C(N)
|
REAL A(N), B(N), C(N)
|
|
|
! Some initializations
|
! Some initializations
|
DO I = 1, N
|
DO I = 1, N
|
A(I) = I * 1.0
|
A(I) = I * 1.0
|
B(I) = A(I)
|
B(I) = A(I)
|
ENDDO
|
ENDDO
|
|
|
!$OMP PARALLEL SHARED(A,B,NTHREADS), PRIVATE(C,I,TID)
|
!$OMP PARALLEL SHARED(A,B,NTHREADS), PRIVATE(C,I,TID)
|
TID = OMP_GET_THREAD_NUM()
|
TID = OMP_GET_THREAD_NUM()
|
IF (TID .EQ. 0) THEN
|
IF (TID .EQ. 0) THEN
|
NTHREADS = OMP_GET_NUM_THREADS()
|
NTHREADS = OMP_GET_NUM_THREADS()
|
PRINT *, 'Number of threads =', NTHREADS
|
PRINT *, 'Number of threads =', NTHREADS
|
END IF
|
END IF
|
PRINT *, 'Thread',TID,' starting...'
|
PRINT *, 'Thread',TID,' starting...'
|
|
|
!$OMP SECTIONS
|
!$OMP SECTIONS
|
|
|
!$OMP SECTION
|
!$OMP SECTION
|
PRINT *, 'Thread',TID,' doing section 1'
|
PRINT *, 'Thread',TID,' doing section 1'
|
DO I = 1, N
|
DO I = 1, N
|
C(I) = A(I) + B(I)
|
C(I) = A(I) + B(I)
|
WRITE(*,100) TID,I,C(I)
|
WRITE(*,100) TID,I,C(I)
|
100 FORMAT(' Thread',I2,': C(',I2,')=',F8.2)
|
100 FORMAT(' Thread',I2,': C(',I2,')=',F8.2)
|
ENDDO
|
ENDDO
|
|
|
!$OMP SECTION
|
!$OMP SECTION
|
PRINT *, 'Thread',TID,' doing section 2'
|
PRINT *, 'Thread',TID,' doing section 2'
|
DO I = 1+N/2, N
|
DO I = 1+N/2, N
|
C(I) = A(I) * B(I)
|
C(I) = A(I) * B(I)
|
WRITE(*,100) TID,I,C(I)
|
WRITE(*,100) TID,I,C(I)
|
ENDDO
|
ENDDO
|
|
|
!$OMP END SECTIONS NOWAIT
|
!$OMP END SECTIONS NOWAIT
|
|
|
PRINT *, 'Thread',TID,' done.'
|
PRINT *, 'Thread',TID,' done.'
|
|
|
!$OMP END PARALLEL
|
!$OMP END PARALLEL
|
|
|
END
|
END
|
|
|
