24小时热门版块排行榜    

查看: 797  |  回复: 8
当前主题已经存档。
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

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 ]
回复此楼

» 猜你喜欢

岂能尽如人意,但求无愧于心。
已阅   回复此楼   关注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
…………
……………
9楼2009-12-19 22:13:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 9 个回答

nono2009

超级版主 (文学泰斗)

No gains, no pains.

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

★ ★
dellus(金币+2,VIP+0):版主帮忙啊 12-17 17:56
以前自己在77下写过一个子程序。不知新版的fortran是否提供了更好的函数?跟着学
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的2000年的问题,不过这个函数计时精度不够 12-18 06:14
我看了一下Compaq Visual Fortran 6.6 Language Reference,里面有date和time两个内部函数,应该可以解决你的问题。
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)'
4楼2009-12-17 18:19:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见