xref: /aosp_15_r20/external/cblas/src/cblas_ssymm.c (revision 1858f9982ea1ad57fb52080c08594e4d1cce4fa1)
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