Skip to content

Commit e2a8c35

Browse files
martin-frbgmartin-frbg
authored
Fixes from netlib PR253
LAPACKE interfaces for Aasen's functions now call ?sytrf_aa and ?hetrf_aa instead of ?sytrf and ?hetrf
1 parent f5959f2 commit e2a8c35

File tree

6 files changed

+18
-18
lines changed

6 files changed

+18
-18
lines changed

‎lapack-netlib/LAPACKE/src/lapacke_chetrf_aa_work.c

Lines changed: 3 additions & 3 deletions
Original file line numberDiff line numberDiff line change
@@ -41,7 +41,7 @@ lapack_int LAPACKE_chetrf_aa_work( int matrix_layout, char uplo, lapack_int n,
4141
lapack_int info = 0;
4242
if( matrix_layout == LAPACK_COL_MAJOR ) {
4343
/* Call LAPACK function and adjust info */
44-
LAPACK_chetrf( &uplo, &n, a, &lda, ipiv, work, &lwork, &info );
44+
LAPACK_chetrf_aa( &uplo, &n, a, &lda, ipiv, work, &lwork, &info );
4545
if( info < 0 ) {
4646
info = info - 1;
4747
}
@@ -56,7 +56,7 @@ lapack_int LAPACKE_chetrf_aa_work( int matrix_layout, char uplo, lapack_int n,
5656
}
5757
/* Query optimal working array(s) size if requested */
5858
if( lwork == -1 ) {
59-
LAPACK_chetrf( &uplo, &n, a, &lda_t, ipiv, work, &lwork, &info );
59+
LAPACK_chetrf_aa( &uplo, &n, a, &lda_t, ipiv, work, &lwork, &info );
6060
return (info < 0) ? (info - 1) : info;
6161
}
6262
/* Allocate memory for temporary array(s) */
@@ -69,7 +69,7 @@ lapack_int LAPACKE_chetrf_aa_work( int matrix_layout, char uplo, lapack_int n,
6969
/* Transpose input matrices */
7070
LAPACKE_che_trans( matrix_layout, uplo, n, a, lda, a_t, lda_t );
7171
/* Call LAPACK function and adjust info */
72-
LAPACK_chetrf( &uplo, &n, a_t, &lda_t, ipiv, work, &lwork, &info );
72+
LAPACK_chetrf_aa( &uplo, &n, a_t, &lda_t, ipiv, work, &lwork, &info );
7373
if( info < 0 ) {
7474
info = info - 1;
7575
}

‎lapack-netlib/LAPACKE/src/lapacke_csytrf_aa_work.c

Lines changed: 3 additions & 3 deletions
Original file line numberDiff line numberDiff line change
@@ -41,7 +41,7 @@ lapack_int LAPACKE_csytrf_aa_work( int matrix_layout, char uplo, lapack_int n,
4141
lapack_int info = 0;
4242
if( matrix_layout == LAPACK_COL_MAJOR ) {
4343
/* Call LAPACK function and adjust info */
44-
LAPACK_csytrf( &uplo, &n, a, &lda, ipiv, work, &lwork, &info );
44+
LAPACK_csytrf_aa( &uplo, &n, a, &lda, ipiv, work, &lwork, &info );
4545
if( info < 0 ) {
4646
info = info - 1;
4747
}
@@ -56,7 +56,7 @@ lapack_int LAPACKE_csytrf_aa_work( int matrix_layout, char uplo, lapack_int n,
5656
}
5757
/* Query optimal working array(s) size if requested */
5858
if( lwork == -1 ) {
59-
LAPACK_csytrf( &uplo, &n, a, &lda_t, ipiv, work, &lwork, &info );
59+
LAPACK_csytrf_aa( &uplo, &n, a, &lda_t, ipiv, work, &lwork, &info );
6060
return (info < 0) ? (info - 1) : info;
6161
}
6262
/* Allocate memory for temporary array(s) */
@@ -69,7 +69,7 @@ lapack_int LAPACKE_csytrf_aa_work( int matrix_layout, char uplo, lapack_int n,
6969
/* Transpose input matrices */
7070
LAPACKE_csy_trans( matrix_layout, uplo, n, a, lda, a_t, lda_t );
7171
/* Call LAPACK function and adjust info */
72-
LAPACK_csytrf( &uplo, &n, a_t, &lda_t, ipiv, work, &lwork, &info );
72+
LAPACK_csytrf_aa( &uplo, &n, a_t, &lda_t, ipiv, work, &lwork, &info );
7373
if( info < 0 ) {
7474
info = info - 1;
7575
}

‎lapack-netlib/LAPACKE/src/lapacke_dsytrf_aa_work.c

Lines changed: 3 additions & 3 deletions
Original file line numberDiff line numberDiff line change
@@ -40,7 +40,7 @@ lapack_int LAPACKE_dsytrf_aa_work( int matrix_layout, char uplo, lapack_int n,
4040
lapack_int info = 0;
4141
if( matrix_layout == LAPACK_COL_MAJOR ) {
4242
/* Call LAPACK function and adjust info */
43-
LAPACK_dsytrf( &uplo, &n, a, &lda, ipiv, work, &lwork, &info );
43+
LAPACK_dsytrf_aa( &uplo, &n, a, &lda, ipiv, work, &lwork, &info );
4444
if( info < 0 ) {
4545
info = info - 1;
4646
}
@@ -55,7 +55,7 @@ lapack_int LAPACKE_dsytrf_aa_work( int matrix_layout, char uplo, lapack_int n,
5555
}
5656
/* Query optimal working array(s) size if requested */
5757
if( lwork == -1 ) {
58-
LAPACK_dsytrf( &uplo, &n, a, &lda_t, ipiv, work, &lwork, &info );
58+
LAPACK_dsytrf_aa( &uplo, &n, a, &lda_t, ipiv, work, &lwork, &info );
5959
return (info < 0) ? (info - 1) : info;
6060
}
6161
/* Allocate memory for temporary array(s) */
@@ -67,7 +67,7 @@ lapack_int LAPACKE_dsytrf_aa_work( int matrix_layout, char uplo, lapack_int n,
6767
/* Transpose input matrices */
6868
LAPACKE_dsy_trans( matrix_layout, uplo, n, a, lda, a_t, lda_t );
6969
/* Call LAPACK function and adjust info */
70-
LAPACK_dsytrf( &uplo, &n, a_t, &lda_t, ipiv, work, &lwork, &info );
70+
LAPACK_dsytrf_aa( &uplo, &n, a_t, &lda_t, ipiv, work, &lwork, &info );
7171
if( info < 0 ) {
7272
info = info - 1;
7373
}

‎lapack-netlib/LAPACKE/src/lapacke_ssytrf_aa_work.c

Lines changed: 3 additions & 3 deletions
Original file line numberDiff line numberDiff line change
@@ -40,7 +40,7 @@ lapack_int LAPACKE_ssytrf_aa_work( int matrix_layout, char uplo, lapack_int n,
4040
lapack_int info = 0;
4141
if( matrix_layout == LAPACK_COL_MAJOR ) {
4242
/* Call LAPACK function and adjust info */
43-
LAPACK_ssytrf( &uplo, &n, a, &lda, ipiv, work, &lwork, &info );
43+
LAPACK_ssytrf_aa( &uplo, &n, a, &lda, ipiv, work, &lwork, &info );
4444
if( info < 0 ) {
4545
info = info - 1;
4646
}
@@ -55,7 +55,7 @@ lapack_int LAPACKE_ssytrf_aa_work( int matrix_layout, char uplo, lapack_int n,
5555
}
5656
/* Query optimal working array(s) size if requested */
5757
if( lwork == -1 ) {
58-
LAPACK_ssytrf( &uplo, &n, a, &lda_t, ipiv, work, &lwork, &info );
58+
LAPACK_ssytrf_aa( &uplo, &n, a, &lda_t, ipiv, work, &lwork, &info );
5959
return (info < 0) ? (info - 1) : info;
6060
}
6161
/* Allocate memory for temporary array(s) */
@@ -67,7 +67,7 @@ lapack_int LAPACKE_ssytrf_aa_work( int matrix_layout, char uplo, lapack_int n,
6767
/* Transpose input matrices */
6868
LAPACKE_ssy_trans( matrix_layout, uplo, n, a, lda, a_t, lda_t );
6969
/* Call LAPACK function and adjust info */
70-
LAPACK_ssytrf( &uplo, &n, a_t, &lda_t, ipiv, work, &lwork, &info );
70+
LAPACK_ssytrf_aa( &uplo, &n, a_t, &lda_t, ipiv, work, &lwork, &info );
7171
if( info < 0 ) {
7272
info = info - 1;
7373
}

‎lapack-netlib/LAPACKE/src/lapacke_zhetrf_aa_work.c

Lines changed: 3 additions & 3 deletions
Original file line numberDiff line numberDiff line change
@@ -41,7 +41,7 @@ lapack_int LAPACKE_zhetrf_aa_work( int matrix_layout, char uplo, lapack_int n,
4141
lapack_int info = 0;
4242
if( matrix_layout == LAPACK_COL_MAJOR ) {
4343
/* Call LAPACK function and adjust info */
44-
LAPACK_zhetrf( &uplo, &n, a, &lda, ipiv, work, &lwork, &info );
44+
LAPACK_zhetrf_aa( &uplo, &n, a, &lda, ipiv, work, &lwork, &info );
4545
if( info < 0 ) {
4646
info = info - 1;
4747
}
@@ -56,7 +56,7 @@ lapack_int LAPACKE_zhetrf_aa_work( int matrix_layout, char uplo, lapack_int n,
5656
}
5757
/* Query optimal working array(s) size if requested */
5858
if( lwork == -1 ) {
59-
LAPACK_zhetrf( &uplo, &n, a, &lda_t, ipiv, work, &lwork, &info );
59+
LAPACK_zhetrf_aa( &uplo, &n, a, &lda_t, ipiv, work, &lwork, &info );
6060
return (info < 0) ? (info - 1) : info;
6161
}
6262
/* Allocate memory for temporary array(s) */
@@ -69,7 +69,7 @@ lapack_int LAPACKE_zhetrf_aa_work( int matrix_layout, char uplo, lapack_int n,
6969
/* Transpose input matrices */
7070
LAPACKE_zhe_trans( matrix_layout, uplo, n, a, lda, a_t, lda_t );
7171
/* Call LAPACK function and adjust info */
72-
LAPACK_zhetrf( &uplo, &n, a_t, &lda_t, ipiv, work, &lwork, &info );
72+
LAPACK_zhetrf_aa( &uplo, &n, a_t, &lda_t, ipiv, work, &lwork, &info );
7373
if( info < 0 ) {
7474
info = info - 1;
7575
}

‎lapack-netlib/LAPACKE/src/lapacke_zsytrf_aa_work.c

Lines changed: 3 additions & 3 deletions
Original file line numberDiff line numberDiff line change
@@ -41,7 +41,7 @@ lapack_int LAPACKE_zsytrf_aa_work( int matrix_layout, char uplo, lapack_int n,
4141
lapack_int info = 0;
4242
if( matrix_layout == LAPACK_COL_MAJOR ) {
4343
/* Call LAPACK function and adjust info */
44-
LAPACK_zsytrf( &uplo, &n, a, &lda, ipiv, work, &lwork, &info );
44+
LAPACK_zsytrf_aa( &uplo, &n, a, &lda, ipiv, work, &lwork, &info );
4545
if( info < 0 ) {
4646
info = info - 1;
4747
}
@@ -56,7 +56,7 @@ lapack_int LAPACKE_zsytrf_aa_work( int matrix_layout, char uplo, lapack_int n,
5656
}
5757
/* Query optimal working array(s) size if requested */
5858
if( lwork == -1 ) {
59-
LAPACK_zsytrf( &uplo, &n, a, &lda_t, ipiv, work, &lwork, &info );
59+
LAPACK_zsytrf_aa( &uplo, &n, a, &lda_t, ipiv, work, &lwork, &info );
6060
return (info < 0) ? (info - 1) : info;
6161
}
6262
/* Allocate memory for temporary array(s) */
@@ -69,7 +69,7 @@ lapack_int LAPACKE_zsytrf_aa_work( int matrix_layout, char uplo, lapack_int n,
6969
/* Transpose input matrices */
7070
LAPACKE_zsy_trans( matrix_layout, uplo, n, a, lda, a_t, lda_t );
7171
/* Call LAPACK function and adjust info */
72-
LAPACK_zsytrf( &uplo, &n, a_t, &lda_t, ipiv, work, &lwork, &info );
72+
LAPACK_zsytrf_aa( &uplo, &n, a_t, &lda_t, ipiv, work, &lwork, &info );
7373
if( info < 0 ) {
7474
info = info - 1;
7575
}

0 commit comments

Comments
 (0)