1 /* @(#)s_copysign.c 5.1 93/09/24 */
2 /*
3  * ====================================================
4  * Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved.
5  *
6  * Developed at SunPro, a Sun Microsystems, Inc. business.
7  * Permission to use, copy, modify, and distribute this
8  * software is freely granted, provided that this notice
9  * is preserved.
10  * ====================================================
11  */
12 
13 /*
14  * copysign(double x, double y)
15  * copysign(x,y) returns a value with the magnitude of x and
16  * with the sign bit of y.
17  */
18 
19 #include <float.h>
20 #include <math.h>
21 
22 #include "math_private.h"
23 
24 double
25 copysign(double x, double y)
26 {
27 	u_int32_t hx,hy;
28 	GET_HIGH_WORD(hx,x);
29 	GET_HIGH_WORD(hy,y);
30 	SET_HIGH_WORD(x,(hx&0x7fffffff)|(hy&0x80000000));
31         return x;
32 }
33 
34 #if	LDBL_MANT_DIG == DBL_MANT_DIG
35 __strong_alias(copysignl, copysign);
36 #endif	/* LDBL_MANT_DIG == DBL_MANT_DIG */
37