From fa3c7a4ff1fb32e82603305aece85e489e405bc7 Mon Sep 17 00:00:00 2001 From: julie Date: Tue, 23 Feb 2016 05:36:56 +0000 Subject: [PATCH] APPLYING INTEL PATCHES sent to Julie on Feb 19th 2016 by Dima from INTEL (dmitry.g.baksheev@intel.com) [PATCH 25/42] Fix lapacke_?gejsv - do not transpose U,V on input - U,V are [out], so no need to transpose on input --- LAPACKE/src/lapacke_cgejsv_work.c | 8 -------- LAPACKE/src/lapacke_dgejsv_work.c | 8 -------- LAPACKE/src/lapacke_sgejsv_work.c | 8 -------- LAPACKE/src/lapacke_zgejsv_work.c | 8 -------- 4 files changed, 32 deletions(-) diff --git a/LAPACKE/src/lapacke_cgejsv_work.c b/LAPACKE/src/lapacke_cgejsv_work.c index 72c3ddce..e9ce6e95 100644 --- 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 379e63d2..ac6eb724 100644 --- 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 fe2fa8fb..6a29d157 100644 --- 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 5d785e9a..16d81b25 100644 --- 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, -- 2.34.1