diff options
author | Thomas Fitzsimmons <fitzsim@redhat.com> | 2002-08-17 01:29:45 +0400 |
---|---|---|
committer | Thomas Fitzsimmons <fitzsim@redhat.com> | 2002-08-17 01:29:45 +0400 |
commit | d3bd3632ac06bdb7f76b13494873b3567edabcbf (patch) | |
tree | f7d4c5083ac5dbdc7b55c198fe3df5651e9c4e02 /newlib/libc/sys/linux/cmath/s_catanhf.c | |
parent | e0f61e607c336e5ad589afd95a13161a202dd1ba (diff) |
* libc/sys/linux/cmath: New directory.
* libc/sys/linux/include/cmathcalls.h: New file.
* libc/sys/linux/include/complex.h: New file.
* libc/sys/linux/machine/i386/huge_val.h: New file
* libm/math/w_sincos.c: New file
* libm/math/wf_sincos.c: New file
* libm/mathfp/s_sincos.c: New file
* libm/mathfp/sf_sincos.c: New file
* Makefile.am (LIBC_OBJECTLISTS): Add cmath/objectlist.awk.in.
* libc/include/math.h: Add sincos and sincosf declarations.
* libc/sys/linux/Makefile.am (SUBDIRS): Add cmath.
(SUBLIBS): Likewise.
* libc/sys/linux/configure.in (AC_OUTPUT): Add cmath.
* libm/math/Makefile.am (src): Add w_sincos.c.
(fsrc): Add wf_sincos.c.
* libm/mathfp/Makefile.am (src): Add s_sincos.c
(fsrc): Add sf_sincos.c.
Diffstat (limited to 'newlib/libc/sys/linux/cmath/s_catanhf.c')
-rw-r--r-- | newlib/libc/sys/linux/cmath/s_catanhf.c | 80 |
1 files changed, 80 insertions, 0 deletions
diff --git a/newlib/libc/sys/linux/cmath/s_catanhf.c b/newlib/libc/sys/linux/cmath/s_catanhf.c new file mode 100644 index 000000000..30a5a1b8a --- /dev/null +++ b/newlib/libc/sys/linux/cmath/s_catanhf.c @@ -0,0 +1,80 @@ +/* Return arc hyperbole tangent for float value. + Copyright (C) 1997 Free Software Foundation, Inc. + This file is part of the GNU C Library. + Contributed by Ulrich Drepper <drepper@cygnus.com>, 1997. + + The GNU C Library is free software; you can redistribute it and/or + modify it under the terms of the GNU Lesser General Public + License as published by the Free Software Foundation; either + version 2.1 of the License, or (at your option) any later version. + + The GNU C Library is distributed in the hope that it will be useful, + but WITHOUT ANY WARRANTY; without even the implied warranty of + MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU + Lesser General Public License for more details. + + You should have received a copy of the GNU Lesser General Public + License along with the GNU C Library; if not, write to the Free + Software Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA + 02111-1307 USA. */ + +#include <complex.h> +#include <math.h> + +#include "math_private.h" + + +__complex__ float +__catanhf (__complex__ float x) +{ + __complex__ float res; + int rcls = fpclassify (__real__ x); + int icls = fpclassify (__imag__ x); + + if (rcls <= FP_INFINITE || icls <= FP_INFINITE) + { + if (icls == FP_INFINITE) + { + __real__ res = __copysignf (0.0, __real__ x); + __imag__ res = __copysignf (M_PI_2, __imag__ x); + } + else if (rcls == FP_INFINITE || rcls == FP_ZERO) + { + __real__ res = __copysignf (0.0, __real__ x); + if (icls >= FP_ZERO) + __imag__ res = __copysignf (M_PI_2, __imag__ x); + else + __imag__ res = __nanf (""); + } + else + { + __real__ res = __nanf (""); + __imag__ res = __nanf (""); + } + } + else if (rcls == FP_ZERO && icls == FP_ZERO) + { + res = x; + } + else + { + float i2, num, den; + + i2 = __imag__ x * __imag__ x; + + num = 1.0 + __real__ x; + num = i2 + num * num; + + den = 1.0 - __real__ x; + den = i2 + den * den; + + __real__ res = 0.25 * (__ieee754_logf (num) - __ieee754_logf (den)); + + den = 1 - __real__ x * __real__ x - i2; + + __imag__ res = 0.5 * __ieee754_atan2f (2.0 * __imag__ x, den); + } + + return res; +} +weak_alias (__catanhf, catanhf) |