Replace FSF snail mail address with URLs.
[glibc.git] / math / s_csinhf.c
blob926bb37990906e0d119e58184f92b6938c0360e7
1 /* Complex sine hyperbole function for float.
2 Copyright (C) 1997, 2011 Free Software Foundation, Inc.
3 This file is part of the GNU C Library.
4 Contributed by Ulrich Drepper <drepper@cygnus.com>, 1997.
6 The GNU C Library is free software; you can redistribute it and/or
7 modify it under the terms of the GNU Lesser General Public
8 License as published by the Free Software Foundation; either
9 version 2.1 of the License, or (at your option) any later version.
11 The GNU C Library is distributed in the hope that it will be useful,
12 but WITHOUT ANY WARRANTY; without even the implied warranty of
13 MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
14 Lesser General Public License for more details.
16 You should have received a copy of the GNU Lesser General Public
17 License along with the GNU C Library; if not, see
18 <http://www.gnu.org/licenses/>. */
20 #include <complex.h>
21 #include <fenv.h>
22 #include <math.h>
24 #include <math_private.h>
27 __complex__ float
28 __csinhf (__complex__ float x)
30 __complex__ float retval;
31 int negate = signbit (__real__ x);
32 int rcls = fpclassify (__real__ x);
33 int icls = fpclassify (__imag__ x);
35 __real__ x = fabsf (__real__ x);
37 if (__builtin_expect (rcls >= FP_ZERO, 1))
39 /* Real part is finite. */
40 if (__builtin_expect (icls >= FP_ZERO, 1))
42 /* Imaginary part is finite. */
43 float sinh_val = __ieee754_sinhf (__real__ x);
44 float cosh_val = __ieee754_coshf (__real__ x);
45 float sinix, cosix;
47 __sincosf (__imag__ x, &sinix, &cosix);
49 __real__ retval = sinh_val * cosix;
50 __imag__ retval = cosh_val * sinix;
52 if (negate)
53 __real__ retval = -__real__ retval;
55 else
57 if (rcls == FP_ZERO)
59 /* Real part is 0.0. */
60 __real__ retval = __copysignf (0.0, negate ? -1.0 : 1.0);
61 __imag__ retval = __nanf ("") + __nanf ("");
63 if (icls == FP_INFINITE)
64 feraiseexcept (FE_INVALID);
66 else
68 __real__ retval = __nanf ("");
69 __imag__ retval = __nanf ("");
71 feraiseexcept (FE_INVALID);
75 else if (__builtin_expect (rcls == FP_INFINITE, 1))
77 /* Real part is infinite. */
78 if (__builtin_expect (icls > FP_ZERO, 1))
80 /* Imaginary part is finite. */
81 float sinix, cosix;
83 __sincosf (__imag__ x, &sinix, &cosix);
85 __real__ retval = __copysignf (HUGE_VALF, cosix);
86 __imag__ retval = __copysignf (HUGE_VALF, sinix);
88 if (negate)
89 __real__ retval = -__real__ retval;
91 else if (icls == FP_ZERO)
93 /* Imaginary part is 0.0. */
94 __real__ retval = negate ? -HUGE_VALF : HUGE_VALF;
95 __imag__ retval = __imag__ x;
97 else
99 /* The addition raises the invalid exception. */
100 __real__ retval = HUGE_VALF;
101 __imag__ retval = __nanf ("") + __nanf ("");
103 #ifdef FE_INVALID
104 if (icls == FP_INFINITE)
105 feraiseexcept (FE_INVALID);
106 #endif
109 else
111 __real__ retval = __nanf ("");
112 __imag__ retval = __imag__ x == 0.0 ? __imag__ x : __nanf ("");
115 return retval;
117 #ifndef __csinhf
118 weak_alias (__csinhf, csinhf)
119 #endif