dgemm in MKL

[meow]

I try the function dgemm in MKL :

call dgemm(AA1,BB1,CC1,transa,transb,1,0)

but I only get zeros in matrix CC1 no matter what goes into AA1 and BB1; matmul gives me the rigth answer; thanks!
Pablo

[Vladimir_Lunev]

Hi Pablo,

Fortran-77 interface for DGEMM is:
call dgemm(transa, transb, m, n, k, alpha, a, lda, b, ldb, beta, c, ldc)

And Fortran-95 interface is:
call gemm(a, b, c [,transa][,transb] [,alpha][,beta])

The default interface is Fortran 77.

Fortran-95 interface library is not included into the original distribution and should be built before using the interface. In the MKL User's Guide you can find how to build Fortran-95 interface library for BLAS if you prefer Fortran-95. Please, see the chapter "Using Language-Specific Interfaces with Intel MKL".

Vladimir

[kiran_nalawade]

Quoting - Vladimir Lunev (Intel)

Hi Pablo,

Fortran-77 interface for DGEMM is:
call dgemm(transa, transb, m, n, k, alpha, a, lda, b, ldb, beta, c, ldc)

And Fortran-95 interface is:
call gemm(a, b, c [,transa][,transb] [,alpha][,beta])

The default interface is Fortran 77.

Fortran-95 interface library is not included into the original distribution and should be built before using the interface. In the MKL User's Guide you can find how to build Fortran-95 interface library for BLAS if you prefer Fortran-95. Please, see the chapter "Using Language-Specific Interfaces with Intel MKL".

Vladimir

Anybody knows this for C version og dgmm() ? ( with all argument info)

[Andrey_Bespalov]

It is not quite clear to medo you mean CBLAS or BLAS but you can find the declarations of the functions in the header files.

void cblas_dgemm(const CBLAS_ORDER Order, const CBLAS_TRANSPOSE TransA,
const CBLAS_TRANSPOSE TransB, const MKL_INT M, const MKL_INT N,
const MKL_INT K, const double alpha, const double *A,
const MKL_INT lda, const double *B, const MKL_INT ldb,

void dgemm(const char *transa, const char *transb, const MKL_INT *m, const MKL_INT *n, const MKL_INT *k,
const double *alpha, const double *a, const MKL_INT *lda, const double *b, const MKL_INT *ldb,
const double *beta, double *c, const MKL_INT *ldc);