Il Linguaggio Fortran 90/95

3.13 Costrutto e istruzione WHERE 165
3.13 Costrutto e istruzione WHERE
Come si è potuto apprendere nelle pagine precedenti, il Fortran 90/95 consente di utilizzare,
nelle istruzioni di assegnazione, singoli elementi di un array, una sua sezione o l’intero array.
Ad esempio, per estrarre il logaritmo degli elementi dell’array bidimensionale matrice di m×n
elementi, si può utilizzare la semplice istruzione:
log_matrice = LOG(matrice)
ben più compatta del doppio ciclo DO innestato:
DO i=1,m
DO j=1,n
log_matrice(i,j) = LOG(matrice(i,j))
END DO
END DO
Tuttavia, supponendo di voler estrarre il logaritmo soltanto di alcuni elementi dell’arraymatrice,
ad esempio (come appare oltremodo chiaro) dei soli elementi positivi, un modo per svolgere questo
compito è quello di applicare l’operatore logaritmo ad un elemento per volta, combinando
i cicli DO con il costrutto IF, al modo seguente:
DO i=1,m
DO j=1,n
IF(matrice(i,j)>0.) THEN
log_matrice(i,j) = LOG(matrice(i,j))
ELSE
log_matrice(i,j) = -HUGE(a)
END IF
END DO
END DO
dove, evidentemente, lo scopo della funzione intrinseca HUGE è quello di produrre un numero
reale di valore assoluto molto alto (per l’esattezza, il valore più grande rappresentabile in
macchina) al fine di simulare numericamente il concetto di ”infinito”. Lo stesso risultato può,
tuttavia, essere ottenuto in maniera più efficace con una forma speciale di istruzione di assegnazione:
la cosiddetta assegnazione con schema o maschera (masked array assignment). Si tratta
di un’istruzione la cui operazione è controllata da un array logico della stessa forma dell’array
specificato nell’assegnazione. L’operazione di assegnazione è svolta soltanto per gli elementi
dell’array che corrispondono ai valori dello schema di assegnazione.
Nel Fortran 90/95 questo tipo di assegnazione è implementato mediante il costrutto WHERE,
la cui forma generale è la seguente:
[nome:] WHERE (espressione_schema)
istruzioni_di_assegnazione_di_array ! primo blocco
ELSEWHERE [nome]
istruzioni_di_assegnazione_di_array ! secondo blocco
END WHERE [nome]