[数学库] fortran代码中dgemm函数调用的疑问

        program TestMKLDgemm
        implicit none
        real*8  :: alpha, beta
        integer :: i, j
        integer, parameter :: m = 5, k = 5, n = 5
        real*8 :: A(m,k), B(k,n), C(m,n)
        !// C = alpha*A*B + beta*C
        alpha = 1.d0
        beta  = 0.d0

        A=5.d0

        B=2.d0

        C = 1.d0

        call DGEMM('N', 'N', M, N, K, alpha, A, M, B, K, beta, C, M)
        print*, B
        PRINT*
        print*, A
        PRINT*
         print*, C
         READ*
        end program TestMKLDgemm