Delete unnecessary work buffer usage of *ktev and *kteqr in lapacke

Author: sh-zheng

LAPACKE/src/lapacke_dkteqr.c

@@ -56,23 +56,8 @@ lapack_int API_SUFFIX(LAPACKE_dkteqr)( int matrix_layout, char compz, lapack_int
         }
     }
 #endif
-    /* Additional scalars initializations for work arrays */
-    if( API_SUFFIX(LAPACKE_lsame)( compz, 'n' ) ) {
-        lwork = 1;
-    } else {
-        lwork = MAX(1,2*n-2);
-    }
-    /* Allocate memory for working array(s) */
-    work = (double*)LAPACKE_malloc( sizeof(double) * lwork );
-    if( work == NULL ) {
-        info = LAPACK_WORK_MEMORY_ERROR;
-        goto exit_level_0;
-    }
     /* Call middle-level interface */
-    info = API_SUFFIX(LAPACKE_dkteqr_work)( matrix_layout, compz, n, e, z, ldz, work );
-    /* Release memory and exit */
-    LAPACKE_free( work );
-exit_level_0:
+    info = API_SUFFIX(LAPACKE_dkteqr_work)( matrix_layout, compz, n, e, z, ldz, NULL );
     if( info == LAPACK_WORK_MEMORY_ERROR ) {
         API_SUFFIX(LAPACKE_xerbla)( "LAPACKE_dkteqr", info );
     }

LAPACKE/src/lapacke_dktev.c

@@ -52,21 +52,8 @@ lapack_int API_SUFFIX(LAPACKE_dktev)( int matrix_layout, char jobz, lapack_int n
         }
     }
 #endif
-    /* Allocate memory for working array(s) */
-    if( API_SUFFIX(LAPACKE_lsame)( jobz, 'v' ) ) {
-        work = (double*)LAPACKE_malloc( sizeof(double) * MAX(1,2*n-2) );
-        if( work == NULL ) {
-            info = LAPACK_WORK_MEMORY_ERROR;
-            goto exit_level_0;
-        }
-    }
     /* Call middle-level interface */
-    info = API_SUFFIX(LAPACKE_dktev_work)( matrix_layout, jobz, n, d, e, z, ldz, work );
-    /* Release memory and exit */
-    if( API_SUFFIX(LAPACKE_lsame)( jobz, 'v' ) ) {
-        LAPACKE_free( work );
-    }
-exit_level_0:
+    info = API_SUFFIX(LAPACKE_dktev_work)( matrix_layout, jobz, n, d, e, z, ldz, NULL );
     if( info == LAPACK_WORK_MEMORY_ERROR ) {
         API_SUFFIX(LAPACKE_xerbla)( "LAPACKE_dktev", info );
     }

LAPACKE/src/lapacke_skteqr.c

@@ -56,23 +56,8 @@ lapack_int API_SUFFIX(LAPACKE_skteqr)( int matrix_layout, char compz, lapack_int
         }
     }
 #endif
-    /* Additional scalars initializations for work arrays */
-    if( API_SUFFIX(LAPACKE_lsame)( compz, 'n' ) ) {
-        lwork = 1;
-    } else {
-        lwork = MAX(1,2*n-2);
-    }
-    /* Allocate memory for working array(s) */
-    work = (float*)LAPACKE_malloc( sizeof(float) * lwork );
-    if( work == NULL ) {
-        info = LAPACK_WORK_MEMORY_ERROR;
-        goto exit_level_0;
-    }
     /* Call middle-level interface */
-    info = API_SUFFIX(LAPACKE_skteqr_work)( matrix_layout, compz, n, e, z, ldz, work );
-    /* Release memory and exit */
-    LAPACKE_free( work );
-exit_level_0:
+    info = API_SUFFIX(LAPACKE_skteqr_work)( matrix_layout, compz, n, e, z, ldz, NULL );
     if( info == LAPACK_WORK_MEMORY_ERROR ) {
         API_SUFFIX(LAPACKE_xerbla)( "LAPACKE_skteqr", info );
     }

LAPACKE/src/lapacke_sktev.c

@@ -52,21 +52,8 @@ lapack_int API_SUFFIX(LAPACKE_sktev)( int matrix_layout, char jobz, lapack_int n
         }
     }
 #endif
-    /* Allocate memory for working array(s) */
-    if( API_SUFFIX(LAPACKE_lsame)( jobz, 'v' ) ) {
-        work = (float*)LAPACKE_malloc( sizeof(float) * MAX(1,2*n-2) );
-        if( work == NULL ) {
-            info = LAPACK_WORK_MEMORY_ERROR;
-            goto exit_level_0;
-        }
-    }
     /* Call middle-level interface */
-    info = API_SUFFIX(LAPACKE_sktev_work)( matrix_layout, jobz, n, d, e, z, ldz, work );
-    /* Release memory and exit */
-    if( API_SUFFIX(LAPACKE_lsame)( jobz, 'v' ) ) {
-        LAPACKE_free( work );
-    }
-exit_level_0:
+    info = API_SUFFIX(LAPACKE_sktev_work)( matrix_layout, jobz, n, d, e, z, ldz, NULL );
     if( info == LAPACK_WORK_MEMORY_ERROR ) {
         API_SUFFIX(LAPACKE_xerbla)( "LAPACKE_sktev", info );
     }

SRC/dkygv.f

@@ -138,8 +138,8 @@
 *> \param[in] LWORK
 *> \verbatim
 *>          LWORK is INTEGER
-*>          The length of the array WORK.  LWORK >= max(1,3*N-1).
-*>          For optimal efficiency, LWORK >= (NB+2)*N,
+*>          The length of the array WORK.  LWORK >= max(1,2*N-1).
+*>          For optimal efficiency, LWORK >= (NB+1)*N,
 *>          where NB is the blocksize for SSYTRD returned by ILAENV.
 *>
 *>          If LWORK = -1, then a workspace query is assumed; the routine

SRC/skygv.f

@@ -138,8 +138,8 @@
 *> \param[in] LWORK
 *> \verbatim
 *>          LWORK is INTEGER
-*>          The length of the array WORK.  LWORK >= max(1,3*N-1).
-*>          For optimal efficiency, LWORK >= (NB+2)*N,
+*>          The length of the array WORK.  LWORK >= max(1,2*N-1).
+*>          For optimal efficiency, LWORK >= (NB+1)*N,
 *>          where NB is the blocksize for SSYTRD returned by ILAENV.
 *>
 *>          If LWORK = -1, then a workspace query is assumed; the routine