Il Linguaggio Fortran 90/95

5.17 Overloading 291
! Esegue la somma di due array 2D di tipo "matrix"
TYPE(matrix), INTENT(IN), DIMENSION(:,:) :: Y
TYPE(matrix), INTENT(IN), DIMENSION(SIZE(Y,1),SIZE(Y,2)) :: Z
TYPE(matrix), DIMENSION(SIZE(Y,1),SIZE(Y,2)) :: X
X(:,:)%elem = Y(:,:)%elem + Z(:,:)%elem
END FUNCTION matrix_add
FUNCTION vector_add(Y,Z) RESULT(X)
! Esegue la somma di due array 1D di tipo "matrix"
TYPE(matrix), INTENT(IN), DIMENSION(:) :: Y
TYPE(matrix), INTENT(IN), DIMENSION(SIZE(Y,1)) :: Z
TYPE(matrix), DIMENSION(SIZE(Y,1)) :: X
X(:)%elem = Y(:)%elem + Z(:)%elem
END FUNCTION vector_add
!
! Operatore "differenza": X = Y - Z
!
FUNCTION matrix_sub(Y,Z) RESULT(X)
! Esegue la differenza fra due array 2D di tipo "matrix"
TYPE(matrix), INTENT(IN), DIMENSION(:,:) :: Y
TYPE(matrix), INTENT(IN), DIMENSION(SIZE(Y,1),SIZE(Y,2)) :: Z
TYPE(matrix), DIMENSION(SIZE(Y,1),SIZE(Y,2)) :: X
X(:,:)%elem = Y(:,:)%elem - Z(:,:)%elem
END FUNCTION matrix_sub
FUNCTION vector_sub(Y,Z) RESULT(X)
! Esegue la differenza fra due array 1D di tipo "matrix"
TYPE(matrix), INTENT(IN), DIMENSION(:) :: Y
TYPE(matrix), INTENT(IN), DIMENSION(SIZE(Y,1)) :: Z
TYPE(matrix), DIMENSION(SIZE(Y,1)) :: X
X(:)%elem = Y(:)%elem - Z(:)%elem
END FUNCTION vector_sub
!
! Operatore "moltiplicazione": X = Y * Z
!
FUNCTION matrix_mul(Y,Z) RESULT(X)
! Esegue il prodotto righe x colonne fra due array 2D di tipo "matrix"
TYPE(matrix), INTENT(IN), DIMENSION(:,:) :: Y
TYPE(matrix), INTENT(IN), DIMENSION(:,:) :: Z
TYPE(matrix), DIMENSION(SIZE(Y,1),SIZE(Z,2)) :: X
X(:,:)%elem = MATMUL(Y(:,:)%elem, Z(:,:)%elem)
END FUNCTION matrix_mul