ifort 有問題!?

[YCTzeng (我有一個大密寶)]

編譯器版本是
Intel(R) Fortran Intel(R) 64 Compiler XE for applications
running on Intel(R) 64, Version 12.1.4.319 Build 20120410
系統是 64位元 Linux fedora core 14
kernel版本 2.6.35.14-106.fc14.x86_64 #1 SMP
4 GB 的記憶體

以下是程式碼：
如果將最後一行註解 write(*,*) v5(100)
就可以編譯，且成功正確執行。
真的很奇怪!


program hmwf
implicit none
integer(4),parameter :: lmax=45000000
real(8),dimension(lmax) :: v1,v2,v3,v4,v5,v6
integer(4) :: max,i,j
real(8) :: su1,su2,su3

max=45000000
do j=1,100
  su1=0.d0
  su2=0.d0
  su3=0.d0
  do i=1,max
    v1(i)=1.d10
    v2(i)=v1(i)
    v3(i)=v1(i)
    v4(i)=v3(i)
    v5(i)=v3(i)
    v6(i)=v3(i)
    su1=su1+v1(i)**2
    su2=su2+v2(i)**2
    su3=su3+v3(i)**2
  enddo
  write(*,*) 'j=',j,v6(max),v4(1)
enddo
write(*,*) v5(100)

end program hmwf

--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 118.168.81.104
→ egg9:2G的array一次就要了這麼多個記憶體當然會不夠用囉
※ 編輯: YCTzeng         來自: 118.168.81.104       (05/23 12:19)
→ egg9:如果資料不是這麼有必要都保存下來就改個寫法吧
→ YCTzeng:v1,v2,v3,v4,v5,v6都只各佔約343MB,加起來2GB,我系統有4GB
→ gilocustom:Ver.12.1.3 build 20120212，可以編譯，不過記憶體不夠
→ gilocustom:所以沒測試可不可以正常運作。
→ YCTzeng:Intel 64,Ver:12.1.1.256 Build 20111011可以編譯正常執行
→ YCTzeng:所以好像真的是這個版本編譯器有問題,大家要小心!
→ YCTzeng:剛剛又測試了20120410的32位元版本,可編譯但不能正確執行,
→ YCTzeng:但是如果用 ifort -fast 就可編譯,正確執行。
推 terryys:有沒有試過用allocate來配置記憶體？可能會跟固定的不同
→ YCTzeng:推樓上,改成allocatable,可編譯,正確執行。真是奇怪。
推 perceval:http://tinyurl.com/87dlfbb large arrays(>2GB memory)
→ perceval:剛查到的解決方法，另外PGI也有同樣問題
→ perceval:不過gfortran 可以編譯，正確執行。
→ YCTzeng:感謝樓上，這樣就解決了。謝謝。
→ YCTzeng:補充： ifort -shared-intel -mcmodel=medium xxx.f90
推 callmei:推一下 原來LINUX下不是直接設STACK的大小@@