Use in-place transform shortcut only if matrix is square

martin-frbg · web-flow · commit 254db9bd7cac · 2017-09-03T09:52:55.000+02:00
diff --git a/interface/imatcopy.c b/interface/imatcopy.c
@@ -121,7 +121,7 @@ void CNAME( enum CBLAS_ORDER CORDER, enum CBLAS_TRANSPOSE CTRANS, blasint crows,
     		return;
   	}
 #ifdef NEW_IMATCOPY
-    if ( *lda == *ldb ) {
+    if ( *lda == *ldb && *cols == *rows ) {
         if ( order == BlasColMajor )
         {
             if ( trans == BlasNoTrans )

Original file line number	Diff line number	Diff line change
`@@ -121,7 +121,7 @@ void CNAME( enum CBLAS_ORDER CORDER, enum CBLAS_TRANSPOSE CTRANS, blasint crows,`
`121`	`121`	`return;`
`122`	`122`	`}`
`123`	`123`	`#ifdef NEW_IMATCOPY`
`124`		`- if ( lda == ldb ) {`
	`124`	`+ if ( lda == ldb && cols == rows ) {`
`125`	`125`	`if ( order == BlasColMajor )`
`126`	`126`	`{`
`127`	`127`	`if ( trans == BlasNoTrans )`