1 /* s_frexpf.c -- float version of s_frexp.c.
5 * ====================================================
6 * Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved.
8 * Developed at SunPro, a Sun Microsystems, Inc. business.
9 * Permission to use, copy, modify, and distribute this
10 * software is freely granted, provided that this notice
12 * ====================================================
15 #if defined(LIBM_SCCS) && !defined(lint)
16 static char rcsid
[] = "$NetBSD: s_frexpf.c,v 1.5 1995/05/10 20:47:26 jtc Exp $";
20 #include <math_private.h>
21 #include <libm-alias-float.h>
24 two25
= 3.3554432000e+07; /* 0x4c000000 */
26 float __frexpf(float x
, int *eptr
)
32 if(ix
>=0x7f800000||(ix
==0)) return x
+ x
; /* 0,inf,nan */
33 if (ix
<0x00800000) { /* subnormal */
39 *eptr
+= (ix
>>23)-126;
40 hx
= (hx
&0x807fffff)|0x3f000000;
44 libm_alias_float (__frexp
, frexp
)