Parallelization Examples/2

• Inner loop strip mining

DO 21 K=1,N

DO 21 J=1,N

DO 21 I=1,N

A(I,J,K)=A(I,J-1,K)+A(I,J,K)

21 CONTINUE

C$OMP PARALLEL PRIVATE(NTH,NPI,K,J,I)

NTH = OMP_NUM_THREADS()

NPI = N/NTH

C$OMP DO

DO M=0,NTH-1

IS = M*NPI

IE = IS + NPI

IS = IS + 1

IF(IE.GT.N) IE = N

DO 21 K=1,N

DO 21 J=1,N

DO 21 I=IS,IE

A(I,J,K)=A(I,J-1,K)+A(I,J,K)

21 CONTINUE

ENDDO

C$OMP END PARALLEL

Previous slide

Next slide

Back to first slide

View graphic version