|
| static const static double double | ln2_hi = 6.93147180369123816490e-01 |
| |
| static const static double double | ln2_lo = 1.90821492927058770002e-10 |
| |
| static const static double double | two54 = 1.80143985094819840000e+16 |
| |
| static const static double double | Lg1 = 6.666666666666735130e-01 |
| |
| static const static double double | Lg2 = 3.999999999940941908e-01 |
| |
| static const static double double | Lg3 = 2.857142874366239149e-01 |
| |
| static const static double double | Lg4 = 2.222219843214978396e-01 |
| |
| static const static double double | Lg5 = 1.818357216161805012e-01 |
| |
| static const static double double | Lg6 = 1.531383769920937332e-01 |
| |
| static const static double double | Lg7 = 1.479819860511658591e-01 |
| |
| static double | zero = 0.0 |
| |
| int | k |
| |
| int | hx = __HI(x) |
| |
| int | i = (hx + 0x95f64) & 0x100000 |
| |
| int | j = 0x6b851 - hx |
| |
| unsigned | lx = __LO(x) |
| |
| | f = x - 1.0 |
| |
| | s = f / (2.0 + f) |
| |
| | dk = (double)k |
| |
| | z = s * s |
| |
| | w = z * z |
| |
| | t1 = w * (Lg2 + w * (Lg4 + w * Lg6)) |
| |
| | t2 = z * (Lg1 + w * (Lg3 + w * (Lg5 + w * Lg7))) |
| |
| | R = t2 + t1 |
| |