lapacke_zgb_trans.c File Reference

#include "lapacke_utils.h"

Go to the source code of this file.

Functions
void	LAPACKE_zgb_trans (int matrix_layout, lapack_int m, lapack_int n, lapack_int kl, lapack_int ku, const lapack_complex_double *in, lapack_int ldin, lapack_complex_double *out, lapack_int ldout)

Function Documentation

LAPACKE_zgb_trans()

void LAPACKE_zgb_trans	(	int	matrix_layout,
		lapack_int	m,
		lapack_int	n,
		lapack_int	kl,
		lapack_int	ku,
		const lapack_complex_double *	in,
		lapack_int	ldin,
		lapack_complex_double *	out,
		lapack_int	ldout )

Definition at line 39 of file lapacke_zgb_trans.c.

{
    lapack_int i, j;
 
    if( in == NULL || out == NULL ) return;
 
    if( matrix_layout == LAPACK_COL_MAJOR ) {
        for( j = 0; j < MIN( ldout, n ); j++ ) {
            for( i = MAX( ku-j, 0 ); i < MIN3( ldin, m+ku-j, kl+ku+1 );
                 i++ ) {
                out[(size_t)i*ldout+j] = in[i+(size_t)j*ldin];
            }
        }
    } else if ( matrix_layout == LAPACK_ROW_MAJOR ) {
        /* TODO: interchange loops for performance.
         * This is just reference impemeltation
         */
        for( j = 0; j < MIN( n, ldin ); j++ ) {
            for( i = MAX( ku-j, 0 ); i < MIN3( ldout, m+ku-j, kl+ku+1 );
                 i++ ) {
                out[i+(size_t)j*ldout] = in[(size_t)i*ldin+j];
            }
        }
    }
}