|
static const static double double | one = 1.0 |
|
static const static double double | halF [2] = {0.5,-0.5,} |
|
static const static double double | huge = 1.0e+300 |
|
static const static double double | twom1000 = 9.33263618503218878990e-302 |
|
static const static double double | o_threshold = 7.09782712893383973096e+02 |
|
static const static double double | u_threshold = -7.45133219101941108420e+02 |
|
static const static double double | ln2HI [2] |
|
static const static double double | ln2LO [2] |
|
static const static double double | invln2 = 1.44269504088896338700e+00 |
|
static const static double double | P1 = 1.66666666666666019037e-01 |
|
static const static double double | P2 = -2.77777777770155933842e-03 |
|
static const static double double | P3 = 6.61375632143793436117e-05 |
|
static const static double double | P4 = -1.65339022054652515390e-06 |
|
static const static double double | P5 = 4.13813679705723846039e-08 |
|
int | k |
|
int | xsb = (hx >> 31) & 1 |
|
unsigned | hx = __HI(x) |
|
| t = x * x |
|
| c = x - t * (P1 + t * (P2 + t * (P3 + t * (P4 + t * P5)))) |
|
| else |
|