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