[444] | 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) |
---|