lapacke_zlantr_work.c File Reference

#include "lapacke_utils.h"

Go to the source code of this file.

Functions
double	LAPACKE_zlantr_work (int matrix_layout, char norm, char uplo, char diag, lapack_int m, lapack_int n, const lapack_complex_double *a, lapack_int lda, double *work)

Function Documentation

LAPACKE_zlantr_work()

double LAPACKE_zlantr_work	(	int	matrix_layout,
		char	norm,
		char	uplo,
		char	diag,
		lapack_int	m,
		lapack_int	n,
		const lapack_complex_double *	a,
		lapack_int	lda,
		double *	work )

Definition at line 35 of file lapacke_zlantr_work.c.

{
    lapack_int info = 0;
    double res = 0.;
    if( matrix_layout == LAPACK_COL_MAJOR ) {
        /* Call LAPACK function */
        res = LAPACK_zlantr( &norm, &uplo, &diag, &m, &n, a, &lda, work );
    } else if( matrix_layout == LAPACK_ROW_MAJOR ) {
        double* work_lapack = NULL;
        char norm_lapack;
        char uplo_lapack;
        /* Check leading dimension(s) */
        if( lda < n ) {
            info = -8;
            LAPACKE_xerbla( "LAPACKE_zlantr_work", info );
            return info;
        }
        if( LAPACKE_lsame( norm, '1' ) || LAPACKE_lsame( norm, 'o' ) ) {
            norm_lapack = 'i';
        } else if( LAPACKE_lsame( norm, 'i' ) ) {
            norm_lapack = '1';
        } else {
            norm_lapack = norm;
        }
        if( LAPACKE_lsame( uplo, 'u' ) ) {
            uplo_lapack = 'l';
        } else {
            uplo_lapack = 'u';
        }
        /* Allocate memory for work array(s) */
        if( LAPACKE_lsame( norm_lapack, 'i' ) ) {
            work_lapack = (double*)LAPACKE_malloc( sizeof(double) * MAX(1,n) );
            if( work_lapack == NULL ) {
                info = LAPACK_WORK_MEMORY_ERROR;
                goto exit_level_0;
            }
        }
        /* Call LAPACK function */
        res = LAPACK_zlantr( &norm_lapack, &uplo_lapack, &diag, &n, &m, a, &lda, work_lapack );
        /* Release memory and exit */
        if( work_lapack ) {
            LAPACKE_free( work_lapack );
        }
exit_level_0:
        if( info == LAPACK_WORK_MEMORY_ERROR ) {
            LAPACKE_xerbla( "LAPACKE_zlantr_work", info );
        }
    } else {
        info = -1;
        LAPACKE_xerbla( "LAPACKE_zlantr_work", info );
    }
    return res;
}