Merge pull request #157 from echeresh/e_lapacke

Fixing NaN-checking in LAPACKE interfaces for sub- and super-diagonal parameters

Author: langou

LAPACKE/src/lapacke_checon_3.c

@@ -40,6 +40,7 @@ lapack_int LAPACKE_checon_3( int matrix_layout, char uplo, lapack_int n,
 {
     lapack_int info = 0;
     lapack_complex_float* work = NULL;
+    lapack_int e_start = LAPACKE_lsame( uplo, 'U' ) ? 1 : 0;
     if( matrix_layout != LAPACK_COL_MAJOR && matrix_layout != LAPACK_ROW_MAJOR ) {
         LAPACKE_xerbla( "LAPACKE_checon_3", -1 );
         return -1;
@@ -49,7 +50,7 @@ lapack_int LAPACKE_checon_3( int matrix_layout, char uplo, lapack_int n,
     if( LAPACKE_che_nancheck( matrix_layout, uplo, n, a, lda ) ) {
         return -4;
     }
-    if( LAPACKE_c_nancheck( n, e, 1 ) ) {
+    if( LAPACKE_c_nancheck( n-1, e + e_start, 1 ) ) {
         return -6;
     }
     if( LAPACKE_s_nancheck( 1, &anorm, 1 ) ) {

LAPACKE/src/lapacke_chesv_rk.c

@@ -52,9 +52,6 @@ lapack_int LAPACKE_chesv_rk( int matrix_layout, char uplo, lapack_int n,
     if( LAPACKE_che_nancheck( matrix_layout, uplo, n, a, lda ) ) {
         return -5;
     }
-    if( LAPACKE_c_nancheck( n, e, 1) ) {
-        return -7;
-    }
     if( LAPACKE_cge_nancheck( matrix_layout, n, nrhs, b, ldb ) ) {
         return -10;
     }

LAPACKE/src/lapacke_chetrf_rk.c

@@ -50,9 +50,6 @@ lapack_int LAPACKE_chetrf_rk( int matrix_layout, char uplo, lapack_int n,
     if( LAPACKE_che_nancheck( matrix_layout, uplo, n, a, lda ) ) {
         return -4;
     }
-    if( LAPACKE_c_nancheck( n, e, 1 ) ) {
-        return -6;
-    }
 #endif
     /* Query optimal working array(s) size */
     info = LAPACKE_chetrf_rk_work( matrix_layout, uplo, n, a, lda, e, ipiv,

LAPACKE/src/lapacke_chetri_3.c

@@ -41,6 +41,7 @@ lapack_int LAPACKE_chetri_3( int matrix_layout, char uplo, lapack_int n,
     lapack_int lwork = -1;
     lapack_complex_float* work = NULL;
     lapack_complex_float work_query;
+    lapack_int e_start = LAPACKE_lsame( uplo, 'U' ) ? 1 : 0;
     if( matrix_layout != LAPACK_COL_MAJOR && matrix_layout != LAPACK_ROW_MAJOR ) {
         LAPACKE_xerbla( "LAPACKE_chetri_3", -1 );
         return -1;
@@ -50,7 +51,7 @@ lapack_int LAPACKE_chetri_3( int matrix_layout, char uplo, lapack_int n,
     if( LAPACKE_che_nancheck( matrix_layout, uplo, n, a, lda ) ) {
         return -4;
     }
-    if( LAPACKE_c_nancheck( n, e, 1 ) ) {
+    if( LAPACKE_c_nancheck( n-1, e + e_start, 1 ) ) {
         return -6;
     }
 #endif

LAPACKE/src/lapacke_cstegr.c

@@ -58,7 +58,7 @@ lapack_int LAPACKE_cstegr( int matrix_layout, char jobz, char range,
     if( LAPACKE_s_nancheck( n, d, 1 ) ) {
         return -5;
     }
-    if( LAPACKE_s_nancheck( n, e, 1 ) ) {
+    if( LAPACKE_s_nancheck( n-1, e, 1 ) ) {
         return -6;
     }
     if( LAPACKE_lsame( range, 'v' ) ) {

LAPACKE/src/lapacke_cstein.c

@@ -51,7 +51,7 @@ lapack_int LAPACKE_cstein( int matrix_layout, lapack_int n, const float* d,
     if( LAPACKE_s_nancheck( n, d, 1 ) ) {
         return -3;
     }
-    if( LAPACKE_s_nancheck( n, e, 1 ) ) {
+    if( LAPACKE_s_nancheck( n-1, e, 1 ) ) {
         return -4;
     }
     if( LAPACKE_s_nancheck( n, w, 1 ) ) {

LAPACKE/src/lapacke_csycon_3.c

@@ -40,6 +40,7 @@ lapack_int LAPACKE_csycon_3( int matrix_layout, char uplo, lapack_int n,
 {
     lapack_int info = 0;
     lapack_complex_float* work = NULL;
+    lapack_int e_start = LAPACKE_lsame( uplo, 'U' ) ? 1 : 0;
     if( matrix_layout != LAPACK_COL_MAJOR && matrix_layout != LAPACK_ROW_MAJOR ) {
         LAPACKE_xerbla( "LAPACKE_csycon_3", -1 );
         return -1;
@@ -49,7 +50,7 @@ lapack_int LAPACKE_csycon_3( int matrix_layout, char uplo, lapack_int n,
     if( LAPACKE_csy_nancheck( matrix_layout, uplo, n, a, lda ) ) {
         return -4;
     }
-    if( LAPACKE_c_nancheck( n, e, 1 ) ) {
+    if( LAPACKE_c_nancheck( n-1, e + e_start, 1 ) ) {
         return -6;
     }
     if( LAPACKE_s_nancheck( 1, &anorm, 1 ) ) {

LAPACKE/src/lapacke_csysv_rk.c

@@ -52,9 +52,6 @@ lapack_int LAPACKE_csysv_rk( int matrix_layout, char uplo, lapack_int n,
     if( LAPACKE_csy_nancheck( matrix_layout, uplo, n, a, lda ) ) {
         return -5;
     }
-    if( LAPACKE_c_nancheck( n, e, 1) ) {
-        return -7;
-    }
     if( LAPACKE_cge_nancheck( matrix_layout, n, nrhs, b, ldb ) ) {
         return -9;
     }

LAPACKE/src/lapacke_csytrf_rk.c

@@ -50,9 +50,6 @@ lapack_int LAPACKE_csytrf_rk( int matrix_layout, char uplo, lapack_int n,
     if( LAPACKE_csy_nancheck( matrix_layout, uplo, n, a, lda ) ) {
         return -4;
     }
-    if( LAPACKE_c_nancheck( n, e, 1 ) ) {
-        return -6;
-    }
 #endif
     /* Query optimal working array(s) size */
     info = LAPACKE_csytrf_rk_work( matrix_layout, uplo, n, a, lda, e, ipiv,

LAPACKE/src/lapacke_csytri_3.c

@@ -41,6 +41,7 @@ lapack_int LAPACKE_csytri_3( int matrix_layout, char uplo, lapack_int n,
     lapack_int lwork = -1;
     lapack_complex_float* work = NULL;
     lapack_complex_float work_query;
+    lapack_int e_start = LAPACKE_lsame( uplo, 'U' ) ? 1 : 0;
     if( matrix_layout != LAPACK_COL_MAJOR && matrix_layout != LAPACK_ROW_MAJOR ) {
         LAPACKE_xerbla( "LAPACKE_csytri_3", -1 );
         return -1;
@@ -50,7 +51,7 @@ lapack_int LAPACKE_csytri_3( int matrix_layout, char uplo, lapack_int n,
     if( LAPACKE_csy_nancheck( matrix_layout, uplo, n, a, lda ) ) {
         return -4;
     }
-    if( LAPACKE_c_nancheck( n, e, 1 ) ) {
+    if( LAPACKE_c_nancheck( n-1, e + e_start, 1 ) ) {
         return -6;
     }
 #endif