[心得] 利用內建副程式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