|
static const static double double | one = 1.00000000000000000000e+00 |
|
static const static double double | pi = 3.14159265358979311600e+00 |
|
static const static double double | pio2_hi = 1.57079632679489655800e+00 |
|
static const static double double | pio2_lo = 6.12323399573676603587e-17 |
|
static const static double double | pS0 = 1.66666666666666657415e-01 |
|
static const static double double | pS1 = -3.25565818622400915405e-01 |
|
static const static double double | pS2 = 2.01212532134862925881e-01 |
|
static const static double double | pS3 = -4.00555345006794114027e-02 |
|
static const static double double | pS4 = 7.91534994289814532176e-04 |
|
static const static double double | pS5 = 3.47933107596021167570e-05 |
|
static const static double double | qS1 = -2.40339491173441421878e+00 |
|
static const static double double | qS2 = 2.02094576023350569471e+00 |
|
static const static double double | qS3 = -6.88283971605453293030e-01 |
|
static const static double double | qS4 = 7.70381505559019352791e-02 |
|
int | hx |
|
int | ix = hx&0x7fffffff |
|
| else |
|
| s = sqrt(z) |
|
| df = s |
|
| c = (z-df*df)/(s+df) |
|
| p = z*(pS0+z*(pS1+z*(pS2+z*(pS3+z*(pS4+z*pS5))))) |
|
| q = one+z*(qS1+z*(qS2+z*(qS3+z*qS4))) |
|
| r = p/q |
|
| w = r*s+c |
|