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