diff options
Diffstat (limited to 'ml/dlib/dlib/external/cblas/cblas_zher2.c')
-rw-r--r-- | ml/dlib/dlib/external/cblas/cblas_zher2.c | 140 |
1 files changed, 0 insertions, 140 deletions
diff --git a/ml/dlib/dlib/external/cblas/cblas_zher2.c b/ml/dlib/dlib/external/cblas/cblas_zher2.c deleted file mode 100644 index 8bf0bd733..000000000 --- a/ml/dlib/dlib/external/cblas/cblas_zher2.c +++ /dev/null @@ -1,140 +0,0 @@ -/* - * cblas_zher2.c - * The program is a C interface to zher2. - * - * Keita Teranishi 3/23/98 - * - */ -#include <stdio.h> -#include <stdlib.h> -#include "cblas.h" -#include "cblas_f77.h" -void cblas_zher2(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; - double *x=(double *)X, *xx=(double *)X, *y=(double *)Y, - *yy=(double *)Y, *tx, *ty, *stx, *sty; - - - if (order == CblasColMajor) - { - if (Uplo == CblasLower) UL = 'L'; - else if (Uplo == CblasUpper) UL = 'U'; - else - { - cblas_xerbla(2, "cblas_zher2", "Illegal Uplo setting, %d\n",Uplo ); - return; - } - #ifdef F77_CHAR - F77_UL = C2F_CHAR(&UL); - #endif - - F77_zher2(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_zher2", "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(double)); - y = malloc(n*sizeof(double)); - 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 = (double *) X; - y = (double *) Y; - } - F77_zher2(F77_UL, &F77_N, alpha, y, &F77_incY, x, - &F77_incX, A, &F77_lda); - } - else - { - cblas_xerbla(1, "cblas_zher2", "Illegal Order setting, %d\n", order); - return; - } - if(X!=x) - free(x); - if(Y!=y) - free(y); - - return; -} |