Commit 80e0818c authored by PRUVOST Florent's avatar PRUVOST Florent

use lapack_const(Magma...) conversion instead of hardcoded parameters like 'N', 'R' etc

parent e0fc3996
...@@ -305,7 +305,7 @@ magma_zparfb_gpu(magma_side_t side, magma_trans_t trans, ...@@ -305,7 +305,7 @@ magma_zparfb_gpu(magma_side_t side, magma_trans_t trans,
transW = storev == MorseColumnwise ? MagmaConjTrans : MagmaNoTrans; transW = storev == MorseColumnwise ? MagmaConjTrans : MagmaNoTrans;
transA2 = storev == MorseColumnwise ? MagmaNoTrans : MagmaConjTrans; transA2 = storev == MorseColumnwise ? MagmaNoTrans : MagmaConjTrans;
cublasZgemm(lapack_const(transW), 'N', cublasZgemm(lapack_const(transW), lapack_const(MagmaNoTrans),
K, N1, M2, K, N1, M2,
zone, zone,
(cuDoubleComplex*)V /* K*M2 */, LDV, (cuDoubleComplex*)V /* K*M2 */, LDV,
...@@ -316,8 +316,8 @@ magma_zparfb_gpu(magma_side_t side, magma_trans_t trans, ...@@ -316,8 +316,8 @@ magma_zparfb_gpu(magma_side_t side, magma_trans_t trans,
WORKC = NULL; WORKC = NULL;
if (WORKC == NULL) { if (WORKC == NULL) {
/* W = op(T) * W */ /* W = op(T) * W */
cublasZtrmm( 'L', 'U', cublasZtrmm( lapack_const(MagmaLeft), lapack_const(MagmaUpper),
lapack_const(trans), 'N', lapack_const(trans), lapack_const(MagmaNonUnit),
K, N2, K, N2,
zone, zone,
(cuDoubleComplex*)T, LDT, (cuDoubleComplex*)T, LDT,
...@@ -332,7 +332,7 @@ magma_zparfb_gpu(magma_side_t side, magma_trans_t trans, ...@@ -332,7 +332,7 @@ magma_zparfb_gpu(magma_side_t side, magma_trans_t trans,
} }
/* A2 = A2 - op(V) * W */ /* A2 = A2 - op(V) * W */
cublasZgemm(lapack_const(transA2), 'N', cublasZgemm(lapack_const(transA2), lapack_const(MagmaNoTrans),
M2, N2, K, M2, N2, K,
mzone, mzone,
(cuDoubleComplex*)V /* M2*K */, LDV, (cuDoubleComplex*)V /* M2*K */, LDV,
...@@ -349,14 +349,14 @@ magma_zparfb_gpu(magma_side_t side, magma_trans_t trans, ...@@ -349,14 +349,14 @@ magma_zparfb_gpu(magma_side_t side, magma_trans_t trans,
zzero, WORKC, LDWORKC ); zzero, WORKC, LDWORKC );
/* A1 = A1 - opt(T) * W */ /* A1 = A1 - opt(T) * W */
cublasZgemm( lapack_const(trans), 'N', cublasZgemm( lapack_const(trans), lapack_const(MagmaNoTrans),
K, N1, K, K, N1, K,
mzone, T, LDT, mzone, T, LDT,
WORK, LDWORK, WORK, LDWORK,
zone, A1, LDA1 ); zone, A1, LDA1 );
/* A2 = A2 - Wc * W */ /* A2 = A2 - Wc * W */
cublasZgemm( 'N', 'N', cublasZgemm( lapack_const(MagmaNoTrans), lapack_const(MagmaNoTrans),
M2, N2, K, M2, N2, K,
mzone, WORKC, LDWORKC, mzone, WORKC, LDWORKC,
WORK, LDWORK, WORK, LDWORK,
...@@ -386,7 +386,7 @@ magma_zparfb_gpu(magma_side_t side, magma_trans_t trans, ...@@ -386,7 +386,7 @@ magma_zparfb_gpu(magma_side_t side, magma_trans_t trans,
transW = storev == MorseColumnwise ? MagmaNoTrans : MagmaConjTrans; transW = storev == MorseColumnwise ? MagmaNoTrans : MagmaConjTrans;
transA2 = storev == MorseColumnwise ? MagmaConjTrans : MagmaNoTrans; transA2 = storev == MorseColumnwise ? MagmaConjTrans : MagmaNoTrans;
cublasZgemm('N', lapack_const(transW), cublasZgemm(lapack_const(MagmaNoTrans), lapack_const(transW),
M1, K, N2, M1, K, N2,
zone, zone,
(cuDoubleComplex*)A2 /* M1*N2 */, LDA2, (cuDoubleComplex*)A2 /* M1*N2 */, LDA2,
...@@ -397,8 +397,8 @@ magma_zparfb_gpu(magma_side_t side, magma_trans_t trans, ...@@ -397,8 +397,8 @@ magma_zparfb_gpu(magma_side_t side, magma_trans_t trans,
WORKC = NULL; WORKC = NULL;
if (WORKC == NULL) { if (WORKC == NULL) {
/* W = W * op(T) */ /* W = W * op(T) */
cublasZtrmm( 'R', 'U', cublasZtrmm( lapack_const(MagmaRight), lapack_const(MagmaUpper),
lapack_const(trans), 'N', lapack_const(trans), lapack_const(MagmaNonUnit),
M2, K, M2, K,
zone, zone,
(cuDoubleComplex*)T, LDT, (cuDoubleComplex*)T, LDT,
...@@ -413,7 +413,7 @@ magma_zparfb_gpu(magma_side_t side, magma_trans_t trans, ...@@ -413,7 +413,7 @@ magma_zparfb_gpu(magma_side_t side, magma_trans_t trans,
} }
/* A2 = A2 - W * op(V) */ /* A2 = A2 - W * op(V) */
cublasZgemm('N', lapack_const(transA2), cublasZgemm(lapack_const(MagmaNoTrans), lapack_const(transA2),
M2, N2, K, M2, N2, K,
mzone, mzone,
(cuDoubleComplex*)WORK /* M2*K */, LDWORK, (cuDoubleComplex*)WORK /* M2*K */, LDWORK,
...@@ -423,7 +423,7 @@ magma_zparfb_gpu(magma_side_t side, magma_trans_t trans, ...@@ -423,7 +423,7 @@ magma_zparfb_gpu(magma_side_t side, magma_trans_t trans,
} else { } else {
/* A1 = A1 - W * opt(T) */ /* A1 = A1 - W * opt(T) */
cublasZgemm( 'N', lapack_const(trans), cublasZgemm( lapack_const(MagmaNoTrans), lapack_const(trans),
M1, K, K, M1, K, K,
mzone, WORK, LDWORK, mzone, WORK, LDWORK,
T, LDT, T, LDT,
...@@ -437,7 +437,7 @@ magma_zparfb_gpu(magma_side_t side, magma_trans_t trans, ...@@ -437,7 +437,7 @@ magma_zparfb_gpu(magma_side_t side, magma_trans_t trans,
zzero, WORKC, LDWORKC ); zzero, WORKC, LDWORKC );
/* A2 = A2 - W * Wc */ /* A2 = A2 - W * Wc */
cublasZgemm( 'N', 'N', cublasZgemm( lapack_const(MagmaNoTrans), lapack_const(MagmaNoTrans),
M2, N2, K, M2, N2, K,
mzone, WORK, LDWORK, mzone, WORK, LDWORK,
WORKC, LDWORKC, WORKC, LDWORKC,
......
Markdown is supported
0%
or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment