Question

我正在尝试使用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必须是＆gt; = MAX（K，1）：lda = 2 K = 3BLAS错误：传递给cblas_sgemm的参数号9具有无效值

Answer 1

你看到的错误对我来说似乎完全正确。

LDA始终是线性存储器中阵列A的间距。如果您使用行主存储顺序，则音高将是列数，而不是行数。所以在这种情况下，LDA应该是3。

vecLib cblas_sgemm文档错了吗？

1 个答案: