[問題] 主程式 call subroutine 會爆掉的原因

看板Fortran作者 (單手挑藏獒)時間8年前 (2015/11/23 18:08), 8年前編輯推噓1(1015)
留言16則, 5人參與, 最新討論串1/1
編譯後錯誤訊息: 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://www.ptt.cc/bbs/Fortran/M.1448273296.A.D92.html

11/23 23:23, , 1F
我不知道為何compile 出錯 但是就算成功兩個也是不一樣的
11/23 23:23, 1F

11/23 23:24, , 2F
啊 懂了 主程式裡宣告了A的大小傳給sub
11/23 23:24, 2F

11/23 23:25, , 3F
但是sub裡又重新宣告了他的大小 而且大小不一樣
11/23 23:25, 3F

11/23 23:25, , 4F
所以你要傳遞部份是可以的 但是需要另外寫一個叫B
11/23 23:25, 4F

11/23 23:26, , 5F
B的大小要跟sub裡宣告的一樣
11/23 23:26, 5F

11/23 23:26, , 6F
如果你是不想寫死array的大小 那就要用allocatable
11/23 23:26, 6F

11/23 23:26, , 7F
但是用allocatable的話 就只能用module傳遞
11/23 23:26, 7F

11/23 23:27, , 8F
事實上這就是為何你去看一堆老程式 老是宣告一個很大的
11/23 23:27, 8F

11/23 23:28, , 9F
array 然後主副程式都用共一個parameter來宣告大小
11/23 23:28, 9F
主程式中 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

11/24 01:26, , 10F
我設A(11,11,11)這樣的矩陣大小 兩種方法都可以執行
11/24 01:26, 10F

11/24 01:26, , 11F
如果設成512 會顯示記憶體不足 所以應該是記憶體不足所造
11/24 01:26, 11F

11/24 01:26, , 12F
成的問題
11/24 01:26, 12F

11/24 17:16, , 13F
sub裡的A陣列記憶體位置應該跟主程式的一樣,所以就……
11/24 17:16, 13F

11/26 21:01, , 14F
副程式的矩陣宣告用A(:,:,:),自動設定維度大小
11/26 21:01, 14F

12/12 22:41, , 15F
8*512**3使用記憶體約1G,程式2的方法副程式又再吃1G,
12/12 22:41, 15F

12/12 22:41, , 16F
你系統如果是32位元就會不夠
12/12 22:41, 16F
文章代碼(AID): #1MKkMGsI (Fortran)
文章代碼(AID): #1MKkMGsI (Fortran)