192 int i,irow,j,jcol,
LDA;
199 A = (
double* )malloc( (*n+*kl)*LDA*
sizeof( double ) );
200 for( i=0; i<*ku; i++ ){
203 for( j=jcol; j<*n; j++ )
204 A[ LDA*(j-jcol)+irow ]=a[ (*lda)*j+i ];
208 for( j=0; j<*n; j++ )
209 A[ LDA*j+irow ]=a[ (*lda)*j+i ];
210 for( i=*ku+1; i<*ku+*kl+1; i++ ){
213 for( j=jcol; j<(*n+*kl); j++ )
214 A[ LDA*j+irow ]=a[ (*lda)*(j-jcol)+i ];
217 A, LDA, x, *incx, *beta, y, *incy );
222 a, *lda, x, *incx, *beta, y, *incy );
void get_transpose_type(char *type, CBLAS_TRANSPOSE *trans)
void cblas_dgbmv(CBLAS_LAYOUT layout, CBLAS_TRANSPOSE TransA, const int M, const int N, const int KL, const int KU, const double alpha, const double *A, const int lda, const double *X, const int incX, const double beta, double *Y, const int incY)