lapack_int lwork, double* rwork );
lapack_int LAPACKE_sgesvdx_work( int matrix_layout, char jobu, char jobvt, char range,
- lapack_int m, lapack_int n, float* a,
- lapack_int lda, float vl, float vu,
- lapack_int il, lapack_int iu, lapack_int* ns,
- float* s, float* u, lapack_int ldu,
- float* vt, lapack_int ldvt,
- float* work, lapack_int lwork, lapack_int* iwork );
+ lapack_int m, lapack_int n, float* a,
+ lapack_int lda, float vl, float vu,
+ lapack_int il, lapack_int iu, lapack_int* ns,
+ float* s, float* u, lapack_int ldu,
+ float* vt, lapack_int ldvt,
+ float* work, lapack_int lwork, lapack_int* iwork );
lapack_int LAPACKE_dgesvdx_work( int matrix_layout, char jobu, char jobvt, char range,
- lapack_int m, lapack_int n, double* a,
- lapack_int lda, double vl, double vu,
- lapack_int il, lapack_int iu, lapack_int* ns,
- double* s, double* u, lapack_int ldu,
- double* vt, lapack_int ldvt,
- double* work, lapack_int lwork, lapack_int* iwork );
+ lapack_int m, lapack_int n, double* a,
+ lapack_int lda, double vl, double vu,
+ lapack_int il, lapack_int iu, lapack_int* ns,
+ double* s, double* u, lapack_int ldu,
+ double* vt, lapack_int ldvt,
+ double* work, lapack_int lwork, lapack_int* iwork );
lapack_int LAPACKE_cgesvdx_work( int matrix_layout, char jobu, char jobvt, char range,
- lapack_int m, lapack_int n, lapack_complex_float* a,
- lapack_int lda, float vl, float vu,
- lapack_int il, lapack_int iu, lapack_int* ns,
- float* s, lapack_complex_float* u, lapack_int ldu,
- lapack_complex_float* vt, lapack_int ldvt,
- lapack_complex_float* work, lapack_int lwork,
- float* rwork, lapack_int* iwork );
+ lapack_int m, lapack_int n, lapack_complex_float* a,
+ lapack_int lda, float vl, float vu,
+ lapack_int il, lapack_int iu, lapack_int* ns,
+ float* s, lapack_complex_float* u, lapack_int ldu,
+ lapack_complex_float* vt, lapack_int ldvt,
+ lapack_complex_float* work, lapack_int lwork,
+ float* rwork, lapack_int* iwork );
lapack_int LAPACKE_zgesvdx_work( int matrix_layout, char jobu, char jobvt, char range,
- lapack_int m, lapack_int n, lapack_complex_double* a,
- lapack_int lda, double vl, double vu,
- lapack_int il, lapack_int iu, lapack_int* ns,
- double* s, lapack_complex_double* u, lapack_int ldu,
- lapack_complex_double* vt, lapack_int ldvt,
- lapack_complex_double* work, lapack_int lwork,
- double* rwork, lapack_int* iwork );
+ lapack_int m, lapack_int n, lapack_complex_double* a,
+ lapack_int lda, double vl, double vu,
+ lapack_int il, lapack_int iu, lapack_int* ns,
+ double* s, lapack_complex_double* u, lapack_int ldu,
+ lapack_complex_double* vt, lapack_int ldvt,
+ lapack_complex_double* work, lapack_int lwork,
+ double* rwork, lapack_int* iwork );
lapack_int LAPACKE_sgesvj_work( int matrix_layout, char joba, char jobu,
char jobv, lapack_int m, lapack_int n, float* a,
lapack_int lda, float* sva, lapack_int mv,
lapack_complex_float* v, lapack_int ldv,
lapack_complex_float* cwork, lapack_int lwork,
- float* rwork,lapack_int lrwork );
+ float* rwork,lapack_int lrwork );
lapack_int LAPACKE_zgesvj_work( int matrix_layout, char joba, char jobu,
char jobv, lapack_int m, lapack_int n,
lapack_complex_double* a, lapack_int lda, double* sva,
lapack_int lda, lapack_complex_double* e, lapack_int* ipiv,
lapack_complex_double* b, lapack_int ldb,
lapack_complex_double* work, lapack_int lwork );
-
+
lapack_int LAPACKE_ssytrf_rk( int matrix_layout, char uplo, lapack_int n, float* a,
lapack_int lda, float* e, lapack_int* ipiv );
lapack_int LAPACKE_dsytrf_rk( int matrix_layout, char uplo, lapack_int n, double* a,
const lapack_int* ipiv, float anorm, float* rcond );
lapack_int LAPACKE_zsycon_3( int matrix_layout, char uplo, lapack_int n,
const lapack_complex_double* a, lapack_int lda,
- const lapack_complex_double* e,
+ const lapack_complex_double* e,
const lapack_int* ipiv, double anorm,
double* rcond );
lapack_int LAPACKE_checon_3( int matrix_layout, char uplo, lapack_int n,
const lapack_complex_float* a, lapack_int lda,
- const lapack_complex_float* e,
+ const lapack_complex_float* e,
const lapack_int* ipiv, float anorm, float* rcond );
lapack_int LAPACKE_zhecon_3( int matrix_layout, char uplo, lapack_int n,
const lapack_complex_double* a, lapack_int lda,
const lapack_complex_double* a, lapack_int lda,
const lapack_complex_double* e,
const lapack_int* ipiv, double anorm,
- double* rcond, lapack_complex_double* work );
+ double* rcond, lapack_complex_double* work );
lapack_int LAPACKE_checon_3_work( int matrix_layout, char uplo, lapack_int n,
const lapack_complex_float* a, lapack_int lda,
const lapack_complex_float* e,
const lapack_int* ipiv,
lapack_complex_double* b, lapack_int* ldb, lapack_int *info );
-void LAPACK_ssytri_3( char* uplo, lapack_int* n, float* a, lapack_int* lda, const float* e,
+void LAPACK_ssytri_3( char* uplo, lapack_int* n, float* a, lapack_int* lda, const float* e,
const lapack_int* ipiv, float* work, lapack_int* lwork, lapack_int *info );
void LAPACK_dsytri_3( char* uplo, lapack_int* n, double* a, lapack_int* lda, const double* e,
const lapack_int* ipiv, double* work, lapack_int* lwork, lapack_int *info );
lapack_int* lda, const lapack_complex_double* e, const lapack_int* ipiv,
lapack_complex_double* work, lapack_int* lwork, lapack_int *info );
-void LAPACK_ssycon_3( char* uplo, lapack_int* n, const float* a, lapack_int* lda, const float* e,
+void LAPACK_ssycon_3( char* uplo, lapack_int* n, const float* a, lapack_int* lda, const float* e,
const lapack_int* ipiv, float* anorm, float* rcond,
float* work, lapack_int* iwork, lapack_int *info );
void LAPACK_dsycon_3( char* uplo, lapack_int* n, const double* a, lapack_int* lda, const double* e,
{
lapack_int info = 0;
lapack_int lwork = (
- // 1.1
- ( LAPACKE_lsame( jobu, 'n' ) && LAPACKE_lsame( jobv, 'n' ) &&
+ // 1.1
+ ( LAPACKE_lsame( jobu, 'n' ) && LAPACKE_lsame( jobv, 'n' ) &&
( LAPACKE_lsame( jobt, 't' ) || LAPACKE_lsame( joba, 'f' ) || LAPACKE_lsame( joba, 'g' ) )) ? 2*n+1 :
//1.2
- ( ( LAPACKE_lsame( jobu, 'n' ) && LAPACKE_lsame( jobv, 'n' ) &&
+ ( ( LAPACKE_lsame( jobu, 'n' ) && LAPACKE_lsame( jobv, 'n' ) &&
!( LAPACKE_lsame( jobt, 't' ) || LAPACKE_lsame( joba, 'f' ) || LAPACKE_lsame( joba, 'g' ) )) ? n*n+3*n :
//2.1
- ( ( ( LAPACKE_lsame( jobv, 'v' ) || LAPACKE_lsame( jobv, 'j' ) ) &&
- !( LAPACKE_lsame( jobu, 'u' ) || LAPACKE_lsame( jobu, 'f' ) ) &&
+ ( ( ( LAPACKE_lsame( jobv, 'v' ) || LAPACKE_lsame( jobv, 'j' ) ) &&
+ !( LAPACKE_lsame( jobu, 'u') || LAPACKE_lsame( jobu, 'f' ) )&&
( LAPACKE_lsame( jobt, 't' ) || LAPACKE_lsame( joba, 'f' ) || LAPACKE_lsame( joba, 'g' ) ))? 3*n :
//2.2
- ( ( ( LAPACKE_lsame( jobv, 'v' ) || LAPACKE_lsame( jobv, 'j' ) ) &&
+ ( ( ( LAPACKE_lsame( jobv, 'v' ) || LAPACKE_lsame( jobv, 'j' ) ) &&
!( LAPACKE_lsame( jobu, 'u' ) || LAPACKE_lsame( jobu, 'f' ) ) &&
!( LAPACKE_lsame( jobt, 't' ) || LAPACKE_lsame( joba, 'f' ) || LAPACKE_lsame( joba, 'g' ) ))? 3*n :
( LAPACKE_lsame( jobt, 't' ) || LAPACKE_lsame( joba, 'f' ) || LAPACKE_lsame( joba, 'g' ) ))? 4*n*n:
1) ) ) ) ) ) ) );
lapack_int lrwork = (
- // 1.1
- ( LAPACKE_lsame( jobu, 'n' ) && LAPACKE_lsame( jobv, 'n' ) &&
+ // 1.1
+ ( LAPACKE_lsame( jobu, 'n' ) && LAPACKE_lsame( jobv, 'n' ) &&
( LAPACKE_lsame( jobt, 't' ) || LAPACKE_lsame( joba, 'f' ) || LAPACKE_lsame( joba, 'g' ) )) ? MAX(7,n+2*m) :
//1.2
!( LAPACKE_lsame( jobt, 't' ) || LAPACKE_lsame( joba, 'f' ) || LAPACKE_lsame( joba, 'g' ) )) ? MAX(7,2*n) :
//2.1
- ( ( ( LAPACKE_lsame( jobv, 'v' ) || LAPACKE_lsame( jobv, 'j' ) ) &&
- !( LAPACKE_lsame( jobu, 'u' ) || LAPACKE_lsame( jobu, 'f' ) ) &&
+ ( ( ( LAPACKE_lsame( jobv, 'v' ) || LAPACKE_lsame( jobv, 'j' ) ) &&
+ !( LAPACKE_lsame( jobu, 'u') || LAPACKE_lsame( jobu, 'f' ) ) &&
( LAPACKE_lsame( jobt, 't' ) || LAPACKE_lsame( joba, 'f' ) || LAPACKE_lsame( joba, 'g' ) ))? MAX( 7, n+ 2*m ) :
//2.2
- ( ( ( LAPACKE_lsame( jobv, 'v' ) || LAPACKE_lsame( jobv, 'j' ) ) &&
+ ( ( ( LAPACKE_lsame( jobv, 'v' ) || LAPACKE_lsame( jobv, 'j' ) ) &&
!( LAPACKE_lsame( jobu, 'u' ) || LAPACKE_lsame( jobu, 'f' ) ) &&
!( LAPACKE_lsame( jobt, 't' ) || LAPACKE_lsame( joba, 'f' ) || LAPACKE_lsame( joba, 'g' ) ))? MAX(7,2*n) :
( ( ( LAPACKE_lsame( jobu, 'u' ) || LAPACKE_lsame( jobu, 'f' ) ) &&
( LAPACKE_lsame( jobv, 'v' ) || LAPACKE_lsame( jobv, 'j' ) ) &&
( LAPACKE_lsame( jobt, 't' ) || LAPACKE_lsame( joba, 'f' ) || LAPACKE_lsame( joba, 'g' ) ))? MAX(7,2*n) :
- 7 ))))))));
+ 7) ) ) ) ) ) ) );
lapack_int* iwork = NULL;
float* rwork = NULL;
lapack_complex_float* cwork = NULL;
#endif
/* Query optimal working array(s) size */
info = LAPACKE_cgesvdx_work( matrix_layout, jobu, jobvt, range,
- m, n, a, lda, vl, vu, il, iu, ns, s, u,
+ m, n, a, lda, vl, vu, il, iu, ns, s, u,
ldu, vt, ldvt, &work_query, lwork, rwork, iwork );
if( info != 0 ) {
goto exit_level_0;
lwork = LAPACK_C2INT (work_query);
/* Allocate memory for work arrays */
work = (lapack_complex_float*)
- LAPACKE_malloc( sizeof(lapack_complex_float) * lwork );
+ LAPACKE_malloc( sizeof(lapack_complex_float) * lwork );
if( work == NULL ) {
info = LAPACK_WORK_MEMORY_ERROR;
goto exit_level_1;
}
/* Call middle-level interface */
info = LAPACKE_cgesvdx_work( matrix_layout, jobu, jobvt, range,
- m, n, a, lda, vl, vu, il, iu, ns, s, u,
- ldu, vt, ldvt, work, lwork, rwork, iwork );
+ m, n, a, lda, vl, vu, il, iu, ns, s, u,
+ ldu, vt, ldvt, work, lwork, rwork, iwork );
/* Backup significant data from working array(s) */
for( i=0; i<12*MIN(m,n)-1; i++ ) {
superb[i] = iwork[i+1];
#include "lapacke_utils.h"
lapack_int LAPACKE_cgesvdx_work( int matrix_layout, char jobu, char jobvt, char range,
- lapack_int m, lapack_int n, lapack_complex_float* a,
- lapack_int lda, float vl, float vu,
- lapack_int il, lapack_int iu, lapack_int* ns,
- float* s, lapack_complex_float* u, lapack_int ldu,
- lapack_complex_float* vt, lapack_int ldvt,
- lapack_complex_float* work, lapack_int lwork,
- float* rwork, lapack_int* iwork )
+ lapack_int m, lapack_int n, lapack_complex_float* a,
+ lapack_int lda, float vl, float vu,
+ lapack_int il, lapack_int iu, lapack_int* ns,
+ float* s, lapack_complex_float* u, lapack_int ldu,
+ lapack_complex_float* vt, lapack_int ldvt,
+ lapack_complex_float* work, lapack_int lwork,
+ float* rwork, lapack_int* iwork )
{
lapack_int info = 0;
if( matrix_layout == LAPACK_COL_MAJOR ) {
/* Call LAPACK function and adjust info */
LAPACK_cgesvdx( &jobu, &jobvt, &range, &m, &n, a, &lda, &vl, &vu,
- &il, &iu, ns, s, u, &ldu, vt, &ldvt,
+ &il, &iu, ns, s, u, &ldu, vt, &ldvt,
work, &lwork, rwork, iwork, &info );
if( info < 0 ) {
info = info - 1;
/* Query optimal working array(s) size if requested */
if( lwork == -1 ) {
LAPACK_cgesvdx( &jobu, &jobvt, &range, &m, &n, a, &lda_t, &vl, &vu,
- &il, &iu, ns, s, u, &ldu_t, vt,
+ &il, &iu, ns, s, u, &ldu_t, vt,
&ldvt_t, work, &lwork, rwork, iwork, &info );
return (info < 0) ? (info - 1) : info;
}
LAPACKE_cge_trans( matrix_layout, m, n, a, lda, a_t, lda_t );
/* Call LAPACK function and adjust info */
LAPACK_cgesvdx( &jobu, &jobvt, &range, &m, &n, a_t, &lda_t, &vl, &vu,
- &il, &iu, ns, s, u_t, &ldu_t, vt_t,
- &ldvt_t, work, &lwork, rwork, iwork, &info );
+ &il, &iu, ns, s, u_t, &ldu_t, vt_t,
+ &ldvt_t, work, &lwork, rwork, iwork, &info );
if( info < 0 ) {
info = info - 1;
}
}
if( LAPACKE_lsame( jobv, 'a' ) || LAPACKE_lsame( jobv, 'v' ) ) {
v_t = (lapack_complex_float*)
- LAPACKE_malloc( sizeof(lapack_complex_float) * ldv_t * MAX(1,n) );
+ LAPACKE_malloc( sizeof(lapack_complex_float) * ldv_t * MAX(1,n) );
if( v_t == NULL ) {
info = LAPACK_TRANSPOSE_MEMORY_ERROR;
goto exit_level_1;
{
lapack_int info = 0;
float* rwork = NULL;
- lapack_int lwork = -1;
+ lapack_int lwork = -1;
lapack_complex_float* work = NULL;
lapack_complex_float work_query;
if( matrix_layout != LAPACK_COL_MAJOR && matrix_layout != LAPACK_ROW_MAJOR ) {
info = LAPACKE_cggsvd3_work( matrix_layout, jobu, jobv, jobq, m, n, p, k, l,
a, lda, b, ldb, alpha, beta, u, ldu, v, ldv, q,
ldq, &work_query, lwork, rwork, iwork );
- if( info != 0 )
- goto exit_level_0;
- lwork = LAPACK_C2INT( work_query );
+ if( info != 0 )
+ goto exit_level_0;
+ lwork = LAPACK_C2INT( work_query );
/* Allocate memory for working array(s) */
rwork = (float*)LAPACKE_malloc( sizeof(float) * MAX(1,2*n) );
if( rwork == NULL ) {
if( matrix_layout == LAPACK_COL_MAJOR ) {
/* Call LAPACK function and adjust info */
LAPACK_chbev_2stage( &jobz, &uplo, &n, &kd, ab, &ldab, w, z, &ldz, work, &lwork,
- rwork, &info );
+ rwork, &info );
if( info < 0 ) {
info = info - 1;
}
lapack_int LAPACKE_chetrf_rk_work( int matrix_layout, char uplo, lapack_int n,
lapack_complex_float* a, lapack_int lda,
- lapack_complex_float* e,
+ lapack_complex_float* e,
lapack_int* ipiv, lapack_complex_float* work,
lapack_int lwork )
{
}
#endif
return LAPACKE_chetrs_3_work( matrix_layout, uplo, n, nrhs, a, lda,
- e, ipiv, b, ldb );
+ e, ipiv, b, ldb );
}
lapack_int LAPACKE_csycon_3_work( int matrix_layout, char uplo, lapack_int n,
const lapack_complex_float* a, lapack_int lda,
- const lapack_complex_float* e,
+ const lapack_complex_float* e,
const lapack_int* ipiv, float anorm,
float* rcond, lapack_complex_float* work )
{
}
/* Call middle-level interface */
info = LAPACKE_dbdsvdx_work( matrix_layout, uplo, jobz, range,
- n, d, e, vl, vu, il, iu, ns, s, z,
- ldz, work, iwork);
+ n, d, e, vl, vu, il, iu, ns, s, z,
+ ldz, work, iwork);
/* Backup significant data from working array(s) */
for( i=0; i<12*n-1; i++ ) {
superb[i] = iwork[i+1];
#include "lapacke_utils.h"
lapack_int LAPACKE_dbdsvdx_work( int matrix_layout, char uplo, char jobz, char range,
- lapack_int n, double* d, double* e,
- double vl, double vu,
- lapack_int il, lapack_int iu, lapack_int* ns,
- double* s, double* z, lapack_int ldz,
- double* work, lapack_int* iwork )
+ lapack_int n, double* d, double* e,
+ double vl, double vu,
+ lapack_int il, lapack_int iu, lapack_int* ns,
+ double* s, double* z, lapack_int ldz,
+ double* work, lapack_int* iwork )
{
lapack_int info = 0;
if( matrix_layout == LAPACK_COL_MAJOR ) {
/* Call LAPACK function and adjust info */
LAPACK_dbdsvdx( &uplo, &jobz, &range, &n, d, e, &vl, &vu,
- &il, &iu, ns, s, z, &ldz,
+ &il, &iu, ns, s, z, &ldz,
work, iwork, &info );
if( info < 0 ) {
info = info - 1;
}
/* Call LAPACK function and adjust info */
LAPACK_dbdsvdx( &uplo, &jobz, &range, &n, d, e, &vl, &vu,
- &il, &iu, ns, s, z_t, &ldz_t, work,
- iwork, &info );
+ &il, &iu, ns, s, z_t, &ldz_t, work,
+ iwork, &info );
if( info < 0 ) {
info = info - 1;
}
#endif
/* Query optimal working array(s) size */
info = LAPACKE_dgesvdx_work( matrix_layout, jobu, jobvt, range,
- m, n, a, lda, vl, vu, il, iu, ns, s, u,
+ m, n, a, lda, vl, vu, il, iu, ns, s, u,
ldu, vt, ldvt, &work_query, lwork, iwork );
if( info != 0 ) {
goto exit_level_0;
}
/* Call middle-level interface */
info = LAPACKE_dgesvdx_work( matrix_layout, jobu, jobvt, range,
- m, n, a, lda, vl, vu, il, iu, ns, s, u,
- ldu, vt, ldvt, work, lwork, iwork );
+ m, n, a, lda, vl, vu, il, iu, ns, s, u,
+ ldu, vt, ldvt, work, lwork, iwork );
/* Backup significant data from working array(s) */
for( i=0; i<12*MIN(m,n)-1; i++ ) {
superb[i] = iwork[i+1];
#include "lapacke_utils.h"
lapack_int LAPACKE_dgesvdx_work( int matrix_layout, char jobu, char jobvt, char range,
- lapack_int m, lapack_int n, double* a,
- lapack_int lda, double vl, double vu,
- lapack_int il, lapack_int iu, lapack_int* ns,
- double* s, double* u, lapack_int ldu,
- double* vt, lapack_int ldvt,
- double* work, lapack_int lwork, lapack_int* iwork )
+ lapack_int m, lapack_int n, double* a,
+ lapack_int lda, double vl, double vu,
+ lapack_int il, lapack_int iu, lapack_int* ns,
+ double* s, double* u, lapack_int ldu,
+ double* vt, lapack_int ldvt,
+ double* work, lapack_int lwork, lapack_int* iwork )
{
lapack_int info = 0;
if( matrix_layout == LAPACK_COL_MAJOR ) {
/* Call LAPACK function and adjust info */
LAPACK_dgesvdx( &jobu, &jobvt, &range, &m, &n, a, &lda, &vl, &vu,
- &il, &iu, ns, s, u, &ldu, vt, &ldvt,
+ &il, &iu, ns, s, u, &ldu, vt, &ldvt,
work, &lwork, iwork, &info );
if( info < 0 ) {
info = info - 1;
/* Query optimal working array(s) size if requested */
if( lwork == -1 ) {
LAPACK_dgesvdx( &jobu, &jobvt, &range, &m, &n, a, &lda_t, &vl, &vu,
- &il, &iu, ns, s, u, &ldu_t, vt,
+ &il, &iu, ns, s, u, &ldu_t, vt,
&ldvt_t, work, &lwork, iwork, &info );
return (info < 0) ? (info - 1) : info;
}
LAPACKE_dge_trans( matrix_layout, m, n, a, lda, a_t, lda_t );
/* Call LAPACK function and adjust info */
LAPACK_dgesvdx( &jobu, &jobvt, &range, &m, &n, a_t, &lda_t, &vl, &vu,
- &il, &iu, ns, s, u_t, &ldu_t, vt_t,
- &ldvt_t, work, &lwork, iwork, &info );
+ &il, &iu, ns, s, u_t, &ldu_t, vt_t,
+ &ldvt_t, work, &lwork, iwork, &info );
if( info < 0 ) {
info = info - 1;
}
lapack_int n, const double* a, lapack_int lda )
{
lapack_int info = 0;
- double res = 0.;
+ double res = 0.;
double* work = NULL;
if( matrix_layout != LAPACK_COL_MAJOR && matrix_layout != LAPACK_ROW_MAJOR ) {
LAPACKE_xerbla( "LAPACKE_dlange", -1 );
double* work )
{
lapack_int info = 0;
- double res = 0.;
+ double res = 0.;
if( matrix_layout == LAPACK_COL_MAJOR ) {
/* Call LAPACK function and adjust info */
res = LAPACK_dlange( &norm, &m, &n, a, &lda, work );
const double* a, lapack_int lda )
{
lapack_int info = 0;
- double res = 0.;
+ double res = 0.;
double* work = NULL;
if( matrix_layout != LAPACK_COL_MAJOR && matrix_layout != LAPACK_ROW_MAJOR ) {
LAPACKE_xerbla( "LAPACKE_dlansy", -1 );
double* work )
{
lapack_int info = 0;
- double res = 0.;
+ double res = 0.;
if( matrix_layout == LAPACK_COL_MAJOR ) {
/* Call LAPACK function and adjust info */
res = LAPACK_dlansy( &norm, &uplo, &n, a, &lda, work );
lapack_int lda )
{
lapack_int info = 0;
- double res = 0.;
+ double res = 0.;
double* work = NULL;
if( matrix_layout != LAPACK_COL_MAJOR && matrix_layout != LAPACK_ROW_MAJOR ) {
LAPACKE_xerbla( "LAPACKE_dlantr", -1 );
}
#endif
return LAPACKE_dsytrs_3_work( matrix_layout, uplo, n, nrhs, a, lda,
- e, ipiv, b, ldb );
+ e, ipiv, b, ldb );
}
if( info != 0 ) {
goto exit_level_0;
}
- lwork = (lapack_int)work_query;
+ lwork = (lapack_int)work_query;
/* Allocate memory for work arrays */
work = (double*)
LAPACKE_malloc( sizeof(double) * lwork );
}
/* Call middle-level interface */
info = LAPACKE_sbdsvdx_work( matrix_layout, uplo, jobz, range,
- n, d, e, vl, vu, il, iu, ns, s, z,
- ldz, work, iwork);
+ n, d, e, vl, vu, il, iu, ns, s, z,
+ ldz, work, iwork);
/* Backup significant data from working array(s) */
for( i=0; i<12*n-1; i++ ) {
superb[i] = iwork[i+1];
#endif
/* Query optimal working array(s) size */
info = LAPACKE_sgesvdx_work( matrix_layout, jobu, jobvt, range,
- m, n, a, lda, vl, vu, il, iu, ns, s, u,
+ m, n, a, lda, vl, vu, il, iu, ns, s, u,
ldu, vt, ldvt, &work_query, lwork, iwork );
if( info != 0 ) {
goto exit_level_0;
}
/* Call middle-level interface */
info = LAPACKE_sgesvdx_work( matrix_layout, jobu, jobvt, range,
- m, n, a, lda, vl, vu, il, iu, ns, s, u,
- ldu, vt, ldvt, work, lwork, iwork );
+ m, n, a, lda, vl, vu, il, iu, ns, s, u,
+ ldu, vt, ldvt, work, lwork, iwork );
/* Backup significant data from working array(s) */
for( i=0; i<12*MIN(m,n)-1; i++ ) {
superb[i] = iwork[i+1];
#include "lapacke_utils.h"
lapack_int LAPACKE_sgesvdx_work( int matrix_layout, char jobu, char jobvt, char range,
- lapack_int m, lapack_int n, float* a,
- lapack_int lda, float vl, float vu,
- lapack_int il, lapack_int iu, lapack_int* ns,
- float* s, float* u, lapack_int ldu,
- float* vt, lapack_int ldvt,
- float* work, lapack_int lwork, lapack_int* iwork )
+ lapack_int m, lapack_int n, float* a,
+ lapack_int lda, float vl, float vu,
+ lapack_int il, lapack_int iu, lapack_int* ns,
+ float* s, float* u, lapack_int ldu,
+ float* vt, lapack_int ldvt,
+ float* work, lapack_int lwork, lapack_int* iwork )
{
lapack_int info = 0;
if( matrix_layout == LAPACK_COL_MAJOR ) {
/* Call LAPACK function and adjust info */
LAPACK_sgesvdx( &jobu, &jobvt, &range, &m, &n, a, &lda, &vl, &vu,
- &il, &iu, ns, s, u, &ldu, vt, &ldvt,
+ &il, &iu, ns, s, u, &ldu, vt, &ldvt,
work, &lwork, iwork, &info );
if( info < 0 ) {
info = info - 1;
/* Query optimal working array(s) size if requested */
if( lwork == -1 ) {
LAPACK_sgesvdx( &jobu, &jobvt, &range, &m, &n, a, &lda_t, &vl, &vu,
- &il, &iu, ns, s, u, &ldu_t, vt,
+ &il, &iu, ns, s, u, &ldu_t, vt,
&ldvt_t, work, &lwork, iwork, &info );
return (info < 0) ? (info - 1) : info;
}
LAPACKE_sge_trans( matrix_layout, m, n, a, lda, a_t, lda_t );
/* Call LAPACK function and adjust info */
LAPACK_sgesvdx( &jobu, &jobvt, &range, &m, &n, a_t, &lda_t, &vl, &vu,
- &il, &iu, ns, s, u_t, &ldu_t, vt_t,
- &ldvt_t, work, &lwork, iwork, &info );
+ &il, &iu, ns, s, u_t, &ldu_t, vt_t,
+ &ldvt_t, work, &lwork, iwork, &info );
if( info < 0 ) {
info = info - 1;
}
lapack_int n, const float* a, lapack_int lda )
{
lapack_int info = 0;
- float res = 0.;
+ float res = 0.;
float* work = NULL;
if( matrix_layout != LAPACK_COL_MAJOR && matrix_layout != LAPACK_ROW_MAJOR ) {
LAPACKE_xerbla( "LAPACKE_slange", -1 );
float* work )
{
lapack_int info = 0;
- float res = 0.;
+ float res = 0.;
if( matrix_layout == LAPACK_COL_MAJOR ) {
/* Call LAPACK function and adjust info */
res = LAPACK_slange( &norm, &m, &n, a, &lda, work );
const float* a, lapack_int lda )
{
lapack_int info = 0;
- float res = 0.;
+ float res = 0.;
float* work = NULL;
if( matrix_layout != LAPACK_COL_MAJOR && matrix_layout != LAPACK_ROW_MAJOR ) {
LAPACKE_xerbla( "LAPACKE_slansy", -1 );
float* work )
{
lapack_int info = 0;
- float res = 0.;
+ float res = 0.;
if( matrix_layout == LAPACK_COL_MAJOR ) {
/* Call LAPACK function and adjust info */
res = LAPACK_slansy( &norm, &uplo, &n, a, &lda, work );
if( info != 0 ) {
goto exit_level_0;
}
- lwork = (lapack_int)work_query;
+ lwork = (lapack_int)work_query;
/* Allocate memory for work arrays */
work = (float*)
LAPACKE_malloc( sizeof(float) * lwork );
{
lapack_int info = 0;
lapack_int lwork = (
- // 1.1
- ( LAPACKE_lsame( jobu, 'n' ) && LAPACKE_lsame( jobv, 'n' ) &&
+ // 1.1
+ ( LAPACKE_lsame( jobu, 'n' ) && LAPACKE_lsame( jobv, 'n' ) &&
( LAPACKE_lsame( jobt, 't' ) || LAPACKE_lsame( joba, 'f' ) || LAPACKE_lsame( joba, 'g' ) )) ? 2*n+1 :
//1.2
- ( ( LAPACKE_lsame( jobu, 'n' ) && LAPACKE_lsame( jobv, 'n' ) &&
+ ( ( LAPACKE_lsame( jobu, 'n' ) && LAPACKE_lsame( jobv, 'n' ) &&
!( LAPACKE_lsame( jobt, 't' ) || LAPACKE_lsame( joba, 'f' ) || LAPACKE_lsame( joba, 'g' ) )) ? n*n+3*n :
//2.1
- ( ( ( LAPACKE_lsame( jobv, 'v' ) || LAPACKE_lsame( jobv, 'j' ) ) &&
- (!( LAPACKE_lsame( jobu, 'u') || LAPACKE_lsame( jobu, 'f' ) )&&
+ ( ( ( LAPACKE_lsame( jobv, 'v' ) || LAPACKE_lsame( jobv, 'j' ) ) &&
+ (!( LAPACKE_lsame( jobu, 'u') || LAPACKE_lsame( jobu, 'f' ) )&&
( LAPACKE_lsame( jobt, 't' ) || LAPACKE_lsame( joba, 'f' ) || LAPACKE_lsame( joba, 'g' ) ))? 3*n :
//2.2
- ( ( ( LAPACKE_lsame( jobv, 'v' ) || LAPACKE_lsame( jobv, 'j' ) ) &&
+ ( ( ( LAPACKE_lsame( jobv, 'v' ) || LAPACKE_lsame( jobv, 'j' ) ) &&
!( LAPACKE_lsame( jobu, 'u' ) || LAPACKE_lsame( jobu, 'f' ) ) &&
!( LAPACKE_lsame( jobt, 't' ) || LAPACKE_lsame( joba, 'f' ) || LAPACKE_lsame( joba, 'g' ) ))? 3*n :
( LAPACKE_lsame( jobt, 't' ) || LAPACKE_lsame( joba, 'f' ) || LAPACKE_lsame( joba, 'g' ) ))? 4*n*n:
1) ) ) ) ) ) ) ) );
lapack_int lrwork = (
- // 1.1
- ( LAPACKE_lsame( jobu, 'n' ) && LAPACKE_lsame( jobv, 'n' ) &&
+ // 1.1
+ ( LAPACKE_lsame( jobu, 'n' ) && LAPACKE_lsame( jobv, 'n' ) &&
( LAPACKE_lsame( jobt, 't' ) || LAPACKE_lsame( joba, 'f' ) || LAPACKE_lsame( joba, 'g' ) )) ? MAX(7,n+2*m) :
//1.2
!( LAPACKE_lsame( jobt, 't' ) || LAPACKE_lsame( joba, 'f' ) || LAPACKE_lsame( joba, 'g' ) )) ? MAX(7,2*n) :
//2.1
- ( ( ( LAPACKE_lsame( jobv, 'v' ) || LAPACKE_lsame( jobv, 'j' ) ) &&
- (!( LAPACKE_lsame( jobu, 'u') || LAPACKE_lsame( jobu, 'f' ) ) &&
+ ( ( ( LAPACKE_lsame( jobv, 'v' ) || LAPACKE_lsame( jobv, 'j' ) ) &&
+ (!( LAPACKE_lsame( jobu, 'u') || LAPACKE_lsame( jobu, 'f' ) ) &&
( LAPACKE_lsame( jobt, 't' ) || LAPACKE_lsame( joba, 'f' ) || LAPACKE_lsame( joba, 'g' ) ))? MAX( 7, n+ 2*m ) :
//2.2
- ( ( ( LAPACKE_lsame( jobv, 'v' ) || LAPACKE_lsame( jobv, 'j' ) ) &&
+ ( ( ( LAPACKE_lsame( jobv, 'v' ) || LAPACKE_lsame( jobv, 'j' ) ) &&
!( LAPACKE_lsame( jobu, 'u' ) || LAPACKE_lsame( jobu, 'f' ) ) &&
!( LAPACKE_lsame( jobt, 't' ) || LAPACKE_lsame( joba, 'f' ) || LAPACKE_lsame( joba, 'g' ) ))? MAX(7,2*n) :
#include "lapacke_utils.h"
lapack_int LAPACKE_zgesvdx( int matrix_layout, char jobu, char jobvt, char range,
- lapack_int m, lapack_int n, lapack_complex_double* a,
- lapack_int lda, double vl, double vu,
- lapack_int il, lapack_int iu, lapack_int* ns,
- double* s, lapack_complex_double* u, lapack_int ldu,
- lapack_complex_double* vt, lapack_int ldvt,
- lapack_int* superb )
+ lapack_int m, lapack_int n, lapack_complex_double* a,
+ lapack_int lda, double vl, double vu,
+ lapack_int il, lapack_int iu, lapack_int* ns,
+ double* s, lapack_complex_double* u, lapack_int ldu,
+ lapack_complex_double* vt, lapack_int ldvt,
+ lapack_int* superb )
{
lapack_int info = 0;
lapack_int lwork = -1;
#endif
/* Query optimal working array(s) size */
info = LAPACKE_zgesvdx_work( matrix_layout, jobu, jobvt, range,
- m, n, a, lda, vl, vu, il, iu, ns, s, u,
+ m, n, a, lda, vl, vu, il, iu, ns, s, u,
ldu, vt, ldvt, &work_query, lwork, rwork, iwork );
if( info != 0 ) {
goto exit_level_0;
lwork = LAPACK_Z2INT (work_query);
/* Allocate memory for work arrays */
work = (lapack_complex_double*)
- LAPACKE_malloc( sizeof(lapack_complex_double) * lwork );
+ LAPACKE_malloc( sizeof(lapack_complex_double) * lwork );
if( work == NULL ) {
info = LAPACK_WORK_MEMORY_ERROR;
goto exit_level_1;
}
/* Call middle-level interface */
info = LAPACKE_zgesvdx_work( matrix_layout, jobu, jobvt, range,
- m, n, a, lda, vl, vu, il, iu, ns, s, u,
- ldu, vt, ldvt, work, lwork, rwork, iwork );
+ m, n, a, lda, vl, vu, il, iu, ns, s, u,
+ ldu, vt, ldvt, work, lwork, rwork, iwork );
/* Backup significant data from working array(s) */
for( i=0; i<12*MIN(m,n)-1; i++ ) {
superb[i] = iwork[i+1];
#include "lapacke_utils.h"
lapack_int LAPACKE_zgesvdx_work( int matrix_layout, char jobu, char jobvt, char range,
- lapack_int m, lapack_int n, lapack_complex_double* a,
- lapack_int lda, double vl, double vu,
- lapack_int il, lapack_int iu, lapack_int* ns,
- double* s, lapack_complex_double* u, lapack_int ldu,
- lapack_complex_double* vt, lapack_int ldvt,
- lapack_complex_double* work, lapack_int lwork,
- double* rwork, lapack_int* iwork )
+ lapack_int m, lapack_int n, lapack_complex_double* a,
+ lapack_int lda, double vl, double vu,
+ lapack_int il, lapack_int iu, lapack_int* ns,
+ double* s, lapack_complex_double* u, lapack_int ldu,
+ lapack_complex_double* vt, lapack_int ldvt,
+ lapack_complex_double* work, lapack_int lwork,
+ double* rwork, lapack_int* iwork )
{
lapack_int info = 0;
if( matrix_layout == LAPACK_COL_MAJOR ) {
/* Call LAPACK function and adjust info */
LAPACK_zgesvdx( &jobu, &jobvt, &range, &m, &n, a, &lda, &vl, &vu,
- &il, &iu, ns, s, u, &ldu, vt, &ldvt,
+ &il, &iu, ns, s, u, &ldu, vt, &ldvt,
work, &lwork, rwork, iwork, &info );
if( info < 0 ) {
info = info - 1;
/* Query optimal working array(s) size if requested */
if( lwork == -1 ) {
LAPACK_zgesvdx( &jobu, &jobvt, &range, &m, &n, a, &lda_t, &vl, &vu,
- &il, &iu, ns, s, u, &ldu_t, vt,
+ &il, &iu, ns, s, u, &ldu_t, vt,
&ldvt_t, work, &lwork, rwork, iwork, &info );
return (info < 0) ? (info - 1) : info;
}
LAPACKE_zge_trans( matrix_layout, m, n, a, lda, a_t, lda_t );
/* Call LAPACK function and adjust info */
LAPACK_zgesvdx( &jobu, &jobvt, &range, &m, &n, a_t, &lda_t, &vl, &vu,
- &il, &iu, ns, s, u_t, &ldu_t, vt_t,
+ &il, &iu, ns, s, u_t, &ldu_t, vt_t,
&ldvt_t, work, &lwork, rwork, iwork, &info );
if( info < 0 ) {
info = info - 1;
}
if( LAPACKE_lsame( jobv, 'a' ) || LAPACKE_lsame( jobv, 'v' ) ) {
v_t = (lapack_complex_double*)
- LAPACKE_malloc( sizeof(lapack_complex_double) * ldv_t * MAX(1,n) );
+ LAPACKE_malloc( sizeof(lapack_complex_double) * ldv_t * MAX(1,n) );
if( v_t == NULL ) {
info = LAPACK_TRANSPOSE_MEMORY_ERROR;
goto exit_level_1;
lapack_int LAPACKE_zhetrf_rk_work( int matrix_layout, char uplo, lapack_int n,
lapack_complex_double* a, lapack_int lda,
- lapack_complex_double* e,
+ lapack_complex_double* e,
lapack_int* ipiv, lapack_complex_double* work,
lapack_int lwork )
{
}
#endif
return LAPACKE_zhetrs_3_work( matrix_layout, uplo, n, nrhs, a, lda,
- e, ipiv, b, ldb );
+ e, ipiv, b, ldb );
}
lapack_int lda )
{
lapack_int info = 0;
- double res = 0.;
+ double res = 0.;
double* work = NULL;
if( matrix_layout != LAPACK_COL_MAJOR && matrix_layout != LAPACK_ROW_MAJOR ) {
LAPACKE_xerbla( "LAPACKE_zlange", -1 );
lapack_int lda, double* work )
{
lapack_int info = 0;
- double res = 0.;
+ double res = 0.;
if( matrix_layout == LAPACK_COL_MAJOR ) {
/* Call LAPACK function and adjust info */
res = LAPACK_zlange( &norm, &m, &n, a, &lda, work );
const lapack_complex_double* a, lapack_int lda )
{
lapack_int info = 0;
- double res = 0.;
+ double res = 0.;
double* work = NULL;
if( matrix_layout != LAPACK_COL_MAJOR && matrix_layout != LAPACK_ROW_MAJOR ) {
LAPACKE_xerbla( "LAPACKE_zlanhe", -1 );
lapack_int lda, double* work )
{
lapack_int info = 0;
- double res = 0.;
+ double res = 0.;
if( matrix_layout == LAPACK_COL_MAJOR ) {
/* Call LAPACK function and adjust info */
res = LAPACK_zlanhe( &norm, &uplo, &n, a, &lda, work );
const lapack_complex_double* a, lapack_int lda )
{
lapack_int info = 0;
- double res = 0.;
+ double res = 0.;
double* work = NULL;
if( matrix_layout != LAPACK_COL_MAJOR && matrix_layout != LAPACK_ROW_MAJOR ) {
LAPACKE_xerbla( "LAPACKE_zlansy", -1 );
lapack_int lda, double* work )
{
lapack_int info = 0;
- double res = 0.;
+ double res = 0.;
if( matrix_layout == LAPACK_COL_MAJOR ) {
/* Call LAPACK function and adjust info */
res = LAPACK_zlansy( &norm, &uplo, &n, a, &lda, work );
const lapack_complex_double* a, lapack_int lda )
{
lapack_int info = 0;
- double res = 0.;
+ double res = 0.;
double* work = NULL;
if( matrix_layout != LAPACK_COL_MAJOR && matrix_layout != LAPACK_ROW_MAJOR ) {
LAPACKE_xerbla( "LAPACKE_zlantr", -1 );
}
#endif
return LAPACKE_zsytrs_3_work( matrix_layout, uplo, n, nrhs, a, lda,
- e, ipiv, b, ldb );
+ e, ipiv, b, ldb );
}