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

看板Fortran作者PTT0207445 (我是鄉民)時間6年前 (2018/03/24 23:35)推噓5(5推 0噓 3→)

留言8則, 5人參與討論串1/1

假設你有一程式不斷產生數據，你要將數據放在一個陣列裡，但是卻不能事先知道一共有幾個數據，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