Wie führen Sie Vektor-Matrix-Multiplikation mit BLAS ?

BLAS definiert die GEMV (Matrix-Vektor-Multiplikation) level-2-Betrieb. Wie eine BLAS-Bibliothek durchführen Vektor-Matrix-Multiplikation ?

Es ist wahrscheinlich offensichtlich, aber ich sehe nicht, wie zu verwenden BLAS Vorgang für diese Multiplikation. Ich hätte erwartet, dass ein GEVM Betrieb.

  • transpose(x) . transpose(A) = A. x
InformationsquelleAutor Baptiste Wicht | 2015-04-06
  1. 7

    Die Matrix-Vektor-Multiplikation einer (M x N) - Matrix mit einer (N x 1) Vektor wird Ergebnis eine (M x 1) - Vektor. In kurzen a*A(MxN)*X(Nx1) + b*Y(Mx1) -> Y(Mx1). Natürlich können Sie INCX und INCY wenn der Vektor in einer matrix.

    Zur Festlegung einer Vektor-Matrix-Multiplikation Vektor umgesetzt werden sollten. also a*X(1xM)*A(MxN) + b*Y(1xN) -> Y(1xN). Im Grunde haben Sie nicht ein Vektor, sondern eine einzelne Zeile der matrix.

    Ab diesem Punkt gibt es zwei Möglichkeiten.

    Entweder level-3 "GEMM"

    ?gemm(transa, transb, m, n, k, alpha, a, lda, b, ldb, beta, c, ldc)

    mit

    ?gemm('N', 'N', 1, N, M, a, X, 1, A, M, b, Y, 1)

    Oder tun Sie etwas mehr Mathematik. Wenn man bedenkt, dass (X*A)^T = A^T * X^T die Zeile matrix X konvertiert zu Vektor X^T(MX1). Auch Y transpose Vektor Y^T(Nx1). Natürlich Speicher-Weise sowohl X und X^T sind-store mit der gleichen Art und Weise, der Reihe nach. Das heißt, Sie können wieder verwenden GEMV mit transpose matrix A

    ?gemv('T', M, N, a, A, M, X, 1, b, Y, 1)
    • Danke für die Erklärung und die genaue Beispiele, es funktioniert gut 🙂
    • Ich habe versucht, sowohl. Aber ?gemm etwa 2~3 mal schneller als ?gemv. Sollte nicht ?gemv immer schneller als ?gemm?
    InformationsquelleAutor ztik
Schreibe einen Kommentar