|
30000*1000个浮点如果是双精度的话,时间是绝对不够的.
time = 3.150000 0.0000000E+00
3464101623798.05
俺的机器还不是太挫.
ifort -O3别的暂时没开.
program test
implicit none
real,dimension(2) :: t1,t2
real :: t0,etime
integer:: i
real*8 :: x
t0=ETIME(t1)
i=1
x=0.0_8
do while(i<=300000*1000)
x=dsqrt(dfloat(i))+x
i=i+1
end do
t0=ETIME(t2)
print*,"time =",t2(1)-t1(1),t2(2)-t1(2)
print*,x
end program test
原帖由 zxjike 于 2007-7-2 23:30 发表 ![]()
我没有跑,但是我觉得应该是全部。 |
|