Mentions légales du service

Skip to content
Snippets Groups Projects

Fix for the TPMQRT routine

Merged Dalal Sukkari requested to merge sukkarde into master
1 file
+ 15
15
Compare changes
  • Side-by-side
  • Inline
@@ -47,7+47,7 @@
m1 = K;
n1 = N;
ldwork = IB;
ldworkc = M;
ws = IB * n1;
}
else {
@@ -61,7+61,7 @@
/* TS case */
if (L == 0) {
CUDA_ztsmqr( side, trans, m1, n1, M, N, K, IB,
A, LDA, B, LDB, V, LDV, T, LDT,
A, LDA, B, LDB, V, LDV, T, LDT,
WORK, ldwork, WORK + ws, ldworkc,
CUBLAS_STREAM_VALUE );
}
CUBLAS_STREAM_VALUE );
}
/* TT case */
else if( L == M ) {
else if( L == M ) {
// CUDA_zttmqr( side, trans, m1, n1, M, N, K, IB,
// A, LDA, B, LDB, V, LDV, T, LDT,
// WORK, ldwork );
CUDA_ztsmqr( side, trans, m1, n1, M, N, K, IB,
A, LDA, B, LDB, V, LDV, T, LDT,
A, LDA, B, LDB, V, LDV, T, LDT,
WORK, ldwork, WORK + ws, ldworkc,
CUBLAS_STREAM_VALUE );
morse_warning("The lower triangular part must be set to 0 for to be correct");
// CUDA_zttmqr( side, trans, m1, n1, M, N, K, IB,
// A, LDA, B, LDB, V, LDV, T, LDT,
// WORK, ldwork );
}
else {
cudablas_error(-6, "TPMQRT not available on GPU yet\n" );
return -6;
CUBLAS_STREAM_VALUE );
morse_warning("TPMQRT is not available on GPU yet: The lower triangular part must be set to 0 for to be correct");
}
else {
//LAPACKE_ztpmqrt_work( LAPACK_COL_MAJOR, M, N, K, L, IB, V, LDV, T, LDT, A, LDA, B, LDB, WORK );
}
cudablas_error(-6, "TPMQRT not available on GPU yet\n" );
return -6;
}
return MORSE_SUCCESS;
}
Loading