360 {
361
367
371
375 *incx);
376 else {
377 LDA = *k+2;
380 for( i=0; i<*k; i++ ){
381 irow=*k-i;
382 jcol=(*k)-i;
383 for( j=jcol; j<*n; j++ ) {
384 A[ LDA*(j-jcol)+irow ].real=a[ (*lda)*j+i ].
real;
385 A[ LDA*(j-jcol)+irow ].imag=a[ (*lda)*j+i ].
imag;
386 }
387 }
388 i=*k;
389 irow=*k-i;
390 for( j=0; j<*n; j++ ) {
391 A[ LDA*j+irow ].
real=a[ (*lda)*j+i ].
real;
392 A[ LDA*j+irow ].
imag=a[ (*lda)*j+i ].
imag;
393 }
394 }
395 else {
396 i=0;
397 irow=*k-i;
398 for( j=0; j<*n; j++ ) {
399 A[ LDA*j+irow ].
real=a[ (*lda)*j+i ].
real;
400 A[ LDA*j+irow ].
imag=a[ (*lda)*j+i ].
imag;
401 }
402 for( i=1; i<*k+1; i++ ){
403 irow=*k-i;
404 jcol=i;
405 for( j=jcol; j<(*n+*k); j++ ) {
406 A[ LDA*j+irow ].
real=a[ (*lda)*(j-jcol)+i ].real;
407 A[ LDA*j+irow ].
imag=a[ (*lda)*(j-jcol)+i ].imag;
408 }
409 }
410 }
412 x, *incx);
413 free(A);
414 }
415 }
418 else
420}
void cblas_ctbsv(CBLAS_LAYOUT layout, CBLAS_UPLO Uplo, CBLAS_TRANSPOSE TransA, CBLAS_DIAG Diag, const CBLAS_INT N, const CBLAS_INT K, const void *A, const CBLAS_INT lda, void *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)