1 #include <math.h>
2 #include "mtest.h"
3 #include "test.h"
4
5 enum {LESS,EQUAL,GREATER,UNORD};
6
7 #define TEST(f,want) do{ \
8 int r, e; \
9 feclearexcept(FE_ALL_EXCEPT); \
10 r = (f); \
11 e = fetestexcept(FE_ALL_EXCEPT); \
12 if (r != (want)) \
13 t_error(#f " failed: got %d want %d\n", r, (want)); \
14 if (e & INVALID) \
15 t_error(#f " raised the invalid exception\n"); \
16 }while(0)
17
18 #undef T
19 #define T(a,b,rel) do{ \
20 TEST(isunordered(a, b), rel == UNORD); \
21 TEST(isless(a, b), rel == LESS); \
22 TEST(islessequal(a, b), rel == LESS || rel == EQUAL); \
23 TEST(islessgreater(a, b), rel == LESS || rel == GREATER); \
24 TEST(isgreater(a, b), rel == GREATER); \
25 TEST(isgreaterequal(a, b), rel == GREATER || rel == EQUAL); \
26 }while(0)
27
main()28 int main()
29 {
30 #pragma STDC FENV_ACCESS ON
31 volatile double huge = DBL_MAX;
32 volatile double tiny = DBL_MIN;
33 volatile double eps = DBL_EPSILON;
34 volatile float hugef = FLT_MAX;
35 volatile float tinyf = FLT_MIN;
36 volatile float epsf = FLT_EPSILON;
37 volatile long double hugel = LDBL_MAX;
38 volatile long double tinyl = LDBL_MIN;
39 volatile long double epsl = LDBL_EPSILON;
40
41 T(nan, 1.0, UNORD);
42 T(1.0, nan, UNORD);
43 T(nan, nan, UNORD);
44 T(nan, nan+1.0, UNORD);
45 T(nan, nan+1.0L, UNORD);
46
47 T(1.0, 1.1, LESS);
48 T(1.0, 1.0+eps, LESS);
49 T(1.0+eps, 1.0, GREATER);
50 T(huge-1, huge, EQUAL);
51 T(huge, huge*huge, LESS);
52 T(-0.0, 0.0, EQUAL);
53 T(-tiny, 0.0, LESS);
54 T(tiny, 2*tiny, LESS);
55 T(tiny*0x1p-9, tiny*0x1p-8, LESS);
56
57 T(1.0f, 1.1f, LESS);
58 T(1.0f, 1.0f+epsf, LESS);
59 T(1.0f+epsf, 1.0f, GREATER);
60 T(hugef-1, hugef, EQUAL);
61 T(hugef, hugef*hugef, LESS);
62 T(-0.0f, 0.0f, EQUAL);
63 T(-tinyf, 0.0f, LESS);
64 T(tinyf, 2*tinyf, LESS);
65 T(tinyf*0x1p-9f, tinyf*0x1p-8f, LESS);
66
67 T(1.0L, 1.1L, LESS);
68 T(1.0L, 1.0L+epsl, LESS);
69 T(1.0L+epsl, 1.0L, GREATER);
70 T(hugel-1, hugel, EQUAL);
71 T(hugel, hugel*hugel, LESS);
72 T(-0.0L, 0.0L, EQUAL);
73 T(-tinyl, 0.0L, LESS);
74 T(tinyl, 2*tinyl, LESS);
75 T(tinyl*0x1p-9L, tinyl*0x1p-8L, LESS);
76
77 #if FLT_EVAL_METHOD == 2
78 T(huge*huge, huge*huge*2, LESS);
79 T(tiny*tiny*0.5, tiny*tiny, LESS);
80 T(-tiny*tiny, 0.0, LESS);
81 T(1.0, 1.0+eps*0.25, LESS);
82 #else
83 T(huge*huge, huge*huge*2, EQUAL);
84 T(tiny*tiny*0.5, tiny*tiny, EQUAL);
85 T(-tiny*tiny, 0.0, EQUAL);
86 T(1.0, 1.0+eps*0.25, EQUAL);
87 #endif
88
89 #if FLT_EVAL_METHOD >= 1
90 T(hugef*hugef, hugef*hugef*2, LESS);
91 T(tinyf*tinyf*0.5f, tinyf*tinyf, LESS);
92 T(-tinyf*tinyf, 0.0f, LESS);
93 T(1.0f, 1.0f+epsf*0.25f, LESS);
94 #else
95 T(hugef*hugef, hugef*hugef*2, EQUAL);
96 T(tinyf*tinyf*0.5f, tinyf*tinyf, EQUAL);
97 T(-tinyf*tinyf, 0.0f, EQUAL);
98 T(1.0f, 1.0f+epsf*0.25f, EQUAL);
99 #endif
100
101 T(hugel*hugel, hugel*hugel*2, EQUAL);
102 T(tinyl*tinyl*0.5L, tinyl*tinyl, EQUAL);
103 T(-tinyl*tinyl, 0.0L, EQUAL);
104 T(1.0L, 1.0L+epsl*0.25L, EQUAL);
105
106 return t_status;
107 }
108