gendat.f

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

      Subroutine gendat( a, b, lda, m, n )
! ----------------------------------------------------------------------
! --- Routine 'gendat' generates matrix 'a' and vector 'b' for the
!     matrix-vector multiplication 'Ab = c'.
! ----------------------------------------------------------------------
      Use         numerics
      Implicit    None
      Integer  :: lda, m, n
      Real(l_) :: a(lda,n), b(n)

      Integer  :: i
! ----------------------------------------------------------------------
      Do i = 1, m
         a(i,1:n) = i
      End Do
      Do i = m+1, lda
         a(i,1:n) = 0.0_l_
      End Do
      b = 1.0_l_
! ----------------------------------------------------------------------
      End Subroutine gendat