Il Linguaggio Fortran 90/95

7.8 Puntatori come parametri formali di procedure 381
Il seguente programma mostra un ulteriore esempio di passaggio di argomenti con attributo
POINTER, anche questa volta ad una procedura interna. Lo stesso programma fornisce, nel
contempo, un nuovo esempio di uso efficiente di puntatori nelle operazioni di swap, molto
frequenti negli algoritmi iterativi. Il programma calcola la radice quadrata degli elementi di un
vettore y di 1000 elementi utilizzando la nota formula iterativa:
xn+1 = xn + y/xn
2
ed utilizza efficacemente i puntatori allo scopo di evitare onerose operazioni di copia. La convergenza
del metodo viene testata esaminando la differenza fra due approssimazioni consecutive
con riferimento al massimo discostamento fra i due passi.
PROGRAM prova_func_punt
IMPLICIT NONE
REAL, PARAMETER :: eps=0.001
INTEGER :: i, j
INTEGER, PARAMETER :: dim=1000
REAL, DIMENSION(dim) :: y=(/(i, i=1,dim)/)
REAL, DIMENSION(dim), target :: guess=1., xold=1., xnew
REAL, DIMENSION(:), POINTER :: ptnew, ptold, swap
LOGICAL :: conv=.FALSE.
ptnew => xnew
ptold => xold
ptold => guess
j = 0
DO WHILE(.NOT.conv)
j = j+1
ptnew = calc(ptold,y)
swap => ptold
ptold => ptnew
ptnew => swap
conv = MAXVAL(ABS(ptnew-ptold))