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

[PATCH 22/42] Fix lapacke_?ggsvp3_work: V is p-by-p

diff --git a/LAPACKE/src/lapacke_cggsvp3_work.c b/LAPACKE/src/lapacke_cggsvp3_work.c
index d532391..77e77b2 100644 (file)
--- a/LAPACKE/src/lapacke_cggsvp3_work.c
+++ b/LAPACKE/src/lapacke_cggsvp3_work.c
@@ -87,7 +87,7 @@ lapack_int LAPACKE_cggsvp3_work( int matrix_layout, char jobu, char jobv,
             LAPACKE_xerbla( "LAPACKE_cggsvp3_work", info );
             return info;
         }
-        if( ldv < m ) {
+        if( ldv < p ) {
             info = -19;
             LAPACKE_xerbla( "LAPACKE_cggsvp3_work", info );
             return info;
@@ -124,7 +124,7 @@ lapack_int LAPACKE_cggsvp3_work( int matrix_layout, char jobu, char jobv,
         if( LAPACKE_lsame( jobv, 'v' ) ) {
             v_t = (lapack_complex_float*)
                 LAPACKE_malloc( sizeof(lapack_complex_float) *
-                                ldv_t * MAX(1,m) );
+                                ldv_t * MAX(1,p) );
             if( v_t == NULL ) {
                 info = LAPACK_TRANSPOSE_MEMORY_ERROR;
                 goto exit_level_3;
@@ -156,7 +156,7 @@ lapack_int LAPACKE_cggsvp3_work( int matrix_layout, char jobu, char jobv,
             LAPACKE_cge_trans( LAPACK_COL_MAJOR, m, m, u_t, ldu_t, u, ldu );
         }
         if( LAPACKE_lsame( jobv, 'v' ) ) {
-            LAPACKE_cge_trans( LAPACK_COL_MAJOR, p, m, v_t, ldv_t, v, ldv );
+            LAPACKE_cge_trans( LAPACK_COL_MAJOR, p, p, v_t, ldv_t, v, ldv );
         }
         if( LAPACKE_lsame( jobq, 'q' ) ) {
             LAPACKE_cge_trans( LAPACK_COL_MAJOR, n, n, q_t, ldq_t, q, ldq );

diff --git a/LAPACKE/src/lapacke_dggsvp3_work.c b/LAPACKE/src/lapacke_dggsvp3_work.c
index e80a755..73acb27 100644 (file)
--- a/LAPACKE/src/lapacke_dggsvp3_work.c
+++ b/LAPACKE/src/lapacke_dggsvp3_work.c
@@ -84,7 +84,7 @@ lapack_int LAPACKE_dggsvp3_work( int matrix_layout, char jobu, char jobv,
             LAPACKE_xerbla( "LAPACKE_dggsvp3_work", info );
             return info;
         }
-        if( ldv < m ) {
+        if( ldv < p ) {
             info = -19;
             LAPACKE_xerbla( "LAPACKE_dggsvp3_work", info );
             return info;
@@ -115,7 +115,7 @@ lapack_int LAPACKE_dggsvp3_work( int matrix_layout, char jobu, char jobv,
             }
         }
         if( LAPACKE_lsame( jobv, 'v' ) ) {
-            v_t = (double*)LAPACKE_malloc( sizeof(double) * ldv_t * MAX(1,m) );
+            v_t = (double*)LAPACKE_malloc( sizeof(double) * ldv_t * MAX(1,p) );
             if( v_t == NULL ) {
                 info = LAPACK_TRANSPOSE_MEMORY_ERROR;
                 goto exit_level_3;
@@ -145,7 +145,7 @@ lapack_int LAPACKE_dggsvp3_work( int matrix_layout, char jobu, char jobv,
             LAPACKE_dge_trans( LAPACK_COL_MAJOR, m, m, u_t, ldu_t, u, ldu );
         }
         if( LAPACKE_lsame( jobv, 'v' ) ) {
-            LAPACKE_dge_trans( LAPACK_COL_MAJOR, p, m, v_t, ldv_t, v, ldv );
+            LAPACKE_dge_trans( LAPACK_COL_MAJOR, p, p, v_t, ldv_t, v, ldv );
         }
         if( LAPACKE_lsame( jobq, 'q' ) ) {
             LAPACKE_dge_trans( LAPACK_COL_MAJOR, n, n, q_t, ldq_t, q, ldq );

diff --git a/LAPACKE/src/lapacke_sggsvp3_work.c b/LAPACKE/src/lapacke_sggsvp3_work.c
index ec200ca..36c74c2 100644 (file)
--- a/LAPACKE/src/lapacke_sggsvp3_work.c
+++ b/LAPACKE/src/lapacke_sggsvp3_work.c
@@ -84,7 +84,7 @@ lapack_int LAPACKE_sggsvp3_work( int matrix_layout, char jobu, char jobv,
             LAPACKE_xerbla( "LAPACKE_sggsvp3_work", info );
             return info;
         }
-        if( ldv < m ) {
+        if( ldv < p ) {
             info = -19;
             LAPACKE_xerbla( "LAPACKE_sggsvp3_work", info );
             return info;
@@ -116,7 +116,7 @@ lapack_int LAPACKE_sggsvp3_work( int matrix_layout, char jobu, char jobv,
             }
         }
         if( LAPACKE_lsame( jobv, 'v' ) ) {
-            v_t = (float*)LAPACKE_malloc( sizeof(float) * ldv_t * MAX(1,m) );
+            v_t = (float*)LAPACKE_malloc( sizeof(float) * ldv_t * MAX(1,p) );
             if( v_t == NULL ) {
                 info = LAPACK_TRANSPOSE_MEMORY_ERROR;
                 goto exit_level_3;
@@ -146,7 +146,7 @@ lapack_int LAPACKE_sggsvp3_work( int matrix_layout, char jobu, char jobv,
             LAPACKE_sge_trans( LAPACK_COL_MAJOR, m, m, u_t, ldu_t, u, ldu );
         }
         if( LAPACKE_lsame( jobv, 'v' ) ) {
-            LAPACKE_sge_trans( LAPACK_COL_MAJOR, p, m, v_t, ldv_t, v, ldv );
+            LAPACKE_sge_trans( LAPACK_COL_MAJOR, p, p, v_t, ldv_t, v, ldv );
         }
         if( LAPACKE_lsame( jobq, 'q' ) ) {
             LAPACKE_sge_trans( LAPACK_COL_MAJOR, n, n, q_t, ldq_t, q, ldq );

diff --git a/LAPACKE/src/lapacke_zggsvp3_work.c b/LAPACKE/src/lapacke_zggsvp3_work.c
index a266beb..74af193 100644 (file)
--- a/LAPACKE/src/lapacke_zggsvp3_work.c
+++ b/LAPACKE/src/lapacke_zggsvp3_work.c
@@ -87,7 +87,7 @@ lapack_int LAPACKE_zggsvp3_work( int matrix_layout, char jobu, char jobv,
             LAPACKE_xerbla( "LAPACKE_zggsvp3_work", info );
             return info;
         }
-        if( ldv < m ) {
+        if( ldv < p ) {
             info = -19;
             LAPACKE_xerbla( "LAPACKE_zggsvp3_work", info );
             return info;
@@ -125,7 +125,7 @@ lapack_int LAPACKE_zggsvp3_work( int matrix_layout, char jobu, char jobv,
         if( LAPACKE_lsame( jobv, 'v' ) ) {
             v_t = (lapack_complex_double*)
                 LAPACKE_malloc( sizeof(lapack_complex_double) *
-                                ldv_t * MAX(1,m) );
+                                ldv_t * MAX(1,p) );
             if( v_t == NULL ) {
                 info = LAPACK_TRANSPOSE_MEMORY_ERROR;
                 goto exit_level_3;
@@ -157,7 +157,7 @@ lapack_int LAPACKE_zggsvp3_work( int matrix_layout, char jobu, char jobv,
             LAPACKE_zge_trans( LAPACK_COL_MAJOR, m, m, u_t, ldu_t, u, ldu );
         }
         if( LAPACKE_lsame( jobv, 'v' ) ) {
-            LAPACKE_zge_trans( LAPACK_COL_MAJOR, p, m, v_t, ldv_t, v, ldv );
+            LAPACKE_zge_trans( LAPACK_COL_MAJOR, p, p, v_t, ldv_t, v, ldv );
         }
         if( LAPACKE_lsame( jobq, 'q' ) ) {
             LAPACKE_zge_trans( LAPACK_COL_MAJOR, n, n, q_t, ldq_t, q, ldq );