1 /* Complex cosine hyperbole function. m68k fpu version 2 Copyright (C) 1997-2022 Free Software Foundation, Inc. 3 This file is part of the GNU C Library. 4 5 The GNU C Library is free software; you can redistribute it and/or 6 modify it under the terms of the GNU Lesser General Public 7 License as published by the Free Software Foundation; either 8 version 2.1 of the License, or (at your option) any later version. 9 10 The GNU C Library is distributed in the hope that it will be useful, 11 but WITHOUT ANY WARRANTY; without even the implied warranty of 12 MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU 13 Lesser General Public License for more details. 14 15 You should have received a copy of the GNU Lesser General Public 16 License along with the GNU C Library. If not, see 17 <https://www.gnu.org/licenses/>. */ 18 19 #include <complex.h> 20 #include <math.h> 21 #include "mathimpl.h" 22 23 #define s(name) M_SUF (name) 24 #define m81(func) __m81_u(s(func)) 25 26 CFLOAT s(__ccosh)27s(__ccosh) (CFLOAT x) 28 { 29 CFLOAT retval; 30 unsigned long ix_cond = __m81_test (__imag__ x); 31 32 if ((ix_cond & (__M81_COND_INF|__M81_COND_NAN)) == 0) 33 { 34 /* Imaginary part is finite. */ 35 FLOAT sin_ix, cos_ix; 36 37 __asm ("fsincos%.x %2,%1:%0" : "=f" (sin_ix), "=f" (cos_ix) 38 : "f" (__imag__ x)); 39 __real__ retval = cos_ix * m81(__ieee754_cosh) (__real__ x); 40 if (ix_cond & __M81_COND_ZERO) 41 __imag__ retval = (signbit (__real__ x) 42 ? -__imag__ x : __imag__ x); 43 else 44 __imag__ retval = sin_ix * m81(__ieee754_sinh) (__real__ x); 45 } 46 else 47 { 48 unsigned long rx_cond = __m81_test (__real__ x); 49 50 if (rx_cond & __M81_COND_ZERO) 51 { 52 __real__ retval = __imag__ x - __imag__ x; 53 __imag__ retval = __real__ x; 54 } 55 else 56 { 57 if (rx_cond & __M81_COND_INF) 58 __real__ retval = s(fabs) (__real__ x); 59 else 60 __real__ retval = s(__nan) (""); 61 __imag__ retval = __imag__ x - __imag__ x; 62 } 63 } 64 65 return retval; 66 } 67 declare_mgen_alias (__ccosh, ccosh) 68