matgen.f

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

      Subroutine matgen( a, lda, n, b, norma )
! ----------------------------------------------------------------------
      Use         numerics
      Integer  :: n
      Real(l_) :: a(lda,n), b(n), norma
! ----------------------------------------------------------------------

      Call random_number( a )
      norma = Maxval( a )

      b = 0.0_l_
      Do j = 1, n
         Do i = 1, n
            b(i) = b(i) + a(i,j)
         End Do
      End Do
! ----------------------------------------------------------------------
      End Subroutine matgen