274 int irow, jcol, i, j,
LDA;
285 A = (
float* )malloc( (*n+*k)*
LDA*
sizeof( float ) );
287 for( i=0; i<*k; i++ ){
290 for( j=jcol; j<*n; j++ )
291 A[
LDA*(j-jcol)+irow ]=a[ (*lda)*j+i ];
295 for( j=0; j<*n; j++ )
296 A[
LDA*j+irow ]=a[ (*lda)*j+i ];
301 for( j=0; j<*n; j++ )
302 A[
LDA*j+irow ]=a[ (*lda)*j+i ];
303 for( i=1; i<*k+1; i++ ){
306 for( j=jcol; j<(*n+*k); j++ )
307 A[
LDA*j+irow ]=a[ (*lda)*(j-jcol)+i ];
void cblas_stbsv(CBLAS_LAYOUT layout, CBLAS_UPLO Uplo, CBLAS_TRANSPOSE TransA, CBLAS_DIAG Diag, const CBLAS_INT N, const CBLAS_INT K, const float *A, const CBLAS_INT lda, float *X, const CBLAS_INT incX)
void get_diag_type(char *type, CBLAS_DIAG *diag)
void get_uplo_type(char *type, CBLAS_UPLO *uplo)
void get_transpose_type(char *type, CBLAS_TRANSPOSE *trans)