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

[PATCH 25/42] Fix lapacke_?gejsv - do not transpose U,V on input

- U,V are [out], so no need to transpose on input

diff --git a/LAPACKE/src/lapacke_cgejsv_work.c b/LAPACKE/src/lapacke_cgejsv_work.c
index 72c3ddce72f36cbefeda5e471dfabd67bb642c92..e9ce6e95322f275a9e1fd54ab8dd99aeebcc05a2 100644 (file)
--- a/LAPACKE/src/lapacke_cgejsv_work.c
+++ b/LAPACKE/src/lapacke_cgejsv_work.c
@@ -103,14 +103,6 @@ lapack_int LAPACKE_cgejsv_work( int matrix_layout, char joba, char jobu,
         }
         /* Transpose input matrices */
         LAPACKE_cge_trans( matrix_layout, m, n, a, lda, a_t, lda_t );
-        if( LAPACKE_lsame( jobu, 'f' ) || LAPACKE_lsame( jobu, 'u' ) ||
-            LAPACKE_lsame( jobu, 'w' ) ) {
-            LAPACKE_cge_trans( matrix_layout, nu, n, u, ldu, u_t, ldu_t );
-        }
-        if( LAPACKE_lsame( jobv, 'j' ) || LAPACKE_lsame( jobv, 'v' ) ||
-            LAPACKE_lsame( jobv, 'w' ) ) {
-            LAPACKE_cge_trans( matrix_layout, nv, n, v, ldv, v_t, ldv_t );
-        }
         /* Call LAPACK function and adjust info */
         LAPACK_cgejsv( &joba, &jobu, &jobv, &jobr, &jobt, &jobp, &m, &n, a_t,
                        &lda_t, sva, u_t, &ldu_t, v_t, &ldv_t, cwork, &lwork,

diff --git a/LAPACKE/src/lapacke_dgejsv_work.c b/LAPACKE/src/lapacke_dgejsv_work.c
index 379e63d2d1fa7cfba7d0f40a4933dfe150578dc5..ac6eb724cbedc3bac0ef04dc85406c73690aa9a7 100644 (file)
--- a/LAPACKE/src/lapacke_dgejsv_work.c
+++ b/LAPACKE/src/lapacke_dgejsv_work.c
@@ -99,14 +99,6 @@ lapack_int LAPACKE_dgejsv_work( int matrix_layout, char joba, char jobu,
         }
         /* Transpose input matrices */
         LAPACKE_dge_trans( matrix_layout, m, n, a, lda, a_t, lda_t );
-        if( LAPACKE_lsame( jobu, 'f' ) || LAPACKE_lsame( jobu, 'u' ) ||
-            LAPACKE_lsame( jobu, 'w' ) ) {
-            LAPACKE_dge_trans( matrix_layout, nu, n, u, ldu, u_t, ldu_t );
-        }
-        if( LAPACKE_lsame( jobv, 'j' ) || LAPACKE_lsame( jobv, 'v' ) ||
-            LAPACKE_lsame( jobv, 'w' ) ) {
-            LAPACKE_dge_trans( matrix_layout, nv, n, v, ldv, v_t, ldv_t );
-        }
         /* Call LAPACK function and adjust info */
         LAPACK_dgejsv( &joba, &jobu, &jobv, &jobr, &jobt, &jobp, &m, &n, a_t,
                        &lda_t, sva, u_t, &ldu_t, v_t, &ldv_t, work, &lwork,

diff --git a/LAPACKE/src/lapacke_sgejsv_work.c b/LAPACKE/src/lapacke_sgejsv_work.c
index fe2fa8fb076aa75161c11e26abaf317df2d3ba69..6a29d157a3703b2cd24a5145e1b715981b03965f 100644 (file)
--- a/LAPACKE/src/lapacke_sgejsv_work.c
+++ b/LAPACKE/src/lapacke_sgejsv_work.c
@@ -99,14 +99,6 @@ lapack_int LAPACKE_sgejsv_work( int matrix_layout, char joba, char jobu,
         }
         /* Transpose input matrices */
         LAPACKE_sge_trans( matrix_layout, m, n, a, lda, a_t, lda_t );
-        if( LAPACKE_lsame( jobu, 'f' ) || LAPACKE_lsame( jobu, 'u' ) ||
-            LAPACKE_lsame( jobu, 'w' ) ) {
-            LAPACKE_sge_trans( matrix_layout, nu, n, u, ldu, u_t, ldu_t );
-        }
-        if( LAPACKE_lsame( jobv, 'j' ) || LAPACKE_lsame( jobv, 'v' ) ||
-            LAPACKE_lsame( jobv, 'w' ) ) {
-            LAPACKE_sge_trans( matrix_layout, nv, n, v, ldv, v_t, ldv_t );
-        }
         /* Call LAPACK function and adjust info */
         LAPACK_sgejsv( &joba, &jobu, &jobv, &jobr, &jobt, &jobp, &m, &n, a_t,
                        &lda_t, sva, u_t, &ldu_t, v_t, &ldv_t, work, &lwork,

diff --git a/LAPACKE/src/lapacke_zgejsv_work.c b/LAPACKE/src/lapacke_zgejsv_work.c
index 5d785e9a23bcf872dd76c658ee5741cb97addc6f..16d81b2552e3dd6f07aa64f195a1654a6a8ec1b7 100644 (file)
--- a/LAPACKE/src/lapacke_zgejsv_work.c
+++ b/LAPACKE/src/lapacke_zgejsv_work.c
@@ -103,14 +103,6 @@ lapack_int LAPACKE_zgejsv_work( int matrix_layout, char joba, char jobu,
         }
         /* Transpose input matrices */
         LAPACKE_zge_trans( matrix_layout, m, n, a, lda, a_t, lda_t );
-        if( LAPACKE_lsame( jobu, 'f' ) || LAPACKE_lsame( jobu, 'u' ) ||
-            LAPACKE_lsame( jobu, 'w' ) ) {
-            LAPACKE_zge_trans( matrix_layout, nu, n, u, ldu, u_t, ldu_t );
-        }
-        if( LAPACKE_lsame( jobv, 'j' ) || LAPACKE_lsame( jobv, 'v' ) ||
-            LAPACKE_lsame( jobv, 'w' ) ) {
-            LAPACKE_zge_trans( matrix_layout, nv, n, v, ldv, v_t, ldv_t );
-        }
         /* Call LAPACK function and adjust info */
         LAPACK_zgejsv( &joba, &jobu, &jobv, &jobr, &jobt, &jobp, &m, &n, a_t,
                        &lda_t, sva, u_t, &ldu_t, v_t, &ldv_t, cwork, &lwork,