diff options
author | Paul E. Murphy <murphyp@linux.vnet.ibm.com> | 2016-06-27 17:11:46 -0500 |
---|---|---|
committer | Paul E. Murphy <murphyp@linux.vnet.ibm.com> | 2016-08-19 11:28:55 -0500 |
commit | 01ee387015a2075c45a4e1ad45d39e50b5a6d40b (patch) | |
tree | c8d2a88b0900c294d5e82fe869908d46b3ca838f /sysdeps/m68k/m680x0/fpu/s_ccosh_template.c | |
parent | 281f5073e5a6d2cd3733acd9c773c8c6340468c4 (diff) | |
download | glibc-01ee387015a2075c45a4e1ad45d39e50b5a6d40b.tar.gz glibc-01ee387015a2075c45a4e1ad45d39e50b5a6d40b.tar.xz glibc-01ee387015a2075c45a4e1ad45d39e50b5a6d40b.zip |
Convert _Complex cosine functions to generated code
This is fairly straight fowards. m68k overrides are updated to use the framework, and thus are simplified a bit.
Diffstat (limited to 'sysdeps/m68k/m680x0/fpu/s_ccosh_template.c')
-rw-r--r-- | sysdeps/m68k/m680x0/fpu/s_ccosh_template.c | 68 |
1 files changed, 68 insertions, 0 deletions
diff --git a/sysdeps/m68k/m680x0/fpu/s_ccosh_template.c b/sysdeps/m68k/m680x0/fpu/s_ccosh_template.c new file mode 100644 index 0000000000..83b449e7a8 --- /dev/null +++ b/sysdeps/m68k/m680x0/fpu/s_ccosh_template.c @@ -0,0 +1,68 @@ +/* Complex cosine hyperbole function. m68k fpu version + Copyright (C) 1997-2016 Free Software Foundation, Inc. + This file is part of the GNU C Library. + Contributed by Andreas Schwab <schwab@issan.informatik.uni-dortmund.de>. + + 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, see + <http://www.gnu.org/licenses/>. */ + +#include <complex.h> +#include <math.h> +#include "mathimpl.h" + +#define s(name) M_SUF (name) +#define m81(func) __m81_u(s(func)) + +CFLOAT +s(__ccosh) (CFLOAT x) +{ + CFLOAT retval; + unsigned long ix_cond = __m81_test (__imag__ x); + + if ((ix_cond & (__M81_COND_INF|__M81_COND_NAN)) == 0) + { + /* Imaginary part is finite. */ + FLOAT sin_ix, cos_ix; + + __asm ("fsincos%.x %2,%1:%0" : "=f" (sin_ix), "=f" (cos_ix) + : "f" (__imag__ x)); + __real__ retval = cos_ix * m81(__ieee754_cosh) (__real__ x); + if (ix_cond & __M81_COND_ZERO) + __imag__ retval = (signbit (__real__ x) + ? -__imag__ x : __imag__ x); + else + __imag__ retval = sin_ix * m81(__ieee754_sinh) (__real__ x); + } + else + { + unsigned long rx_cond = __m81_test (__real__ x); + + if (rx_cond & __M81_COND_ZERO) + { + __real__ retval = __imag__ x - __imag__ x; + __imag__ retval = __real__ x; + } + else + { + if (rx_cond & __M81_COND_INF) + __real__ retval = s(fabs) (__real__ x); + else + __real__ retval = s(__nan) (""); + __imag__ retval = __imag__ x - __imag__ x; + } + } + + return retval; +} +weak_alias (s(__ccosh), s(ccosh)) |