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