1 /* $OpenBSD: nextafter.c,v 1.2 2011/07/09 03:33:07 martynas Exp $ */
2
3 /*
4 * Written by Martynas Venckus. Public domain
5 */
6
7 #include <assert.h>
8 #include <math.h>
9
10 #define test(f, r) ( \
11 ((f) == (r) && signbit((f)) == signbit((r))) || \
12 (isnan((f)) && isnan((r))) \
13 )
14
15 int
main(int argc,char * argv[])16 main(int argc, char *argv[])
17 {
18 assert(test(nextafter(0.0, 0.0), 0.0));
19 assert(test(nextafter(-0.0, 0.0), 0.0));
20 assert(test(nextafter(0.0, -0.0), -0.0));
21 assert(test(nextafter(-0.0, -0.0), -0.0));
22
23 assert(test(nextafterf(0.0F, 0.0F), 0.0F));
24 assert(test(nextafterf(-0.0F, 0.0F), 0.0F));
25 assert(test(nextafterf(0.0F, -0.0F), -0.0F));
26 assert(test(nextafterf(-0.0F, -0.0F), -0.0F));
27
28 assert(test(nextafterl(0.0L, 0.0L), 0.0L));
29 assert(test(nextafterl(-0.0L, 0.0L), 0.0L));
30 assert(test(nextafterl(0.0L, -0.0L), -0.0L));
31 assert(test(nextafterf(-0.0L, -0.0L), -0.0L));
32
33 assert(test(nextafter(NAN, 1.0), NAN));
34 assert(test(nextafter(1.0, NAN), NAN));
35 assert(test(nextafter(NAN, NAN), NAN));
36
37 assert(test(nextafterf(NAN, 1.0F), NAN));
38 assert(test(nextafterf(1.0F, NAN), NAN));
39 assert(test(nextafterf(NAN, NAN), NAN));
40
41 assert(test(nextafterl(NAN, 1.0L), NAN));
42 assert(test(nextafterl(1.0L, NAN), NAN));
43 assert(test(nextafterl(NAN, NAN), NAN));
44
45 assert(test(nextafter(0x1.fffffffffffffp+0, INFINITY), 0x1p1));
46 assert(test(nextafter(0x1p1, -INFINITY), 0x1.fffffffffffffp+0));
47
48 assert(test(nextafterf(0x1.fffffep+0f, INFINITY), 0x1p1f));
49 assert(test(nextafterf(0x1p1f, -INFINITY), 0x1.fffffep+0f));
50
51 return (0);
52 }
53