Re: [問題] 效率

[janhon (海牛)]

: 我再貼一次完整的CODE
:  do iz=0,ngrid(3)-1 ; do iy=0,ngrid(2)-1 ; do ix=0,ngrid(1)-1
:   xyz(1) = dble(ix)/dble(ngrid(1))
:   xyz(2) = dble(iy)/dble(ngrid(2))
:   xyz(3) = dble(iz)/dble(ngrid(3))
: atmp = pi2 * (wkiG1*xyz(1) + wkiG2*xyz(2) + wkiG3*xyz(3))
     ^^(?)這到底是函數還是變數?

: csumout(ix,iy,iz) = sum(coeff(1:nplane)*cdexp(atmp(1:nplane)))/dsqrt(Vol)
                               (:)                  ^^


--
※ 發信站: 批踢踢實業坊(ptt.cc) 
◆ From: 218.161.81.224
→ janhon:open(1) 又close(1) 又open(1) 數字還是用100以上不要重複.
→ Cypresslin:#1IWGajWZ 有你要的答案
→ janhon:有看到 但是atmp=沒指定位置,幹嘛全陣列都指定成同數值