testcases for rint double rounding issue
authorSzabolcs Nagy <nsz@port70.net>
Tue, 28 Oct 2014 01:24:16 +0000 (02:24 +0100)
committerSzabolcs Nagy <nsz@port70.net>
Tue, 28 Oct 2014 01:24:16 +0000 (02:24 +0100)
src/math/special/rint.h
src/math/special/rintf.h

index f36f319..d6a917e 100644 (file)
@@ -43,3 +43,5 @@ T(RU,             0x1.ffffp-1,                  0x1p+0,          0x0p+0, INEXACT
 T(RU,            -0x1.ffffp-1,                 -0x0p+0,          0x0p+0, INEXACT)
 T(RU,                0x1p-100,                  0x1p+0,          0x0p+0, INEXACT)
 T(RU,               -0x1p-100,                 -0x0p+0,          0x0p+0, INEXACT)
 T(RU,            -0x1.ffffp-1,                 -0x0p+0,          0x0p+0, INEXACT)
 T(RU,                0x1p-100,                  0x1p+0,          0x0p+0, INEXACT)
 T(RU,               -0x1p-100,                 -0x0p+0,          0x0p+0, INEXACT)
+T(RN,    0x1.0800000000001p+4,                0x1.1p+4,          0x0p+0, INEXACT)
+T(RN,    0x1.17fffffffffffp+4,                0x1.1p+4,          0x0p+0, INEXACT)
index e7b6047..09ac456 100644 (file)
@@ -43,3 +43,5 @@ T(RU,     0x1.ffffp-1,          0x1p+0,          0x0p+0, INEXACT)
 T(RU,    -0x1.ffffp-1,         -0x0p+0,          0x0p+0, INEXACT)
 T(RU,        0x1p-100,          0x1p+0,          0x0p+0, INEXACT)
 T(RU,       -0x1p-100,         -0x0p+0,          0x0p+0, INEXACT)
 T(RU,    -0x1.ffffp-1,         -0x0p+0,          0x0p+0, INEXACT)
 T(RU,        0x1p-100,          0x1p+0,          0x0p+0, INEXACT)
 T(RU,       -0x1p-100,         -0x0p+0,          0x0p+0, INEXACT)
+T(RN,   0x1.400002p+1,        0x1.8p+1,          0x0p+0, INEXACT)
+T(RN,   0x1.7ffffep+0,          0x1p+0,          0x0p+0, INEXACT)