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
53
54
55
|
#include <assert.h>
#include <float.h>
#include <math.h>
int main()
{
double xxx;
// Visual Studio needs to be 2013 onwards
#if defined(_MSC_VER) && !defined(__CYGWIN__) && _MSC_VER < 1800
// see http://www.johndcook.com/math_h.html
#else
assert(fpclassify(DBL_MAX + DBL_MAX) == FP_INFINITE);
assert(fpclassify(0 * (DBL_MAX + DBL_MAX)) == FP_NAN);
assert(fpclassify(1.0) == FP_NORMAL);
assert(fpclassify(DBL_MIN) == FP_NORMAL);
assert(fpclassify(DBL_MIN / 2) == FP_SUBNORMAL);
assert(fpclassify(-0.0) == FP_ZERO);
#endif
#if !defined(__clang__) && defined(__GNUC__)
assert(__builtin_fpclassify(0, 1, 2, 3, 4, DBL_MAX + DBL_MAX) == 1);
assert(__builtin_fpclassify(0, 1, 2, 3, 4, 0 * (DBL_MAX + DBL_MAX)) == 0);
assert(__builtin_fpclassify(0, 1, 2, 3, 4, 1.0) == 2);
assert(__builtin_fpclassify(0, 1, 2, 3, 4, DBL_MIN) == 2);
assert(__builtin_fpclassify(0, 1, 2, 3, 4, DBL_MIN / 2) == 3);
assert(__builtin_fpclassify(0, 1, 2, 3, 4, -0.0) == 4);
assert(__builtin_isinf(DBL_MAX + DBL_MAX) == 1);
assert(__builtin_isinf(0.0) == 0);
assert(__builtin_isinf(-(DBL_MAX + DBL_MAX)) == 1);
assert(__builtin_isinf_sign(DBL_MAX + DBL_MAX) == 1);
assert(__builtin_isinf_sign(0.0) == 0);
assert(__builtin_isinf_sign(-(DBL_MAX + DBL_MAX)) == -1);
// these are compile-time
_Static_assert(
__builtin_fpclassify(0, 1, 2, 3, 4, -0.0) == 4,
"__builtin_fpclassify is constant");
_Static_assert(__builtin_isnormal(DBL_MIN), "__builtin_isnormal is constant");
_Static_assert(!__builtin_isinf(0.0), "__builtin_isinf is constant");
_Static_assert(
__builtin_isinf_sign(0.0) == 0, "__builtin_isinf_sign is constant");
#endif
assert(signbit(-1.0) != 0);
assert(signbit(1.0) == 0);
}
|
