Il Linguaggio Fortran 90/95

7.7 Puntatori a sezioni di array 377
INTEGER, DIMENSION(16), TARGET :: vet=(/(i, i=1,16)/)
INTEGER, DIMENSION(16), POINTER :: ptr1, ptr2, ptr3, ptr4, ptr5
ptr1 => vet
ptr2 => ptr1(2::2)
ptr3 => ptr2(2::2)
ptr4 => ptr3(2::2)
ptr5 => ptr4(2::2)
WRITE(*,’(A,16I3)’) " ptr1 = ", ptr1
WRITE(*,’(A,16I3)’) " ptr2 = ", ptr2
WRITE(*,’(A,16I3)’) " ptr3 = ", ptr3
WRITE(*,’(A,16I3)’) " ptr4 = ", ptr4
WRITE(*,’(A,16I3)’) " ptr5 = ", ptr5
END PROGRAM point_sez_array
Il risultato fornito da questo programma è, chiaramente, il seguente:
ptr1 = 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16
ptr2 = 2 4 6 8 10 12 14 16
ptr3 = 4 8 12 16
ptr4 = 8 16
ptr5 = 16
Si noti, però, che sebbene i puntatori lavorino con sezioni di array definite da triplette di
indici, essi non funzionano con sezioni di array definiti da indici vettoriali. Ad esempio, il
programma seguente fornirà un errore in fase di compilazione:
PROGRAM bad_sez_array
IMPLICIT NONE
INTEGER :: i
INTEGER, DIMENSION(3) :: indvet=(1,3,11)
INTEGER, DIMENSION(16), TARGET :: vet=(/(i, i=1,16)/)
INTEGER, DIMENSION(:), POINTER :: punt
punt => vet(indvet) ! errore: sezione di array con indice vettoriale
WRITE(*,’(1X,A,16I3)’) ’punt = ’, punt
END PROGRAM bad_sez_array
Come ulteriore esempio di puntatori a sezioni di array si potrebbe pensare di riscrivere il
programma di risoluzione dell’equazione di Laplace già incontrato al capitolo 3. Ancora una
volta si perviene ad un codice eccezionalmente compatto ed efficiente, magari appena un pò
meno ”immediato” da interpretare:
PROGRAM laplace_con_punt
! *** Sezione dichiarativa ***
IMPLICIT NONE
REAL,DIMENSION(10,10),TARGET :: Told