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
|
#include <stdio.h>
#include <fenv.h>
#include <float.h>
#pragma STDC FENV_ACCESS ON
#if defined(__K1C__) && !defined(__COMPCERT__)
int fetestexcept(int excepts) {
int mask = (K1_SFR_CS_IO_MASK | K1_SFR_CS_DZ_MASK | K1_SFR_CS_OV_MASK | K1_SFR_CS_UN_MASK | K1_SFR_CS_IN_MASK) & excepts;
unsigned long long cs = __builtin_k1_get(K1_SFR_CS);
return cs & mask;
}
int feclearexcept(int excepts) {
int mask = (K1_SFR_CS_IO_MASK | K1_SFR_CS_DZ_MASK | K1_SFR_CS_OV_MASK | K1_SFR_CS_UN_MASK | K1_SFR_CS_IN_MASK) & excepts;
__builtin_k1_wfxl(K1_SFR_CS, mask);
return 0;
}
#endif
double add(double x, double y) {
return x+y;
}
double mul(double x, double y) {
return x*y;
}
float double2float(double x) {
return x;
}
float uint2float(unsigned x) {
return x;
}
double ulong2double(unsigned long x) {
return x;
}
unsigned double2uint(double x) {
return x;
}
int main() {
printf("%x\n", fetestexcept(FE_ALL_EXCEPT));
double v1 = add(3.0, 0.1);
printf("%x\n", fetestexcept(FE_ALL_EXCEPT));
feclearexcept(FE_INEXACT);
printf("%x\n", fetestexcept(FE_ALL_EXCEPT));
double v2 = mul(DBL_MAX, DBL_MAX);
printf("%g %x\n", v2, fetestexcept(FE_ALL_EXCEPT));
feclearexcept(FE_ALL_EXCEPT);
double v3 = mul(DBL_MIN, DBL_MIN);
printf("%g %x\n", v3, fetestexcept(FE_ALL_EXCEPT));
feclearexcept(FE_ALL_EXCEPT);
double v4 = double2float(DBL_MAX);
printf("%g %x\n", v4, fetestexcept(FE_ALL_EXCEPT));
feclearexcept(FE_ALL_EXCEPT);
float v5 = uint2float(0xC07FDFFFU);
printf("%g %x\n", v5, fetestexcept(FE_ALL_EXCEPT)); // BUG 0 should have INEXACT
feclearexcept(FE_ALL_EXCEPT);
double v6 = ulong2double(0x11217763AFF77C7CUL);
printf("%g %x\n", v6, fetestexcept(FE_ALL_EXCEPT)); // BUG 0 should have INEXACT
feclearexcept(FE_ALL_EXCEPT);
unsigned v7 = double2uint(-0.25); // softfloat says "0 and inexact" but here we have "0 and overflow" (due to negative input for unsigned?)
printf("%u %x\n", v7, fetestexcept(FE_ALL_EXCEPT));
feclearexcept(FE_ALL_EXCEPT);
// +41F.307672C5496EF
double d8 = 0x1.307672C5496EFp32;
unsigned v8 = double2uint(d8);
printf("%g %x %x\n", d8, v8, fetestexcept(FE_ALL_EXCEPT));
// BUG reports 307672C5 and inexact, but should report overflow
// bug comes from x86 https://gcc.gnu.org/bugzilla/show_bug.cgi?id=89175
feclearexcept(FE_ALL_EXCEPT);
}
|