• Home
  • Line#
  • Scopes#
  • Navigate#
  • Raw
  • Download
1 /*
2  * cblas_cher2.c
3  * The program is a C interface to cher2.
4  *
5  * Keita Teranishi  3/23/98
6  *
7  */
8 #include <stdio.h>
9 #include <stdlib.h>
10 #include "cblas.h"
11 #include "cblas_f77.h"
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)12 void cblas_cher2(const enum CBLAS_ORDER order, const enum CBLAS_UPLO Uplo,
13                  const int N, const void *alpha, const void *X, const int incX,
14                  const void *Y, const int incY, void *A, const int lda)
15 {
16    char UL;
17 #ifdef F77_CHAR
18    F77_CHAR F77_UL;
19 #else
20    #define F77_UL &UL
21 #endif
22 
23 #ifdef F77_INT
24    F77_INT F77_N=N, F77_lda=lda, F77_incX=incX, F77_incY=incY;
25 #else
26    #define F77_N N
27    #define F77_lda lda
28    #define F77_incX incx
29    #define F77_incY incy
30 #endif
31    int n, i, j, tincx, tincy, incx=incX, incy=incY;
32    float *x=(float *)X, *xx=(float *)X, *y=(float *)Y,
33          *yy=(float *)Y, *tx, *ty, *stx, *sty;
34 
35    extern int CBLAS_CallFromC;
36    extern int RowMajorStrg;
37    RowMajorStrg = 0;
38 
39    CBLAS_CallFromC = 1;
40    if (order == CblasColMajor)
41    {
42       if (Uplo == CblasLower) UL = 'L';
43       else if (Uplo == CblasUpper) UL = 'U';
44       else
45       {
46          cblas_xerbla(2, "cblas_cher2","Illegal Uplo setting, %d\n",Uplo );
47          CBLAS_CallFromC = 0;
48          RowMajorStrg = 0;
49          return;
50       }
51       #ifdef F77_CHAR
52          F77_UL = C2F_CHAR(&UL);
53       #endif
54 
55       F77_cher2(F77_UL, &F77_N, alpha, X, &F77_incX,
56                                             Y, &F77_incY, A, &F77_lda);
57 
58    }  else if (order == CblasRowMajor)
59    {
60       RowMajorStrg = 1;
61       if (Uplo == CblasUpper) UL = 'L';
62       else if (Uplo == CblasLower) UL = 'U';
63       else
64       {
65          cblas_xerbla(2, "cblas_cher2","Illegal Uplo setting, %d\n", Uplo);
66          CBLAS_CallFromC = 0;
67          RowMajorStrg = 0;
68          return;
69       }
70       #ifdef F77_CHAR
71          F77_UL = C2F_CHAR(&UL);
72       #endif
73       if (N > 0)
74       {
75          n = N << 1;
76          x = malloc(n*sizeof(float));
77          y = malloc(n*sizeof(float));
78          tx = x;
79          ty = y;
80          if( incX > 0 ) {
81             i = incX << 1 ;
82             tincx = 2;
83             stx= x+n;
84          } else {
85             i = incX *(-2);
86             tincx = -2;
87             stx = x-2;
88             x +=(n-2);
89          }
90 
91          if( incY > 0 ) {
92             j = incY << 1;
93             tincy = 2;
94             sty= y+n;
95          } else {
96             j = incY *(-2);
97             tincy = -2;
98             sty = y-2;
99             y +=(n-2);
100          }
101 
102          do
103          {
104             *x = *xx;
105             x[1] = -xx[1];
106             x += tincx ;
107             xx += i;
108          }
109          while (x != stx);
110 
111          do
112          {
113             *y = *yy;
114             y[1] = -yy[1];
115             y += tincy ;
116             yy += j;
117          }
118          while (y != sty);
119 
120          x=tx;
121          y=ty;
122 
123          #ifdef F77_INT
124             F77_incX = 1;
125             F77_incY = 1;
126          #else
127             incx = 1;
128             incy = 1;
129          #endif
130       }  else
131       {
132          x = (float *) X;
133          y = (float *) Y;
134       }
135       F77_cher2(F77_UL, &F77_N, alpha, y, &F77_incY, x,
136                                       &F77_incX, A, &F77_lda);
137    } else
138    {
139       cblas_xerbla(1, "cblas_cher2","Illegal Order setting, %d\n", order);
140       CBLAS_CallFromC = 0;
141       RowMajorStrg = 0;
142       return;
143    }
144    if(X!=x)
145       free(x);
146    if(Y!=y)
147       free(y);
148 
149    CBLAS_CallFromC = 0;
150    RowMajorStrg = 0;
151    return;
152 }
153