Il Linguaggio Fortran 90/95

2.3 Istruzioni di ripetizione 117
IF (j500) EXIT outer
IF (sum>100) EXIT inner
END DO inner
sum = sum+i
i = i+10
END DO outer
PRINT*, ’sum = ’,sum
END PROGRAM sommatoria
2.3.4 Osservazioni sull’uso dei cicli
E’ errore assai comune (e spesso molto meno ovvio di quanto si crede) quello di eseguire uno
o più calcoli non ”indiciati” all’interno di un ciclo. Il seguente loop, ad esempio, è molto poco
efficiente in quanto esegue il medesimo calcolo (la divisione per 4 della variabile pi e la radice
quadrata del risultato) 10 volte:
DO i = 1,10
x(i) = x(i)*SQRT(pi/4.0) ! NB: x è un array: gli array saranno
END DO ! trattati nel prossimo capitolo
mentre sarebbe molto più efficiente calcolare il termine SQRT(pi/4.0) una sola volta all’esterno
del ciclo dal momento che esso non dipende dall’indice i del loop (si parla in tal caso di loop
invariant code motion). Così, ad esempio, è senz’altro preferibile sostituire il ciclo precedente
con questo:
factor = SQRT(pi/4.0)
DO i=1,10
x(i) = x(i)*factor
END DO
o, meglio ancora, con un’assegnazione in forma globale (quest’ultima soluzione, come si vedrà
nel capitolo 3, è senz’altro da preferirsi ogni qualvolta le circostanze lo permettano):
factor = SQRT(pi/4.0)
x = x*factor
C’è da dire che molti compilatori hanno la possibilità di eseguire questo tipo di ottimizzazione
sul codice sorgente ogni qualvolta la sostituzione appaia ”ovvia”. Esistono, però, dei casi
in cui l’interpretazione dell’invarianza di un fattore non è immediata, soprattutto in presenza
di funzioni definite dall’utente. Il prossimo ciclo ne è un esempio: