Changeset 1130Tweet

Timestamp:

Jan 12, 2012, 12:36:29 AM (10 years ago)

Author:

sam

Message:

math: implement ulp() for reals, which returns the smallest real y > 0 such
that x + y != x, and nextafter() which behaves like the C function.

Location:

trunk

Files:

• src/lol/math/real.h (modified) (1 diff)
• src/real.cpp (modified) (1 diff)
• test/unit/real.cpp (modified) (1 diff)

trunk/src/lol/math/real.h

@@ -89 +89 @@
     friend real ldexp(real const &x, int exp);
     friend real modf(real const &x, real *iptr);
+    friend real ulp(real const &x);
+    friend real nextafter(real const &x, real const &y);
 
     /* Power functions */

trunk/src/real.cpp

@@ -939 +939 @@
 }
 
+real ulp(real const &x)
+{
+    real ret = real::R_1;
+    if (x)
+        ret.m_signexp = x.m_signexp + 1 - real::BIGITS * real::BIGIT_BITS;
+    else
+        ret.m_signexp = 0;
+    return ret;
+}
+
+real nextafter(real const &x, real const &y)
+{
+    if (x == y)
+        return x;
+    else if (x < y)
+        return x + ulp(x);
+    else
+        return x - ulp(x);
+}
+
 real copysign(real const &x, real const &y)
 {

trunk/test/unit/real.cpp

@@ -241 +241 @@
         LOLUNIT_ASSERT_EQUAL((double)ldexp(a3, 7), 0.0);
         LOLUNIT_ASSERT_EQUAL((double)ldexp(a3, -7), 0.0);
+    }
+
+    LOLUNIT_TEST(Ulp)
+    {
+        real a1 = real::R_PI;
+
+        LOLUNIT_ASSERT_NOT_EQUAL((double)(a1 + ulp(a1) - a1), 0.0);
+        LOLUNIT_ASSERT_EQUAL((double)(a1 + ulp(a1) / 2 - a1), 0.0);
     }