Modul BLAS

!
! simple version of dgemm
!
module blas
  use kinds
  implicit none

contains

  subroutine gemm(a, b, c)
    real(kind=REAL8), dimension(:,:), intent(in)    :: a, b
    real(kind=REAL8), dimension(:,:), intent(inout) :: c
    
    integer          :: n, m, k
    
    n = size(a,1)
    k = size(a,2)
    m = size(b,2)

    call dgemm('n', 'n', n, k, m, 1.0D0, a, n, b, k, 0.0D0, c, n)
    
  end subroutine gemm

end module blas