应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » Fortran » 【求助】Fortran中有什么好方法计算CPU时间？【已完结】
91/1返回列表
查看: 796  |  回复: 8
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前主题已经存档。

dellus

木虫 (职业作家)

[交流] 【求助】Fortran中有什么好方法计算CPU时间？【已完结】

Fortran中SECOND子程序(or function)是怎么用的？
大家有什么好的办法计算CPU时间？
我试过secnds和CPU_Time，好象都有BUG

比如把一个运算放在一个do循环里运行，分别在循环前和循环后记录时间，最后写出时间
如果最后的write语句中加上循环中运行的结果，需要7、8秒，如果只输出时间，只有10（－5）秒

[ Last edited by 余泽成 on 2010-1-27 at 22:20 ]
回复此楼

» 猜你喜欢
岂能尽如人意，但求无愧于心。
1楼 2009-12-17 17:21:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

nono2009

超级版主 (文学泰斗)

No gains, no pains.

优秀区长优秀区长优秀区长优秀区长优秀版主
★ ★
dellus(金币+2,VIP+0):版主帮忙啊 12-17 17:56
以前自己在77下写过一个子程序。不知新版的fortran是否提供了更好的函数？跟着学
赞一下(6人) 回复此楼
2楼2009-12-17 17:25:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

nono2009

超级版主 (文学泰斗)

No gains, no pains.

优秀区长优秀区长优秀区长优秀区长优秀版主
★ ★ ★
波不动(金币+1,VIP+0):辛苦了。 12-17 19:30
dellus(金币+2,VIP+0):谢谢!F95里有DATE_AND_TIME,包括date和time的作用,而且不存在date的２０００年的问题,不过这个函数计时精度不够 12-18 06:14
我看了一下Compaq Visual Fortran 6.6 Language Reference，里面有date和time两个内部函数，应该可以解决你的问题。
赞一下(12人) 回复此楼
3楼2009-12-17 18:14:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

nono2009

超级版主 (文学泰斗)

No gains, no pains.

优秀区长优秀区长优秀区长优秀区长优秀版主

从我N年以前的程序中摘了点给你看，也不知行不行

★ ★ ★ ★ ★ ★ ★ ★ ★
波不动(金币+2,VIP+0):感谢应助！ 12-17 19:30
dellus(金币+3,VIP+0):谢谢!很热心的版主 12-18 06:16
dellus(金币+4,VIP+0): 12-23 14:05
...
      INTEGER(4) int_time0,int_time1
      character*8 char_time

!  show the starting time of this calculation
        call TIME(char_time)
        print *, 'Starting time on the PC is ', char_time
...
        int_time0=TIME( )
...
        int_time1=TIME( )
        print *, 'This calculation lasts ',int_time1-int_time0,'(s)'
赞一下(13人) 回复此楼
4楼2009-12-17 18:19:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

holmescn

金虫 (正式写手)
★ ★ ★ ★ ★ ★ ★ ★
jjdg(金币+2,VIP+0):感谢应助！欢迎新虫！ 12-17 18:59
dellus(金币+3,VIP+0):谢谢!试过了,也有BUG 12-18 06:18
dellus(金币+1,VIP+0):3 12-23 14:06
dellus(金币+2,VIP+0): 12-23 14:07
http://www.programfan.com/club/showtxt.asp?id=305205

从上面这个地方看到一种方法，用的是system_clock函数

call system_clock(t1)
call system_clock(t2)
t=t2-t1
单位是微秒。
赞一下(17人) 回复此楼
5楼2009-12-17 18:41:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

holmescn

金虫 (正式写手)
★ ★ ★ ★ ★
波不动(金币+2,VIP+0):感谢应助！ 12-17 19:30
dellus(金币+1,VIP+0):谢谢! 12-18 06:19
dellus(金币+2,VIP+0): 12-23 14:07
话说，如果是*unix系统，直接用time命令就好了。我在win下实现了一个，凑合能用。
赞一下(11人) 回复此楼
6楼2009-12-17 18:44:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dellus

木虫 (职业作家)
senlia(金币+0,VIP+0):我算出来都是5.156左右 每次运行时间都不完全一样 但都差不过十几毫秒 12-18 08:43
两段程序几乎一样，除了输出语句略有差别，但t１的运行结果一个是1.999999767576810E-006，一个是7.15873999999985。有谁知道怎么回事吗？还是只有我用的系统有这样的问题？

program main
      integer :: i
      double precision :: f0, f1, t0, t1
      f0 = 0
      f1 = 0
      t0 = secnds (0.0d0)
      do i = 1, 100000000
            f1 = f1 + dcos(f1)
      end do
      t1 = secnds(t0)
      write(*,*) t1
end program main


program main
      integer :: i
      double precision :: f0, f1, t0, t1
      f0 = 0
      f1 = 0
      t0 = secnds (0.0d0)
      do i = 1, 100000000
            f1 = f1 + dcos(f1)
      end do
      t1 = secnds(t0)
      write(*,*) t1, f1
end program main
赞一下(6人) 回复此楼
岂能尽如人意，但求无愧于心。
7楼2009-12-18 06:27:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

holmescn

金虫 (正式写手)
★ ★ ★ ★ ★ ★ ★ ★ ★ ★
senlia(金币+2,VIP+0):感谢应助 12-18 15:38
dellus(金币+5,VIP+0):谢谢！我用的是Linux，gcc。我试过real的精度，结果和双精度是一样的。我同一个程序不同次运行之间的差距也不止百分之一秒。 12-18 22:02
dellus(金币+3,VIP+0): 12-23 14:06
首先，不知道你用的什么编译器，什么系统。
我在winxp下，用gfortran 4.4.2编译了这两段程序。有如下问题

1、seconds的参数必须是kind 4的。而你的是double，也就是kind 8的。这个我改了一下，然后编译通过。
2、我的运行时间是：13.218750和13.171875，不知道为什么第二个会快百分之一秒。

另弱的问一下nono老大，write和print在输出到console的时候有什么区别吗？
赞一下(11人) 回复此楼
8楼2009-12-18 15:27:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

15890734984

木虫 (小有名气)
★ ★ ★ ★
jjdg(金币+2,VIP+0):辛苦了！感谢参与！ 12-20 02:33
dellus(金币+2,VIP+0):这个程序是做什么的？ 12-20 15:01
real K(1:8922,1:8922),a(1:89103,1:3)
  open(unit=1,file='filestiffness.txt',status='old')
  open(unit=2,file='kk.txt',status='new')
        read(1,*)((a(i,j),j=1,3),i=1,89103)
        do 10 i=1,8922
            do 20 j=1,8922
               do 30 n=1,89103
                   if (a(n,1).eq.i.and.a(n,2).eq.j) then
                             K(i,j)=a(n,3)
       end if
30      continue
20      continue
10      continue
!100     format(1x,8922(e19.12,1x))
  write(2,*) ((K(i,j),j=1,8922),i=1,8922)
        end


这个程序的功能是将数组a的第一列和第二列作为行号和列号给第三列在k中定位.
filestiffness.txt内容如下:
10   9   2558
6    4    5285
3    3  5582
…………
……………
赞一下(13人) 回复此楼
9楼2009-12-19 22:13:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 dellus 的主题更新
91/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
信息提示
关闭
请填处理意见
关闭