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
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
|
// gcc -Wall -W -lm double-to-float-with-rounding-modes.c -o double-to-float-with-rounding-modes -frounding-math -fsignaling-nans -ffp-contract=off -msse2 -mfpmath=sse
#ifdef __GNUC__
#include <assert.h>
#include <math.h>
#include <fenv.h>
float setRoundingModeAndCast(int mode, double d) {
fesetround(mode);
return (float)d;
}
void test (int mode, double d, float result) {
float f = setRoundingModeAndCast(mode,d);
assert(f == result);
return;
}
#endif
int main (void)
{
#ifdef __GNUC__
// 0x1.fffffep+127f is the largest binary32 float so
// these should all be representable as floats...
test(FE_TONEAREST, 0x1.fffffep+127, 0x1.fffffep+127);
#ifdef FE_UPWARD
test(FE_UPWARD, 0x1.fffffep+127, 0x1.fffffep+127);
#endif
#ifdef FE_DOWNWARD
test(FE_DOWNWARD, 0x1.fffffep+127, 0x1.fffffep+127);
#endif
test(FE_TOWARDZERO, 0x1.fffffep+127, 0x1.fffffep+127);
// Likewise, this should be an obvious sanity check
test(FE_TONEAREST, +INFINITY, +INFINITY);
#ifdef FE_UPWARD
test(FE_UPWARD, +INFINITY, +INFINITY);
#endif
#ifdef FE_DOWNWARD
test(FE_DOWNWARD, +INFINITY, +INFINITY);
#endif
test(FE_TOWARDZERO, +INFINITY, +INFINITY);
// Nearer to 0x1.fffffep+127 than to 0x1.000000p+128
test(FE_TONEAREST, 0x1.fffffe0000001p+127, 0x1.fffffep+127);
#ifdef FE_UPWARD
test(FE_UPWARD, 0x1.fffffe0000001p+127, +INFINITY);
#endif
#ifdef FE_DOWNWARD
test(FE_DOWNWARD, 0x1.fffffe0000001p+127, 0x1.fffffep+127);
#endif
test(FE_TOWARDZERO, 0x1.fffffe0000001p+127, 0x1.fffffep+127);
// 0x1.fffffefffffffp+127 is immediately below half way
test(FE_TONEAREST, 0x1.fffffefffffffp+127, 0x1.fffffep+127);
#ifdef FE_UPWARD
test(FE_UPWARD, 0x1.fffffefffffffp+127, +INFINITY);
#endif
#ifdef FE_DOWNWARD
test(FE_DOWNWARD, 0x1.fffffefffffffp+127, 0x1.fffffep+127);
#endif
test(FE_TOWARDZERO, 0x1.fffffefffffffp+127, 0x1.fffffep+127);
// Half way
test(FE_TONEAREST, 0x1.ffffffp+127, +INFINITY);
#ifdef FE_UPWARD
test(FE_UPWARD, 0x1.ffffffp+127, +INFINITY);
#endif
#ifdef FE_DOWNWARD
test(FE_DOWNWARD, 0x1.ffffffp+127, 0x1.fffffep+127);
#endif
test(FE_TOWARDZERO, 0x1.ffffffp+127, 0x1.fffffep+127);
// Larger
test(FE_TONEAREST, 0x1.0p+128, +INFINITY);
#ifdef FE_UPWARD
test(FE_UPWARD, 0x1.0p+128, +INFINITY);
#endif
#ifdef FE_DOWNWARD
test(FE_DOWNWARD, 0x1.0p+128, 0x1.fffffep+127);
#endif
test(FE_TOWARDZERO, 0x1.0p+128, 0x1.fffffep+127);
// Huge
test(FE_TONEAREST, 0x1.fffffffffffffp+1023, +INFINITY);
#ifdef FE_UPWARD
test(FE_UPWARD, 0x1.fffffffffffffp+1023, +INFINITY);
#endif
#ifdef FE_DOWNWARD
test(FE_DOWNWARD, 0x1.fffffffffffffp+1023, 0x1.fffffep+127);
#endif
test(FE_TOWARDZERO, 0x1.fffffffffffffp+1023, 0x1.fffffep+127);
// Same again but negative
test(FE_TONEAREST, -0x1.fffffep+127, -0x1.fffffep+127);
#ifdef FE_UPWARD
test(FE_UPWARD, -0x1.fffffep+127, -0x1.fffffep+127);
#endif
#ifdef FE_DOWNWARD
test(FE_DOWNWARD, -0x1.fffffep+127, -0x1.fffffep+127);
#endif
test(FE_TOWARDZERO, -0x1.fffffep+127, -0x1.fffffep+127);
test(FE_TONEAREST, -INFINITY, -INFINITY);
#ifdef FE_UPWARD
test(FE_UPWARD, -INFINITY, -INFINITY);
#endif
#ifdef FE_DOWNWARD
test(FE_DOWNWARD, -INFINITY, -INFINITY);
#endif
test(FE_TOWARDZERO, -INFINITY, -INFINITY);
test(FE_TONEAREST, -0x1.fffffe0000001p+127, -0x1.fffffep+127);
#ifdef FE_UPWARD
test(FE_UPWARD, -0x1.fffffe0000001p+127, -0x1.fffffep+127);
#endif
#ifdef FE_DOWNWARD
test(FE_DOWNWARD, -0x1.fffffe0000001p+127, -INFINITY);
#endif
test(FE_TOWARDZERO, -0x1.fffffe0000001p+127, -0x1.fffffep+127);
test(FE_TONEAREST, -0x1.fffffefffffffp+127, -0x1.fffffep+127);
#ifdef FE_UPWARD
test(FE_UPWARD, -0x1.fffffefffffffp+127, -0x1.fffffep+127);
#endif
#ifdef FE_DOWNWARD
test(FE_DOWNWARD, -0x1.fffffefffffffp+127, -INFINITY);
#endif
test(FE_TOWARDZERO, -0x1.fffffefffffffp+127, -0x1.fffffep+127);
test(FE_TONEAREST, -0x1.ffffffp+127, -INFINITY);
#ifdef FE_UPWARD
test(FE_UPWARD, -0x1.ffffffp+127, -0x1.fffffep+127);
#endif
#ifdef FE_DOWNWARD
test(FE_DOWNWARD, -0x1.ffffffp+127, -INFINITY);
#endif
test(FE_TOWARDZERO, -0x1.ffffffp+127, -0x1.fffffep+127);
test(FE_TONEAREST, -0x1.0p+128, -INFINITY);
#ifdef FE_UPWARD
test(FE_UPWARD, -0x1.0p+128, -0x1.fffffep+127);
#endif
#ifdef FE_DOWNWARD
test(FE_DOWNWARD, -0x1.0p+128, -INFINITY);
#endif
test(FE_TOWARDZERO, -0x1.0p+128, -0x1.fffffep+127);
test(FE_TONEAREST, -0x1.fffffffffffffp+1023, -INFINITY);
#ifdef FE_UPWARD
test(FE_UPWARD, -0x1.fffffffffffffp+1023, -0x1.fffffep+127);
#endif
#ifdef FE_DOWNWARD
test(FE_DOWNWARD, -0x1.fffffffffffffp+1023, -INFINITY);
#endif
test(FE_TOWARDZERO, -0x1.fffffffffffffp+1023, -0x1.fffffep+127);
#endif
return 1;
}
|
