Apply patch from Nadezhda Mozartova from Intel. (Sent to Julien on November

20th.) Thanks Nadezhda! Thanks Intel!

diff --git a/lapacke/src/lapacke_cgesdd_work.c b/lapacke/src/lapacke_cgesdd_work.c
index bad2fd9..ad27f69 100644 (file)
--- a/lapacke/src/lapacke_cgesdd_work.c
+++ b/lapacke/src/lapacke_cgesdd_work.c
@@ -57,7 +57,7 @@ lapack_int LAPACKE_cgesdd_work( int matrix_order, char jobz, lapack_int m,
                              ( LAPACKE_lsame( jobz, 'o' ) && m<n) ) ? m :
                              ( LAPACKE_lsame( jobz, 's' ) ? MIN(m,n) : 1);
         lapack_int nrows_vt = ( LAPACKE_lsame( jobz, 'a' ) ||
-                              ( LAPACKE_lsame( jobz, 'o' ) && m<n) ) ? n :
+                              ( LAPACKE_lsame( jobz, 'o' ) && m>=n) ) ? n :
                               ( LAPACKE_lsame( jobz, 's' ) ? MIN(m,n) : 1);
         lapack_int lda_t = MAX(1,m);
         lapack_int ldu_t = MAX(1,nrows_u);

diff --git a/lapacke/src/lapacke_dgesdd_work.c b/lapacke/src/lapacke_dgesdd_work.c
index 1d8bb4c..60b8fd5 100644 (file)
--- a/lapacke/src/lapacke_dgesdd_work.c
+++ b/lapacke/src/lapacke_dgesdd_work.c
@@ -55,7 +55,7 @@ lapack_int LAPACKE_dgesdd_work( int matrix_order, char jobz, lapack_int m,
                              ( LAPACKE_lsame( jobz, 'o' ) && m<n) ) ? m :
                              ( LAPACKE_lsame( jobz, 's' ) ? MIN(m,n) : 1);
         lapack_int nrows_vt = ( LAPACKE_lsame( jobz, 'a' ) ||
-                              ( LAPACKE_lsame( jobz, 'o' ) && m<n) ) ? n :
+                              ( LAPACKE_lsame( jobz, 'o' ) && m>=n) ) ? n :
                               ( LAPACKE_lsame( jobz, 's' ) ? MIN(m,n) : 1);
         lapack_int lda_t = MAX(1,m);
         lapack_int ldu_t = MAX(1,nrows_u);

diff --git a/lapacke/src/lapacke_sgesdd_work.c b/lapacke/src/lapacke_sgesdd_work.c
index a002dd6..3fc0ec4 100644 (file)
--- a/lapacke/src/lapacke_sgesdd_work.c
+++ b/lapacke/src/lapacke_sgesdd_work.c
@@ -55,7 +55,7 @@ lapack_int LAPACKE_sgesdd_work( int matrix_order, char jobz, lapack_int m,
                              ( LAPACKE_lsame( jobz, 'o' ) && m<n) ) ? m :
                              ( LAPACKE_lsame( jobz, 's' ) ? MIN(m,n) : 1);
         lapack_int nrows_vt = ( LAPACKE_lsame( jobz, 'a' ) ||
-                              ( LAPACKE_lsame( jobz, 'o' ) && m<n) ) ? n :
+                              ( LAPACKE_lsame( jobz, 'o' ) && m>=n) ) ? n :
                               ( LAPACKE_lsame( jobz, 's' ) ? MIN(m,n) : 1);
         lapack_int lda_t = MAX(1,m);
         lapack_int ldu_t = MAX(1,nrows_u);

diff --git a/lapacke/src/lapacke_zgesdd_work.c b/lapacke/src/lapacke_zgesdd_work.c
index 0b0d524..566c652 100644 (file)
--- a/lapacke/src/lapacke_zgesdd_work.c
+++ b/lapacke/src/lapacke_zgesdd_work.c
@@ -57,7 +57,7 @@ lapack_int LAPACKE_zgesdd_work( int matrix_order, char jobz, lapack_int m,
                              ( LAPACKE_lsame( jobz, 'o' ) && m<n) ) ? m :
                              ( LAPACKE_lsame( jobz, 's' ) ? MIN(m,n) : 1);
         lapack_int nrows_vt = ( LAPACKE_lsame( jobz, 'a' ) ||
-                              ( LAPACKE_lsame( jobz, 'o' ) && m<n) ) ? n :
+                              ( LAPACKE_lsame( jobz, 'o' ) && m>=n) ) ? n :
                               ( LAPACKE_lsame( jobz, 's' ) ? MIN(m,n) : 1);
         lapack_int lda_t = MAX(1,m);
         lapack_int ldu_t = MAX(1,nrows_u);