diff options
Diffstat (limited to 'lib/libm/j1.c')
| -rw-r--r-- | lib/libm/j1.c | 25 |
1 files changed, 11 insertions, 14 deletions
diff --git a/lib/libm/j1.c b/lib/libm/j1.c index 839160c4..211c78bf 100644 --- a/lib/libm/j1.c +++ b/lib/libm/j1.c @@ -63,7 +63,7 @@ static double pone(double), qone(double); static const double invsqrtpi = 5.64189583547756279280e-01, /* 0x3FE20DD7, 0x50429B6D */ - tpi = 6.36619772367581382433e-01; /* 0x3FE45F30, 0x6DC9C883 */ + tpi = 6.36619772367581382433e-01; /* 0x3FE45F30, 0x6DC9C883 */ static double common(uint32_t ix, double x, int y1, int sign) { @@ -104,14 +104,14 @@ static double common(uint32_t ix, double x, int y1, int sign) /* R0/S0 on [0,2] */ static const double r00 = -6.25000000000000000000e-02, /* 0xBFB00000, 0x00000000 */ - r01 = 1.40705666955189706048e-03, /* 0x3F570D9F, 0x98472C61 */ - r02 = -1.59955631084035597520e-05, /* 0xBEF0C5C6, 0xBA169668 */ - r03 = 4.96727999609584448412e-08, /* 0x3E6AAAFA, 0x46CA0BD9 */ - s01 = 1.91537599538363460805e-02, /* 0x3F939D0B, 0x12637E53 */ - s02 = 1.85946785588630915560e-04, /* 0x3F285F56, 0xB9CDF664 */ - s03 = 1.17718464042623683263e-06, /* 0x3EB3BFF8, 0x333F8498 */ - s04 = 5.04636257076217042715e-09, /* 0x3E35AC88, 0xC97DFF2C */ - s05 = 1.23542274426137913908e-11; /* 0x3DAB2ACF, 0xCFB97ED8 */ + r01 = 1.40705666955189706048e-03, /* 0x3F570D9F, 0x98472C61 */ + r02 = -1.59955631084035597520e-05, /* 0xBEF0C5C6, 0xBA169668 */ + r03 = 4.96727999609584448412e-08, /* 0x3E6AAAFA, 0x46CA0BD9 */ + s01 = 1.91537599538363460805e-02, /* 0x3F939D0B, 0x12637E53 */ + s02 = 1.85946785588630915560e-04, /* 0x3F285F56, 0xB9CDF664 */ + s03 = 1.17718464042623683263e-06, /* 0x3EB3BFF8, 0x333F8498 */ + s04 = 5.04636257076217042715e-09, /* 0x3E35AC88, 0xC97DFF2C */ + s05 = 1.23542274426137913908e-11; /* 0x3DAB2ACF, 0xCFB97ED8 */ double j1(double x) { @@ -172,8 +172,7 @@ double y1(double x) return -tpi / x; z = x * x; u = U0[0] + z * (U0[1] + z * (U0[2] + z * (U0[3] + z * U0[4]))); - v = 1 + - z * (V0[0] + z * (V0[1] + z * (V0[2] + z * (V0[3] + z * V0[4])))); + v = 1 + z * (V0[0] + z * (V0[1] + z * (V0[2] + z * (V0[3] + z * V0[4])))); return x * (u / v) + tpi * (j1(x) * log(x) - 1 / x); } @@ -385,8 +384,6 @@ static double qone(double x) } z = 1.0 / (x * x); r = p[0] + z * (p[1] + z * (p[2] + z * (p[3] + z * (p[4] + z * p[5])))); - s = 1.0 + - z * (q[0] + - z * (q[1] + z * (q[2] + z * (q[3] + z * (q[4] + z * q[5]))))); + s = 1.0 + z * (q[0] + z * (q[1] + z * (q[2] + z * (q[3] + z * (q[4] + z * q[5]))))); return (.375 + r / s) / x; } |