假设你有一程式不断产生数据,你要将数据放在一个阵列里,但是却不能事先知道
一共有几个数据,Fortran 77的做法是宣告一个足够大的阵列,执行的时候万一
数据的个数超过阵列的长度,则程式会出现内存错误而终止程式。
Fortran 90以后可以宣告一个可变大小的阵列,例如:
integer(4) :: n
real(8),dimension(:),allocatable :: a
此时阵列a并不占用内存空间,使用阵列a之前必须先把阵列的长度确定好,
利用allocate定义长度,例如:
n=1024
allocate(a(n))
此时程式在内存配置了一段空间给a阵列使用,a阵列的长度是1024。
由于一个倍精确度浮点数占用8个byte,所以a阵列使用了8KB的记忆空间。
等到阵列a不需要再使用的时候,可以利用deallocate把空间释放出来,像这样:
deallocate(a)
但是如果资料的长度超过allocate时的长度,事情变得有点麻烦。
必须先宣告另一个可变大小的阵列暂时存放a的所有值,然后
反复利用allocate和deallocate是可以做到不断将a的长度增长又保留资料,
但是此时可以利用内建副程式move_alloc简化一些程式码,像这样:
program demo
integer(4) :: n,m
real(8),dimension(:),allocatable :: a,temp
n=1024
allocate(a(n)) !此时a在内存中占用8KB。
a(:)=1.d0 !计算过程中a阵列记录了很多浮点数。
m=n+1024 !发现其实有更多浮点数必须记录。
allocate(temp(m)) !定义temp阵列的长度,此时temp阵列占用16KB的内存空间。
temp(1:n)=a !把阵列a暂存在temp,
temp(n+1:m)=2.d0 !然后把更多的数据紧接着也存在temp。
call move_alloc(temp,a) !重新命名原本的阵列temp现在叫做a,此时a的长度是2048,
print*,allocated(temp) !而temp阵列的状态变为deallocated。
end program demo
以上简单分享,希望对大家有帮助。