diff options
| author | Vladimir Chalupecky <vladimir.chalupecky@gmail.com> | 2017-02-21 06:32:02 +0100 |
|---|---|---|
| committer | Vladimir Chalupecky <vladimir.chalupecky@gmail.com> | 2017-02-21 06:32:02 +0100 |
| commit | 463e1f5a97d82781bea940b3081457af277e75c2 (patch) | |
| tree | 573ebd1e3b4aa8954cd68e9a11401a32a8767d70 /LAPACKE | |
| parent | 4b03bc241909f81f9b30fee4c1feab250f474809 (diff) | |
| download | lapack-463e1f5a97d82781bea940b3081457af277e75c2.tar.gz lapack-463e1f5a97d82781bea940b3081457af277e75c2.tar.bz2 lapack-463e1f5a97d82781bea940b3081457af277e75c2.zip | |
LAPACKE: fix wrong matrix size in lapacke_?ormbr
Diffstat (limited to 'LAPACKE')
| -rw-r--r-- | LAPACKE/src/lapacke_dormbr.c | 7 | ||||
| -rw-r--r-- | LAPACKE/src/lapacke_sormbr.c | 7 |
2 files changed, 8 insertions, 6 deletions
diff --git a/LAPACKE/src/lapacke_dormbr.c b/LAPACKE/src/lapacke_dormbr.c index 5f1de948..bee4b1ae 100644 --- a/LAPACKE/src/lapacke_dormbr.c +++ b/LAPACKE/src/lapacke_dormbr.c @@ -42,7 +42,7 @@ lapack_int LAPACKE_dormbr( int matrix_layout, char vect, char side, char trans, lapack_int lwork = -1; double* work = NULL; double work_query; - lapack_int nq, r; + lapack_int nq, ar, ac; if( matrix_layout != LAPACK_COL_MAJOR && matrix_layout != LAPACK_ROW_MAJOR ) { LAPACKE_xerbla( "LAPACKE_dormbr", -1 ); return -1; @@ -50,8 +50,9 @@ lapack_int LAPACKE_dormbr( int matrix_layout, char vect, char side, char trans, #ifndef LAPACK_DISABLE_NAN_CHECK /* Optionally check input matrices for NaNs */ nq = LAPACKE_lsame( side, 'l' ) ? m : n; - r = LAPACKE_lsame( vect, 'q' ) ? nq : MIN(nq,k); - if( LAPACKE_dge_nancheck( matrix_layout, r, MIN(nq,k), a, lda ) ) { + ar = LAPACKE_lsame( vect, 'q' ) ? nq : MIN(nq,k); + ac = LAPACKE_lsame( vect, 'q' ) ? MIN(nq,k) : nq; + if( LAPACKE_dge_nancheck( matrix_layout, ar, ac, a, lda ) ) { return -8; } if( LAPACKE_dge_nancheck( matrix_layout, m, n, c, ldc ) ) { diff --git a/LAPACKE/src/lapacke_sormbr.c b/LAPACKE/src/lapacke_sormbr.c index c27896dc..4df425ef 100644 --- a/LAPACKE/src/lapacke_sormbr.c +++ b/LAPACKE/src/lapacke_sormbr.c @@ -42,7 +42,7 @@ lapack_int LAPACKE_sormbr( int matrix_layout, char vect, char side, char trans, lapack_int lwork = -1; float* work = NULL; float work_query; - lapack_int nq, r; + lapack_int nq, ar, ac; if( matrix_layout != LAPACK_COL_MAJOR && matrix_layout != LAPACK_ROW_MAJOR ) { LAPACKE_xerbla( "LAPACKE_sormbr", -1 ); return -1; @@ -50,8 +50,9 @@ lapack_int LAPACKE_sormbr( int matrix_layout, char vect, char side, char trans, #ifndef LAPACK_DISABLE_NAN_CHECK /* Optionally check input matrices for NaNs */ nq = LAPACKE_lsame( side, 'l' ) ? m : n; - r = LAPACKE_lsame( vect, 'q' ) ? nq : MIN(nq,k); - if( LAPACKE_sge_nancheck( matrix_layout, r, MIN(nq,k), a, lda ) ) { + ar = LAPACKE_lsame( vect, 'q' ) ? nq : MIN(nq,k); + ac = LAPACKE_lsame( vect, 'q' ) ? MIN(nq,k) : nq; + if( LAPACKE_sge_nancheck( matrix_layout, ar, ac, a, lda ) ) { return -8; } if( LAPACKE_sge_nancheck( matrix_layout, m, n, c, ldc ) ) { |