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
|
/* { dg-do run { target { powerpc*-*-linux* } } } */
/* { dg-skip-if "" { powerpc*-*-darwin* } { "*" } { "" } } */
/* { dg-skip-if "" { powerpc*-*-*spe* } { "*" } { "" } } */
/* { dg-require-effective-target powerpc_fprs } */
/* { dg-require-effective-target longdouble128 } */
/* { dg-options "-O2 -mhard-float" } */
#include <stddef.h>
#include <stdlib.h>
#include <math.h>
#ifdef DEBUG
#include <stdio.h>
#endif
int
main (void)
{
double high = pow (2.0, 60);
double low = 2.0;
long double a = ((long double)high) + ((long double)low);
double x0 = __builtin_unpack_longdouble (a, 0);
double x1 = __builtin_unpack_longdouble (a, 1);
long double b = __builtin_pack_longdouble (x0, x1);
#ifdef DEBUG
{
size_t i;
union {
long double ld;
double d;
unsigned char uc[sizeof (long double)];
char c[sizeof (long double)];
} u;
printf ("a = 0x");
u.ld = a;
for (i = 0; i < sizeof (long double); i++)
printf ("%.2x", u.uc[i]);
printf (", %Lg\n", a);
printf ("b = 0x");
u.ld = b;
for (i = 0; i < sizeof (long double); i++)
printf ("%.2x", u.uc[i]);
printf (", %Lg\n", b);
printf ("hi = 0x");
u.d = high;
for (i = 0; i < sizeof (double); i++)
printf ("%.2x", u.uc[i]);
printf (",%*s %g\n", (int)(2 * (sizeof (long double) - sizeof (double))), "", high);
printf ("lo = 0x");
u.d = low;
for (i = 0; i < sizeof (double); i++)
printf ("%.2x", u.uc[i]);
printf (",%*s %g\n", (int)(2 * (sizeof (long double) - sizeof (double))), "", low);
printf ("x0 = 0x");
u.d = x0;
for (i = 0; i < sizeof (double); i++)
printf ("%.2x", u.uc[i]);
printf (",%*s %g\n", (int)(2 * (sizeof (long double) - sizeof (double))), "", x0);
printf ("x1 = 0x");
u.d = x1;
for (i = 0; i < sizeof (double); i++)
printf ("%.2x", u.uc[i]);
printf (",%*s %g\n", (int)(2 * (sizeof (long double) - sizeof (double))), "", x1);
}
#endif
if (high != x0)
abort ();
if (low != x1)
abort ();
if (a != b)
abort ();
if (x0 != high)
abort ();
if (x1 != low)
abort ();
return 0;
}
|
