1 #include <fenv.h>
2 #include <float.h>
3 #include <math.h>
4
5 #undef RN
6 #undef RZ
7 #undef RD
8 #undef RU
9 #ifdef FE_TONEAREST
10 #define RN FE_TONEAREST
11 #else
12 #define RN 0
13 #endif
14 #ifdef FE_TOWARDZERO
15 #define RZ FE_TOWARDZERO
16 #else
17 #define RZ -1
18 #endif
19 #ifdef FE_DOWNWARD
20 #define RD FE_DOWNWARD
21 #else
22 #define RD -1
23 #endif
24 #ifdef FE_UPWARD
25 #define RU FE_UPWARD
26 #else
27 #define RU -1
28 #endif
29
30 #undef INEXACT
31 #undef INVALID
32 #undef DIVBYZERO
33 #undef UNDERFLOW
34 #undef OVERFLOW
35 #ifdef FE_INEXACT
36 #define INEXACT FE_INEXACT
37 #else
38 #define INEXACT 0
39 #endif
40 #ifdef FE_INVALID
41 #define INVALID FE_INVALID
42 #else
43 #define INVALID 0
44 #endif
45 #ifdef FE_DIVBYZERO
46 #define DIVBYZERO FE_DIVBYZERO
47 #else
48 #define DIVBYZERO 0
49 #endif
50 #ifdef FE_UNDERFLOW
51 #define UNDERFLOW FE_UNDERFLOW
52 #else
53 #define UNDERFLOW 0
54 #endif
55 #ifdef FE_OVERFLOW
56 #define OVERFLOW FE_OVERFLOW
57 #else
58 #define OVERFLOW 0
59 #endif
60
61 #undef inf
62 #undef nan
63 #define inf INFINITY
64 #define nan NAN
65
66 #define T(...) {__FILE__, __LINE__, __VA_ARGS__},
67
68 #define POS char *file; int line;
69 struct d_d {POS int r; double x; double y; float dy; int e; };
70 struct f_f {POS int r; float x; float y; float dy; int e; };
71 struct l_l {POS int r; long double x; long double y; float dy; int e; };
72 struct ff_f {POS int r; float x; float x2; float y; float dy; int e; };
73 struct dd_d {POS int r; double x; double x2; double y; float dy; int e; };
74 struct ll_l {POS int r; long double x; long double x2; long double y; float dy; int e; };
75 struct d_di {POS int r; double x; double y; float dy; long long i; int e; };
76 struct f_fi {POS int r; float x; float y; float dy; long long i; int e; };
77 struct l_li {POS int r; long double x; long double y; float dy; long long i; int e; };
78 struct di_d {POS int r; double x; long long i; double y; float dy; int e; };
79 struct fi_f {POS int r; float x; long long i; float y; float dy; int e; };
80 struct li_l {POS int r; long double x; long long i; long double y; float dy; int e; };
81 struct d_i {POS int r; double x; long long i; int e; };
82 struct f_i {POS int r; float x; long long i; int e; };
83 struct l_i {POS int r; long double x; long long i; int e; };
84 struct d_dd {POS int r; double x; double y; float dy; double y2; float dy2; int e; };
85 struct f_ff {POS int r; float x; float y; float dy; float y2; float dy2; int e; };
86 struct l_ll {POS int r; long double x; long double y; float dy; long double y2; float dy2; int e; };
87 struct ff_fi {POS int r; float x; float x2; float y; float dy; long long i; int e; };
88 struct dd_di {POS int r; double x; double x2; double y; float dy; long long i; int e; };
89 struct ll_li {POS int r; long double x; long double x2; long double y; float dy; long long i; int e; };
90 struct fff_f {POS int r; float x; float x2; float x3; float y; float dy; int e; };
91 struct ddd_d {POS int r; double x; double x2; double x3; double y; float dy; int e; };
92 struct lll_l {POS int r; long double x; long double x2; long double x3; long double y; float dy; int e; };
93 #undef POS
94
95 char *estr(int);
96 char *rstr(int);
97
98 float ulperr(double got, double want, float dwant);
99 float ulperrf(float got, float want, float dwant);
100 float ulperrl(long double got, long double want, float dwant);
101
checkexcept(int got,int want,int r)102 static int checkexcept(int got, int want, int r)
103 {
104 if (r == RN)
105 #if defined CHECK_INEXACT
106 return got == want;
107 #elif defined CHECK_INEXACT_OMISSION
108 return got == want || got == (want|INEXACT);
109 #else
110 return (got|INEXACT) == (want|INEXACT);
111 #endif
112 return (got|INEXACT|UNDERFLOW) == (want|INEXACT|UNDERFLOW);
113 }
114
checkexceptall(int got,int want,int r)115 static int checkexceptall(int got, int want, int r)
116 {
117 return got == want;
118 }
119
checkulp(float d,int r)120 static int checkulp(float d, int r)
121 {
122 // TODO: we only care about >=1.5 ulp errors for now, should be 1.0
123 if (r == RN)
124 return fabsf(d) < 1.5;
125 // accept larger error in non-nearest rounding mode
126 return fabsf(d) < 3.0;
127 }
128
checkcr(long double y,long double ywant,int r)129 static int checkcr(long double y, long double ywant, int r)
130 {
131 if (isnan(ywant))
132 return isnan(y);
133 return y == ywant && signbit(y) == signbit(ywant);
134 }
135
136