lapacke_stb_trans.c File Reference

#include "lapacke_utils.h"

Go to the source code of this file.

Functions
void	LAPACKE_stb_trans (int matrix_layout, char uplo, char diag, lapack_int n, lapack_int kd, const float *in, lapack_int ldin, float *out, lapack_int ldout)

Function Documentation

LAPACKE_stb_trans()

void LAPACKE_stb_trans	(	int	matrix_layout,
		char	uplo,
		char	diag,
		lapack_int	n,
		lapack_int	kd,
		const float *	in,
		lapack_int	ldin,
		float *	out,
		lapack_int	ldout )

Definition at line 39 of file lapacke_stb_trans.c.

{
    lapack_logical colmaj, upper, unit;
 
    if( in == NULL || out == NULL ) return;
 
    colmaj = ( matrix_layout == LAPACK_COL_MAJOR );
    upper  = LAPACKE_lsame( uplo, 'u' );
    unit   = LAPACKE_lsame( diag, 'u' );
 
    if( ( !colmaj && ( matrix_layout != LAPACK_ROW_MAJOR ) ) ||
        ( !upper  && !LAPACKE_lsame( uplo, 'l' ) ) ||
        ( !unit   && !LAPACKE_lsame( diag, 'n' ) ) ) {
        /* Just exit if any of input parameters are wrong */
        return;
    }
 
    if( unit ) {
        /* Unit case, diagonal excluded from transposition */
        if( colmaj ) {
            if( upper ) {
                LAPACKE_sgb_trans( matrix_layout, n-1, n-1, 0, kd-1,
                                   &in[ldin], ldin, &out[1], ldout );
            } else {
                LAPACKE_sgb_trans( matrix_layout, n-1, n-1, kd-1, 0,
                                   &in[1], ldin, &out[ldout], ldout );
            }
        } else {
            if( upper ) {
                LAPACKE_sgb_trans( matrix_layout, n-1, n-1, 0, kd-1,
                                   &in[1], ldin, &out[ldout], ldout );
            } else {
                LAPACKE_sgb_trans( matrix_layout, n-1, n-1, kd-1, 0,
                                   &in[ldin], ldin, &out[1], ldout );
            }
        }
    } else {
        /* Non-unit case */
        if( upper ) {
            LAPACKE_sgb_trans( matrix_layout, n, n, 0, kd, in, ldin, out,
                               ldout );
        } else {
            LAPACKE_sgb_trans( matrix_layout, n, n, kd, 0, in, ldin, out,
                               ldout );
        }
    }
}