作者PTT0207445 (我是乡民)
看板Fortran
标题[心得] 利用内建副程式move_alloc调整可变阵列
时间Sat Mar 24 23:35:23 2018
假设你有一程式不断产生数据，你要将数据放在一个阵列里，但是却不能事先知道
一共有几个数据，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
以上简单分享，希望对大家有帮助。
--
※ 发信站: 批踢踢实业坊(ptt.cc), 来自: 129.12.197.57
※ 文章网址: https://webptt.com/cn.aspx?n=bbs/Fortran/M.1521905725.A.691.html
1F：推 alen84204: 推个 03/25 22:30
2F：推 espresso1: 推 03/30 21:07
3F：推 roythepigmy: 谢谢分享 04/09 22:25
4F：推 jubilee2: push 05/05 05:13
5F：推 fragmentwing: 推个 我终於看明白这个的用处了 04/18 11:20
6F：→ fragmentwing: 上上次程式作业别人写了50行而已 虽说有自找麻烦的 04/18 11:20
7F：→ fragmentwing: 成份在 不过为了做出以不停成长的矩阵来记录数据的 04/18 11:20
8F：→ fragmentwing: 方法 我总共打了300多行 04/18 11:20