| Current Path : /compat/linux/proc/self/root/usr/src/lib/msun/src/ |
FreeBSD hs32.drive.ne.jp 9.1-RELEASE FreeBSD 9.1-RELEASE #1: Wed Jan 14 12:18:08 JST 2015 root@hs32.drive.ne.jp:/sys/amd64/compile/hs32 amd64 |
| Current File : //compat/linux/proc/self/root/usr/src/lib/msun/src/s_logbf.c |
/* s_logbf.c -- float version of s_logb.c.
* Conversion to float by Ian Lance Taylor, Cygnus Support, ian@cygnus.com.
*/
/*
* ====================================================
* Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved.
*
* Developed at SunPro, a Sun Microsystems, Inc. business.
* Permission to use, copy, modify, and distribute this
* software is freely granted, provided that this notice
* is preserved.
* ====================================================
*/
#include <sys/cdefs.h>
__FBSDID("$FreeBSD: release/9.1.0/lib/msun/src/s_logbf.c 176451 2008-02-22 02:30:36Z das $");
#include "math.h"
#include "math_private.h"
static const float
two25 = 3.355443200e+07; /* 0x4c000000 */
float
logbf(float x)
{
int32_t ix;
GET_FLOAT_WORD(ix,x);
ix &= 0x7fffffff; /* high |x| */
if(ix==0) return (float)-1.0/fabsf(x);
if(ix>=0x7f800000) return x*x;
if(ix<0x00800000) {
x *= two25; /* convert subnormal x to normal */
GET_FLOAT_WORD(ix,x);
ix &= 0x7fffffff;
return (float) ((ix>>23)-127-25);
} else
return (float) ((ix>>23)-127);
}