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