diff options
Diffstat (limited to 'libquadmath/math/j0q.c')
| -rw-r--r-- | libquadmath/math/j0q.c | 3 |
1 files changed, 3 insertions, 0 deletions
diff --git a/libquadmath/math/j0q.c b/libquadmath/math/j0q.c index 3ef9201a8a4..8c6f811125e 100644 --- a/libquadmath/math/j0q.c +++ b/libquadmath/math/j0q.c @@ -807,6 +807,7 @@ static __float128 Y0_2D[NY0_2D + 1] = { /* 1.000000000000000000000000000000000000000E0 */ }; +static const long double U0 = -7.3804295108687225274343927948483016310862e-02Q; /* Bessel function of the second kind, order zero. */ @@ -829,6 +830,8 @@ y0q (__float128 x) return -HUGE_VALQ + x; } xx = fabsq (x); + if (xx <= 0x1p-57) + return U0 + TWOOPI * logq (x); if (xx <= 2.0Q) { /* 0 <= x <= 2 */ |