lib/glibc
1
0
Fork 0
mirror of https://sourceware.org/git/glibc.git synced 2025-07-29 11:41:21 +03:00
Code Activity

Fix sign errr in some cacosh results

Browse Source
This commit is contained in:
Richard B. Kreckel
2011-12-21 21:01:29 -05:00
committed by Ulrich Drepper
parent 60e8585f2a
commit b27e24b874
5 changed files with 13 additions and 32 deletions
Download Patch File Download Diff File Expand all files Collapse all files

22
math/s_cacoshf.c
View File

@ -1,5 +1,5 @@
/* Return arc hyperbole cosine for float value.
Copyright (C) 1997, 2006 Free Software Foundation, Inc.
Copyright (C) 1997, 2006, 2011 Free Software Foundation, Inc.
This file is part of the GNU C Library.
Contributed by Ulrich Drepper <drepper@cygnus.com>, 1997.
@ -67,7 +67,6 @@ __cacoshf (__complex__ float x)
}
else
{
#if 1
__complex__ float y;
__real__ y = (__real__ x - __imag__ x) * (__real__ x + __imag__ x) - 1.0;
@ -75,30 +74,13 @@ __cacoshf (__complex__ float x)
y = __csqrtf (y);
if (__real__ x < 0.0)
if (signbit (__real__ x))
y = -y;
__real__ y += __real__ x;
__imag__ y += __imag__ x;
res = __clogf (y);
#else
float re2 = __real__ x * __real__ x;
float im2 = __imag__ x * __imag__ x;
float sq = re2 - im2 - 1.0;
float ro = __ieee754_sqrtf (sq * sq + 4 * re2 * im2);
float a = __ieee754_sqrtf ((sq + ro) / 2.0);
float b = __ieee754_sqrtf ((-sq + ro) / 2.0);
__real__ res = 0.5 * __ieee754_logf (re2 + __real__ x * 2 * a
+ im2 + __imag__ x * 2 * b
+ ro);
__imag__ res = __ieee754_atan2f (__imag__ x + b, __real__ x + a);
#endif
/* We have to use the positive branch. */
if (__real__ res < 0.0)
res = -res;
}
return res;