321 {
327
331
335 x, *incx);
336 else {
337 LDA = *k+2;
340 for( i=0; i<*k; i++ ){
341 irow=*k-i;
342 jcol=(*k)-i;
343 for( j=jcol; j<*n; j++ ) {
344 A[ LDA*(j-jcol)+irow ].real=a[ (*lda)*j+i ].
real;
345 A[ LDA*(j-jcol)+irow ].imag=a[ (*lda)*j+i ].
imag;
346 }
347 }
348 i=*k;
349 irow=*k-i;
350 for( j=0; j<*n; j++ ) {
351 A[ LDA*j+irow ].
real=a[ (*lda)*j+i ].
real;
352 A[ LDA*j+irow ].
imag=a[ (*lda)*j+i ].
imag;
353 }
354 }
355 else {
356 i=0;
357 irow=*k-i;
358 for( j=0; j<*n; j++ ) {
359 A[ LDA*j+irow ].
real=a[ (*lda)*j+i ].
real;
360 A[ LDA*j+irow ].
imag=a[ (*lda)*j+i ].
imag;
361 }
362 for( i=1; i<*k+1; i++ ){
363 irow=*k-i;
364 jcol=i;
365 for( j=jcol; j<(*n+*k); j++ ) {
366 A[ LDA*j+irow ].
real=a[ (*lda)*(j-jcol)+i ].real;
367 A[ LDA*j+irow ].
imag=a[ (*lda)*(j-jcol)+i ].imag;
368 }
369 }
370 }
372 *incx);
373 free(A);
374 }
375 }
378 else
380}
void cblas_ztbmv(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)