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
|
#include<stdio.h>
#include<stdlib.h>
#define STR_EXPAND(tok) #tok
#define STR(tok) STR_EXPAND(tok)
#if defined(__ppc__) || defined(__PPC__)
#define ARCH_BIG_ENDIAN
#elif defined(__i386__) || defined(__x86_64__) || defined(__ARMEL__)
#undef ARCH_BIG_ENDIAN
#else
#error "unknown endianness"
#endif
union converter64 {
double dbl;
struct {
#ifdef ARCH_BIG_ENDIAN
unsigned h, l;
#else
unsigned l, h;
#endif
} u64;
};
int num_errors = 0;
void comp64(double d, unsigned high, unsigned low, char* s) {
union converter64 c;
c.dbl = d;
if((c.u64.h & 0x7FF00000) == 0x7FF00000 && (c.u64.l | (c.u64.h & 0xFFFFF)) != 0) {
c.u64.l = 0xFFFFFFFF;
c.u64.h = 0x7FFFFFFF;
}
if((high & 0x7FF00000) == 0x7FF00000 && (low | (high & 0xFFFFF)) != 0) {
low = 0xFFFFFFFF;
high = 0x7FFFFFFF;
}
if(low != c.u64.l || high != c.u64.h) {
printf("%s : %08x %08x (%a)\n", s, c.u64.h, c.u64.l, d);
num_errors++;
}
}
void compd(double d1, double d2, char* s) {
union converter64 c1, c2;
c1.dbl = d1;
c2.dbl = d2;
if(c1.u64.l != c2.u64.l || c1.u64.h != c2.u64.h) {
printf("%s : %a %a\n", s, d1, d2);
num_errors++;
}
}
int main(void) {
comp64(3.518437208883201171875E+013, 0x42c00000, 0x00000002, STR(__LINE__));
compd(1.50000000000000011102230246251565404236316680908203125, 0x1.8p+0, STR(__LINE__));
compd(0.500000000000000166533453693773481063544750213623046875, 0x1.0000000000002p-1, STR(__LINE__));
compd(1.2, 1.20000000000000005, STR(__LINE__));
compd(1.2f, 1.2000001f, STR(__LINE__));
compd(1, 1., STR(__LINE__));
compd(0.5, .5, STR(__LINE__));
compd(1, 1E0, STR(__LINE__));
compd(0.5, 0x1p-1, STR(__LINE__));
compd(0.5, 0x1.p-1, STR(__LINE__));
compd(0.5, 0x.1p3, STR(__LINE__));
compd(0.5, 0x.8p0, STR(__LINE__));
compd(15./16, 0x.fp0, STR(__LINE__));
compd(15./16, 0x.Fp0, STR(__LINE__));
compd(15./16, 0x.fP0, STR(__LINE__));
compd(15./16, 0X.fp0, STR(__LINE__));
printf("%d error(s) detected.\n", num_errors);
return 0;
}
|
