Il Linguaggio Fortran 90/95

146 Array
! Sezione dichiarativa
IMPLICIT NONE
INTEGER , PARAMETER :: long=SELECTED_INT_KIND(9)
INTEGER(KIND=long) :: j
INTEGER :: i
REAL(KIND(1.D0)), DIMENSION(5) :: x=1.0
REAL :: inizio, fine, somma
! Sezione esecutiva
somma=0.
CALL CPU_TIME(inizio)
DO j=1_long,1000000_long
DO i=1,5
somma=somma+x(i)
END DO
somma=0.
END DO
CALL CPU_TIME(fine)
WRITE(*,’(1X,A17,G14.7)’) "DO loop: ", fine-inizio
somma=0.
CALL CPU_TIME(inizio)
DO j=1_long,1000000_long
somma=x(1)+x(2)+x(3)+x(4)+x(5)
somma=0.
END DO
CALL CPU_TIME(fine)
WRITE(*,’(1X,A17,G14.7)’) "Loop unrolling: ", fine-inizio
STOP
END PROGRAM test_loop_unrolling
compilato con il compilatore Digital Visual Fortran 6.0 e mandato in esecuzione su un PC con
processore Intel Pentium III, produce il seguente output:
DO loop: 0.1301872
Loop unrolling: 0.2002880E-01
che evidenzia i benefici associati al meccanismo di loop unrolling. Naturalmente l’operazione
di srotolamento degli array può essere estesa anche ai casi in cui le dimensioni dell’array non
consentono di esplicitarne tutti gli elementi. Allo scopo basta eseguire l’unrolling un certo
numero di volte su un numero ben preciso di elementi. Ad esempio, dato il costrutto:
n = 500
somma = 0.
DO i = 1,n
somma = somma+x(i)
END DO