1 #include <string.h>
2 #include <stdlib.h>
3 #include <float.h>
4
5 #define fpsizeoff sizeof(float)
6 #define fpsizeof sizeof(double)
7 #define fpsizeofl sizeof(long double)
8
9 /* Work around the fact that with the Intel double-extended precision,
10 we've got a 10 byte type stuffed into some amount of padding. And
11 the fact that -ffloat-store is going to stuff this value temporarily
12 into some bit of stack frame that we've no control over and can't zero. */
13 #if LDBL_MANT_DIG == 64
14 # if defined(__i386__) || defined(__x86_64__) || defined (__ia64__)
15 # undef fpsizeofl
16 # define fpsizeofl 10
17 # endif
18 #endif
19
20 /* Work around the fact that the sign of the second double in the IBM
21 double-double format is not strictly specified when it contains a zero.
22 For instance, -0.0L can be represented with either (-0.0, +0.0) or
23 (-0.0, -0.0). The former is what we'll get from the compiler when it
24 builds constants; the later is what we'll get from the negation operator
25 at runtime. */
26 /* ??? This hack only works for big-endian, which is fortunately true for
27 AIX and, Darwin. */
28 #if LDBL_MANT_DIG == 106
29 # undef fpsizeofl
30 # define fpsizeofl sizeof(double)
31 #endif
32
33
34 #define TEST(TYPE, EXT) \
35 TYPE c##EXT (TYPE x, TYPE y) \
36 { \
37 return __builtin_copysign##EXT (x, y); \
38 } \
39 \
40 struct D##EXT { TYPE x, y, z; }; \
41 \
42 static const struct D##EXT T##EXT[] = { \
43 { 1.0, 2.0, 1.0 }, \
44 { 1.0, -2.0, -1.0 }, \
45 { -1.0, -2.0, -1.0 }, \
46 { 0.0, -2.0, -0.0 }, \
47 { -0.0, -2.0, -0.0 }, \
48 { -0.0, 2.0, 0.0 }, \
49 { __builtin_inf##EXT (), -0.0, -__builtin_inf##EXT () }, \
50 { -__builtin_nan##EXT (""), __builtin_inf##EXT (), \
51 __builtin_nan##EXT ("") } \
52 }; \
53 \
54 void test##EXT (void) \
55 { \
56 int i, n = sizeof (T##EXT) / sizeof (T##EXT[0]); \
57 TYPE r; \
58 for (i = 0; i < n; ++i) \
59 { \
60 r = c##EXT (T##EXT[i].x, T##EXT[i].y); \
61 if (memcmp (&r, &T##EXT[i].z, fpsizeof##EXT) != 0) \
62 abort (); \
63 } \
64 }
65
TEST(float,f)66 TEST(float, f)
67 TEST(double, )
68 TEST(long double, l)
69
70 int main()
71 {
72 testf();
73 test();
74 testl();
75 return 0;
76 }
77