Il Linguaggio Fortran 90/95

294 Unità di Programma
definite. In particolare, definite due matrici 3×3 A e B ed un vettore c lo scopo del programma
è quello di eseguire le operazioni:
A×B
A×c
A −1
A/c = A −1 ×c
PROGRAM prova_operatori
!
USE operatori
IMPLICIT NONE
TYPE(matrix),DIMENSION(3,3) :: A, B, Unity
TYPE(matrix),DIMENSION(3) :: c
TYPE(matrix),DIMENSION(3,3) :: AperB, Inv_A
TYPE(matrix),DIMENSION(3) :: Aperc, cdivA
INTEGER :: i, j, dim
REAL :: h
! Inizializzazione delle matrici
DO i=1,3
DO j=1,3
Unity%elem=0.0
CALL RANDOM_NUMBER(HARVEST=h)
A(i,j)%elem = h+0.5
CALL RANDOM_NUMBER(HARVEST=h)
B(i,j)%elem = h-0.5
END DO
CALL RANDOM_NUMBER(HARVEST=h)
c(i)%elem = h
END DO
FORALL (i=1:3)
A(i,i)%elem = 2.0*SUM(ABS(A(i,:)%elem)) ! assicura la dominanza
! diagonale forte
Unity(i,i)%elem = 1.0 ! Unity e’ la matrice identica 3x3
END FORALL
AperB = A*B ! ris: A*B righe x colonne
Aperc = A*c ! ris. A*c righe x colonne
Inv_A = Unity/A ! ris: INV(A)*Unity = INV(A)
cdivA = c/A ! ris: INV(A)*c
dim=10
PRINT 100, "Matrice A:"
DO i=1,3