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

source code of glibc/math/bug-nexttoward.c