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

Definition at line 40 of file lapacke_cgb_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];
            }
        }
    }
}

Here is the caller graph for this function: