Fortran Coder

查看: 1206|回复: 1
打印 上一主题 下一主题

[求助] fortran中CPU_TIME计时程序

[复制链接]

10

帖子

5

主题

0

精华

入门

F 币
60 元
贡献
30 点
跳转到指定楼层
楼主
发表于 2023-9-19 12:21:24 | 只看该作者 回帖奖励 |倒序浏览 |阅读模式
      
[Fortran] 纯文本查看 复制代码
program main
      implicit none
      real(kind=4)::t1_real, t2_real
      integer i, j ,k

      call CPU_TIME(t1_real)


      do i = 1,10000
          do j = 1,10000
              k = j
          end do
      end do


      call CPU_TIME(t2_real)

      write(*,*) "count1",t1_real
      write(*,*) "count2",t2_real
      write(*,*) "时间",t2_real - t1_real

      end program

以上程序的输出结果为: count1  0.0000000E+00
                                     count2  3.1250000E-02
                                     时间  3.1250000E-02

当我把以上程序写入abaqus子程序中:
      
[Fortran] 纯文本查看 复制代码
subroutine umat(stress,statev,ddsdde,sse,spd,scd,......)
...
...
...
...
real(kind=4) :: t1_real, t2_real
integer  i,j,k
...
...
...
...
      call CPU_TIME(t1_real)
      do i = 1,10000
          do j = 1,10000
              k = j
          end do
      end do
      call CPU_TIME(t2_real)
      write(*,*) "count1",t1_real                               !默认输出在log文件里
      write(*,*) "count2",t2_real                               !默认输出在log文件里
      write(*,*) "时间",t2_real - t1_real                     !默认输出在log文件里
...
...
...
...
end

在log文件中查看输出结果:
count1   6.250000   
count2   6.250000  
时间  0.0000000E+00
count1   6.250000   
count2   6.250000      
时间  0.0000000E+00
....
....
...

所有的 count1和count2都相等,时间为0.  是由于某些原因CPU_TIME函数在abaqus子程序中失效了吗?

分享到:  微信微信
收藏收藏 点赞点赞 点踩点踩

2033

帖子

12

主题

5

精华

论坛跑堂

臭石头雪球

F 币
1641 元
贡献
709 点

美女勋章热心勋章星光勋章新人勋章贡献勋章管理勋章帅哥勋章爱心勋章规矩勋章元老勋章水王勋章

沙发
发表于 2023-9-19 14:08:33 | 只看该作者
可能你的umat是用release编译的,开了优化。而主程序代码是debug编译的,没有开优化。

一旦你开启了优化,你的代码
      do i = 1,10000
          do j = 1,10000
              k = j
          end do
      end do
就可能自动被优化为 k=10000 这句代码。
您需要登录后才可以回帖 登录 | 极速注册

本版积分规则

捐赠本站|Archiver|关于我们 About Us|小黑屋|Fcode ( 京ICP备18005632-2号 )

GMT+8, 2024-12-23 22:47

Powered by Tencent X3.4

© 2013-2024 Tencent

快速回复 返回顶部 返回列表