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

[PATCH 34/42] Fix lapacke_?sytri2 - work is real, not complex

diff --git a/LAPACKE/include/lapacke.h b/LAPACKE/include/lapacke.h
index 31b8ed2..45cf358 100644 (file)
--- a/LAPACKE/include/lapacke.h
+++ b/LAPACKE/include/lapacke.h
@@ -10682,7 +10682,7 @@ lapack_int LAPACKE_dsytri2( int matrix_layout, char uplo, lapack_int n,
 lapack_int LAPACKE_dsytri2_work( int matrix_layout, char uplo, lapack_int n,
                                  double* a, lapack_int lda,
                                  const lapack_int* ipiv,
-                                 lapack_complex_double* work, lapack_int lwork );
+                                 double* work, lapack_int lwork );
 lapack_int LAPACKE_dsytri2x( int matrix_layout, char uplo, lapack_int n,
                              double* a, lapack_int lda, const lapack_int* ipiv,
                              lapack_int nb );
@@ -10775,7 +10775,7 @@ lapack_int LAPACKE_ssytri2( int matrix_layout, char uplo, lapack_int n, float* a
 lapack_int LAPACKE_ssytri2_work( int matrix_layout, char uplo, lapack_int n,
                                  float* a, lapack_int lda,
                                  const lapack_int* ipiv,
-                                 lapack_complex_float* work, lapack_int lwork );
+                                 float* work, lapack_int lwork );
 lapack_int LAPACKE_ssytri2x( int matrix_layout, char uplo, lapack_int n,
                              float* a, lapack_int lda, const lapack_int* ipiv,
                              lapack_int nb );
@@ -17179,7 +17179,7 @@ void LAPACK_dsyswapr( char* uplo, lapack_int* n,
 void LAPACK_dsytri2( char* uplo, lapack_int* n,
                      double* a, lapack_int* lda,
                      const lapack_int* ipiv,
-                     lapack_complex_double* work, lapack_int* lwork , lapack_int *info );
+                     double* work, lapack_int* lwork , lapack_int *info );
 void LAPACK_dsytri2x( char* uplo, lapack_int* n,
                       double* a, lapack_int* lda,
                       const lapack_int* ipiv, double* work,
@@ -17233,7 +17233,7 @@ void LAPACK_ssyswapr( char* uplo, lapack_int* n,
 void LAPACK_ssytri2( char* uplo, lapack_int* n,
                      float* a, lapack_int* lda,
                      const lapack_int* ipiv,
-                     lapack_complex_float* work, lapack_int* lwork , lapack_int *info );
+                     float* work, lapack_int* lwork , lapack_int *info );
 void LAPACK_ssytri2x( char* uplo, lapack_int* n,
                       float* a, lapack_int* lda,
                       const lapack_int* ipiv, float* work,

diff --git a/LAPACKE/src/lapacke_dsytri2.c b/LAPACKE/src/lapacke_dsytri2.c
index 4797491..3fa583a 100644 (file)
--- a/LAPACKE/src/lapacke_dsytri2.c
+++ b/LAPACKE/src/lapacke_dsytri2.c
@@ -38,8 +38,8 @@ lapack_int LAPACKE_dsytri2( int matrix_layout, char uplo, lapack_int n,
 {
     lapack_int info = 0;
     lapack_int lwork = -1;
-    lapack_complex_double* work = NULL;
-    lapack_complex_double work_query;
+    double* work = NULL;
+    double work_query;
     if( matrix_layout != LAPACK_COL_MAJOR && matrix_layout != LAPACK_ROW_MAJOR ) {
         LAPACKE_xerbla( "LAPACKE_dsytri2", -1 );
         return -1;
@@ -58,8 +58,8 @@ lapack_int LAPACKE_dsytri2( int matrix_layout, char uplo, lapack_int n,
     }
     lwork = LAPACK_Z2INT( work_query );
     /* Allocate memory for work arrays */
-    work = (lapack_complex_double*)
-        LAPACKE_malloc( sizeof(lapack_complex_double) * lwork );
+    work = (double*)
+        LAPACKE_malloc( sizeof(double) * lwork );
     if( work == NULL ) {
         info = LAPACK_WORK_MEMORY_ERROR;
         goto exit_level_0;

diff --git a/LAPACKE/src/lapacke_dsytri2_work.c b/LAPACKE/src/lapacke_dsytri2_work.c
index 2168baf..dba46c1 100644 (file)
--- a/LAPACKE/src/lapacke_dsytri2_work.c
+++ b/LAPACKE/src/lapacke_dsytri2_work.c
@@ -36,7 +36,7 @@
 lapack_int LAPACKE_dsytri2_work( int matrix_layout, char uplo, lapack_int n,
                                  double* a, lapack_int lda,
                                  const lapack_int* ipiv,
-                                 lapack_complex_double* work, lapack_int lwork )
+                                 double* work, lapack_int lwork )
 {
     lapack_int info = 0;
     if( matrix_layout == LAPACK_COL_MAJOR ) {

diff --git a/LAPACKE/src/lapacke_ssytri2.c b/LAPACKE/src/lapacke_ssytri2.c
index c817a8b..fe77c13 100644 (file)
--- a/LAPACKE/src/lapacke_ssytri2.c
+++ b/LAPACKE/src/lapacke_ssytri2.c
@@ -38,8 +38,8 @@ lapack_int LAPACKE_ssytri2( int matrix_layout, char uplo, lapack_int n, float* a
 {
     lapack_int info = 0;
     lapack_int lwork = -1;
-    lapack_complex_float* work = NULL;
-    lapack_complex_float work_query;
+    float* work = NULL;
+    float work_query;
     if( matrix_layout != LAPACK_COL_MAJOR && matrix_layout != LAPACK_ROW_MAJOR ) {
         LAPACKE_xerbla( "LAPACKE_ssytri2", -1 );
         return -1;
@@ -58,8 +58,8 @@ lapack_int LAPACKE_ssytri2( int matrix_layout, char uplo, lapack_int n, float* a
     }
     lwork = LAPACK_C2INT( work_query );
     /* Allocate memory for work arrays */
-    work = (lapack_complex_float*)
-        LAPACKE_malloc( sizeof(lapack_complex_float) * lwork );
+    work = (float*)
+        LAPACKE_malloc( sizeof(float) * lwork );
     if( work == NULL ) {
         info = LAPACK_WORK_MEMORY_ERROR;
         goto exit_level_0;

diff --git a/LAPACKE/src/lapacke_ssytri2_work.c b/LAPACKE/src/lapacke_ssytri2_work.c
index 1084194..b678bec 100644 (file)
--- a/LAPACKE/src/lapacke_ssytri2_work.c
+++ b/LAPACKE/src/lapacke_ssytri2_work.c
@@ -36,7 +36,7 @@
 lapack_int LAPACKE_ssytri2_work( int matrix_layout, char uplo, lapack_int n,
                                  float* a, lapack_int lda,
                                  const lapack_int* ipiv,
-                                 lapack_complex_float* work, lapack_int lwork )
+                                 float* work, lapack_int lwork )
 {
     lapack_int info = 0;
     if( matrix_layout == LAPACK_COL_MAJOR ) {