181253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray/* 281253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray * cblas_zhpr2.c 381253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray * The program is a C interface to zhpr2. 481253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray * 581253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray * Keita Teranishi 5/20/98 681253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray * 781253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray */ 881253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray#include <stdio.h> 981253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray#include <stdlib.h> 1081253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray#include "cblas.h" 1181253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray#include "cblas_f77.h" 1281253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murrayvoid cblas_zhpr2(const enum CBLAS_ORDER order, const enum CBLAS_UPLO Uplo, 1381253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray const int N,const void *alpha, const void *X, 1481253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray const int incX,const void *Y, const int incY, void *Ap) 1581253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray 1681253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray{ 1781253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray char UL; 1881253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray#ifdef F77_CHAR 1981253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray F77_CHAR F77_UL; 2081253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray#else 2181253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray #define F77_UL &UL 2281253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray#endif 2381253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray 2481253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray#ifdef F77_INT 2581253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray F77_INT F77_N=N, F77_incX=incX, F77_incY=incY; 2681253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray#else 2781253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray #define F77_N N 2881253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray #define F77_incX incx 2981253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray #define F77_incY incy 3081253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray#endif 3181253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray int n, i, j, incx=incX, incy=incY; 3281253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray double *x=(double *)X, *xx=(double *)X, *y=(double *)Y, 3381253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray *yy=(double *)Y, *stx, *sty; 3481253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray 3581253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray extern int CBLAS_CallFromC; 3681253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray extern int RowMajorStrg; 3781253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray RowMajorStrg = 0; 3881253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray 3981253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray CBLAS_CallFromC = 1; 4081253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray if (order == CblasColMajor) 4181253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray { 4281253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray if (Uplo == CblasLower) UL = 'L'; 4381253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray else if (Uplo == CblasUpper) UL = 'U'; 4481253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray else 4581253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray { 4681253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray cblas_xerbla(2, "cblas_zhpr2","Illegal Uplo setting, %d\n",Uplo ); 4781253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray CBLAS_CallFromC = 0; 4881253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray RowMajorStrg = 0; 4981253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray return; 5081253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray } 5181253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray #ifdef F77_CHAR 5281253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray F77_UL = C2F_CHAR(&UL); 5381253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray #endif 5481253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray 5581253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray F77_zhpr2(F77_UL, &F77_N, alpha, X, &F77_incX, Y, &F77_incY, Ap); 5681253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray 5781253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray } else if (order == CblasRowMajor) 5881253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray { 5981253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray RowMajorStrg = 1; 6081253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray if (Uplo == CblasUpper) UL = 'L'; 6181253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray else if (Uplo == CblasLower) UL = 'U'; 6281253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray else 6381253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray { 6481253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray cblas_xerbla(2, "cblas_zhpr2","Illegal Uplo setting, %d\n", Uplo); 6581253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray CBLAS_CallFromC = 0; 6681253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray RowMajorStrg = 0; 6781253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray return; 6881253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray } 6981253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray #ifdef F77_CHAR 7081253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray F77_UL = C2F_CHAR(&UL); 7181253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray #endif 7281253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray if (N > 0) 7381253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray { 7481253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray n = N << 1; 7581253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray x = malloc(n*sizeof(double)); 7681253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray y = malloc(n*sizeof(double)); 7781253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray stx = x + n; 7881253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray sty = y + n; 7981253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray if( incX > 0 ) 8081253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray i = incX << 1; 8181253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray else 8281253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray i = incX *(-2); 8381253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray 8481253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray if( incY > 0 ) 8581253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray j = incY << 1; 8681253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray else 8781253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray j = incY *(-2); 8881253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray do 8981253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray { 9081253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray *x = *xx; 9181253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray x[1] = -xx[1]; 9281253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray x += 2; 9381253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray xx += i; 9481253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray } while (x != stx); 9581253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray do 9681253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray { 9781253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray *y = *yy; 9881253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray y[1] = -yy[1]; 9981253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray y += 2; 10081253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray yy += j; 10181253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray } 10281253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray while (y != sty); 10381253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray x -= n; 10481253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray y -= n; 10581253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray 10681253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray #ifdef F77_INT 10781253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray if(incX > 0 ) 10881253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray F77_incX = 1; 10981253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray else 11081253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray F77_incX = -1; 11181253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray 11281253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray if(incY > 0 ) 11381253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray F77_incY = 1; 11481253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray else 11581253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray F77_incY = -1; 11681253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray 11781253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray #else 11881253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray if(incX > 0 ) 11981253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray incx = 1; 12081253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray else 12181253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray incx = -1; 12281253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray 12381253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray if(incY > 0 ) 12481253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray incy = 1; 12581253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray else 12681253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray incy = -1; 12781253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray #endif 12881253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray 12981253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray } else 13081253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray { 13181253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray x = (double *) X; 13281253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray y = (void *) Y; 13381253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray } 13481253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray F77_zhpr2(F77_UL, &F77_N, alpha, y, &F77_incY, x, &F77_incX, Ap); 13581253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray } 13681253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray else 13781253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray { 13881253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray cblas_xerbla(1, "cblas_zhpr2","Illegal Order setting, %d\n", order); 13981253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray CBLAS_CallFromC = 0; 14081253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray RowMajorStrg = 0; 14181253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray return; 14281253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray } 14381253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray if(X!=x) 14481253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray free(x); 14581253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray if(Y!=y) 14681253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray free(y); 14781253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray CBLAS_CallFromC = 0; 14881253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray RowMajorStrg = 0; 14981253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray return; 15081253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray} 151