Il Linguaggio Fortran 90/95

3.4 Operazioni globali su array 133
• Soluzione elemento×elemento:
REAL a(20), b(20), c(20)
...
DO i=1,20
a(i) = 0.
END DO
...
DO i=1,20
a(i) = a(i)/3.1 + b(i)*SQRT(c(i))
END DO
• Soluzione globale:
REAL, DIMENSION(20) :: a, b, c
...
a = 0.
...
a = a/3.1 + b*SQRT(c)
Si noti che, nell’esempio in esame, la funzione intrinseca SQRT opera su ciascun elemento
dell’array c.
L’utilità di lavorare con array in forma globale è resa tanto più evidente quanto maggiore
è il rango e meno ”concordi” sono i limiti delle estensioni degli array operandi. Ad esempio,
dichiarati tre array di rango quattro al modo seguente:
REAL, DIMENSION(10,10,21,21) :: x
REAL, DIMENSION(0:9,0:9,-10:10,-10:10) :: y
REAL, DIMENSION(11:20,-9:0,0:20,-20:0) :: z
la semplice istruzione:
x = y+z
è perfettamente equivalente al seguente insieme di cicli DO innestati:
DO i=1,10
DO j=1,10
DO k=1,21
DO l=1,21
x(i,j,k,l) = y(i-1,j-1,k-11,l-11) + z(i+10,j-10,k-1,l-21)
END DO
END DO
END DO
END DO