※ 引述《Yagyu (近在眼前)》之铭言：
: ※ 引述《DMFC (sole)》之铭言：
: : 贴一部份程式码，fortran90语法
: do i=1,nx ; do j=1,ny ; do k=1,nz
: csum=sum(coeff(1:np)*cdexp)/Vol
: csumout(i,j,k)=csum
: end do ; end do ; end do
: sum(coeff(1:np)*cdexp)
: 请问这边是要将整个 coeff(1:23000) 乘上 cdexp 再做 sum 吗?
: 会这么问是因为我不确定你未贴上程式码部分是否还有 np 的 loop
: 如果没有 np loop, 同时 cdexp 只是个常数,不会随 nx,ny,nz 变动
: 那可以试着将这部份移出循环外, 这样省得每次 loop 都要重算一遍
: 如果有 np loop, 那请无视我的想法
: 另外输出的部份
: do i=1,nx
: do j=1,ny
: do k=1,nz
: csumout(i,j,k)=sum(coeff(1:np)*cdexp)/Vol
: end do
: end do
: end do
: do k=1,nz
: do j=1,ny
: do i=1,nx
: csumout(i,j,k)=sum(coeff(1:np)*cdexp)/Vol
: end do
: end do
: end do
: 两者差异 请参考彭国伦先生的fortran工具书 应该是在阵列章节中的多维阵列那边
: 没记错的话 是跟内存存放资料方式有关 这边变动我想绝对有帮助
谢谢
没错~顺序影响很大
不过这是因为我手误
由于我没法直接贴SOURCE CODE
所以是用手KEYIN
LOOP顺序我贴错了
且~很无奈的
我那个 cdexp 是与 i,j,k 有关
cdexp 非 常数~无法提出
我再贴一次完整的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)
end do ; end do ; end do
ngrid(1) ngrid(2) ngrid(3) 都是常数
atmp, wkiG1, wkiG2, wkiG3 都是维度23000(nplane)的大矩阵
Vol, pi2 是常数
cdexp 是 fortran 内有默认的function
意思是对 double precision 的 complex 取 exp
(atmp 是double precision 的 complex)
: : 请问一下
: : 由于资料量不算小
: : 目前测试 nx=ny=nz=47, np=23000
: : 上面贴的部份是需要最久计算的部份
: : 普通执行需要约两分钟(测试用参数)
: : 若是把第三行去掉
: : 则可以10秒内结束
: : 可是我就是要把这个资料输出
: : 才用了一个三维 ARRAY 把资料存起来
: : 然后离开这个 LOOP 后再写出
: : 如果把 WRITE 写在 LOOP 内
: : 则情况会更糟
: : 请问有什么方法可以处理这窘境

cdexp不是标准的内建函数，是你用的编译器有额外支援。

喔喔喔~~了解

同样的运算在回圈内做了很多次阿wkiG1*xyz(1)如果循环执行中没有变动 应该包在第一层

拉丁兄~你说的没错~但那不是重点呀包在那一层对整体效率的帮助不大因为我目前的问题是卡在写入 csumout忽然有个问题gilocustom 请问为什么他不是标准的内建函数呢我用 gfortran ifort 都有 cdexp或是这样问~那里可以查到"标准内建函数"?因为我又做了很多测试~我觉得问题是出在 cdexp请问~如果不这样写~要怎么写才能处理具有 double precision 的 complex 的 exponential 呢?

听起来有点矛盾，所以如果不要把第六行的值给到阵列里直接给到一个变量，e.g. csum，就会变快囉？

google "fortran standard document"就有很多标准书可查阅。从欧拉公式着手，用dexp与cmplx就可以兜出cdexp。

有一个可能，如果你没有开csumout，compiler说不定直接跳过所有计算过程 因为你没有输出 所以把他最佳化掉了!所以产生了速度提高的错觉.....没有看到完整的code很拿抓药 只能猜测如果不介意 可以私信讨论 以前做最佳化很有心得

coeff & atmp both are arry....cdexp & sum are a function..你在函数理面再作函数运算....要丢哪里暂存得看COMPILAR而定会慢会慢在这儿吧.....~"~a.....

Cypresslin ~ 是的~所以我理解不能......