From 35d5105922445adeec359d42cf5972df88e213af Mon Sep 17 00:00:00 2001 From: Martin Kroeker Date: Sat, 26 Feb 2022 19:23:40 +0100 Subject: [PATCH 1/2] Enable xGEMMT functions --- relapack/config.h | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/relapack/config.h b/relapack/config.h index e4fab0a12..9d6919463 100644 --- a/relapack/config.h +++ b/relapack/config.h @@ -115,7 +115,7 @@ #define INCLUDE_CTGSYL INCLUDE_XTGSYL #define INCLUDE_ZTGSYL INCLUDE_XTGSYL -#define INCLUDE_XGEMMT 0 +#define INCLUDE_XGEMMT 1 #define INCLUDE_SGEMMT INCLUDE_XGEMMT #define INCLUDE_DGEMMT INCLUDE_XGEMMT #define INCLUDE_CGEMMT INCLUDE_XGEMMT From 4058f324923cd293b2117eafb65bc758a9a34a19 Mon Sep 17 00:00:00 2001 From: Martin Kroeker Date: Sat, 26 Feb 2022 19:24:27 +0100 Subject: [PATCH 2/2] Fix xGEMMT argument lists --- relapack/src/lapack_wrappers.c | 12 ++++++++---- 1 file changed, 8 insertions(+), 4 deletions(-) diff --git a/relapack/src/lapack_wrappers.c b/relapack/src/lapack_wrappers.c index 0252f3d92..fc3dbc11e 100644 --- a/relapack/src/lapack_wrappers.c +++ b/relapack/src/lapack_wrappers.c @@ -566,7 +566,8 @@ void LAPACK(sgemmt)( const float *B, const blasint *ldB, const float *beta, float *C, const blasint *ldC ) { - RELAPACK_sgemmt(uplo, n, A, ldA, info); + blasint info; + RELAPACK_sgemmt(uplo, transA, transB, n, k, alpha, A, ldA, B, ldB, beta, C, info); } #endif @@ -578,7 +579,8 @@ void LAPACK(dgemmt)( const double *B, const blasint *ldB, const double *beta, double *C, const blasint *ldC ) { - RELAPACK_dgemmt(uplo, n, A, ldA, info); + blasint info; + RELAPACK_dgemmt(uplo, transA, transB, n, k, alpha, A, ldA, B, ldB, beta, C, info); } #endif @@ -590,7 +592,8 @@ void LAPACK(cgemmt)( const float *B, const blasint *ldB, const float *beta, float *C, const blasint *ldC ) { - RELAPACK_cgemmt(uplo, n, A, ldA, info); + blasint info; + RELAPACK_cgemmt(uplo, transA, transB, n, k, alpha, A, ldA, B, ldB, beta, C, info); } #endif @@ -602,6 +605,7 @@ void LAPACK(zgemmt)( const double *B, const blasint *ldB, const double *beta, double *C, const blasint *ldC ) { - RELAPACK_zgemmt(uplo, n, A, ldA, info); + blasint info; + RELAPACK_zgemmt(uplo, transA, transB, n, k, alpha, A, ldA, B, ldB, beta, C, info); } #endif