Il Linguaggio Fortran 90/95

326 Array Processing
CONTAINS
SUBROUTINE sub(a,res)
IMPLICIT NONE
! Parametri formali
REAL, INTENT(OUT) :: res
REAL, DIMENSION(:,:), INTENT(IN) :: a ! assumed shape array
! Variabili locali
REAL, DIMENSION (SIZE(a,1),SIZE(a,2)) :: work ! automatic array
...
res = a(...)
...
END SUBROUTINE sub
END PROGRAM array
Un ulteriore esempio di uso combinato delle diverse tipologie di array è fornito dal seguente
programma che valuta il determinante di una matrice quadrata A. Detta matrice viene definita
come allocabile nel programma principale, viene quindi allocata in base alla dimensione n il cui
valore è di input. Una volta che la lettura della matrice abbia avuto termine, il main invoca
la funzione ricorsiva det che valuta il determinante secondo la regola di Laplace applicata alla
prima riga ossia come la somma dei prodotti degli elementi della suddetta riga per i rispettivi
complementi algebrici:
n
detA =
j=1
(−1) (1+j) a1j det A1j
dove, si ricorda, A1j è la sotto-matrice ottenuta da A cancellando la prima riga e la j_ma
colonna. Il determinante detA1j è detto minore complementare dell’elemento a1j, il prodotto
(−1) (1+j) detA1j è detto complemento algebrico di a1j.
Si noti come la matrice venga passata alla funzione come array correntemente allocato e
venga, invece, letta come un comune array di forma presunta. La stessa function, poi, crea ad
ogni invocazione ricorsiva un array automatico (di dimensioni sempre diverse) che rappresenta
una estratta di A il cui determinante è il minore relativo al generico elemento di A.
PROGRAM determinante
! Scopo: calcolare il determinante di una matrice quadrata.
! Gli elementi della matrice vengono letti dal file matrix.dat
IMPLICIT NONE
INTEGER :: n ! dimensione della matrice
REAL,DIMENSION(:,:),ALLOCATABLE :: A
INTEGER :: stato
!
INTERFACE
RECURSIVE FUNCTION det(n,A) RESULT(detA)
IMPLICIT NONE
REAL :: detA