Fortran - Matrices dinámicas
UNA dynamic array es una matriz, cuyo tamaño no se conoce en el momento de la compilación, pero se conocerá en el momento de la ejecución.
Las matrices dinámicas se declaran con el atributo allocatable.
Por ejemplo,
real, dimension (:,:), allocatable :: darrayEl rango de la matriz, es decir, las dimensiones deben mencionarse; sin embargo, para asignar memoria a dicha matriz, use el allocate función.
allocate ( darray(s1,s2) )Después de usar la matriz, en el programa, la memoria creada debe liberarse usando el deallocate función
deallocate (darray)Ejemplo
El siguiente ejemplo demuestra los conceptos discutidos anteriormente.
program dynamic_array
implicit none
!rank is 2, but size not known
real, dimension (:,:), allocatable :: darray
integer :: s1, s2
integer :: i, j
print*, "Enter the size of the array:"
read*, s1, s2
! allocate memory
allocate ( darray(s1,s2) )
do i = 1, s1
do j = 1, s2
darray(i,j) = i*j
print*, "darray(",i,",",j,") = ", darray(i,j)
end do
end do
deallocate (darray)
end program dynamic_arrayCuando el código anterior se compila y ejecuta, produce el siguiente resultado:
Enter the size of the array: 3,4
darray( 1 , 1 ) = 1.00000000
darray( 1 , 2 ) = 2.00000000
darray( 1 , 3 ) = 3.00000000
darray( 1 , 4 ) = 4.00000000
darray( 2 , 1 ) = 2.00000000
darray( 2 , 2 ) = 4.00000000
darray( 2 , 3 ) = 6.00000000
darray( 2 , 4 ) = 8.00000000
darray( 3 , 1 ) = 3.00000000
darray( 3 , 2 ) = 6.00000000
darray( 3 , 3 ) = 9.00000000
darray( 3 , 4 ) = 12.0000000Declaración de uso de datos
los data La instrucción se puede utilizar para inicializar más de una matriz o para la inicialización de la sección de la matriz.
La sintaxis de la declaración de datos es:
data variable / list / ...Ejemplo
El siguiente ejemplo demuestra el concepto:
program dataStatement
implicit none
integer :: a(5), b(3,3), c(10),i, j
data a /7,8,9,10,11/
data b(1,:) /1,1,1/
data b(2,:)/2,2,2/
data b(3,:)/3,3,3/
data (c(i),i = 1,10,2) /4,5,6,7,8/
data (c(i),i = 2,10,2)/5*2/
Print *, 'The A array:'
do j = 1, 5
print*, a(j)
end do
Print *, 'The B array:'
do i = lbound(b,1), ubound(b,1)
write(*,*) (b(i,j), j = lbound(b,2), ubound(b,2))
end do
Print *, 'The C array:'
do j = 1, 10
print*, c(j)
end do
end program dataStatementCuando el código anterior se compila y ejecuta, produce el siguiente resultado:
The A array:
7
8
9
10
11
The B array:
1 1 1
2 2 2
3 3 3
The C array:
4
2
5
2
6
2
7
2
8
2Uso de la declaración Where
los wheredeclaración le permite usar algunos elementos de una matriz en una expresión, dependiendo del resultado de alguna condición lógica. Permite la ejecución de la expresión, en un elemento, si la condición dada es verdadera.
Ejemplo
El siguiente ejemplo demuestra el concepto:
program whereStatement
implicit none
integer :: a(3,5), i , j
do i = 1,3
do j = 1, 5
a(i,j) = j-i
end do
end do
Print *, 'The A array:'
do i = lbound(a,1), ubound(a,1)
write(*,*) (a(i,j), j = lbound(a,2), ubound(a,2))
end do
where( a<0 )
a = 1
elsewhere
a = 5
end where
Print *, 'The A array:'
do i = lbound(a,1), ubound(a,1)
write(*,*) (a(i,j), j = lbound(a,2), ubound(a,2))
end do
end program whereStatementCuando el código anterior se compila y ejecuta, produce el siguiente resultado:
The A array:
0 1 2 3 4
-1 0 1 2 3
-2 -1 0 1 2
The A array:
5 5 5 5 5
1 5 5 5 5
1 1 5 5 5