[心得] 利用内建副程式move_alloc调整可变阵列

楼主: PTT0207445 (我是乡民)   2018-03-24 23:35:23
假设你有一程式不断产生数据,你要将数据放在一个阵列里,但是却不能事先知道
一共有几个数据,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
以上简单分享,希望对大家有帮助。
作者: alen84204 (Dana)   2018-03-25 22:30:00
推个
作者: espresso1   2018-03-30 21:07:00
作者: roythepigmy (金色三倍速)   2018-04-09 22:25:00
谢谢分享
作者: jubilee2 (3321)   2018-05-05 05:13:00
push
作者: fragmentwing (片翼碎梦)   2018-04-18 11:20:00
推个 我终于看明白这个的用处了上上次程式作业别人写了50行而已 虽说有自找麻烦的成份在 不过为了做出以不停成长的矩阵来记录数据的方法 我总共打了300多行

Links booklink

Contact Us: admin [ a t ] ucptt.com