Il Linguaggio Fortran 90/95

342 Array Processing
0 8 0
0 0 0
7 0 0
allora l’istruzione:
PACK(mat,MASK=mat/=0,VECTOR=(/1,3,5,9,11,13/))
produrrà il risultato: (/7,8,5,9,11,13/), mentre l’istruzione:
PACK(mat,MASK=mat/=0)
fornirà come risultato il vettore (/7,8/).
Un più interessante impiego della funzione in oggetto è offerto dalla seguente subroutine il
cui compito è quello di estrarre gli eleemnti della diagonale principale da una matrice quadrata
di input:
SUBROUTINE Matrix_Diagonal(A,diag,n)
! Scopo: estrarre la diagonale principale
! da una matrice quadrata A di n elementi
IMPLICIT NONE
REAL, INTENT(IN), DIMENSION(:,:) :: A
REAL, INTENT(OUT), DIMENSION(:) :: diag
INTEGER, INTENT(IN) :: n
REAL, DIMENSION (1:SIZE(A,1)*SIZE(A,1)) :: temp
!
IF(SIZE(A,1)==N .AND. SIZE(A,2)==N) THEN
temp = PACK(A,.TRUE.)
diag = temp(1:n*n:n+1)
ELSE
PRINT*, "La matrice A non e’ quadrata"
PRINT*, "Impossibile estrarre la diagonale"
END IF
RETURN
END SUBROUTINE Matrix_Diagonal
SPREAD(SOURCE,DIM,NCOPIES)
Crea un array ”replicato” con una dimensione aggiuntiva eseguendo copie degli elementi esistenti
lungo una dimensione specificata.
SOURCE è uno scalare o un array di tipo qualsiasi (nel caso di array, il rango deve essere
minore di 7).
DIM è uno scalare di tipo INTEGER con valore compreso tra 1 ed n +1, essendo n il rango di
SOURCE.