1 | /* Complex sine hyperbole function for long double. |
---|
2 | Copyright (C) 1997 Free Software Foundation, Inc. |
---|
3 | This file is part of the GNU C Library. |
---|
4 | Contributed by Ulrich Drepper <drepper@cygnus.com>, 1997. |
---|
5 | |
---|
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. |
---|
10 | |
---|
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. |
---|
15 | |
---|
16 | You should have received a copy of the GNU Lesser General Public |
---|
17 | License along with the GNU C Library; if not, write to the Free |
---|
18 | Software Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA |
---|
19 | 02111-1307 USA. */ |
---|
20 | |
---|
21 | #include <complex.h> |
---|
22 | #include <fenv.h> |
---|
23 | #include <math.h> |
---|
24 | |
---|
25 | #include "math_private.h" |
---|
26 | |
---|
27 | |
---|
28 | __complex__ long double |
---|
29 | __csinhl (__complex__ long double x) |
---|
30 | { |
---|
31 | __complex__ long double retval; |
---|
32 | int negate = signbit (__real__ x); |
---|
33 | int rcls = fpclassify (__real__ x); |
---|
34 | int icls = fpclassify (__imag__ x); |
---|
35 | |
---|
36 | __real__ x = fabsl (__real__ x); |
---|
37 | |
---|
38 | if (rcls >= FP_ZERO) |
---|
39 | { |
---|
40 | /* Real part is finite. */ |
---|
41 | if (icls >= FP_ZERO) |
---|
42 | { |
---|
43 | /* Imaginary part is finite. */ |
---|
44 | long double sinh_val = __ieee754_sinhl (__real__ x); |
---|
45 | long double cosh_val = __ieee754_coshl (__real__ x); |
---|
46 | long double sinix, cosix; |
---|
47 | |
---|
48 | __sincosl (__imag__ x, &sinix, &cosix); |
---|
49 | |
---|
50 | __real__ retval = sinh_val * cosix; |
---|
51 | __imag__ retval = cosh_val * sinix; |
---|
52 | |
---|
53 | if (negate) |
---|
54 | __real__ retval = -__real__ retval; |
---|
55 | } |
---|
56 | else |
---|
57 | { |
---|
58 | if (rcls == FP_ZERO) |
---|
59 | { |
---|
60 | /* Real part is 0.0. */ |
---|
61 | __real__ retval = __copysignl (0.0, negate ? -1.0 : 1.0); |
---|
62 | __imag__ retval = __nanl ("") + __nanl (""); |
---|
63 | |
---|
64 | #ifdef FE_INVALID |
---|
65 | if (icls == FP_INFINITE) |
---|
66 | feraiseexcept (FE_INVALID); |
---|
67 | #endif |
---|
68 | } |
---|
69 | else |
---|
70 | { |
---|
71 | __real__ retval = __nanl (""); |
---|
72 | __imag__ retval = __nanl (""); |
---|
73 | |
---|
74 | #ifdef FE_INVALID |
---|
75 | feraiseexcept (FE_INVALID); |
---|
76 | #endif |
---|
77 | } |
---|
78 | } |
---|
79 | } |
---|
80 | else if (rcls == FP_INFINITE) |
---|
81 | { |
---|
82 | /* Real part is infinite. */ |
---|
83 | if (icls == FP_ZERO) |
---|
84 | { |
---|
85 | /* Imaginary part is 0.0. */ |
---|
86 | __real__ retval = negate ? -HUGE_VALL : HUGE_VALL; |
---|
87 | __imag__ retval = __imag__ x; |
---|
88 | } |
---|
89 | else if (icls > FP_ZERO) |
---|
90 | { |
---|
91 | /* Imaginary part is finite. */ |
---|
92 | long double sinix, cosix; |
---|
93 | |
---|
94 | __sincosl (__imag__ x, &sinix, &cosix); |
---|
95 | |
---|
96 | __real__ retval = __copysignl (HUGE_VALL, cosix); |
---|
97 | __imag__ retval = __copysignl (HUGE_VALL, sinix); |
---|
98 | |
---|
99 | if (negate) |
---|
100 | __real__ retval = -__real__ retval; |
---|
101 | } |
---|
102 | else |
---|
103 | { |
---|
104 | /* The addition raises the invalid exception. */ |
---|
105 | __real__ retval = HUGE_VALL; |
---|
106 | __imag__ retval = __nanl ("") + __nanl (""); |
---|
107 | |
---|
108 | #ifdef FE_INVALID |
---|
109 | if (icls == FP_INFINITE) |
---|
110 | feraiseexcept (FE_INVALID); |
---|
111 | #endif |
---|
112 | } |
---|
113 | } |
---|
114 | else |
---|
115 | { |
---|
116 | __real__ retval = __nanl (""); |
---|
117 | __imag__ retval = __imag__ x == 0.0 ? __imag__ x : __nanl (""); |
---|
118 | } |
---|
119 | |
---|
120 | return retval; |
---|
121 | } |
---|
122 | weak_alias (__csinhl, csinhl) |
---|