47 int i,j,irow,jcol,
LDA;
54 for( i=0; i<*ku; i++ ){
57 for( j=jcol; j<*n; j++ ){
58 A[ LDA*(j-jcol)+irow ].real=a[ (*lda)*j+i ].
real;
59 A[ LDA*(j-jcol)+irow ].imag=a[ (*lda)*j+i ].
imag;
64 for( j=0; j<*n; j++ ){
65 A[ LDA*j+irow ].
real=a[ (*lda)*j+i ].
real;
66 A[ LDA*j+irow ].
imag=a[ (*lda)*j+i ].
imag;
68 for( i=*ku+1; i<*ku+*kl+1; i++ ){
71 for( j=jcol; j<(*n+*kl); j++ ){
72 A[ LDA*j+irow ].
real=a[ (*lda)*(j-jcol)+i ].real;
73 A[ LDA*j+irow ].
imag=a[ (*lda)*(j-jcol)+i ].imag;
76 cblas_zgbmv(
CblasRowMajor, trans, *m, *n, *kl, *ku, alpha, A, LDA, x,
77 *incx, beta, y, *incy );
81 cblas_zgbmv(
CblasColMajor, trans, *m, *n, *kl, *ku, alpha, a, *lda, x,
82 *incx, beta, y, *incy );
84 cblas_zgbmv(
UNDEFINED, trans, *m, *n, *kl, *ku, alpha, a, *lda, x,
85 *incx, beta, y, *incy );
void get_transpose_type(char *type, CBLAS_TRANSPOSE *trans)
void cblas_zgbmv(CBLAS_LAYOUT layout, CBLAS_TRANSPOSE TransA, const int M, const int N, const int KL, const int KU, const void *alpha, const void *A, const int lda, const void *X, const int incX, const void *beta, void *Y, const int incY)