fmaf test case that's hard to round in nearest rounding mode

[libc-test] / src / math / special / y0.h

diff --git a/src/math/special/y0.h b/src/math/special/y0.h
index f9267e2..a47c03d 100644 (file)
--- a/src/math/special/y0.h
+++ b/src/math/special/y0.h
@@ -4,3 +4,10 @@ T(RN,                 -0x1p+0,                     nan,          0x0p+0, INVALID
 T(RN,                     inf,                  0x0p+0,          0x0p+0, 0)
 T(RN,                    -inf,                     nan,          0x0p+0, INVALID)
 T(RN,                     nan,                     nan,          0x0p+0, 0)
+// hard cases
+T(RN,    0x1.c982eb8d417eap-1,  -0x1.af74bfa0f1304p-56,   0x1.11721cp-2, INEXACT)
+T(RN,    0x1.c982eb8d417ebp-1,   0x1.5666419c0f3c9p-54,   0x1.dc02a8p-2, INEXACT)
+T(RN,    0x1.fa9534d98569bp+1,   0x1.384a000f3fcecp-53,    0x1.58e17p-6, INEXACT)
+T(RN,    0x1.fa9534d98569cp+1,  -0x1.8fa8956b4b481p-55,   0x1.54cfbep-2, INEXACT)
+T(RN,    0x1.c581dc4e72102p+2,  -0x1.14bb186dc408dp-52,   -0x1.163cfp-2, INEXACT)
+T(RN,    0x1.c581dc4e72103p+2,   0x1.e91b198d39ce2p-56,  -0x1.ee1b6ep-2, INEXACT)