AuroraMiddleware/glibc
1
0
Fork 0
mirror of https://sourceware.org/git/glibc.git synced 2024-11-10 07:10:06 +00:00
Code Issues Packages Projects Releases Wiki Activity

* sysdeps/generic/s_ctanh.c (__ctanh): Handle case of zero den better.

Browse Source 
* sysdeps/generic/s_ctanhf.c (__ctanhf): Likewise.
	* sysdeps/generic/s_ctanhl.c (__ctanhl): Likewise.
This commit is contained in:
Ulrich Drepper 2005-07-07 23:31:46 +00:00
parent 531680aebf
commit 96191db217
3 changed files with 32 additions and 6 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

6
ChangeLog
View File

@ -1,3 +1,9 @@
2005-07-07 Ulrich Drepper <drepper@redhat.com>
* sysdeps/generic/s_ctanh.c (__ctanh): Handle case of zero den better.
* sysdeps/generic/s_ctanhf.c (__ctanhf): Likewise.
* sysdeps/generic/s_ctanhl.c (__ctanhl): Likewise.
2005-04-13 H.J. Lu <hongjiu.lu@intel.com>
[BZ #974]

16
sysdeps/generic/s_ctanh.c
View File

@ -1,5 +1,5 @@
/* Complex hyperbole tangent for double.
Copyright (C) 1997 Free Software Foundation, Inc.
Copyright (C) 1997, 2005 Free Software Foundation, Inc.
This file is part of the GNU C Library.
Contributed by Ulrich Drepper <drepper@cygnus.com>, 1997.
@ -61,8 +61,18 @@ __ctanh (__complex__ double x)
den = (__ieee754_cosh (2.0 * __real__ x) + cos2ix);
__real__ res = __ieee754_sinh (2.0 * __real__ x) / den;
__imag__ res = sin2ix / den;
if (den == 0.0)
{
__complex__ double ez = __cexp (x);
__complex__ double emz = __cexp (-x);
res = (ez - emz) / (ez + emz);
}
else
{
__real__ res = __ieee754_sinh (2.0 * __real__ x) / den;
__imag__ res = sin2ix / den;
}
}
return res;

16
sysdeps/generic/s_ctanhf.c
View File

@ -1,5 +1,5 @@
/* Complex hyperbole tangent for float.
Copyright (C) 1997 Free Software Foundation, Inc.
Copyright (C) 1997, 2005 Free Software Foundation, Inc.
This file is part of the GNU C Library.
Contributed by Ulrich Drepper <drepper@cygnus.com>, 1997.
@ -61,8 +61,18 @@ __ctanhf (__complex__ float x)
den = (__ieee754_coshf (2.0 * __real__ x) + cos2ix);
__real__ res = __ieee754_sinhf (2.0 * __real__ x) / den;
__imag__ res = sin2ix / den;
if (den == 0.0f)
{
__complex__ float ez = __cexpf (x);
__complex__ float emz = __cexpf (-x);
res = (ez - emz) / (ez + emz);
}
else
{
__real__ res = __ieee754_sinhf (2.0 * __real__ x) / den;
__imag__ res = sin2ix / den;
}
}
return res;