1 /* s_copysignl.c -- long double version of s_copysign.c.
2 * Conversion to long double by Jakub Jelinek, jj@ultra.linux.cz.
3 */
4
5 /*
6 * ====================================================
7 * Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved.
8 *
9 * Developed at SunPro, a Sun Microsystems, Inc. business.
10 * Permission to use, copy, modify, and distribute this
11 * software is freely granted, provided that this notice
12 * is preserved.
13 * ====================================================
14 */
15
16 #if defined(LIBM_SCCS) && !defined(lint)
17 static char rcsid[] = "NetBSD: ";
18 #endif
19
20 /*
21 * copysignq(long double x, long double y)
22 * copysignq(x,y) returns a value with the magnitude of x and
23 * with the sign bit of y.
24 */
25
26 #define NO_MATH_REDIRECT
27
28 #include "quadmath-imp.h"
29
copysignq(__float128 x,__float128 y)30 __float128 copysignq(__float128 x, __float128 y)
31 {
32 uint64_t hx,hy;
33 GET_FLT128_MSW64(hx,x);
34 GET_FLT128_MSW64(hy,y);
35 SET_FLT128_MSW64(x,(hx&0x7fffffffffffffffULL)
36 |(hy&0x8000000000000000ULL));
37 return x;
38 }
39