| ! { dg-do run } |
| SUBROUTINE SUBDOMAIN(X, ISTART, IPOINTS) |
| INTEGER ISTART, IPOINTS |
| REAL X(*) |
| INTEGER I |
| DO 100 I=1,IPOINTS |
| X(ISTART+I) = 123.456 |
| 100 CONTINUE |
| END SUBROUTINE SUBDOMAIN |
| SUBROUTINE SUB(X, NPOINTS) |
| INCLUDE "omp_lib.h" ! or USE OMP_LIB |
| REAL X(*) |
| INTEGER NPOINTS |
| INTEGER IAM, NT, IPOINTS, ISTART |
| !$OMP PARALLEL DEFAULT(PRIVATE) SHARED(X,NPOINTS) |
| IAM = OMP_GET_THREAD_NUM() |
| NT = OMP_GET_NUM_THREADS() |
| IPOINTS = NPOINTS/NT |
| ISTART = IAM * IPOINTS |
| IF (IAM .EQ. NT-1) THEN |
| IPOINTS = NPOINTS - ISTART |
| ENDIF |
| CALL SUBDOMAIN(X,ISTART,IPOINTS) |
| !$OMP END PARALLEL |
| END SUBROUTINE SUB |
| PROGRAM A4 |
| REAL ARRAY(10000) |
| CALL SUB(ARRAY, 10000) |
| END PROGRAM A4 |
