floraachy/OpenBLAS
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

Fix cols/rows mixup in omatcopy 2nd step for BlasTrans cases

Browse Source 
Equivalent of #1244 (issue #899) for the non-complex cases. Fixes #1289
This commit is contained in:
Martin Kroeker
2017-09-14 19:59:05 +02:00
committed by GitHub
parent 00740c0e34
commit d674fbb4c7
1 changed files with 3 additions and 3 deletions
Download Patch File Download Diff File Expand all files Collapse all files

6
interface/imatcopy.c
View File

@@ -121,7 +121,7 @@ void CNAME( enum CBLAS_ORDER CORDER, enum CBLAS_TRANSPOSE CTRANS, blasint crows,
return; return;
} }
#ifdef NEW_IMATCOPY #ifdef NEW_IMATCOPY
if ( *lda == *ldb && *cols == *rows ) { if ( *lda == *ldb && *rows == *cols) {
if ( order == BlasColMajor ) if ( order == BlasColMajor )
{ {
if ( trans == BlasNoTrans ) if ( trans == BlasNoTrans )
@@ -171,7 +171,7 @@ void CNAME( enum CBLAS_ORDER CORDER, enum CBLAS_TRANSPOSE CTRANS, blasint crows,
else else
{ {
OMATCOPY_K_CT(*rows, *cols, *alpha, a, *lda, b, *ldb ); OMATCOPY_K_CT(*rows, *cols, *alpha, a, *lda, b, *ldb );
OMATCOPY_K_CN(*rows, *cols, (FLOAT) 1.0, b, *ldb, a, *ldb ); OMATCOPY_K_CN(*cols, *rows, (FLOAT) 1.0, b, *ldb, a, *ldb );
} }
} }
else else
@@ -184,7 +184,7 @@ void CNAME( enum CBLAS_ORDER CORDER, enum CBLAS_TRANSPOSE CTRANS, blasint crows,
else else
{ {
OMATCOPY_K_RT(*rows, *cols, *alpha, a, *lda, b, *ldb ); OMATCOPY_K_RT(*rows, *cols, *alpha, a, *lda, b, *ldb );
OMATCOPY_K_RN(*rows, *cols, (FLOAT) 1.0, b, *ldb, a, *ldb ); OMATCOPY_K_RN(*cols, *rows, (FLOAT) 1.0, b, *ldb, a, *ldb );
} }
} }