4

vecLibs の cblas を使用して 2 つの行列を乗算しようとしています。

#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;  
}

ドキュメントによると、すべての引数が一致しているように見えますが、このエラーが発生します:

lda は >= MAX(K,1) でなければなりません: lda=2 K=3BLAS エラー: cblas_sgemm に渡されたパラメータ番号 9 に無効な値がありました

4

1 に答える 1

6

あなたが見ているエラーは私の目には完全に正しいようです。

LDAは、常に線形メモリ内の配列Aのピッチです。行の主保管順序を使用している場合、ピッチは行数ではなく列数になります。したがって、この場合、LDAは3である必要があります。

于 2012-03-20T18:31:00.133 に答える