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
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
|
/* Automatic switching between software and hardware IEEE 128-bit
floating-point emulation for PowerPC.
Copyright (C) 2016-2019 Free Software Foundation, Inc.
This file is part of the GNU C Library.
Contributed by Michael Meissner (meissner@linux.vnet.ibm.com)
Code is based on the main soft-fp library written by:
Richard Henderson (rth@cygnus.com) and
Jakub Jelinek (jj@ultra.linux.cz).
The GNU C Library is free software; you can redistribute it and/or
modify it under the terms of the GNU Lesser General Public
License as published by the Free Software Foundation; either
version 2.1 of the License, or (at your option) any later version.
In addition to the permissions in the GNU Lesser General Public
License, the Free Software Foundation gives you unlimited
permission to link the compiled version of this file into
combinations with other programs, and to distribute those
combinations without any restriction coming from the use of this
file. (The Lesser General Public License restrictions do apply in
other respects; for example, they cover modification of the file,
and distribution when not linked into a combine executable.)
The GNU C Library is distributed in the hope that it will be useful,
but WITHOUT ANY WARRANTY; without even the implied warranty of
MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
Lesser General Public License for more details.
You should have received a copy of the GNU Lesser General Public
License along with the GNU C Library; if not, see
<http://www.gnu.org/licenses/>. */
#include <soft-fp.h>
#include <quad-float128.h>
#ifndef __FLOAT128_HARDWARE__
#error "This module must be compiled with IEEE 128-bit hardware support"
#endif
TFtype
__addkf3_hw (TFtype a, TFtype b)
{
return a + b;
}
TFtype
__subkf3_hw (TFtype a, TFtype b)
{
return a - b;
}
TFtype
__mulkf3_hw (TFtype a, TFtype b)
{
return a * b;
}
TFtype
__divkf3_hw (TFtype a, TFtype b)
{
return a / b;
}
TFtype
__negkf2_hw (TFtype a)
{
return -a;
}
TFtype
__floatsikf_hw (SItype_ppc a)
{
return (TFtype) a;
}
TFtype
__floatunsikf_hw (USItype_ppc a)
{
return (TFtype) a;
}
TFtype
__floatdikf_hw (DItype_ppc a)
{
return (TFtype) a;
}
TFtype
__floatundikf_hw (UDItype_ppc a)
{
return (TFtype) a;
}
SItype_ppc
__fixkfsi_hw (TFtype a)
{
return (SItype_ppc) a;
}
USItype_ppc
__fixunskfsi_hw (TFtype a)
{
return (USItype_ppc) a;
}
DItype_ppc
__fixkfdi_hw (TFtype a)
{
return (DItype_ppc) a;
}
UDItype_ppc
__fixunskfdi_hw (TFtype a)
{
return (UDItype_ppc) a;
}
TFtype
__extendsfkf2_hw (float a)
{
return (TFtype) a;
}
TFtype
__extenddfkf2_hw (double a)
{
return (TFtype) a;
}
float
__trunckfsf2_hw (TFtype a)
{
return (float) a;
}
double
__trunckfdf2_hw (TFtype a)
{
return (double) a;
}
/* __eqkf2 returns 0 if equal, or 1 if not equal or NaN. */
CMPtype
__eqkf2_hw (TFtype a, TFtype b)
{
return (a != b);
}
/* __gekf2 returns -1 if a < b, 0 if a == b, +1 if a > b, or -2 if NaN. */
CMPtype
__gekf2_hw (TFtype a, TFtype b)
{
if (a < b)
return -1;
else if (__builtin_isunordered (a, b))
return -2;
else if (a == b)
return 0;
return 1;
}
/* __lekf2 returns -1 if a < b, 0 if a == b, +1 if a > b, or +2 if NaN. */
CMPtype
__lekf2_hw (TFtype a, TFtype b)
{
if (a < b)
return -1;
else if (__builtin_isunordered (a, b))
return 2;
else if (a == b)
return 0;
return 1;
}
/* __unordkf2 returns 1 if NaN or 0 otherwise. */
CMPtype
__unordkf2_hw (TFtype a, TFtype b)
{
return (__builtin_isunordered (a, b)) ? 1 : 0;
}
/* Convert __float128 to __ibm128. */
IBM128_TYPE
__extendkftf2_hw (TFtype value)
{
IBM128_TYPE ret;
CVT_FLOAT128_TO_IBM128 (ret, value);
return ret;
}
/* Convert __ibm128 to __float128. */
TFtype
__trunctfkf2_hw (IBM128_TYPE value)
{
TFtype ret;
CVT_IBM128_TO_FLOAT128 (ret, value);
return ret;
}
|
