summaryrefslogtreecommitdiffstats
path: root/ml/dlib/dlib/external/cblas/cblas_ctbmv.c
diff options
context:
space:
mode:
Diffstat (limited to 'ml/dlib/dlib/external/cblas/cblas_ctbmv.c')
-rw-r--r--ml/dlib/dlib/external/cblas/cblas_ctbmv.c139
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;
-}