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