1 /* 2 * 3 * cblas_sspmv.c 4 * This program is a C interface to sspmv. 5 * Written by Keita Teranishi 6 * 4/6/1998 7 * 8 */ 9 #include "cblas.h" 10 #include "cblas_f77.h" 11 void cblas_sspmv(const enum CBLAS_ORDER order, 12 const enum CBLAS_UPLO Uplo, const int N, 13 const float alpha, const float *AP, 14 const float *X, const int incX, const float beta, 15 float *Y, const int incY) 16 { 17 char UL; 18 #ifdef F77_CHAR 19 F77_CHAR F77_UL; 20 #else 21 #define F77_UL &UL 22 #endif 23 #ifdef F77_INT 24 F77_INT F77_N=N, F77_incX=incX, F77_incY=incY; 25 #else 26 #define F77_N N 27 #define F77_incX incX 28 #define F77_incY incY 29 #endif 30 extern int CBLAS_CallFromC; 31 extern int RowMajorStrg; 32 RowMajorStrg = 0; 33 34 CBLAS_CallFromC = 1; 35 if (order == CblasColMajor) 36 { 37 if (Uplo == CblasUpper) UL = 'U'; 38 else if (Uplo == CblasLower) UL = 'L'; 39 else 40 { 41 cblas_xerbla(2, "cblas_sspmv","Illegal Uplo setting, %d\n",Uplo ); 42 CBLAS_CallFromC = 0; 43 RowMajorStrg = 0; 44 return; 45 } 46 #ifdef F77_CHAR 47 F77_UL = C2F_CHAR(&UL); 48 #endif 49 F77_sspmv(F77_UL, &F77_N, &alpha, AP, X, 50 &F77_incX, &beta, Y, &F77_incY); 51 } 52 else if (order == CblasRowMajor) 53 { 54 RowMajorStrg = 1; 55 if (Uplo == CblasUpper) UL = 'L'; 56 else if (Uplo == CblasLower) UL = 'U'; 57 else 58 { 59 cblas_xerbla(2, "cblas_sspmv","Illegal Uplo setting, %d\n", Uplo); 60 CBLAS_CallFromC = 0; 61 RowMajorStrg = 0; 62 return; 63 } 64 #ifdef F77_CHAR 65 F77_UL = C2F_CHAR(&UL); 66 #endif 67 F77_sspmv(F77_UL, &F77_N, &alpha, 68 AP, X,&F77_incX, &beta, Y, &F77_incY); 69 } 70 else cblas_xerbla(1, "cblas_sspmv", "Illegal Order setting, %d\n", order); 71 CBLAS_CallFromC = 0; 72 RowMajorStrg = 0; 73 } 74
