1*1858f998SYi Kong /*
2*1858f998SYi Kong *
3*1858f998SYi Kong * cblas_ssymm.c
4*1858f998SYi Kong * This program is a C interface to ssymm.
5*1858f998SYi Kong * Written by Keita Teranishi
6*1858f998SYi Kong * 4/8/1998
7*1858f998SYi Kong *
8*1858f998SYi Kong */
9*1858f998SYi Kong
10*1858f998SYi Kong #include "cblas.h"
11*1858f998SYi Kong #include "cblas_f77.h"
cblas_ssymm(const enum CBLAS_ORDER Order,const enum CBLAS_SIDE Side,const enum CBLAS_UPLO Uplo,const int M,const int N,const float alpha,const float * A,const int lda,const float * B,const int ldb,const float beta,float * C,const int ldc)12*1858f998SYi Kong void cblas_ssymm(const enum CBLAS_ORDER Order, const enum CBLAS_SIDE Side,
13*1858f998SYi Kong const enum CBLAS_UPLO Uplo, const int M, const int N,
14*1858f998SYi Kong const float alpha, const float *A, const int lda,
15*1858f998SYi Kong const float *B, const int ldb, const float beta,
16*1858f998SYi Kong float *C, const int ldc)
17*1858f998SYi Kong {
18*1858f998SYi Kong char SD, UL;
19*1858f998SYi Kong #ifdef F77_CHAR
20*1858f998SYi Kong F77_CHAR F77_SD, F77_UL;
21*1858f998SYi Kong #else
22*1858f998SYi Kong #define F77_SD &SD
23*1858f998SYi Kong #define F77_UL &UL
24*1858f998SYi Kong #endif
25*1858f998SYi Kong
26*1858f998SYi Kong #ifdef F77_INT
27*1858f998SYi Kong F77_INT F77_M=M, F77_N=N, F77_lda=lda, F77_ldb=ldb;
28*1858f998SYi Kong F77_INT F77_ldc=ldc;
29*1858f998SYi Kong #else
30*1858f998SYi Kong #define F77_M M
31*1858f998SYi Kong #define F77_N N
32*1858f998SYi Kong #define F77_lda lda
33*1858f998SYi Kong #define F77_ldb ldb
34*1858f998SYi Kong #define F77_ldc ldc
35*1858f998SYi Kong #endif
36*1858f998SYi Kong
37*1858f998SYi Kong extern int CBLAS_CallFromC;
38*1858f998SYi Kong extern int RowMajorStrg;
39*1858f998SYi Kong RowMajorStrg = 0;
40*1858f998SYi Kong CBLAS_CallFromC = 1;
41*1858f998SYi Kong
42*1858f998SYi Kong if( Order == CblasColMajor )
43*1858f998SYi Kong {
44*1858f998SYi Kong if( Side == CblasRight) SD='R';
45*1858f998SYi Kong else if ( Side == CblasLeft ) SD='L';
46*1858f998SYi Kong else
47*1858f998SYi Kong {
48*1858f998SYi Kong cblas_xerbla(2, "cblas_ssymm",
49*1858f998SYi Kong "Illegal Side setting, %d\n", Side);
50*1858f998SYi Kong CBLAS_CallFromC = 0;
51*1858f998SYi Kong RowMajorStrg = 0;
52*1858f998SYi Kong return;
53*1858f998SYi Kong }
54*1858f998SYi Kong
55*1858f998SYi Kong if( Uplo == CblasUpper) UL='U';
56*1858f998SYi Kong else if ( Uplo == CblasLower ) UL='L';
57*1858f998SYi Kong else
58*1858f998SYi Kong {
59*1858f998SYi Kong cblas_xerbla(3, "cblas_ssymm",
60*1858f998SYi Kong "Illegal Uplo setting, %d\n", Uplo);
61*1858f998SYi Kong CBLAS_CallFromC = 0;
62*1858f998SYi Kong RowMajorStrg = 0;
63*1858f998SYi Kong return;
64*1858f998SYi Kong }
65*1858f998SYi Kong
66*1858f998SYi Kong #ifdef F77_CHAR
67*1858f998SYi Kong F77_UL = C2F_CHAR(&UL);
68*1858f998SYi Kong F77_SD = C2F_CHAR(&SD);
69*1858f998SYi Kong #endif
70*1858f998SYi Kong
71*1858f998SYi Kong F77_ssymm(F77_SD, F77_UL, &F77_M, &F77_N, &alpha, A, &F77_lda, B, &F77_ldb, &beta, C, &F77_ldc);
72*1858f998SYi Kong } else if (Order == CblasRowMajor)
73*1858f998SYi Kong {
74*1858f998SYi Kong RowMajorStrg = 1;
75*1858f998SYi Kong if( Side == CblasRight) SD='L';
76*1858f998SYi Kong else if ( Side == CblasLeft ) SD='R';
77*1858f998SYi Kong else
78*1858f998SYi Kong {
79*1858f998SYi Kong cblas_xerbla(2, "cblas_ssymm",
80*1858f998SYi Kong "Illegal Side setting, %d\n", Side);
81*1858f998SYi Kong CBLAS_CallFromC = 0;
82*1858f998SYi Kong RowMajorStrg = 0;
83*1858f998SYi Kong return;
84*1858f998SYi Kong }
85*1858f998SYi Kong
86*1858f998SYi Kong if( Uplo == CblasUpper) UL='L';
87*1858f998SYi Kong else if ( Uplo == CblasLower ) UL='U';
88*1858f998SYi Kong else
89*1858f998SYi Kong {
90*1858f998SYi Kong cblas_xerbla(3, "cblas_ssymm",
91*1858f998SYi Kong "Illegal Uplo setting, %d\n", Uplo);
92*1858f998SYi Kong CBLAS_CallFromC = 0;
93*1858f998SYi Kong RowMajorStrg = 0;
94*1858f998SYi Kong return;
95*1858f998SYi Kong }
96*1858f998SYi Kong
97*1858f998SYi Kong #ifdef F77_CHAR
98*1858f998SYi Kong F77_UL = C2F_CHAR(&UL);
99*1858f998SYi Kong F77_SD = C2F_CHAR(&SD);
100*1858f998SYi Kong #endif
101*1858f998SYi Kong
102*1858f998SYi Kong F77_ssymm(F77_SD, F77_UL, &F77_N, &F77_M, &alpha, A, &F77_lda, B, &F77_ldb, &beta, C, &F77_ldc);
103*1858f998SYi Kong } else cblas_xerbla(1, "cblas_ssymm",
104*1858f998SYi Kong "Illegal Order setting, %d\n", Order);
105*1858f998SYi Kong CBLAS_CallFromC = 0;
106*1858f998SYi Kong RowMajorStrg = 0;
107*1858f998SYi Kong return;
108*1858f998SYi Kong }
109