作者janhon (海牛)
看板Fortran
标题Re: [问题] 效率
时间Sat Nov 16 20:11:47 2013
: 我再贴一次完整的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
1F：→ janhon:open(1) 又close(1) 又open(1) 数字还是用100以上不要重复. 11/16 20:17
3F：→ janhon:有看到 但是atmp=没指定位置,干嘛全阵列都指定成同数值 11/26 12:17