Il Linguaggio Fortran 90/95

6.6 Procedure intrinseche per gli array 345
INTEGER,PARAMETER :: last_number = 100
INTEGER,DIMENSION(last_number) :: numbers
INTEGER :: i, number_of_primes, ac
! Start
! Inizializza l’array con i valori (/0, 2, 3, ..., last_number/)
! NB: c’e’ 0 al posto di 1 perché 1 rappresenta un caso particolare.
numbers = (/0,(ac,ac=2,last_number)/)
DO i=2,last_number
IF (numbers(i) /= 0) THEN ! if questo numero e’ primo...
numbers(2*i : last_number : i) = 0 ! ...se ne eliminano i multipli
END IF
END DO
! Conta i primi
number_of_primes = COUNT(numbers/=0)
! Raccoglie i numeri primi nelle prime posizioni dell’array
numbers(1:number_of_primes) = PACK(numbers, numbers /= 0)
! Li stampa
PRINT "(A,I3,A)", "Numeri primi fra 1 e ",last_number,": "
PRINT "(5I7)", numbers(1:number_of_primes)
END PROGRAM Crivello_di_Eratostene
Si riporta, per completezza, anche l’output del programma:
Numeri primi fra 1 e 100:
2 3 5 7 11
13 17 19 23 29
31 37 41 43 47
53 59 61 67 71
73 79 83 89 97
6.6.4 Procedure di trasformazione
RESHAPE(SOURCE,SHAPE[,PAD][,ORDER])
Costruisce un array di forma differente a partire da un array di input.
L’argomento SOURCE è un array di tipo qualsiasi. Esso fornisce gli elementi per l’array
risultante. La sua ampiezza deve essere maggiore o uguale a PRODUCT(SHAPE) se PAD è omesso
oppure se ha ampiezza nulla.
L’argomento SHAPE è un array di sette elementi al massimo, rango unitario e dimensione
costante. Esso definisce la forma dell’array risultante. Non può avere ampiezza nulla ed i suoi
elementi non possono avere valori negativi.
L’argomento opzionalePAD è un array avente stesso tipo e stessi parametri di kind diSOURCE.
Il suo compito è quello di fornire valori ”di riserva” nel caso in cui l’array risultante avesse
ampiezza maggiore di SOURCE.