summaryrefslogtreecommitdiffstats
path: root/ml/dlib/dlib/external/cblas/cblas_zher.c
diff options
context:
space:
mode:
Diffstat (limited to 'ml/dlib/dlib/external/cblas/cblas_zher.c')
-rw-r--r--ml/dlib/dlib/external/cblas/cblas_zher.c99
1 files changed, 0 insertions, 99 deletions
diff --git a/ml/dlib/dlib/external/cblas/cblas_zher.c b/ml/dlib/dlib/external/cblas/cblas_zher.c
deleted file mode 100644
index 30453737c..000000000
--- a/ml/dlib/dlib/external/cblas/cblas_zher.c
+++ /dev/null
@@ -1,99 +0,0 @@
-/*
- * cblas_zher.c
- * The program is a C interface to zher.
- *
- * Keita Teranishi 5/20/98
- *
- */
-#include <stdio.h>
-#include <stdlib.h>
-#include "cblas.h"
-#include "cblas_f77.h"
-void cblas_zher(const enum CBLAS_ORDER order, const enum CBLAS_UPLO Uplo,
- const int N, const double alpha, const void *X, const int incX
- ,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;
-#else
- #define F77_N N
- #define F77_lda lda
- #define F77_incX incx
-#endif
- int n, i, tincx, incx=incX;
- double *x=(double *)X, *xx=(double *)X, *tx, *st;
-
-
- if (order == CblasColMajor)
- {
- if (Uplo == CblasLower) UL = 'L';
- else if (Uplo == CblasUpper) UL = 'U';
- else
- {
- cblas_xerbla(2, "cblas_zher","Illegal Uplo setting, %d\n",Uplo );
- return;
- }
- #ifdef F77_CHAR
- F77_UL = C2F_CHAR(&UL);
- #endif
-
- F77_zher(F77_UL, &F77_N, &alpha, X, &F77_incX, A, &F77_lda);
-
- } else if (order == CblasRowMajor)
- {
- if (Uplo == CblasUpper) UL = 'L';
- else if (Uplo == CblasLower) UL = 'U';
- else
- {
- cblas_xerbla(2, "cblas_zher","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));
- tx = x;
- if( incX > 0 ) {
- i = incX << 1 ;
- tincx = 2;
- st= x+n;
- } else {
- i = incX *(-2);
- tincx = -2;
- st = x-2;
- x +=(n-2);
- }
- do
- {
- *x = *xx;
- x[1] = -xx[1];
- x += tincx ;
- xx += i;
- }
- while (x != st);
- x=tx;
-
- #ifdef F77_INT
- F77_incX = 1;
- #else
- incx = 1;
- #endif
- }
- else x = (double *) X;
- F77_zher(F77_UL, &F77_N, &alpha, x, &F77_incX, A, &F77_lda);
- } else cblas_xerbla(1, "cblas_zher", "Illegal Order setting, %d\n", order);
- if(X!=x)
- free(x);
-
- return;
-}