diff options
Diffstat (limited to 'ml/dlib/dlib/external/cblas/cblas_ctbmv.c')
-rw-r--r-- | ml/dlib/dlib/external/cblas/cblas_ctbmv.c | 139 |
1 files changed, 0 insertions, 139 deletions
diff --git a/ml/dlib/dlib/external/cblas/cblas_ctbmv.c b/ml/dlib/dlib/external/cblas/cblas_ctbmv.c deleted file mode 100644 index 0b313d858..000000000 --- a/ml/dlib/dlib/external/cblas/cblas_ctbmv.c +++ /dev/null @@ -1,139 +0,0 @@ -/* - * cblas_ctbmv.c - * The program is a C interface to ctbmv. - * - * Keita Teranishi 5/20/98 - * - */ -#include "cblas.h" -#include "cblas_f77.h" -void cblas_ctbmv(const enum CBLAS_ORDER order, const enum CBLAS_UPLO Uplo, - const enum CBLAS_TRANSPOSE TransA, const enum CBLAS_DIAG Diag, - const int N, const int K, const void *A, const int lda, - void *X, const int incX) -{ - char TA; - char UL; - char DI; -#ifdef F77_CHAR - F77_CHAR F77_TA, F77_UL, F77_DI; -#else - #define F77_TA &TA - #define F77_UL &UL - #define F77_DI &DI -#endif -#ifdef F77_INT - F77_INT F77_N=N, F77_lda=lda, F77_K=K, F77_incX=incX; -#else - #define F77_N N - #define F77_K K - #define F77_lda lda - #define F77_incX incX -#endif - int n, i=0, tincX; - float *st=0, *x=(float *)X; - - if (order == CblasColMajor) - { - if (Uplo == CblasUpper) UL = 'U'; - else if (Uplo == CblasLower) UL = 'L'; - else - { - cblas_xerbla(2, "cblas_ctbmv","Illegal Uplo setting, %d\n", Uplo); - return; - } - if (TransA == CblasNoTrans) TA = 'N'; - else if (TransA == CblasTrans) TA = 'T'; - else if (TransA == CblasConjTrans) TA = 'C'; - else - { - cblas_xerbla(3, "cblas_ctbmv","Illegal TransA setting, %d\n", TransA); - return; - } - if (Diag == CblasUnit) DI = 'U'; - else if (Diag == CblasNonUnit) DI = 'N'; - else - { - cblas_xerbla(4, "cblas_ctbmv","Illegal Diag setting, %d\n", Diag); - return; - } - #ifdef F77_CHAR - F77_UL = C2F_CHAR(&UL); - F77_TA = C2F_CHAR(&TA); - F77_DI = C2F_CHAR(&DI); - #endif - F77_ctbmv( F77_UL, F77_TA, F77_DI, &F77_N, &F77_K, A, &F77_lda, X, - &F77_incX); - } - else if (order == CblasRowMajor) - { - if (Uplo == CblasUpper) UL = 'L'; - else if (Uplo == CblasLower) UL = 'U'; - else - { - cblas_xerbla(2, "cblas_ctbmv","Illegal Uplo setting, %d\n", Uplo); - return; - } - - if (TransA == CblasNoTrans) TA = 'T'; - else if (TransA == CblasTrans) TA = 'N'; - else if (TransA == CblasConjTrans) - { - TA = 'N'; - if ( N > 0) - { - if(incX > 0) - tincX = incX; - else - tincX = -incX; - i = tincX << 1; - n = i * N; - x++; - st = x + n; - do - { - *x = -(*x); - x+= i; - } - while (x != st); - x -= n; - } - } - else - { - cblas_xerbla(3, "cblas_ctbmv","Illegal TransA setting, %d\n", TransA); - return; - } - - if (Diag == CblasUnit) DI = 'U'; - else if (Diag == CblasNonUnit) DI = 'N'; - else - { - cblas_xerbla(4, "cblas_ctbmv","Illegal Uplo setting, %d\n", Uplo); - return; - } - #ifdef F77_CHAR - F77_UL = C2F_CHAR(&UL); - F77_TA = C2F_CHAR(&TA); - F77_DI = C2F_CHAR(&DI); - #endif - - F77_ctbmv( F77_UL, F77_TA, F77_DI, &F77_N, &F77_K, A, &F77_lda, X, - &F77_incX); - - if (TransA == CblasConjTrans) - { - if (N > 0) - { - do - { - *x = -(*x); - x += i; - } - while (x != st); - } - } - } - else cblas_xerbla(1, "cblas_ctbmv", "Illegal Order setting, %d\n", order); - return; -} |