Я пытался использовать мкл 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
что не так с моим кодом?
Я был неправ. Порядок аргументов m,n,k
как указано выше
Других решений пока нет …