mkl cblas_dgemm错误参数

Question

我尝试用mkl cblas_dgemm来计算矩阵乘法。

如我所知，lda、ldb、ldc应该是简单行主矩阵的列数。我试着做以下几件事：

double a[3 * 2] = { 1,2,3,4,5,6 };      // 3 x 2 matrix
double b[2 * 4] = { 1,2,3,4,5,6,7,8 };  // 2 x 4 matrix
double c[3 * 4] = { 0, };               // 3 x 4 matrix

// c <- 1.0 * (a*b) + 0.0 * c
cblas_dgemm(CblasRowMajor, CblasNoTrans, CblasNoTrans, 3, 2, 4, 1.0, a, 2, b, 4, 0.0, c, 4);
//            rowmajor       no trans      no trans    m  k  n  alph A lda B ldb beta C ldc

但是它给了我一条错误消息，c没有更新。

Intel MKL ERROR: Parameter 9 was incorrect on entry to cblas_dgemm.

参数9是lda

我的密码怎么了？

Answer 1

我错了。参数排序为m,n,k，如上文注释所示