summaryrefslogtreecommitdiff
path: root/LAPACKE
diff options
context:
space:
mode:
authoreugene.chereshnev <eugenechereshnev@gmail.com>2017-01-31 18:40:37 -0800
committereugene.chereshnev <eugenechereshnev@gmail.com>2017-02-03 18:47:01 -0800
commit9da5c0b890bfa58da8dffec251fcd5361e0f73ab (patch)
treed1b2169bf44d4028a24d3ebd0fcbe01d531dab98 /LAPACKE
parent4f3d96e44fba63295680ed43679578409bb5aeea (diff)
downloadlapack-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.c2
-rw-r--r--LAPACKE/src/lapacke_sbdsvdx_work.c18
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;
}