Il Linguaggio Fortran 90/95

5.17 Overloading 293
DO k = n,1,-1
X(k,:)%elem = X(k,:)%elem / W(k,k)
DO i=1,k-1
X(i,:)%elem = X(i,:)%elem - W(i,k) * X(k,:)%elem
END DO
END DO
END FUNCTION matrix_div
FUNCTION matrix_vector_div(Y,Z) RESULT(X)
! "Divide" un array 1D per un array 2D di tipo "matrix"
TYPE(matrix), INTENT(IN), DIMENSION(:) :: Y
TYPE(matrix), INTENT(IN), DIMENSION(:,:) :: Z
TYPE(matrix), DIMENSION(SIZE(Y,1)) :: X
REAL, DIMENSION(SIZE(Z,1),SIZE(Z,2)) :: W
INTEGER :: i, j, k, n
! effettua una copia di riserva degli argomenti
W(:,:) = Z(:,:)%elem
X(:)%elem = Y(:)%elem
! eliminazione di Gauss sulla matrice argomento (W|X)
n = SIZE(Z,2)
DO k = 1,n-1
DO i=k+1,n
W(i,k) = W(i,k)/W(k,k)
X(i)%elem = X(i)%elem - W(i,k) * X(k)%elem
END DO
DO j=k+1,n
DO i=k+1,n
W(i,j) = W(i,j) - W(i,k) * W(k,j)
END DO
END DO
END DO
! "back substitution" su X
DO k = n,1,-1
X(k)%elem = X(k)%elem / W(k,k)
DO i=1,k-1
X(i)%elem = X(i)%elem - W(i,k) * X(k)%elem
END DO
END DO
END FUNCTION matrix_vector_div
END MODULE operatori
Il programma che segue mostra un semplice utilizzo di questo modulo e delle operazioni da esso