Il Linguaggio Fortran 90/95

3.8 Ottimizzazione delle operazioni con array 151
• S3 assegna un valore a c(i) che sarà usato dalla stessa S3 all’iterazione successiva.
• S3 assegna un valore a c(i) che sarà usato da S4 all’iterazione corrente.
Da quanto detto si può subito notare che il ciclo originale può essere ”splittato” nelle tre
componenti S1, (S2, S3) e S4. Naturalmente, poiché le operazioni eseguite da S1 e da S4 sono
subordinate a quelle eseguite dal blocco (S2, S3) è necessario che nella versione ”ottimizzata”
del ciclo il blocco (S2, S3) preceda S1 e S4.
! dopo il loop fission
DO i=1,n
b(i) = c(i-1)*x + y ! S2
c(i) = 1./b(i) ! S3
END DO
DO i=1,n
a(i) = a(i)+b(i-1) ! S1
END DO
DO i=1,n
d(i) = SQRT(c(i)) ! S4
END DO
3.8.8 Loop fusion
Se due cicli vicini operano sugli stessi dati, conviene di solito combinarli in un unico ciclo
operando in tal modo l’operazione nota come loop fusion. Questo procedimento, che rappresenta
un pò l’inverso del loop fission, non solo rende un codice più compatto e leggibile ma, al
tempo stesso, consente di utilizzare più efficientemente la cache memory limitando al massimo
il numero delle operazioni di lettura dei dati in memoria. Inoltre esso può aumentare la velocità
di esecuzione dei cicli in quanto, fondendo insieme diversi cicli controllati dal medesimo indice,
riduce il numero di condizioni logiche da verificare per la chiusura del loop. Un esempio di
applicazione dell’operazione di fusione è riportata di seguito:
! ciclo originale
DO i=1,n
a(i) = a(i) + k
END DO
DO i=1,n
b(i+1) = b(i) + a(i)
END DO
! dopo il loop fusion
DO i=1,n
a(i) = a(i) + k