APPLYING INTEL PATCHES sent to Julie on Feb 19th 2016 by Dima from INTEL (dmitry...

[PATCH 24/42] Fix lapacke_?gejsv - remove NaN check of U,V

- U,V are [out] and should not be checked for NaNs

diff --git a/LAPACKE/src/lapacke_cgejsv.c b/LAPACKE/src/lapacke_cgejsv.c
index 32143ae..7a79805 100644 (file)
--- a/LAPACKE/src/lapacke_cgejsv.c
+++ b/LAPACKE/src/lapacke_cgejsv.c
@@ -136,18 +136,6 @@ lapack_int LAPACKE_cgejsv( int matrix_layout, char joba, char jobu, char jobv,
     if( LAPACKE_cge_nancheck( matrix_layout, m, n, a, lda ) ) {
         return -10;
     }
-    if( LAPACKE_lsame( jobu, 'f' ) || LAPACKE_lsame( jobu, 'u' ) ||
-        LAPACKE_lsame( jobu, 'w' ) ) {
-        if( LAPACKE_cge_nancheck( matrix_layout, nu, n, u, ldu ) ) {
-            return -13;
-        }
-    }
-    if( LAPACKE_lsame( jobv, 'j' ) || LAPACKE_lsame( jobv, 'v' ) ||
-        LAPACKE_lsame( jobv, 'w' ) ) {
-        if( LAPACKE_cge_nancheck( matrix_layout, nv, n, v, ldv ) ) {
-            return -15;
-        }
-    }
 #endif
     /* Allocate memory for working array(s) */
     iwork = (lapack_int*)LAPACKE_malloc( sizeof(lapack_int) * MAX(1,m+3*n) );

diff --git a/LAPACKE/src/lapacke_dgejsv.c b/LAPACKE/src/lapacke_dgejsv.c
index 093afee..e2eee75 100644 (file)
--- a/LAPACKE/src/lapacke_dgejsv.c
+++ b/LAPACKE/src/lapacke_dgejsv.c
@@ -86,18 +86,6 @@ lapack_int LAPACKE_dgejsv( int matrix_layout, char joba, char jobu, char jobv,
     if( LAPACKE_dge_nancheck( matrix_layout, m, n, a, lda ) ) {
         return -10;
     }
-    if( LAPACKE_lsame( jobu, 'f' ) || LAPACKE_lsame( jobu, 'u' ) ||
-        LAPACKE_lsame( jobu, 'w' ) ) {
-        if( LAPACKE_dge_nancheck( matrix_layout, nu, n, u, ldu ) ) {
-            return -13;
-        }
-    }
-    if( LAPACKE_lsame( jobv, 'j' ) || LAPACKE_lsame( jobv, 'v' ) ||
-        LAPACKE_lsame( jobv, 'w' ) ) {
-        if( LAPACKE_dge_nancheck( matrix_layout, nv, n, v, ldv ) ) {
-            return -15;
-        }
-    }
 #endif
     /* Allocate memory for working array(s) */
     iwork = (lapack_int*)LAPACKE_malloc( sizeof(lapack_int) * MAX(1,m+3*n) );

diff --git a/LAPACKE/src/lapacke_sgejsv.c b/LAPACKE/src/lapacke_sgejsv.c
index 7af9901..62ea03e 100644 (file)
--- a/LAPACKE/src/lapacke_sgejsv.c
+++ b/LAPACKE/src/lapacke_sgejsv.c
@@ -86,18 +86,6 @@ lapack_int LAPACKE_sgejsv( int matrix_layout, char joba, char jobu, char jobv,
     if( LAPACKE_sge_nancheck( matrix_layout, m, n, a, lda ) ) {
         return -10;
     }
-    if( LAPACKE_lsame( jobu, 'f' ) || LAPACKE_lsame( jobu, 'u' ) ||
-        LAPACKE_lsame( jobu, 'w' ) ) {
-        if( LAPACKE_sge_nancheck( matrix_layout, nu, n, u, ldu ) ) {
-            return -13;
-        }
-    }
-    if( LAPACKE_lsame( jobv, 'j' ) || LAPACKE_lsame( jobv, 'v' ) ||
-        LAPACKE_lsame( jobv, 'w' ) ) {
-        if( LAPACKE_sge_nancheck( matrix_layout, nv, n, v, ldv ) ) {
-            return -15;
-        }
-    }
 #endif
     /* Allocate memory for working array(s) */
     iwork = (lapack_int*)LAPACKE_malloc( sizeof(lapack_int) * MAX(1,m+3*n) );

diff --git a/LAPACKE/src/lapacke_zgejsv.c b/LAPACKE/src/lapacke_zgejsv.c
index 7a93886..046704c 100644 (file)
--- a/LAPACKE/src/lapacke_zgejsv.c
+++ b/LAPACKE/src/lapacke_zgejsv.c
@@ -136,18 +136,6 @@ lapack_int LAPACKE_zgejsv( int matrix_layout, char joba, char jobu, char jobv,
     if( LAPACKE_zge_nancheck( matrix_layout, m, n, a, lda ) ) {
         return -10;
     }
-    if( LAPACKE_lsame( jobu, 'f' ) || LAPACKE_lsame( jobu, 'u' ) ||
-        LAPACKE_lsame( jobu, 'w' ) ) {
-        if( LAPACKE_zge_nancheck( matrix_layout, nu, n, u, ldu ) ) {
-            return -13;
-        }
-    }
-    if( LAPACKE_lsame( jobv, 'j' ) || LAPACKE_lsame( jobv, 'v' ) ||
-        LAPACKE_lsame( jobv, 'w' ) ) {
-        if( LAPACKE_zge_nancheck( matrix_layout, nv, n, v, ldv ) ) {
-            return -15;
-        }
-    }
 #endif
     /* Allocate memory for working array(s) */
     iwork = (lapack_int*)LAPACKE_malloc( sizeof(lapack_int) * MAX(1,m+3*n) );