1 /*
2 *
3 * cblas_csyrk.c
4 * This program is a C interface to csyrk.
5 * Written by Keita Teranishi
6 * 4/8/1998
7 *
8 */
9
10 #include "cblas.h"
11 #include "cblas_f77.h"
cblas_csyrk(const enum CBLAS_ORDER Order,const enum CBLAS_UPLO Uplo,const enum CBLAS_TRANSPOSE Trans,const CBLAS_INT_TYPE N,const CBLAS_INT_TYPE K,const void * alpha,const void * A,const CBLAS_INT_TYPE lda,const void * beta,void * C,const CBLAS_INT_TYPE ldc)12 void cblas_csyrk(const enum CBLAS_ORDER Order, const enum CBLAS_UPLO Uplo,
13 const enum CBLAS_TRANSPOSE Trans, const CBLAS_INT_TYPE N, const CBLAS_INT_TYPE K,
14 const void *alpha, const void *A, const CBLAS_INT_TYPE lda,
15 const void *beta, void *C, const CBLAS_INT_TYPE ldc)
16 {
17 char UL, TR;
18 #ifdef F77_CHAR
19 F77_CHAR F77_TR, F77_UL;
20 #else
21 #define F77_TR &TR
22 #define F77_UL &UL
23 #endif
24
25 #ifdef F77_INT
26 F77_INT F77_N=N, F77_K=K, F77_lda=lda;
27 F77_INT F77_ldc=ldc;
28 #else
29 #define F77_N N
30 #define F77_K K
31 #define F77_lda lda
32 #define F77_ldc ldc
33 #endif
34
35
36 if( Order == CblasColMajor )
37 {
38
39 if( Uplo == CblasUpper) UL='U';
40 else if ( Uplo == CblasLower ) UL='L';
41 else
42 {
43 cblas_xerbla(2, "cblas_csyrk", "Illegal Uplo setting, %d\n", Uplo);
44 return;
45 }
46
47 if( Trans == CblasTrans) TR ='T';
48 else if ( Trans == CblasConjTrans ) TR='C';
49 else if ( Trans == CblasNoTrans ) TR='N';
50 else
51 {
52 cblas_xerbla(3, "cblas_csyrk", "Illegal Trans setting, %d\n", Trans);
53 return;
54 }
55
56
57 #ifdef F77_CHAR
58 F77_UL = C2F_CHAR(&UL);
59 F77_TR = C2F_CHAR(&TR);
60 #endif
61
62 F77_csyrk(F77_UL, F77_TR, &F77_N, &F77_K, alpha, A, &F77_lda,
63 beta, C, &F77_ldc);
64 } else if (Order == CblasRowMajor)
65 {
66 if( Uplo == CblasUpper) UL='L';
67 else if ( Uplo == CblasLower ) UL='U';
68 else
69 {
70 cblas_xerbla(3, "cblas_csyrk", "Illegal Uplo setting, %d\n", Uplo);
71 return;
72 }
73 if( Trans == CblasTrans) TR ='N';
74 else if ( Trans == CblasConjTrans ) TR='N';
75 else if ( Trans == CblasNoTrans ) TR='T';
76 else
77 {
78 cblas_xerbla(3, "cblas_csyrk", "Illegal Trans setting, %d\n", Trans);
79 return;
80 }
81
82 #ifdef F77_CHAR
83 F77_UL = C2F_CHAR(&UL);
84 F77_TR = C2F_CHAR(&TR);
85 #endif
86
87 F77_csyrk(F77_UL, F77_TR, &F77_N, &F77_K, alpha, A, &F77_lda,
88 beta, C, &F77_ldc);
89 }
90 else cblas_xerbla(1, "cblas_csyrk", "Illegal Order setting, %d\n", Order);
91 return;
92 }
93
94