<< linalg_dsyev blaslapack linalg_gesv >>

linalg >> linalg > blaslapack > linalg_gemm

linalg_gemm

Computes D = Alpha*A*B+ Beta*C for a real or complex matrix.

Calling Sequence

D = linalg_gemm ( Alpha , A , B , Beta , C )

Parameters

Alpha :

a 1-by-1 matrix of doubles (real or complex)

A :

a m-by-n matrix of doubles (real or complex)

B :

a n-by-p matrix of doubles (real or complex)

Beta :

a 1-by-1 matrix of doubles (real or complex)

C :

a m-by-p matrix of doubles (real or complex)

D :

a m-by-p matrix of doubles (real or complex), D = Alpha*A*B+ Beta*C

Description

Computes D = Alpha*A*B+ Beta*C for a real or complex matrix. Uses BLAS/DGEMM if all inputs are real, uses ZGEMM if any input is complex. Convert to complex matrices when necessary to call ZGEMM.

Examples

// A real matrix
Alpha = 0.1234;
Beta = 5.678;
m = 3;
n = 4;
k = 2;
C = matrix(1:m*n,m,n);
A = matrix(1:m*k,m,k);
B = matrix(1:k*n,k,n);
D = linalg_gemm(Alpha, A, B, Beta, C)

// See with matrices of mixed types
Alpha = 1 + %i;
Beta = Alpha;
A = matrix(1:4,2,2);
B = A;
C = A;
D = linalg_gemm(Alpha, A, B, Beta, C)

Authors


Report an issue
<< linalg_dsyev blaslapack linalg_gesv >>