]>
git.proxmox.com Git - mirror_edk2.git/blob - StdLib/LibC/Math/s_modf.c
1 /* @(#)s_modf.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 * ====================================================
12 #include <LibConfig.h>
13 #include <sys/EfiCdefs.h>
14 #if defined(LIBM_SCCS) && !defined(lint)
15 __RCSID("$NetBSD: s_modf.c,v 1.11 2002/05/26 22:01:57 wiz Exp $");
19 * modf(double x, double *iptr)
20 * return fraction part of x, and return x's integral part in *iptr.
29 #include "math_private.h"
31 static const double one
= 1.0;
34 modf(double x
, double *iptr
)
38 EXTRACT_WORDS(i0
,i1
,x
);
39 j0
= ((i0
>>20)&0x7ff)-0x3ff; /* exponent of x */
40 if(j0
<20) { /* integer part in high x */
41 if(j0
<0) { /* |x|<1 */
42 INSERT_WORDS(*iptr
,i0
&0x80000000,0); /* *iptr = +-0 */
46 if(((i0
&i
)|i1
)==0) { /* x is integral */
49 GET_HIGH_WORD(high
,x
);
50 INSERT_WORDS(x
,high
&0x80000000,0); /* return +-0 */
53 INSERT_WORDS(*iptr
,i0
&(~i
),0);
57 } else if (j0
>51) { /* no fraction part */
60 GET_HIGH_WORD(high
,x
);
61 INSERT_WORDS(x
,high
&0x80000000,0); /* return +-0 */
63 } else { /* fraction part in low x */
64 i
= ((u_int32_t
)(0xffffffff))>>(j0
-20);
65 if((i1
&i
)==0) { /* x is integral */
68 GET_HIGH_WORD(high
,x
);
69 INSERT_WORDS(x
,high
&0x80000000,0); /* return +-0 */
72 INSERT_WORDS(*iptr
,i0
,i1
&(~i
));