lapacke_ctf_trans.c File Reference

#include "lapacke_utils.h"

Go to the source code of this file.

Functions
void	LAPACKE_ctf_trans (int matrix_layout, char transr, char uplo, char diag, lapack_int n, const lapack_complex_float *in, lapack_complex_float *out)

Function Documentation

LAPACKE_ctf_trans()

void LAPACKE_ctf_trans	(	int	matrix_layout,
		char	transr,
		char	uplo,
		char	diag,
		lapack_int	n,
		const lapack_complex_float *	in,
		lapack_complex_float *	out )

Definition at line 40 of file lapacke_ctf_trans.c.

{
    lapack_int row, col;
    lapack_logical rowmaj, ntr, lower, unit;
 
    if( in == NULL || out == NULL ) return ;
 
    rowmaj = (matrix_layout == LAPACK_ROW_MAJOR);
    ntr    = LAPACKE_lsame( transr, 'n' );
    lower  = LAPACKE_lsame( uplo,   'l' );
    unit   = LAPACKE_lsame( diag,   'u' );
 
    if( ( !rowmaj && ( matrix_layout != LAPACK_COL_MAJOR ) ) ||
        ( !ntr    && !LAPACKE_lsame( transr, 't' ) &&
                     !LAPACKE_lsame( transr, 'c' ) ) ||
        ( !lower  && !LAPACKE_lsame( uplo,   'u' ) ) ||
        ( !unit   && !LAPACKE_lsame( diag,   'n' ) ) ) {
        /* Just exit if input parameters are wrong */
        return;
    }
 
    /* Determine parameters of array representing RFP */
    if( ntr ) {
        if( n%2 == 0 ) {
            row = n + 1;
            col = n / 2;
        } else {
            row = n;
            col = (n + 1) / 2;
        }
    } else {
        if( n%2 == 0 ) {
            row = n / 2;
            col = n + 1;
        } else {
            row = (n + 1) / 2;
            col = n;
        }
    }
 
    /* Perform conversion: */
    if( rowmaj ) {
        LAPACKE_cge_trans( LAPACK_ROW_MAJOR, row, col, in, col, out, row );
    } else {
        LAPACKE_cge_trans( LAPACK_COL_MAJOR, row, col, in, row, out, col );
    }
}