作者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/m.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