<< Chapter < Page | Chapter >> Page > |
Tal como sucede con las variables locales de un programa en C, los arreglos automáticos de FORTRAN 90 sólo tienen asignado espacio de almacenamiento durante la vida de la subrutina o función que los contiene. Esto es diferente del almacenamiento local de arreglos tradicional de FORTRAN, donde algo de espacio se apartaba a tiempo de compilación o enlace. El tamaño y forma de los arreglos automáticos puede esculpirse mediante ua combinación de constantes y argumentos. Por ejemplo, he aquí una declaración de un arreglo automático,
B
, usando la nueva especificación de sintaxis de FORTRAN 90:
SUBROUTINE RELAX(N,A)
INTEGER NREAL, DIMENSION (N) :: A, B
Se declararon dos arreglos:
A
, el argumento de prueba, y
B
, un arreglo automático de tamaño explícito. Cuando se regresa de la subrutina,
B
deja de existir. Observe que el tamaño de
B
se toma de uno de los argumentos,
N
.
Los arreglos asignables le dan a usted la flexibilidad de elegir el tamaño de un arreglo después de examinar otras variables en el programa. Por ejemplo, puede que desee determinar la cantidad de datos de entrada antes de asignar los arreglos. Este pequeño programa pregunta al usuario el tamaño de la matriz antes de asignarle espacio de almacenamiento:
INTEGER M,N
REAL, ALLOCATABLE, DIMENSION (:,:) :: X...
WRITE (*,*) 'INTRODUZCA LAS DIMENSIONES DE X'READ (*,*) M,N
ALLOCATE (X(M,N))...
hacer algo con X...
DEALLOCATE (X)...
La sentencia
ALLOCATE
crea un arreglo de
M × N
que posteriormente se libera mediante la sentencia
DEALLOCATE
. Tal como sucede con los programas en
C
, es importante devolver la memoria asignada cuando se ha terminado de usar; de otra forma, su programa pudiera consumir toda la memoria virtual de almacenamiento disponible.
El problema de flujo calórico es un programa ideal para demostrar cuán agradablemente puede expresar FORTRAN 90 los problemas que emplean arreglos regulares:
PROGRAM HEATROD
PARAMETER(MAXTIME=200)INTEGER TICKS,I,MAXTIME
REAL*4 ROD(10)ROD(1) = 100.0
DO I=2,9ROD(I) = 0.0
ENDDOROD(10) = 0.0
DO TICKS=1,MAXTIMEIF ( MOD(TICKS,20) .EQ. 1 ) PRINT 100,TICKS,(ROD(I),I=1,10)
ROD(2:9) = (ROD(1:8) + ROD(3:10) ) / 2ENDDO
100 FORMAT(I4,10F7.2)END
El programa es idéntico, excepto por el hecho de que el bucle interno ha sido reemplazado por una sola sentencia que calcula la "nueva" sección, al promediar una tira de los elementos "izquierdos" y una tira de los elementos "derechos".
La salida de este programa luce así:
E6000: f90 heat90.f
E6000:a.out1 100.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00
21 100.00 82.38 66.34 50.30 38.18 26.06 18.20 10.35 5.18 0.0041 100.00 87.04 74.52 61.99 50.56 39.13 28.94 18.75 9.38 0.00
61 100.00 88.36 76.84 65.32 54.12 42.91 32.07 21.22 10.61 0.0081 100.00 88.74 77.51 66.28 55.14 44.00 32.97 21.93 10.97 0.00
101 100.00 88.84 77.70 66.55 55.44 44.32 33.23 22.14 11.07 0.00121 100.00 88.88 77.76 66.63 55.52 44.41 33.30 22.20 11.10 0.00
141 100.00 88.89 77.77 66.66 55.55 44.43 33.32 22.22 11.11 0.00161 100.00 88.89 77.78 66.66 55.55 44.44 33.33 22.22 11.11 0.00
181 100.00 88.89 77.78 66.67 55.55 44.44 33.33 22.22 11.11 0.00E6000:
Notification Switch
Would you like to follow the 'Cómputo de alto rendimiento' conversation and receive update notifications?