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_csqrtf.c | |
| parent | bc62c2fb152d6ffec63975d88fd8f1bc3d3b7c01 (diff) | |
| download | glibc-77425c63e72bc0c01d81bc7e9ba4bb41d11679e6.tar.gz glibc-77425c63e72bc0c01d81bc7e9ba4bb41d11679e6.tar.xz glibc-77425c63e72bc0c01d81bc7e9ba4bb41d11679e6.zip | |
Add branch predictions to complex math code
Diffstat (limited to 'math/s_csqrtf.c')
| -rw-r--r-- | math/s_csqrtf.c | 9 |
1 files changed, 4 insertions, 5 deletions
diff --git a/math/s_csqrtf.c b/math/s_csqrtf.c index 1cf3b79c31..1613192a5e 100644 --- a/math/s_csqrtf.c +++ b/math/s_csqrtf.c @@ -1,5 +1,5 @@ /* Complex square root of float value. - Copyright (C) 1997, 1998, 2005 Free Software Foundation, Inc. + Copyright (C) 1997, 1998, 2005, 2011 Free Software Foundation, Inc. This file is part of the GNU C Library. Based on an algorithm by Stephen L. Moshier <moshier@world.std.com>. Contributed by Ulrich Drepper <drepper@cygnus.com>, 1997. @@ -21,7 +21,6 @@ #include <complex.h> #include <math.h> - #include <math_private.h> @@ -32,7 +31,7 @@ __csqrtf (__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) { @@ -61,7 +60,7 @@ __csqrtf (__complex__ float x) } else { - if (icls == FP_ZERO) + if (__builtin_expect (icls == FP_ZERO, 0)) { if (__real__ x < 0.0) { @@ -75,7 +74,7 @@ __csqrtf (__complex__ float x) __imag__ res = __copysignf (0.0, __imag__ x); } } - else if (rcls == FP_ZERO) + else if (__builtin_expect (rcls == FP_ZERO, 0)) { float r = __ieee754_sqrtf (0.5 * fabsf (__imag__ x)); |