dotpr.f

来自「网络带宽测试工具」· F 代码 · 共 20 行

      Function dotpr( n, x, y )                           Result( sum )
      Use        numerics
      Use        floptime
      Implicit   None

      Integer  :: n
      Real(l_) :: x(n), y(n)
 
      Integer  :: i
      Real(l_) :: sum
! ---------------------------------------------------------------------
      sum = 0.0_l_
!$omp parallel do reduction(+:sum)
      Do i = 1, n
         sum = sum + x(i)*y(i)
      End Do
      flops = flops + 2*n                    ! --- Keep track of flops.
! ---------------------------------------------------------------------
      End Function dotpr