您所在位置：主题：如何提高速度！急，如何改下面的代码使其计算更快.
burger
[专家分：0] 发布于 2011-06-22 22:45:00
x=1.0/M
do i=1,M
pxmx(i)=x
aymx(i)=pxmx(i)/LceLL*pi
axmx(i)=pi/LceLL
x=x+1.0/M
end do
x=0.0
do i=1,M+1
pxtx(i)=x
pxmt(i)=x
axtx(i)=pxtx(i)/LceLL*pi
aytx(i)=axtx(i)*0
axmt(i)=pi/LceLL-pxmt(i)/LceLL*pi
aymt(i)=axmt(i)
x=x+1.0/M
end do
do i=1,M+1
arphx(i)=axtx(i)
arphy(i)=aytx(i)
px(i)=pxtx(i)*pi
end do
do i=1,M
arphx(i+M+1)=axmx(i)
arphy(i+m+1)=aymx(i)
px(i+m+1)=pxmx(i)*pi+1*pi
end do
do i=1,m+1
arphx(i+2*M+1)=axmt(i)
arphy(i+2*M+1)=aymt(i)
px(i+2*M+1)=pxmt(i)*2**0.5*pi+2*pi
end do
do i=1,M
pxmx(i)=x
aymx(i)=pxmx(i)/LceLL*pi
axmx(i)=pi/LceLL
x=x+1.0/M
end do
x=0.0
do i=1,M+1
pxtx(i)=x
pxmt(i)=x
axtx(i)=pxtx(i)/LceLL*pi
aytx(i)=axtx(i)*0
axmt(i)=pi/LceLL-pxmt(i)/LceLL*pi
aymt(i)=axmt(i)
x=x+1.0/M
end do
do i=1,M+1
arphx(i)=axtx(i)
arphy(i)=aytx(i)
px(i)=pxtx(i)*pi
end do
do i=1,M
arphx(i+M+1)=axmx(i)
arphy(i+m+1)=aymx(i)
px(i+m+1)=pxmx(i)*pi+1*pi
end do
do i=1,m+1
arphx(i+2*M+1)=axmt(i)
arphy(i+2*M+1)=aymt(i)
px(i+2*M+1)=pxmt(i)*2**0.5*pi+2*pi
end do