Il Linguaggio Fortran 90/95

5.13 Effetti collaterali nelle funzioni 275
PURE FUNCTION lunghezza(x,y)
IMPLICIT NONE
REAL, INTENT(IN) :: x, y
REAL :: lunghezza
lunghezza = SQRT(x**2+y**2)
END FUNCTION lunghezza
Analogamente, le pure subroutine si definiscono aggiungendo il prefisso PURE all’istruzione di
dichiarazione di subroutine. I loro vincoli sono esattamente uguali a quelli delle pure function
tranne per il fatto che esse possono modificare i valori degli argomenti dichiarati con l’attributo
INTENT(OUT) o INTENT(INOUT).
Si noti che le funzioni PURE sono le uniche funzioni, a parte quelle intrinseche, che possono
essere invocate all’interno di un costrutto FORALL. Così, ad esempio, il seguente costrutto:
FORALL(i=1:n,j=1:m)
x(i,j) = my_func(y(i,j))
END FORALL
è valido soltanto se la funzione my_func viene definita con l’attributo PURE, come ad esempio:
PURE FUNCTION my_func(var1) RESULT(var2)
REAL(KIND(1.D0)), INTENT(IN) :: var1
REAL(KIND(1.D0)) :: var2
...
END FUNCTION my_func
In tutti gli altri casi l’istruzione di chiamata di funzione produce un errore di compilazione.
Una particolare classe di procedure PURE, sono le elemental procedure, le quali si definiscono
sostituendo al prefisso PURE il prefisso ELEMENTAL nell’istruzione di dichiarazione di procedura.
Le elemental function non sono altro che pure function aventi come parametri formali soltanto
variabili scalari (quindi non sono ammessi come parametri formali né puntatori né procedure)
e che forniscono un risultato scalare. Si noti che il prefisso RECURSIVE è incompatibile con il
prefisso ELEMENTAL.
Le elemental subroutine sono definite allo stesso modo. L’unica differenza rispetto alle corrispondenti
function è il fatto che è permesso modificare i valori dei parametri formali specificati
con gli attributi INTENT(OUT) o INTENT(INOUT).
Il vantaggio sostanziale offerto dalle procedure ELEMENTAL rispetto alle PURE consiste in una
più efficace parallelizzazione delle operazioni. A parte le differenze testé citate, le elemental
procedure si comportano esattamente come le comuni procedure, come dimostrato dal seguente
esempio in cui si confrontano i valori della funzione esponenziale e x ottenuti con la funzione
intrinseca EXP e con i primi venti termini dello sviluppo in serie di Taylor.
PROGRAM prova_user_exp
!
IMPLICIT NONE