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)						\
35TYPE c##EXT (TYPE x, TYPE y)					\
36{								\
37  return __builtin_copysign##EXT (x, y);			\
38}								\
39								\
40struct D##EXT { TYPE x, y, z; };				\
41								\
42static 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								\
54void 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
66TEST(float, f)
67TEST(double, )
68TEST(long double, l)
69
70int main()
71{
72  testf();
73  test();
74  testl();
75  return 0;
76}
77