mxm.f

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

      Subroutine mxm( a, b, c, m, l, n )
      Use          numerics
      Implicit     None

      Integer   :: l, m, n
      Real(l_)  :: a(m,l), b(l,n), c(m,n)

      Integer   :: i, j, lf, ll
! ---------------------------------------------------------------------
      c = 0.0_l_                                     ! --- Initialise c.   
      lf = l - Mod( l, 8 )
      Do j = 1, n
         Do ll = 1, lf, 8
            Do i = 1, m
               c(i,j) = c(i,j)
     &                   + a(i,ll)*b(ll,j)     + a(i,ll+1)*b(ll+1,j)
     &                   + a(i,ll+2)*b(ll+2,j) + a(i,ll+3)*b(ll+3,j)
     &                   + a(i,ll+4)*b(ll+4,j) + a(i,ll+5)*b(ll+5,j)
     &                   + a(i,ll+6)*b(ll+6,j) + a(i,ll+7)*b(ll+7,j)
            End Do
         End Do
         Do ll = lf + 1, l
            Do i = 1, m
               c(i,j) = c(i,j) + a(i,ll)*b(ll,j)
            End Do
         End Do
      End Do
! ---------------------------------------------------------------------
      End Subroutine mxm