Skip to content

Commit 9b24dee

Browse files
sergey-v-kuznetsovsergey-v-kuznetsov
committed
Fixes from MKL team in LAPACKE interfaces
1 parent 799ef93 commit 9b24dee

Some content is hidden

Large Commits have some content hidden by default. Use the searchbox below for content that may be hidden.

74 files changed

+161
-103
lines changed

‎LAPACKE/src/lapacke_cgesvd_work.c

Lines changed: 4 additions & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -28,6 +28,7 @@
2828
*****************************************************************************
2929
* Contents: Native middle-level C interface to LAPACK function cgesvd
3030
* Author: Intel Corporation
31+
* Generated November 2015
3132
*****************************************************************************/
3233

3334
#include "lapacke_utils.h"
@@ -55,6 +56,8 @@ lapack_int LAPACKE_cgesvd_work( int matrix_layout, char jobu, char jobvt,
5556
( LAPACKE_lsame( jobu, 's' ) ? MIN(m,n) : 1);
5657
lapack_int nrows_vt = LAPACKE_lsame( jobvt, 'a' ) ? n :
5758
( LAPACKE_lsame( jobvt, 's' ) ? MIN(m,n) : 1);
59+
lapack_int ncols_vt = ( LAPACKE_lsame( jobvt, 'a' ) ||
60+
LAPACKE_lsame( jobvt, 's' ) ) ? n : 1;
5861
lapack_int lda_t = MAX(1,m);
5962
lapack_int ldu_t = MAX(1,nrows_u);
6063
lapack_int ldvt_t = MAX(1,nrows_vt);
@@ -72,7 +75,7 @@ lapack_int LAPACKE_cgesvd_work( int matrix_layout, char jobu, char jobvt,
7275
LAPACKE_xerbla( "LAPACKE_cgesvd_work", info );
7376
return info;
7477
}
75-
if( ldvt < n ) {
78+
if( ldvt < ncols_vt ) {
7679
info = -12;
7780
LAPACKE_xerbla( "LAPACKE_cgesvd_work", info );
7881
return info;

‎LAPACKE/src/lapacke_cheev_work.c

Lines changed: 2 additions & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -28,6 +28,7 @@
2828
*****************************************************************************
2929
* Contents: Native middle-level C interface to LAPACK function cheev
3030
* Author: Intel Corporation
31+
* Generated November 2015
3132
*****************************************************************************/
3233

3334
#include "lapacke_utils.h"
@@ -77,7 +78,7 @@ lapack_int LAPACKE_cheev_work( int matrix_layout, char jobz, char uplo,
7778
info = info - 1;
7879
}
7980
/* Transpose output matrices */
80-
if ( jobz == 'V') {
81+
if ( jobz == 'V' || jobz == 'v' ) {
8182
LAPACKE_cge_trans( LAPACK_COL_MAJOR, n, n, a_t, lda_t, a, lda );
8283
} else {
8384
LAPACKE_che_trans( LAPACK_COL_MAJOR, uplo, n, a_t, lda_t, a, lda );

‎LAPACKE/src/lapacke_cheevd_2stage_work.c

Lines changed: 2 additions & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -28,6 +28,7 @@
2828
*****************************************************************************
2929
* Contents: Native middle-level C interface to LAPACK function cheevd_2stage
3030
* Author: Intel Corporation
31+
* Generated December 2016
3132
*****************************************************************************/
3233

3334
#include "lapacke_utils.h"
@@ -78,7 +79,7 @@ lapack_int LAPACKE_cheevd_2stage_work( int matrix_layout, char jobz, char uplo,
7879
info = info - 1;
7980
}
8081
/* Transpose output matrices */
81-
if ( jobz == 'V') {
82+
if ( jobz == 'V' || jobz == 'v' ) {
8283
LAPACKE_cge_trans( LAPACK_COL_MAJOR, n, n, a_t, lda_t, a, lda );
8384
} else {
8485
LAPACKE_che_trans( LAPACK_COL_MAJOR, uplo, n, a_t, lda_t, a, lda );

‎LAPACKE/src/lapacke_cheevd_work.c

Lines changed: 2 additions & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -28,6 +28,7 @@
2828
*****************************************************************************
2929
* Contents: Native middle-level C interface to LAPACK function cheevd
3030
* Author: Intel Corporation
31+
* Generated November 2015
3132
*****************************************************************************/
3233

3334
#include "lapacke_utils.h"
@@ -78,7 +79,7 @@ lapack_int LAPACKE_cheevd_work( int matrix_layout, char jobz, char uplo,
7879
info = info - 1;
7980
}
8081
/* Transpose output matrices */
81-
if ( jobz == 'V') {
82+
if ( jobz == 'V' || jobz == 'v' ) {
8283
LAPACKE_cge_trans( LAPACK_COL_MAJOR, n, n, a_t, lda_t, a, lda );
8384
} else {
8485
LAPACKE_che_trans( LAPACK_COL_MAJOR, uplo, n, a_t, lda_t, a, lda );

‎LAPACKE/src/lapacke_chegst.c

Lines changed: 3 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -28,13 +28,14 @@
2828
*****************************************************************************
2929
* Contents: Native high-level C interface to LAPACK function chegst
3030
* Author: Intel Corporation
31+
* Generated November 2015
3132
*****************************************************************************/
3233

3334
#include "lapacke_utils.h"
3435

3536
lapack_int LAPACKE_chegst( int matrix_layout, lapack_int itype, char uplo,
3637
lapack_int n, lapack_complex_float* a,
37-
lapack_int lda, lapack_complex_float* b,
38+
lapack_int lda, const lapack_complex_float* b,
3839
lapack_int ldb )
3940
{
4041
if( matrix_layout != LAPACK_COL_MAJOR && matrix_layout != LAPACK_ROW_MAJOR ) {
@@ -47,7 +48,7 @@ lapack_int LAPACKE_chegst( int matrix_layout, lapack_int itype, char uplo,
4748
if( LAPACKE_che_nancheck( matrix_layout, uplo, n, a, lda ) ) {
4849
return -5;
4950
}
50-
if( LAPACKE_cge_nancheck( matrix_layout, n, n, b, ldb ) ) {
51+
if( LAPACKE_che_nancheck( matrix_layout, uplo, n, b, ldb ) ) {
5152
return -7;
5253
}
5354
}

‎LAPACKE/src/lapacke_chegst_work.c

Lines changed: 2 additions & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -28,13 +28,14 @@
2828
*****************************************************************************
2929
* Contents: Native middle-level C interface to LAPACK function chegst
3030
* Author: Intel Corporation
31+
* Generated November 2015
3132
*****************************************************************************/
3233

3334
#include "lapacke_utils.h"
3435

3536
lapack_int LAPACKE_chegst_work( int matrix_layout, lapack_int itype, char uplo,
3637
lapack_int n, lapack_complex_float* a,
37-
lapack_int lda, lapack_complex_float* b,
38+
lapack_int lda, const lapack_complex_float* b,
3839
lapack_int ldb )
3940
{
4041
lapack_int info = 0;

‎LAPACKE/src/lapacke_chegv.c

Lines changed: 3 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -28,6 +28,7 @@
2828
*****************************************************************************
2929
* Contents: Native high-level C interface to LAPACK function chegv
3030
* Author: Intel Corporation
31+
* Generated November 2015
3132
*****************************************************************************/
3233

3334
#include "lapacke_utils.h"
@@ -49,10 +50,10 @@ lapack_int LAPACKE_chegv( int matrix_layout, lapack_int itype, char jobz,
4950
#ifndef LAPACK_DISABLE_NAN_CHECK
5051
if( LAPACKE_get_nancheck() ) {
5152
/* Optionally check input matrices for NaNs */
52-
if( LAPACKE_cge_nancheck( matrix_layout, n, n, a, lda ) ) {
53+
if( LAPACKE_che_nancheck( matrix_layout, uplo, n, a, lda ) ) {
5354
return -6;
5455
}
55-
if( LAPACKE_cge_nancheck( matrix_layout, n, n, b, ldb ) ) {
56+
if( LAPACKE_che_nancheck( matrix_layout, uplo, n, b, ldb ) ) {
5657
return -8;
5758
}
5859
}

‎LAPACKE/src/lapacke_chegv_2stage.c

Lines changed: 3 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -28,6 +28,7 @@
2828
*****************************************************************************
2929
* Contents: Native high-level C interface to LAPACK function chegv_2stage
3030
* Author: Intel Corporation
31+
* Generated December 2016
3132
*****************************************************************************/
3233

3334
#include "lapacke_utils.h"
@@ -49,10 +50,10 @@ lapack_int LAPACKE_chegv_2stage( int matrix_layout, lapack_int itype, char jobz,
4950
#ifndef LAPACK_DISABLE_NAN_CHECK
5051
if( LAPACKE_get_nancheck() ) {
5152
/* Optionally check input matrices for NaNs */
52-
if( LAPACKE_cge_nancheck( matrix_layout, n, n, a, lda ) ) {
53+
if( LAPACKE_che_nancheck( matrix_layout, uplo, n, a, lda ) ) {
5354
return -6;
5455
}
55-
if( LAPACKE_cge_nancheck( matrix_layout, n, n, b, ldb ) ) {
56+
if( LAPACKE_che_nancheck( matrix_layout, uplo, n, b, ldb ) ) {
5657
return -8;
5758
}
5859
}

‎LAPACKE/src/lapacke_chegvd.c

Lines changed: 3 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -28,6 +28,7 @@
2828
*****************************************************************************
2929
* Contents: Native high-level C interface to LAPACK function chegvd
3030
* Author: Intel Corporation
31+
* Generated November 2015
3132
*****************************************************************************/
3233

3334
#include "lapacke_utils.h"
@@ -54,10 +55,10 @@ lapack_int LAPACKE_chegvd( int matrix_layout, lapack_int itype, char jobz,
5455
#ifndef LAPACK_DISABLE_NAN_CHECK
5556
if( LAPACKE_get_nancheck() ) {
5657
/* Optionally check input matrices for NaNs */
57-
if( LAPACKE_cge_nancheck( matrix_layout, n, n, a, lda ) ) {
58+
if( LAPACKE_che_nancheck( matrix_layout, uplo, n, a, lda ) ) {
5859
return -6;
5960
}
60-
if( LAPACKE_cge_nancheck( matrix_layout, n, n, b, ldb ) ) {
61+
if( LAPACKE_che_nancheck( matrix_layout, uplo, n, b, ldb ) ) {
6162
return -8;
6263
}
6364
}

‎LAPACKE/src/lapacke_chegvx.c

Lines changed: 2 additions & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -28,6 +28,7 @@
2828
*****************************************************************************
2929
* Contents: Native high-level C interface to LAPACK function chegvx
3030
* Author: Intel Corporation
31+
* Generated November 2015
3132
*****************************************************************************/
3233

3334
#include "lapacke_utils.h"
@@ -59,7 +60,7 @@ lapack_int LAPACKE_chegvx( int matrix_layout, lapack_int itype, char jobz,
5960
if( LAPACKE_s_nancheck( 1, &abstol, 1 ) ) {
6061
return -15;
6162
}
62-
if( LAPACKE_cge_nancheck( matrix_layout, n, n, b, ldb ) ) {
63+
if( LAPACKE_che_nancheck( matrix_layout, uplo, n, b, ldb ) ) {
6364
return -9;
6465
}
6566
if( LAPACKE_lsame( range, 'v' ) ) {

0 commit comments

Comments
 (0)