OpenBLAS/interface/rotg.c at v0.3.22 · OpenMathLib/OpenBLAS · GitHub

Name: OpenBLAS/interface/rotg.c at v0.3.22 · OpenMathLib/OpenBLAS · GitHub
Rating: 4.9 (8738 reviews)
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
107
108
#include<math.h>
#include"common.h"
#ifdefFUNCTION_PROFILE
#include"functable.h"
#endif

#ifndefCBLAS

voidNAME(FLOAT*DA, FLOAT*DB, FLOAT*C, FLOAT*S){

#else

voidCNAME(FLOAT*DA, FLOAT*DB, FLOAT*C, FLOAT*S){

#endif

#if defined(__i386__) || defined(__x86_64__) || defined(__ia64__) || defined(_M_X64) || defined(_M_IX86)

long doubleda=*DA;
long doubledb=*DB;
long doublec;
long doubles;
long doubler, roe, z;

long doubleada=fabsl(da);
long doubleadb=fabsl(db);
long doublescale=ada+adb;

#ifndefCBLAS
PRINT_DEBUG_NAME;
#else
PRINT_DEBUG_CNAME;
#endif

roe=db;
if (ada>adb) roe=da;

if (scale==ZERO) {
*C=ONE;
*S=ZERO;
*DA=ZERO;
*DB=ZERO;
 } else {
r=sqrt(da*da+db*db);
if (roe<0) r=-r;
c=da / r;
s=db / r;
z=ONE;
if (da!=ZERO) {
if (ada>adb){
z=s;
 } else {
z=ONE / c;
 }
 }

*C=c;
*S=s;
*DA=r;
*DB=z;
 }

#else
FLOATda=*DA;
FLOATdb=*DB;
FLOATc=*C;
FLOATs=*S;
FLOATr, roe, z;

FLOATada=fabs(da);
FLOATadb=fabs(db);
FLOATscale=ada+adb;

#ifndefCBLAS
PRINT_DEBUG_NAME;
#else
PRINT_DEBUG_CNAME;
#endif

roe=db;
if (ada>adb) roe=da;

if (scale==ZERO) {
*C=ONE;
*S=ZERO;
*DA=ZERO;
*DB=ZERO;
 } else {
FLOATaa=da / scale;
FLOATbb=db / scale;

r=scale*sqrt(aa*aa+bb*bb);
if (roe<0) r=-r;
c=da / r;
s=db / r;
z=ONE;
if (ada>adb) z=s;
if ((ada <= adb) && (c!=ZERO)) z=ONE / c;

*C=c;
*S=s;
*DA=r;
*DB=z;
 }
#endif

return;
}