Data Initialization: Example

• The F90 allocatable arrays should be initialized explicitly to place the pages:

PARAMETER (PAGESIZE=16384, SIZEDBLE=8)

REAL, ALLOCATABLE:: WORK(:)

READ *, MXRVAR

ALLOCATE(WORK(MXRVAR))

!$OMP PARALLEL DEFAULT(PRIVATE) SHARED(WORK,MXRVAR)

ID = OMP_GET_THREAD_NUM()

NP = OMP_GET_NUM_THREADS()

NBLOCK = MXRVAR/NP

NPAGES = MXRVAR*SIZEDBLE/PAGESIZE

NLOWER = ID*NBLOCK

NUPPER = NLOWER + NBLOCK

NLOWER = NLOWER + 1

NUPPER = MIN(NUPPER,MXRVAR)

c round-robin distribution of pages:

DO I=ID,NPAGES-NP+1,NP

WORK(1+I*PAGESIZE/SIZEDBLE) = 0.0

ENDDO

!$OMP END PARALLEL

Previous slide

Next slide

Back to first slide

View graphic version