Il Linguaggio Fortran 90/95

7.4 Allocazione dinamica della memoria con i puntatori 371
TYPE(matrix) :: scale
scale%row = A%row
scale%col = A%col
ALLOCATE(scale%array(A%row,A%col))
scale%array = c*A%array
END FUNCTION scale
SUBROUTINE setm(row,column,array,A)
TYPE(matrix), INTENT(OUT) :: A
INTEGER, INTENT(IN) :: row, column
REAL, DIMENSION(:,:), INTENT(IN) :: array
A%row = row
A%col = column
ALLOCATE(A%array(A%row,A%col))
A%array = array
END SUBROUTINE setm
SUBROUTINE remove(A)
TYPE(matrix), INTENT(INOUT) :: A
DEALLOCATE(A%array)
END SUBROUTINE remove
SUBROUTINE printm(A)
TYPE(matrix), INTENT(IN) :: A
INTEGER :: i, j
DO i = 1,A%row
PRINT*, (A%array(i,j), j=1,A%col)
END DO
END SUBROUTINE printm
END MODULE matrici
Naturalmente è possibile estendere questo modulo al fine di prevedere altre operazioni come
il calcolo della trasposta, dell’inversa oppure del determinante. A tal fine si può fare uso, con
lievi modifiche, di alcune procedure già riportate nei capitoli precedenti. Un programma utile
a testare il modulo precedente potrebbe essere il seguente:
PROGRAM test
USE matrici
IMPLICIT NONE
TYPE(matrix) :: A, B, C
REAL, DIMENSION(2,4) :: tmp
tmp(1,:) = (/1.0, 2.0, 3.0, 4.0 /)
tmp(2,:) = (/5.0, 6.0, 7.0, 8.0 /)