]> git.ipfire.org Git - thirdparty/glibc.git/blame - math/bug-nexttoward.c
Improve the accuracy of tgamma (BZ #26983)
[thirdparty/glibc.git] / math / bug-nexttoward.c
CommitLineData
07449987
UD
1#include <fenv.h>
2#include <math.h>
3#include <float.h>
4#include <stdlib.h>
5#include <stdio.h>
94f2c076 6#include <math-tests.h>
07449987 7
0af797de
CM
8#if !defined(FE_OVERFLOW) && !defined(FE_UNDERFLOW)
9/* If there's no support for the exceptions this test is checking,
10 then just return success and allow the test to be compiled. */
11# define fetestexcept(e) 1
12#endif
13
3e336a87
UD
14float zero = 0.0;
15float inf = INFINITY;
16
07449987
UD
17int
18main (void)
19{
20 int result = 0;
21
22 long double tl = (long double) FLT_MAX + 0x1.0p128L;
23 float fi = INFINITY;
24 float m = FLT_MAX;
25 feclearexcept (FE_ALL_EXCEPT);
26 if (nexttowardf (m, tl) != fi)
27 {
28 puts ("nexttowardf+ failed");
29 ++result;
30 }
94f2c076 31 if (EXCEPTION_TESTS (float) && fetestexcept (FE_OVERFLOW) == 0)
07449987
UD
32 {
33 puts ("nexttowardf+ did not overflow");
34 ++result;
35 }
36 feclearexcept (FE_ALL_EXCEPT);
37 if (nexttowardf (-m, -tl) != -fi)
38 {
39 puts ("nexttowardf- failed");
40 ++result;
41 }
94f2c076 42 if (EXCEPTION_TESTS (float) && fetestexcept (FE_OVERFLOW) == 0)
07449987
UD
43 {
44 puts ("nexttowardf- did not overflow");
45 ++result;
46 }
47
3e336a87
UD
48 fi = 0;
49 m = FLT_MIN;
50 feclearexcept (FE_ALL_EXCEPT);
51 fi = nexttowardf (m, fi);
52 if (fi < 0 || fi >= FLT_MIN)
53 {
54 puts ("nexttowardf+ failed");
55 ++result;
56 }
94f2c076 57 if (EXCEPTION_TESTS (float) && fetestexcept (FE_UNDERFLOW) == 0)
3e336a87
UD
58 {
59 puts ("nexttowardf+ did not underflow");
60 ++result;
61 }
62 fi = 0;
63 feclearexcept (FE_ALL_EXCEPT);
64 fi = nexttowardf (-m, -fi);
65 if (fi > 0 || fi <= -FLT_MIN)
66 {
67 puts ("nexttowardf- failed");
68 ++result;
69 }
94f2c076 70 if (EXCEPTION_TESTS (float) && fetestexcept (FE_UNDERFLOW) == 0)
3e336a87
UD
71 {
72 puts ("nexttowardf- did not underflow");
73 ++result;
74 }
75 fi = -INFINITY;
76 feclearexcept (FE_ALL_EXCEPT);
77 m = nexttowardf (zero, inf);
78 if (m < 0.0 || m >= FLT_MIN)
79 {
80 puts ("nexttowardf+ failed");
81 ++result;
82 }
94f2c076 83 if (EXCEPTION_TESTS (float) && fetestexcept (FE_UNDERFLOW) == 0)
3e336a87
UD
84 {
85 puts ("nexttowardf+ did not underflow");
86 ++result;
87 }
88 feclearexcept (FE_ALL_EXCEPT);
89 if (nexttowardf (m, fi) != 0.0)
90 {
91 puts ("nexttowardf+ failed");
92 ++result;
93 }
94f2c076 94 if (EXCEPTION_TESTS (float) && fetestexcept (FE_UNDERFLOW) == 0)
3e336a87
UD
95 {
96 puts ("nexttowardf+ did not underflow");
97 ++result;
98 }
99 feclearexcept (FE_ALL_EXCEPT);
100 m = nexttowardf (copysignf (zero, -1.0), -inf);
101 if (m > 0.0 || m <= -FLT_MIN)
102 {
103 puts ("nexttowardf- failed");
104 ++result;
105 }
94f2c076 106 if (EXCEPTION_TESTS (float) && fetestexcept (FE_UNDERFLOW) == 0)
3e336a87
UD
107 {
108 puts ("nexttowardf- did not underflow");
109 ++result;
110 }
111 feclearexcept (FE_ALL_EXCEPT);
112 if (nexttowardf (m, -fi) != 0.0)
113 {
114 puts ("nexttowardf- failed");
115 ++result;
116 }
94f2c076 117 if (EXCEPTION_TESTS (float) && fetestexcept (FE_UNDERFLOW) == 0)
3e336a87
UD
118 {
119 puts ("nexttowardf- did not underflow");
120 ++result;
121 }
122
07449987
UD
123 tl = (long double) DBL_MAX + 1.0e305L;
124 double di = INFINITY;
125 double dm = DBL_MAX;
126 feclearexcept (FE_ALL_EXCEPT);
127 if (nexttoward (dm, tl) != di)
128 {
129 puts ("nexttoward+ failed");
130 ++result;
131 }
94f2c076 132 if (EXCEPTION_TESTS (double) && fetestexcept (FE_OVERFLOW) == 0)
07449987
UD
133 {
134 puts ("nexttoward+ did not overflow");
135 ++result;
136 }
137 feclearexcept (FE_ALL_EXCEPT);
138 if (nexttoward (-dm, -tl) != -di)
139 {
140 puts ("nexttoward- failed");
141 ++result;
142 }
94f2c076 143 if (EXCEPTION_TESTS (double) && fetestexcept (FE_OVERFLOW) == 0)
07449987
UD
144 {
145 puts ("nexttoward- did not overflow");
146 ++result;
147 }
148
3e336a87
UD
149 di = 0;
150 dm = DBL_MIN;
151 feclearexcept (FE_ALL_EXCEPT);
152 di = nexttoward (dm, di);
153 if (di < 0 || di >= DBL_MIN)
154 {
155 puts ("nexttoward+ failed");
156 ++result;
157 }
94f2c076 158 if (EXCEPTION_TESTS (double) && fetestexcept (FE_UNDERFLOW) == 0)
3e336a87
UD
159 {
160 puts ("nexttoward+ did not underflow");
161 ++result;
162 }
163 di = 0;
164 feclearexcept (FE_ALL_EXCEPT);
165 di = nexttoward (-dm, -di);
166 if (di > 0 || di <= -DBL_MIN)
167 {
168 puts ("nexttoward- failed");
169 ++result;
170 }
94f2c076 171 if (EXCEPTION_TESTS (double) && fetestexcept (FE_UNDERFLOW) == 0)
3e336a87
UD
172 {
173 puts ("nexttoward- did not underflow");
174 ++result;
175 }
176 di = -INFINITY;
177 feclearexcept (FE_ALL_EXCEPT);
178 dm = nexttoward (zero, inf);
179 if (dm < 0.0 || dm >= DBL_MIN)
180 {
181 puts ("nexttoward+ failed");
182 ++result;
183 }
94f2c076 184 if (EXCEPTION_TESTS (double) && fetestexcept (FE_UNDERFLOW) == 0)
3e336a87
UD
185 {
186 puts ("nexttoward+ did not underflow");
187 ++result;
188 }
189 feclearexcept (FE_ALL_EXCEPT);
190 if (nexttoward (dm, di) != 0.0)
191 {
192 puts ("nexttoward+ failed");
193 ++result;
194 }
94f2c076 195 if (EXCEPTION_TESTS (double) && fetestexcept (FE_UNDERFLOW) == 0)
3e336a87
UD
196 {
197 puts ("nexttoward+ did not underflow");
198 ++result;
199 }
200 feclearexcept (FE_ALL_EXCEPT);
201 dm = nexttoward (copysign (zero, -1.0), -inf);
202 if (dm > 0.0 || dm <= -DBL_MIN)
203 {
204 puts ("nexttoward- failed");
205 ++result;
206 }
94f2c076 207 if (EXCEPTION_TESTS (double) && fetestexcept (FE_UNDERFLOW) == 0)
3e336a87
UD
208 {
209 puts ("nexttoward- did not underflow");
210 ++result;
211 }
212 feclearexcept (FE_ALL_EXCEPT);
213 if (nexttoward (dm, -di) != 0.0)
214 {
215 puts ("nexttoward- failed");
216 ++result;
217 }
94f2c076 218 if (EXCEPTION_TESTS (double) && fetestexcept (FE_UNDERFLOW) == 0)
3e336a87
UD
219 {
220 puts ("nexttoward- did not underflow");
221 ++result;
222 }
223
3e336a87
UD
224 long double li = INFINITY;
225 long double lm = LDBL_MAX;
226 feclearexcept (FE_ALL_EXCEPT);
227 if (nexttowardl (lm, li) != li)
228 {
229 puts ("nexttowardl+ failed");
230 ++result;
231 }
94f2c076 232 if (EXCEPTION_TESTS (long double) && fetestexcept (FE_OVERFLOW) == 0)
3e336a87
UD
233 {
234 puts ("nexttowardl+ did not overflow");
235 ++result;
236 }
237 feclearexcept (FE_ALL_EXCEPT);
238 if (nexttowardl (-lm, -li) != -li)
239 {
240 puts ("nexttowardl failed");
241 ++result;
242 }
94f2c076 243 if (EXCEPTION_TESTS (long double) && fetestexcept (FE_OVERFLOW) == 0)
3e336a87
UD
244 {
245 puts ("nexttowardl- did not overflow");
246 ++result;
247 }
248
249 li = 0;
250 lm = LDBL_MIN;
251 feclearexcept (FE_ALL_EXCEPT);
252 li = nexttowardl (lm, li);
253 if (li < 0 || li >= LDBL_MIN)
254 {
255 puts ("nexttowardl+ failed");
256 ++result;
257 }
94f2c076 258 if (EXCEPTION_TESTS (long double) && fetestexcept (FE_UNDERFLOW) == 0)
3e336a87
UD
259 {
260 puts ("nexttowardl+ did not underflow");
261 ++result;
262 }
263 li = 0;
264 feclearexcept (FE_ALL_EXCEPT);
265 li = nexttowardl (-lm, -li);
266 if (li > 0 || li <= -LDBL_MIN)
267 {
268 puts ("nexttowardl- failed");
269 ++result;
270 }
94f2c076 271 if (EXCEPTION_TESTS (long double) && fetestexcept (FE_UNDERFLOW) == 0)
3e336a87
UD
272 {
273 puts ("nexttowardl- did not underflow");
274 ++result;
275 }
276 li = -INFINITY;
277 feclearexcept (FE_ALL_EXCEPT);
278 lm = nexttowardl (zero, inf);
279 if (lm < 0.0 || lm >= LDBL_MIN)
280 {
281 puts ("nexttowardl+ failed");
282 ++result;
283 }
94f2c076 284 if (EXCEPTION_TESTS (long double) && fetestexcept (FE_UNDERFLOW) == 0)
3e336a87
UD
285 {
286 puts ("nexttowardl+ did not underflow");
287 ++result;
288 }
289 feclearexcept (FE_ALL_EXCEPT);
290 if (nexttowardl (lm, li) != 0.0)
291 {
292 puts ("nexttowardl+ failed");
293 ++result;
294 }
94f2c076 295 if (EXCEPTION_TESTS (long double) && fetestexcept (FE_UNDERFLOW) == 0)
3e336a87
UD
296 {
297 puts ("nexttowardl+ did not underflow");
298 ++result;
299 }
300 feclearexcept (FE_ALL_EXCEPT);
301 lm = nexttowardl (copysign (zero, -1.0), -inf);
302 if (lm > 0.0 || lm <= -LDBL_MIN)
303 {
304 puts ("nexttowardl- failed");
305 ++result;
306 }
94f2c076 307 if (EXCEPTION_TESTS (long double) && fetestexcept (FE_UNDERFLOW) == 0)
3e336a87
UD
308 {
309 puts ("nexttowardl- did not underflow");
310 ++result;
311 }
312 feclearexcept (FE_ALL_EXCEPT);
313 if (nexttowardl (lm, -li) != 0.0)
314 {
315 puts ("nexttowardl- failed");
316 ++result;
317 }
94f2c076 318 if (EXCEPTION_TESTS (long double) && fetestexcept (FE_UNDERFLOW) == 0)
3e336a87
UD
319 {
320 puts ("nexttowardl- did not underflow");
321 ++result;
322 }
3e336a87 323
07449987
UD
324 return result;
325}