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

宇都宫紫苑

代码如下：

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

查看代码

打印?

01	program TestMKLDgemm

02	implicit none

03	real*8  :: alpha, beta

04	integer :: i, j

05	integer, parameter :: m = 5, k = 5, n = 5

06	real*8 :: A(m,k), B(k,n), C(m,n)

07	!// C = alpha*A*B + beta*C

08	alpha = 1.d0

09	beta  = 0.d0

10

11

A=5.d0

12

13

B=2.d0

14

15

C = 1.d0

16

17	call DGEMM('N', 'N', M, N, K, alpha, A, M, B, K, beta, C, M)

18

print*, B

19

PRINT*

20

print*, A

21

PRINT*

22

print*, C

23

READ*

24	end program TestMKLDgemm

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

宇都宫紫苑

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

fcode

哪儿不对？
【2】【5】不就是 【50】 吗？2*5+2*5+2*5+2*5+2*5 = 50

宇都宫紫苑

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

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

fcode

建议你用95接口，调用更简单

use blas95
call GEMM(a, b, c, 'N','N', alpha, beta)

宇都宫紫苑

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

use blas95

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

fcode

条件就是你的编译器提供95接口，很显然ivf是满足这个条件的。

代码的语法是不是95没关系。而且，只要编译器支持，一段代码可以77,90,95,2003,2008,2015混合着用（这些规范之间本来也有很大的交集）。没有人能明确的指出一段代码到底是什么标准。

zxc123

fcode 发表于 2019-7-12 07:59
条件就是你的编译器提供95接口，很显然ivf是满足这个条件的。

代码的语法是不是95没关系。而且，只要编译 ...

大佬，在Linux下调用gemm函数，需要在makefile文件里怎么加啊，我用make进行编译的，MKL库我也下载了，但这个gemm函数调用不过来