diff options
author | eugene.chereshnev <eugenechereshnev@gmail.com> | 2017-01-31 18:40:37 -0800 |
---|---|---|
committer | eugene.chereshnev <eugenechereshnev@gmail.com> | 2017-02-03 18:47:01 -0800 |
commit | 9da5c0b890bfa58da8dffec251fcd5361e0f73ab (patch) | |
tree | d1b2169bf44d4028a24d3ebd0fcbe01d531dab98 /LAPACKE | |
parent | 4f3d96e44fba63295680ed43679578409bb5aeea (diff) | |
download | lapack-9da5c0b890bfa58da8dffec251fcd5361e0f73ab.tar.gz lapack-9da5c0b890bfa58da8dffec251fcd5361e0f73ab.tar.bz2 lapack-9da5c0b890bfa58da8dffec251fcd5361e0f73ab.zip |
lapacke_*bdsvdx: corrected # of columns used for allocation of z_t
Diffstat (limited to 'LAPACKE')
-rw-r--r-- | LAPACKE/src/lapacke_dbdsvdx_work.c | 2 | ||||
-rw-r--r-- | LAPACKE/src/lapacke_sbdsvdx_work.c | 18 |
2 files changed, 10 insertions, 10 deletions
diff --git a/LAPACKE/src/lapacke_dbdsvdx_work.c b/LAPACKE/src/lapacke_dbdsvdx_work.c index a795110f..90a520fe 100644 --- a/LAPACKE/src/lapacke_dbdsvdx_work.c +++ b/LAPACKE/src/lapacke_dbdsvdx_work.c @@ -64,7 +64,7 @@ lapack_int LAPACKE_dbdsvdx_work( int matrix_layout, char uplo, char jobz, char r /* Allocate memory for temporary array(s) */ if( LAPACKE_lsame( jobz, 'v' ) ) { z_t = (double*) - LAPACKE_malloc( sizeof(double) * ldz_t * MAX(2*n,1) ); + LAPACKE_malloc( sizeof(double) * ldz_t * MAX(ncols_z,1) ); if( z_t == NULL ) { info = LAPACK_TRANSPOSE_MEMORY_ERROR; goto exit_level_0; diff --git a/LAPACKE/src/lapacke_sbdsvdx_work.c b/LAPACKE/src/lapacke_sbdsvdx_work.c index f632b3ea..69d970b7 100644 --- a/LAPACKE/src/lapacke_sbdsvdx_work.c +++ b/LAPACKE/src/lapacke_sbdsvdx_work.c @@ -34,17 +34,17 @@ #include "lapacke_utils.h" lapack_int LAPACKE_sbdsvdx_work( int matrix_layout, char uplo, char jobz, char range, - lapack_int n, float* d, float* e, - float vl, float vu, - lapack_int il, lapack_int iu, lapack_int* ns, - float* s, float* z, lapack_int ldz, - float* work, lapack_int* iwork ) + lapack_int n, float* d, float* e, + float vl, float vu, + lapack_int il, lapack_int iu, lapack_int* ns, + float* s, float* z, lapack_int ldz, + float* work, lapack_int* iwork ) { lapack_int info = 0; if( matrix_layout == LAPACK_COL_MAJOR ) { /* Call LAPACK function and adjust info */ LAPACK_sbdsvdx( &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; @@ -64,7 +64,7 @@ lapack_int LAPACKE_sbdsvdx_work( int matrix_layout, char uplo, char jobz, char r /* Allocate memory for temporary array(s) */ if( LAPACKE_lsame( jobz, 'v' ) ) { z_t = (float*) - LAPACKE_malloc( sizeof(float) * ldz_t * MAX(2*n,1) ); + LAPACKE_malloc( sizeof(float) * ldz_t * MAX(ncols_z,1) ); if( z_t == NULL ) { info = LAPACK_TRANSPOSE_MEMORY_ERROR; goto exit_level_0; @@ -72,8 +72,8 @@ lapack_int LAPACKE_sbdsvdx_work( int matrix_layout, char uplo, char jobz, char r } /* Call LAPACK function and adjust info */ LAPACK_sbdsvdx( &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; } |