感謝大家對這問題有興趣
但我沒法貼全部程式碼
因為我是用 chrome 在上 ptt
所以我每次複製貼上
chrome 都傻傻的給我"一行"貼上
它不懂什麼叫換行
因此我去找了一下網路空間
並上傳SOURCE CODE
http://ge.tt/7QJGSAz/v/0
檔名是 trans-WAVECAR.f90
WAVECAR 是一個 binary file
所以這個程式目的是要把 WAVECAR 轉換成 ASCII 並轉換成別種格式
這樣我就可以用繪圖軟體把WAVECAR畫出來
我是在 linux server 上作業
用的是 intel compiler 12.0
指令是用 ifort -assume byterecl trans-WAVECAR.f90
我所貼上來的部份是在 line 185~196 間
我有試過將 csumout 拿掉
程式執行通常都是秒殺
加上去後就不得了的慢了
: : xyz(1) = dble(ix)/dble(ngrid(1))
: : xyz(2) = dble(iy)/dble(ngrid(2))
: : xyz(3) = dble(iz)/dble(ngrid(3))
: 這邊浪費很多加減乘除的處理時間...推文中講過...不再累述....
: 程式碼無全數貼出故無法得知中間少了啥....囧
: 理論上atmp & coeff這兩個array要先填完...接下去作才有意義...
: 建議你在atmp = ......這一行用Euler formular拆解re part and im part
: 並先算完cdexp後...
: 再到下一行直接用dot函數處理並填到csumout裡頭....
: 最後3層loop結束後再加一行 csumout=csumout/dsqrt(Vol)
Euler formular 是指什麼
是這個嬤??
http://zh.wikipedia.org/wiki/%E6%AC%A7%E6%8B%89%E5%85%AC%E5%BC%8F
令我難以接受的是
我試過
csum = sum(coeff(1:nplane)*cdexp(atmp(1:nplane)))/dsqrt(Vol) ! csum is complex
如果只是用一個參數去存取
程式效率很好
但只要再加入
csumout(ix,iy,iz)=csum
把 csum 存到三維矩陣
那效率就差很多很多很多 .....