diff options
Diffstat (limited to 'ml/dlib/dlib/external/cblas/cblas_cher2.c')
-rw-r--r-- | ml/dlib/dlib/external/cblas/cblas_cher2.c | 139 |
1 files changed, 0 insertions, 139 deletions
diff --git a/ml/dlib/dlib/external/cblas/cblas_cher2.c b/ml/dlib/dlib/external/cblas/cblas_cher2.c deleted file mode 100644 index 89d36a0ee..000000000 --- a/ml/dlib/dlib/external/cblas/cblas_cher2.c +++ /dev/null @@ -1,139 +0,0 @@ -/* - * cblas_cher2.c - * The program is a C interface to cher2. - * - * Keita Teranishi 3/23/98 - * - */ -#include <stdio.h> -#include <stdlib.h> -#include "cblas.h" -#include "cblas_f77.h" -void cblas_cher2(const enum CBLAS_ORDER order, const enum CBLAS_UPLO Uplo, - const int N, const void *alpha, const void *X, const int incX, - const void *Y, const int incY, void *A, const int lda) -{ - char UL; -#ifdef F77_CHAR - F77_CHAR F77_UL; -#else - #define F77_UL &UL -#endif - -#ifdef F77_INT - F77_INT F77_N=N, F77_lda=lda, F77_incX=incX, F77_incY=incY; -#else - #define F77_N N - #define F77_lda lda - #define F77_incX incx - #define F77_incY incy -#endif - int n, i, j, tincx, tincy, incx=incX, incy=incY; - float *x=(float *)X, *xx=(float *)X, *y=(float *)Y, - *yy=(float *)Y, *tx, *ty, *stx, *sty; - - - if (order == CblasColMajor) - { - if (Uplo == CblasLower) UL = 'L'; - else if (Uplo == CblasUpper) UL = 'U'; - else - { - cblas_xerbla(2, "cblas_cher2","Illegal Uplo setting, %d\n",Uplo ); - return; - } - #ifdef F77_CHAR - F77_UL = C2F_CHAR(&UL); - #endif - - F77_cher2(F77_UL, &F77_N, alpha, X, &F77_incX, - Y, &F77_incY, A, &F77_lda); - - } else if (order == CblasRowMajor) - { - if (Uplo == CblasUpper) UL = 'L'; - else if (Uplo == CblasLower) UL = 'U'; - else - { - cblas_xerbla(2, "cblas_cher2","Illegal Uplo setting, %d\n", Uplo); - return; - } - #ifdef F77_CHAR - F77_UL = C2F_CHAR(&UL); - #endif - if (N > 0) - { - n = N << 1; - x = malloc(n*sizeof(float)); - y = malloc(n*sizeof(float)); - tx = x; - ty = y; - if( incX > 0 ) { - i = incX << 1 ; - tincx = 2; - stx= x+n; - } else { - i = incX *(-2); - tincx = -2; - stx = x-2; - x +=(n-2); - } - - if( incY > 0 ) { - j = incY << 1; - tincy = 2; - sty= y+n; - } else { - j = incY *(-2); - tincy = -2; - sty = y-2; - y +=(n-2); - } - - do - { - *x = *xx; - x[1] = -xx[1]; - x += tincx ; - xx += i; - } - while (x != stx); - - do - { - *y = *yy; - y[1] = -yy[1]; - y += tincy ; - yy += j; - } - while (y != sty); - - x=tx; - y=ty; - - #ifdef F77_INT - F77_incX = 1; - F77_incY = 1; - #else - incx = 1; - incy = 1; - #endif - } else - { - x = (float *) X; - y = (float *) Y; - } - F77_cher2(F77_UL, &F77_N, alpha, y, &F77_incY, x, - &F77_incX, A, &F77_lda); - } else - { - cblas_xerbla(1, "cblas_cher2","Illegal Order setting, %d\n", order); - return; - } - if(X!=x) - free(x); - if(Y!=y) - free(y); - - return; -} |