作者sin55688 (單手挑藏獒)
看板Fortran
標題[問題] 主程式 call subroutine 會爆掉的原因
時間Mon Nov 23 18:08:13 2015
編譯後錯誤訊息:
forrtl: severe (157): Program Exception - access violation
除錯的結果，我認為跟記憶體使用限制有關，但不清楚原因，想請教板上的前輩
下面是我碰到的情形，我用簡單的語法表示
=======================================
程式1 可執行成功
=======================================
real(8) :: A(512,512,512)
A = 1d0
call sub(512,512,512,A)
A = A+A
=======================================
程式2 執行失敗
=======================================
real(8) :: A(512,512,512)
A = 1d0
call sub(511,511,511,A(1:511,1:511,1:511)
A = A+A <-------在這步會出現錯誤訊息
=======================================
subroutine sub(n1,n2,n3,A)
implict none
integer,intent(in) :: n1,n2,n3
real(8),intent(in) :: A(n1,n2,n3)
end subroutine
副程式本身並無任何作用
所以認為是呼叫副程式時，傳入資料導致主程式記憶體空間出現問題
但不了解傳入整個陣列A 與 部分陣列A(1:511,1:511,1:511)
因為我的理解，subroutine 是傳送記憶體位置，並非傳送值
為何後者會導致錯誤發生，希望有經驗的前輩可以指教，謝謝
--
※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 140.113.22.208
※ 文章網址: https://webptt.com/m.aspx?n=bbs/Fortran/M.1448273296.A.D92.html
1F：→ noonee: 我不知道為何compile 出錯 但是就算成功兩個也是不一樣的 11/23 23:23
2F：→ noonee: 啊 懂了 主程式裡宣告了A的大小傳給sub 11/23 23:24
3F：→ noonee: 但是sub裡又重新宣告了他的大小 而且大小不一樣 11/23 23:25
4F：→ noonee: 所以你要傳遞部份是可以的 但是需要另外寫一個叫B 11/23 23:25
5F：→ noonee: B的大小要跟sub裡宣告的一樣 11/23 23:26
6F：→ noonee: 如果你是不想寫死array的大小 那就要用allocatable 11/23 23:26
7F：→ noonee: 但是用allocatable的話 就只能用module傳遞 11/23 23:26
8F：→ noonee: 事實上這就是為何你去看一堆老程式 老是宣告一個很大的 11/23 23:27
9F：→ noonee: array 然後主副程式都用共一個parameter來宣告大小 11/23 23:28
主程式中 A 陣列其實是用 allocate 配置空間的
而 sub 的 A 陣列應該只是告知傳入陣列的大小，不具備配置記憶體功能
不知我的理解是否錯誤?
在我測試例子當中，A陣列非常的大，如果縮小成 256^3 則不會有問題
所以我在猜想是不是 A(1:511,1:511,1:511)傳入副程式時
這種部分陣列會而外再吃記憶體，導致記憶體不足?
※ 編輯: sin55688 (140.113.125.200), 11/24/2015 01:16:17
10F：→ rex0707: 我設A(11,11,11)這樣的矩陣大小 兩種方法都可以執行 11/24 01:26
11F：→ rex0707: 如果設成512 會顯示記憶體不足 所以應該是記憶體不足所造 11/24 01:26
12F：→ rex0707: 成的問題 11/24 01:26
13F：→ blc: sub裡的A陣列記憶體位置應該跟主程式的一樣，所以就…… 11/24 17:16
14F：推 hiyiyi: 副程式的矩陣宣告用A(:,:,:)，自動設定維度大小 11/26 21:01
15F：→ crazy10167: 8*512**3使用記憶體約1G,程式2的方法副程式又再吃1G, 12/12 22:41
16F：→ crazy10167: 你系統如果是32位元就會不夠 12/12 22:41