Il Linguaggio Fortran 90/95

370 Puntatori
END TYPE matrix
INTERFACE OPERATOR (+)
MODULE PROCEDURE sum
END INTERFACE
INTERFACE OPERATOR (*)
MODULE PROCEDURE prod, scale
END INTERFACE
PRIVATE sum, prod, scale
CONTAINS
FUNCTION sum(A,B)
TYPE(matrix), INTENT(IN) :: A, B
TYPE(matrix) :: sum
IF (A%row /= B%row .OR. A%col /= B%col) THEN
WRITE(*,*) " Impossibile effettuare la somma "
ELSE
sum%row = A%row
sum%col = A%col
ALLOCATE(sum%array(A%row, A%col))
sum%array = A%array + B%array
END IF
END FUNCTION sum
FUNCTION prod(A,B)
TYPE(matrix), INTENT(IN) :: A, B
TYPE(matrix) :: prod
IF (A%col /= B%row ) THEN
WRITE(*,*) " Impossibile effettuare il prodotto "
ELSE
prod%row = A%row
prod%col = B%col
ALLOCATE(prod%array(A%row,B%col))
prod%array = MATMUL(A%array, B%array)
END IF
END FUNCTION prod
FUNCTION scale(c,A)
REAL, INTENT(IN) :: c
TYPE(matrix), INTENT(IN) :: A