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