321 int i,j,irow,jcol,
LDA;
328 A = (
float* )malloc( (*n+*k)*LDA*
sizeof( float ) );
330 for( i=0; i<*k; i++ ){
333 for( j=jcol; j<*n; j++ )
334 A[ LDA*(j-jcol)+irow ]=a[ (*lda)*j+i ];
338 for( j=0; j<*n; j++ )
339 A[ LDA*j+irow ]=a[ (*lda)*j+i ];
344 for( j=0; j<*n; j++ )
345 A[ LDA*j+irow ]=a[ (*lda)*j+i ];
346 for( i=1; i<*k+1; i++ ){
349 for( j=jcol; j<(*n+*k); j++ )
350 A[ LDA*j+irow ]=a[ (*lda)*(j-jcol)+i ];
void cblas_ssbmv(CBLAS_LAYOUT layout, CBLAS_UPLO Uplo, const int N, const int K, const float alpha, const float *A, const int lda, const float *X, const int incX, const float beta, float *Y, const int incY)
void get_uplo_type(char *type, CBLAS_UPLO *uplo)