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 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