fortran代码中dgemm函数调用的疑问

宇都宫紫苑

代码如下：

[Fortran] 纯文本查看 复制代码

        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

生成解决方案没有问题，但是得出的结果如图，根本不对啊，请教是什么问题呀

宇都宫紫苑

添加的库是mkl_blas95_ilp64.lib， mkl_intel_lp64_dll.lib， mkl_intel_thread_dll.lib， mkl_core_dll.lib， libiomp5md.lib

宇都宫紫苑

fcode 发表于 2019-7-11 17:46
哪儿不对？
【2】【5】不就是 【50】 吗？2*5+2*5+2*5+2*5+2*5 = 50

对不起对不起，脑子懵了，还以为是10。。。

宇都宫紫苑

fcode 发表于 2019-7-11 18:23
建议你用95接口，调用更简单

use blas95

老师我想请教下，95接口的使用条件是什么呢？是跟代码的语法是否是for95有关系吗？