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 | |
14 | float zero = 0.0; |
15 | float inf = INFINITY; |
16 | |
17 | int |
18 | main (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 | |