version 1.7, 1999/07/02 15:37:40 |
version 1.8, 2002/05/26 22:01:51 |
Line 21 __RCSID("$NetBSD$"); |
|
Line 21 __RCSID("$NetBSD$"); |
|
#include "math.h" |
#include "math.h" |
#include "math_private.h" |
#include "math_private.h" |
|
|
#ifdef __STDC__ |
|
static const float |
static const float |
#else |
|
static float |
|
#endif |
|
ln2_hi = 6.9313812256e-01, /* 0x3f317180 */ |
ln2_hi = 6.9313812256e-01, /* 0x3f317180 */ |
ln2_lo = 9.0580006145e-06, /* 0x3717f7d1 */ |
ln2_lo = 9.0580006145e-06, /* 0x3717f7d1 */ |
two25 = 3.355443200e+07, /* 0x4c000000 */ |
two25 = 3.355443200e+07, /* 0x4c000000 */ |
Line 37 Lg5 = 1.8183572590e-01, /* 3E3A3325 */ |
|
Line 33 Lg5 = 1.8183572590e-01, /* 3E3A3325 */ |
|
Lg6 = 1.5313838422e-01, /* 3E1CD04F */ |
Lg6 = 1.5313838422e-01, /* 3E1CD04F */ |
Lg7 = 1.4798198640e-01; /* 3E178897 */ |
Lg7 = 1.4798198640e-01; /* 3E178897 */ |
|
|
#ifdef __STDC__ |
|
static const float zero = 0.0; |
static const float zero = 0.0; |
#else |
|
static float zero = 0.0; |
|
#endif |
|
|
|
#ifdef __STDC__ |
float |
float __ieee754_logf(float x) |
__ieee754_logf(float x) |
#else |
|
float __ieee754_logf(x) |
|
float x; |
|
#endif |
|
{ |
{ |
float hfsq,f,s,z,R,w,t1,t2,dk; |
float hfsq,f,s,z,R,w,t1,t2,dk; |
int32_t k,ix,i,j; |
int32_t k,ix,i,j; |