diff options
author | Paul E. Murphy <murphyp@linux.vnet.ibm.com> | 2016-06-28 08:49:23 -0500 |
---|---|---|
committer | Paul E. Murphy <murphyp@linux.vnet.ibm.com> | 2016-08-19 16:46:41 -0500 |
commit | c50eee19c447d3f2c182dc3a22f2b01a053dca41 (patch) | |
tree | 3b5f0d5c832bad20fce31502026f27fd6915ea8f /sysdeps/m68k/m680x0/fpu/s_csinh_template.c | |
parent | ffb84f5e197aaa9d46a35df84689c75d689d73cb (diff) | |
download | glibc-c50eee19c447d3f2c182dc3a22f2b01a053dca41.tar glibc-c50eee19c447d3f2c182dc3a22f2b01a053dca41.tar.gz glibc-c50eee19c447d3f2c182dc3a22f2b01a053dca41.tar.bz2 glibc-c50eee19c447d3f2c182dc3a22f2b01a053dca41.zip |
Convert _Complex sine functions to generated code
Refactor s_c{,a}sin{,h}{f,,l} into a single templated
macro.
Diffstat (limited to 'sysdeps/m68k/m680x0/fpu/s_csinh_template.c')
-rw-r--r-- | sysdeps/m68k/m680x0/fpu/s_csinh_template.c | 62 |
1 files changed, 62 insertions, 0 deletions
diff --git a/sysdeps/m68k/m680x0/fpu/s_csinh_template.c b/sysdeps/m68k/m680x0/fpu/s_csinh_template.c new file mode 100644 index 0000000000..56f3babd95 --- /dev/null +++ b/sysdeps/m68k/m680x0/fpu/s_csinh_template.c @@ -0,0 +1,62 @@ +/* Complex sine 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 CONCATX(a,b) __CONCAT(a,b) +#define s(name) M_SUF (name) +#define m81(func) __m81_u(s(func)) + +CFLOAT +s(__csinh) (CFLOAT x) +{ + CFLOAT retval; + unsigned long ix_cond; + + 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_sinh) (__real__ x); + if (ix_cond & __M81_COND_ZERO) + __imag__ retval = __imag__ x; + else + __imag__ retval = sin_ix * m81(__ieee754_cosh) (__real__ x); + } + else + { + unsigned long rx_cond = __m81_test (__real__ x); + + __imag__ retval = __imag__ x - __imag__ x; + if (rx_cond & (__M81_COND_ZERO|__M81_COND_INF|__M81_COND_NAN)) + __real__ retval = __real__ x; + else + __real__ retval = __imag__ retval; + } + + return retval; +} +weak_alias (s(__csinh), s(csinh)) |