diff options
author | Ulrich Drepper <drepper@gmail.com> | 2011-10-22 13:17:30 -0400 |
---|---|---|
committer | Ulrich Drepper <drepper@gmail.com> | 2011-10-22 13:17:30 -0400 |
commit | 77425c63e72bc0c01d81bc7e9ba4bb41d11679e6 (patch) | |
tree | c4fea5ed8070fc251140eba66b009bca1ff0e3b5 /math/s_catanhf.c | |
parent | bc62c2fb152d6ffec63975d88fd8f1bc3d3b7c01 (diff) | |
download | glibc-77425c63e72bc0c01d81bc7e9ba4bb41d11679e6.tar.gz |
Add branch predictions to complex math code
Diffstat (limited to 'math/s_catanhf.c')
-rw-r--r-- | math/s_catanhf.c | 15 |
1 files changed, 6 insertions, 9 deletions
diff --git a/math/s_catanhf.c b/math/s_catanhf.c index 7424bda9b9..f3d07f2354 100644 --- a/math/s_catanhf.c +++ b/math/s_catanhf.c @@ -1,5 +1,5 @@ /* Return arc hyperbole tangent for float value. - Copyright (C) 1997 Free Software Foundation, Inc. + Copyright (C) 1997, 2011 Free Software Foundation, Inc. This file is part of the GNU C Library. Contributed by Ulrich Drepper <drepper@cygnus.com>, 1997. @@ -20,7 +20,6 @@ #include <complex.h> #include <math.h> - #include <math_private.h> @@ -31,7 +30,7 @@ __catanhf (__complex__ float x) int rcls = fpclassify (__real__ x); int icls = fpclassify (__imag__ x); - if (rcls <= FP_INFINITE || icls <= FP_INFINITE) + if (__builtin_expect (rcls <= FP_INFINITE || icls <= FP_INFINITE, 0)) { if (icls == FP_INFINITE) { @@ -52,20 +51,18 @@ __catanhf (__complex__ float x) __imag__ res = __nanf (""); } } - else if (rcls == FP_ZERO && icls == FP_ZERO) + else if (__builtin_expect (rcls == FP_ZERO && icls == FP_ZERO, 0)) { res = x; } else { - float i2, num, den; - - i2 = __imag__ x * __imag__ x; + float i2 = __imag__ x * __imag__ x; - num = 1.0 + __real__ x; + float num = 1.0 + __real__ x; num = i2 + num * num; - den = 1.0 - __real__ x; + float den = 1.0 - __real__ x; den = i2 + den * den; __real__ res = 0.25 * (__ieee754_logf (num) - __ieee754_logf (den)); |