- Notifications
You must be signed in to change notification settings - Fork 1.6k
/
Copy pathcblas_chemm.c
106 lines (96 loc) · 2.58 KB
/
cblas_chemm.c
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
/*
*
* cblas_chemm.c
* This program is a C interface to chemm.
* Written by Keita Teranishi
* 4/8/1998
*
*/
#include"cblas.h"
#include"cblas_f77.h"
voidcblas_chemm(constCBLAS_LAYOUTlayout, constCBLAS_SIDESide,
constCBLAS_UPLOUplo, constintM, constintN,
constvoid*alpha, constvoid*A, constintlda,
constvoid*B, constintldb, constvoid*beta,
void*C, constintldc)
{
charSD, UL;
#ifdefF77_CHAR
F77_CHARF77_SD, F77_UL;
#else
#defineF77_SD &SD
#defineF77_UL &UL
#endif
#ifdefF77_INT
F77_INTF77_M=M, F77_N=N, F77_lda=lda, F77_ldb=ldb;
F77_INTF77_ldc=ldc;
#else
#defineF77_M M
#defineF77_N N
#defineF77_lda lda
#defineF77_ldb ldb
#defineF77_ldc ldc
#endif
externintCBLAS_CallFromC;
externintRowMajorStrg;
RowMajorStrg=0;
CBLAS_CallFromC=1;
if( layout==CblasColMajor )
{
if( Side==CblasRight) SD='R';
elseif ( Side==CblasLeft ) SD='L';
else
{
cblas_xerbla(2, "cblas_chemm", "Illegal Side setting, %d\n", Side);
CBLAS_CallFromC=0;
RowMajorStrg=0;
return;
}
if( Uplo==CblasUpper) UL='U';
elseif ( Uplo==CblasLower ) UL='L';
else
{
cblas_xerbla(3, "cblas_chemm", "Illegal Uplo setting, %d\n", Uplo);
CBLAS_CallFromC=0;
RowMajorStrg=0;
return;
}
#ifdefF77_CHAR
F77_UL=C2F_CHAR(&UL);
F77_SD=C2F_CHAR(&SD);
#endif
F77_chemm(F77_SD, F77_UL, &F77_M, &F77_N, alpha, A, &F77_lda,
B, &F77_ldb, beta, C, &F77_ldc);
} elseif (layout==CblasRowMajor)
{
RowMajorStrg=1;
if( Side==CblasRight) SD='L';
elseif ( Side==CblasLeft ) SD='R';
else
{
cblas_xerbla(2, "cblas_chemm", "Illegal Side setting, %d\n", Side);
CBLAS_CallFromC=0;
RowMajorStrg=0;
return;
}
if( Uplo==CblasUpper) UL='L';
elseif ( Uplo==CblasLower ) UL='U';
else
{
cblas_xerbla(3, "cblas_chemm", "Illegal Uplo setting, %d\n", Uplo);
CBLAS_CallFromC=0;
RowMajorStrg=0;
return;
}
#ifdefF77_CHAR
F77_UL=C2F_CHAR(&UL);
F77_SD=C2F_CHAR(&SD);
#endif
F77_chemm(F77_SD, F77_UL, &F77_N, &F77_M, alpha, A,
&F77_lda, B, &F77_ldb, beta, C, &F77_ldc);
}
elsecblas_xerbla(1, "cblas_chemm", "Illegal layout setting, %d\n", layout);
CBLAS_CallFromC=0;
RowMajorStrg=0;
return;
}