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
|
/* { dg-do run { target powerpc*-*-darwin* } } */
/* { dg-options "" } */
/* No options so 'long long' can be used. */
#include <stdio.h>
#include <stdlib.h>
#include <string.h>
typedef unsigned long long uint64_t;
typedef uint64_t ldbits[2];
union ldu
{
ldbits lb;
long double ld;
};
static const struct {
ldbits a;
ldbits b;
ldbits result;
} single_tests[] = {
/* Test of values that add to near +Inf. */
{ { 0x7FEFFFFFFFFFFFFFLL, 0xFC88000000000000LL },
{ 0x7C94000000000000LL, 0x0000000000000000LL },
{ 0x7FEFFFFFFFFFFFFFLL, 0x7C80000000000000LL } },
{ { 0x7FEFFFFFFFFFFFFFLL, 0x7C8FFFFFFFFFFFFFLL },
{ 0x792FFFFFFFFFFFFFLL, 0x0000000000000000LL },
{ 0x7FEFFFFFFFFFFFFFLL, 0x7C8FFFFFFFFFFFFFLL } },
{ { 0x7FEFFFFFFFFFFFFFLL, 0x7C8FFFFFFFFFFFFFLL },
{ 0x7930000000000000LL, 0xF5DFFFFFFFFFFFFFLL },
/* correct result is: { 0x7FEFFFFFFFFFFFFFLL, 0x7C8FFFFFFFFFFFFFLL } */
{ 0x7FF0000000000000LL, 0x0000000000000000LL } },
/* Test of values that add to +Inf. */
{ { 0x7FEFFFFFFFFFFFFFLL, 0x7C8FFFFFFFFFFFFFLL },
{ 0x7930000000000000LL, 0x0000000000000000LL },
{ 0x7FF0000000000000LL, 0x0000000000000000LL } },
/* Tests of Inf addition. */
{ { 0x7FF0000000000000LL, 0x0000000000000000LL },
{ 0x0000000000000000LL, 0x0000000000000000LL },
{ 0x7FF0000000000000LL, 0x0000000000000000LL } },
{ { 0x7FF0000000000000LL, 0x0000000000000000LL },
{ 0x7FF0000000000000LL, 0x0000000000000000LL },
{ 0x7FF0000000000000LL, 0x0000000000000000LL } },
/* Test of Inf addition producing NaN. */
{ { 0x7FF0000000000000LL, 0x0000000000000000LL },
{ 0xFFF0000000000000LL, 0x0000000000000000LL },
{ 0x7FF8000000000000LL, 0x0000000000000000LL } },
/* Tests of NaN addition. */
{ { 0x7FF8000000000000LL, 0x0000000000000000LL },
{ 0x0000000000000000LL, 0x0000000000000000LL },
{ 0x7FF8000000000000LL, 0x7FF8000000000000LL } },
{ { 0x7FF8000000000000LL, 0x0000000000000000LL },
{ 0x7FF0000000000000LL, 0x0000000000000000LL },
{ 0x7FF8000000000000LL, 0x7FF8000000000000LL } },
/* Addition of positive integers, with interesting rounding properties. */
{ { 0x4690000000000000LL, 0x4330000000000000LL },
{ 0x4650000000000009LL, 0xC2FFFFFFFFFFFFF2LL },
/* correct result is: { 0x4691000000000001LL, 0xC32C000000000000LL } */
{ 0x4691000000000001LL, 0xc32bfffffffffffeLL } },
{ { 0x4690000000000000LL, 0x4330000000000000LL },
{ 0x4650000000000008LL, 0x42F0000000000010LL },
{ 0x4691000000000001LL, 0xC32E000000000000LL } },
{ { 0x469FFFFFFFFFFFFFLL, 0x433FFFFFFFFFFFFFLL },
{ 0x4340000000000000LL, 0x3FF0000000000000LL },
{ 0x46A0000000000000LL, 0x0000000000000000LL } },
{ { 0x469FFFFFFFFFFFFFLL, 0x433FFFFFFFFFFFFFLL },
{ 0x4340000000000000LL, 0x0000000000000000LL },
{ 0x46A0000000000000LL, 0xBFF0000000000000LL } },
/* Subtraction of integers, with cancellation. */
{ { 0x4690000000000000LL, 0x4330000000000000LL },
{ 0xC690000000000000LL, 0xC330000000000000LL },
{ 0x0000000000000000LL, 0x0000000000000000LL } },
{ { 0x4690000000000000LL, 0x4330000000000000LL },
{ 0xC330000000000000LL, 0x0000000000000000LL },
{ 0x4690000000000000LL, 0x0000000000000000LL } },
{ { 0x4690000000000000LL, 0x4330000000000000LL },
{ 0xC330000000000000LL, 0x3FA0000000000000LL },
{ 0x4690000000000000LL, 0x3FA0000000000000LL } },
{ { 0x4690000000000000LL, 0x4330000000000000LL },
{ 0xC690000000000000LL, 0x3FA0000000000000LL },
/* correct result is: { 0x4330000000000000LL, 0x3FA0000000000000LL } */
{ 0x4330000000000000LL, 0x0000000000000000LL } }
};
static int fail = 0;
static void
run_single_tests (void)
{
size_t i;
for (i = 0; i < sizeof (single_tests) / sizeof (single_tests[0]); i++)
{
union ldu a, b, result, expected;
memcpy (a.lb, single_tests[i].a, sizeof (ldbits));
memcpy (b.lb, single_tests[i].b, sizeof (ldbits));
memcpy (expected.lb, single_tests[i].result, sizeof (ldbits));
result.ld = a.ld + b.ld;
if (memcmp (result.lb, expected.lb,
result.ld == result.ld ? sizeof (ldbits) : sizeof (double))
!= 0)
{
printf ("FAIL: %016llx %016llx + %016llx %016llx\n",
a.lb[0], a.lb[1], b.lb[0], b.lb[1]);
printf (" = %016llx %016llx not %016llx %016llx\n",
result.lb[0], result.lb[1], expected.lb[0], expected.lb[1]);
fail = 1;
}
}
}
int main(void)
{
run_single_tests();
if (fail)
abort ();
else
exit (0);
}
|
