Il Linguaggio Fortran 90/95

150 Array
stesso spazio di iterazione ma ciascuno racchiudente soltanto un sottoinsieme delle operazioni
del ciclo originale. Tale operazione può avere diversi effetti benefici sul codice:
• Può dare luogo a cicli perfettamente innestati.
• Può creare sub-loop privi (o con ridotto numero) di dipendenze fra le relative istruzioni.
• Consentendo, ad esempio, di operare con array più piccoli che possono essere ospitati per
intero nella cache può migliorare l’utilizzo della cache memory.
• Diminuendo la pressione di registrazione può aumentare la frequenza di riutilizzo degli
stessi registri di memoria.
L’esempio che segue mostra come il loop fission possa rimuovere le dipendenze e consentire a
una parte di ciclo di essere eseguita in parallelo:
! ciclo originale
DO i=1,n
a(i) = a(i)*c
x(i+1) = x(i)*7+x(i+1)+a(i)
END DO
! dopo il loop fission
DO i=1,n
a(i) = a(i)*c
END DO
DO i=1,n
x(i+1) = x(i)*7+x(i+1)+a(i)
END DO
Un ulteriore esempio, leggermente più complesso, può servire a chiarire ulteriormente il concetto.
Si consideri il seguente ciclo:
! ciclo originale
DO i=1,n
a(i) = a(i)+b(i-1) ! S1
b(i) = c(i-1)*x + y ! S2
c(i) = 1./b(i) ! S3
d(i) = SQRT(c(i)) ! S4
END DO
Per poter comprendere le dipendenze fra le varie istruzioni sarà bene riportare in successione
l’effetto di ciascuna istruzione ad ogni iterazione:
• S2 assegna un valore a b(i) che sarà usato da S1 all’iterazione successiva.
• S2 assegna un valore a b(i) che sarà usato da S3 all’iterazione corrente.