1 /* @(#)s_copysign.c 5.1 93/09/24 */
3 * ====================================================
4 * Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved.
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
10 * ====================================================
13 #if defined(LIBM_SCCS) && !defined(lint)
14 static char rcsid
[] = "$NetBSD: s_copysign.c,v 1.8 1995/05/10 20:46:57 jtc Exp $";
18 * copysign(double x, double y)
19 * copysign(x,y) returns a value with the magnitude of x and
20 * with the sign bit of y.
24 #include <math_private.h>
25 #include <libm-alias-double.h>
28 __copysign (double x
, double y
)
31 GET_HIGH_WORD (hx
, x
);
32 GET_HIGH_WORD (hy
, y
);
33 SET_HIGH_WORD (x
, (hx
& 0x7fffffff) | (hy
& 0x80000000));
36 libm_alias_double (__copysign
, copysign
)