228 int irow, jcol, i, j,
LDA;
239 A = (
double* )malloc( (*n+*k)*LDA*
sizeof( double ) );
241 for( i=0; i<*k; i++ ){
244 for( j=jcol; j<*n; j++ )
245 A[ LDA*(j-jcol)+irow ]=a[ (*lda)*j+i ];
249 for( j=0; j<*n; j++ )
250 A[ LDA*j+irow ]=a[ (*lda)*j+i ];
255 for( j=0; j<*n; j++ )
256 A[ LDA*j+irow ]=a[ (*lda)*j+i ];
257 for( i=1; i<*k+1; i++ ){
260 for( j=jcol; j<(*n+*k); j++ )
261 A[ LDA*j+irow ]=a[ (*lda)*(j-jcol)+i ];
void get_transpose_type(char *type, CBLAS_TRANSPOSE *trans)
void cblas_dtbmv(CBLAS_LAYOUT layout, CBLAS_UPLO Uplo, CBLAS_TRANSPOSE TransA, CBLAS_DIAG Diag, const int N, const int K, const double *A, const int lda, double *X, const int incX)
void get_diag_type(char *type, CBLAS_DIAG *diag)
void get_uplo_type(char *type, CBLAS_UPLO *uplo)