vecLib cblas_sgemm 文档错误？

Posted 2023-02-25

【中文标题】vecLib cblas_sgemm 文档错误？

【英文标题】：vecLib cblas_sgemm documentation wrong?

【发布时间】：2012-03-20 18:01:03

【问题描述】：

我正在尝试使用 vecLibs 的 cblas 将两个矩阵相乘：

#include <stdlib.h>
#include <string.h>
#include <stdio.h>

#include <vecLib/cblas.h>

int main (void) 

    float *A = malloc(sizeof(float) *  2 * 3);
    float *B = malloc(sizeof(float) *  3 * 1);
    float *C = malloc(sizeof(float) *  2 * 1);

    cblas_sgemm(CblasRowMajor, 
                CblasNoTrans,
                CblasNoTrans, 
                2,
                1,
                3,
                1.0,
                A, 2,
                B, 3,
                0.0,
                C, 2);

    printf ("[ %f, %f]\n", C[0], C[1]);

    return 0;  

根据docs，每个参数似乎都匹配但我收到此错误：

lda 必须 >= MAX(K,1)：lda=2 K=3BLAS 错误：传递给 cblas_sgemm 的参数号 9 具有无效值

【参考方案1】：

您看到的错误在我看来完全正确。

LDA 总是线性内存中数组 A 的间距。如果您使用 row major 存储顺序，间距将是列数，而不是行数。所以在这种情况下，LDA 应该是 3。

【讨论】：

文档说 lda The size of the first dimention of matrix A; if you are passing a matrix A[m][n], the value should be m. 但这个 m 是矩阵行数。

文档说的是正确的 - lda 是矩阵 A 在线性内存中 的第一维。冒号后面的内容适用于列的主要顺序大小写。