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