Il Linguaggio Fortran 90/95

3.8 Ottimizzazione delle operazioni con array 141
PROGRAM prova_ordinamento
IMPLICIT NONE
INTEGER :: i,j
REAL,DIMENSION(5000,5000) :: array1=1.0,array2=2.0,array3=0.0
REAL :: inizio,fine
CALL CPU_TIME(inizio)
DO j=1,5000
DO i=1,5000
array3(i,j) = array1(i,j)+array2(i,j)
END DO
END DO
CALL CPU_TIME(fine)
WRITE(*,*) "Tempo trascorso: ", fine-inizio
STOP
END PROGRAM prova_ordinamento
si può subito notare che esso accede agli elementi dei array per colonna, ossia secondo il meccanismo
ottimale che rispetta l’ordinamento in memoria degli elementi stessi. Compilato con
l’Essential Lahey Fortran 90 Compiler (Rel. 4.00b) e mandato in esecuzione su un PC con
processore Intel Pentium III, l’output prodotto è:
Tempo trascorso: 3.20461
Se, poi, si sostituisce il nido di DO con una più compatta ed elegante assegnazione in forma
globale:
array3 = array1+array2
il risultato è molto prossimo al precedente:
Tempo trascorso: 3.10446
Ben altro risultato, invece, si ottiene quando si ritorna al doppio ciclo e si invertono gli indici
i e j:
DO i=1,5000
DO j=1,5000
array3(i,j) = array1(i,j)+array2(i,j)
END DO
END DO
In questo caso, infatti, l’output è:
Tempo trascorso: 24.6555