?gemm_alloc

Syntax

dest = sgemm_alloc (identifier, m, n, k)

dest = dgemm_alloc (identifier, m, n, k)

Include Files

• mkl.fi

Description

The ?gemm_alloc routine is one of a set of related routines that enable use of an internal packed storage. Call the ?gemm_alloc routine first to allocate storage for a packed matrix structure to be used in subsequent calls, ultimately to compute

C := alpha*op(A)*op(B) + beta*C,

where:

• op(X) is one of the operations op(X) = X, op(X) = X^T, or op(X) = X^H,
• alpha and beta are scalars,
• A , B, and C are matrices:
• op(A) is an m-by-k matrix,
• op(B) is a k-by-n matrix,
• C is an m-by-n matrix.

Input Parameters

identifier

CHARACTER*1. Specifies which matrix is to be packed:

If identifier = 'A' or 'a', the routine allocates storage to pack matrix A.

If identifier = 'B' or 'b', the routine allocates storage to pack matrix B.

m

INTEGER. Specifies the number of rows of matrix op(A) and of the matrix C. The value of m must be at least zero.

n

INTEGER. Specifies the number of columns of matrix op(B) and the number of columns of matrix C. The value of n must be at least zero.

k

INTEGER. Specifies the number of columns of matrix op(A) and the number of rows of matrix op(B). The value of k must be at least zero.

Output Parameters