求助！求问用oneMKL的BLAS库求矩阵乘法

Proteus

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

查看代码

打印?

01	program main

02	implicit none

03	integer :: n, m, k, lda, ldb, ldc, info

04	real(kind=8), dimension(:,:), allocatable :: A, B, C

05	real(kind=8) :: alpha, beta

06	external dgemm

07

! 设置矩阵的大小

08

m = 3

09

n = 3

10

k = 3

11	! 设置矩阵 A、B、C 的值

12	allocate(A(m,k))

13	allocate(B(k,n))

14	allocate(C(m,n))

15	A = reshape([1., 2., 3., 4., 5., 6., 7., 8., 9.], [m, k])

16	B = reshape([1., 2., 3., 4., 5., 6., 7., 8., 9.], [k, n])

17

C = 0.0

18	! 设置矩阵乘法的参数

19	alpha = 1.0

20	beta = 0.0

21

lda = m

22

ldb = k

23

ldc = m

24	! 调用 BLAS 库的 DGEMM 子程序进行矩阵乘法

25	call dgemm('N', 'N', m, n, k, alpha, A, lda, B, ldb, beta, C, ldc)

26

! 输出结果

27	write(*, *) 'Result:'

28	write(*, *) C

29	end program main

问题：提示“Intel MKL ERROR: Parameter 3 was incorrect on entry to DGEMM”，计算出C的结果全是0.0

编译器：Win10 + Visual Studio 2019 + oneAPI
已在项目属性中配置好oneMKL库了，之前运行过调用oneMKL库中的pardiso函数，是可以运行的，且计算结果是没问题的

不知道这个为什么出来的结果不对
求问大佬，求大佬相助！！！

Proteus

是我配置函数库的时候出了问题，已解决。呜呜呜呜呜。