version 1.10, 1999/07/02 15:37:42 |
version 1.11, 2002/05/26 22:01:55 |
Line 113 __RCSID("$NetBSD$"); |
|
Line 113 __RCSID("$NetBSD$"); |
|
#include "math.h" |
#include "math.h" |
#include "math_private.h" |
#include "math_private.h" |
|
|
#ifdef __STDC__ |
|
static const double |
static const double |
#else |
|
static double |
|
#endif |
|
tiny = 1e-300, |
tiny = 1e-300, |
half= 5.00000000000000000000e-01, /* 0x3FE00000, 0x00000000 */ |
half= 5.00000000000000000000e-01, /* 0x3FE00000, 0x00000000 */ |
one = 1.00000000000000000000e+00, /* 0x3FF00000, 0x00000000 */ |
one = 1.00000000000000000000e+00, /* 0x3FF00000, 0x00000000 */ |
Line 192 sb5 = 2.55305040643316442583e+03, /* 0 |
|
Line 188 sb5 = 2.55305040643316442583e+03, /* 0 |
|
sb6 = 4.74528541206955367215e+02, /* 0x407DA874, 0xE79FE763 */ |
sb6 = 4.74528541206955367215e+02, /* 0x407DA874, 0xE79FE763 */ |
sb7 = -2.24409524465858183362e+01; /* 0xC03670E2, 0x42712D62 */ |
sb7 = -2.24409524465858183362e+01; /* 0xC03670E2, 0x42712D62 */ |
|
|
#ifdef __STDC__ |
double |
double erf(double x) |
erf(double x) |
#else |
|
double erf(x) |
|
double x; |
|
#endif |
|
{ |
{ |
int32_t hx,ix,i; |
int32_t hx,ix,i; |
double R,S,P,Q,s,y,z,r; |
double R,S,P,Q,s,y,z,r; |
Line 248 sb7 = -2.24409524465858183362e+01; /* 0 |
|
Line 240 sb7 = -2.24409524465858183362e+01; /* 0 |
|
if(hx>=0) return one-r/x; else return r/x-one; |
if(hx>=0) return one-r/x; else return r/x-one; |
} |
} |
|
|
#ifdef __STDC__ |
double |
double erfc(double x) |
erfc(double x) |
#else |
|
double erfc(x) |
|
double x; |
|
#endif |
|
{ |
{ |
int32_t hx,ix; |
int32_t hx,ix; |
double R,S,P,Q,s,y,z,r; |
double R,S,P,Q,s,y,z,r; |