[心得] 利用內建副程式move_alloc調整可變陣列

[PTT0207445 (我是鄉民)]

假設你有一程式不斷產生數據，你要將數據放在一個陣列裡，但是卻不能事先知道
一共有幾個數據，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://www.ptt.cc/bbs/Fortran/M.1521905725.A.691.html
推 alen84204: 推個
推 espresso1: 推
推 roythepigmy: 謝謝分享
推 jubilee2: push
推 fragmentwing: 推個 我終於看明白這個的用處了
→ fragmentwing: 上上次程式作業別人寫了50行而已 雖說有自找麻煩的
→ fragmentwing: 成份在 不過為了做出以不停成長的矩陣來記錄數據的
→ fragmentwing: 方法 我總共打了300多行