Changeset 1006Tweet

Timestamp:

Oct 4, 2011, 3:35:55 PM (12 years ago)

Author:

sam

Message:

core: replace the cos(real) implementation with a simple call to
sin(pi/2-x). No loss in precision in theory.

File:

• trunk/src/real.cpp (modified) (1 diff)

trunk/src/real.cpp

@@ -739 +739 @@
 real cos(real const &x)
 {
-    bool switch_sign = false;
-    real absx = fmod(fabs(x), real::R_PI << 1);
-
-    if (absx > real::R_PI)
-        absx = (real::R_PI << 1) - absx;
-
-    if (absx > real::R_PI_2)
-    {
-        absx = real::R_PI - absx;
-        switch_sign = true;
-    }
-
-    real ret = 0.0, fact = 1.0, xn = 1.0, x2 = absx * absx;
-    for (int i = 1; ; i += 2)
-    {
-        real newret = ret + xn / fact;
-        if (ret == newret)
-            break;
-        ret = newret;
-        xn *= x2;
-        fact *= (real)(-i * (i + 1));
-    }
-
-    /* Propagate sign */
-    if (switch_sign)
-        ret.m_signexp ^= 0x80000000u;
-    return ret;
+    return sin(real::R_PI_2 - x);
 }