diff options
Diffstat (limited to 'sysdeps/generic/math_private_calls.h')
| -rw-r--r-- | sysdeps/generic/math_private_calls.h | 124 |
1 files changed, 124 insertions, 0 deletions
diff --git a/sysdeps/generic/math_private_calls.h b/sysdeps/generic/math_private_calls.h new file mode 100644 index 0000000000..16b4ea3736 --- /dev/null +++ b/sysdeps/generic/math_private_calls.h @@ -0,0 +1,124 @@ +/* Private function declarations for libm. + Copyright (C) 2011-2018 Free Software Foundation, Inc. + This file is part of the GNU C Library. + + 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/>. */ + +#define __MSUF_X(x, suffix) x ## suffix +#define __MSUF_S(...) __MSUF_X (__VA_ARGS__) +#define __MSUF(x) __MSUF_S (x, _MSUF_) + +#define __MSUF_R_X(x, suffix) x ## suffix ## _r +#define __MSUF_R_S(...) __MSUF_R_X (__VA_ARGS__) +#define __MSUF_R(x) __MSUF_R_S (x, _MSUF_) + +/* IEEE style elementary functions. */ +extern _Mdouble_ __MSUF (__ieee754_acos) (_Mdouble_); +extern _Mdouble_ __MSUF (__ieee754_acosh) (_Mdouble_); +extern _Mdouble_ __MSUF (__ieee754_asin) (_Mdouble_); +extern _Mdouble_ __MSUF (__ieee754_atan2) (_Mdouble_, _Mdouble_); +extern _Mdouble_ __MSUF (__ieee754_atanh) (_Mdouble_); +extern _Mdouble_ __MSUF (__ieee754_cosh) (_Mdouble_); +extern _Mdouble_ __MSUF (__ieee754_exp) (_Mdouble_); +extern _Mdouble_ __MSUF (__ieee754_exp10) (_Mdouble_); +extern _Mdouble_ __MSUF (__ieee754_exp2) (_Mdouble_); +extern _Mdouble_ __MSUF (__ieee754_fmod) (_Mdouble_, _Mdouble_); +extern _Mdouble_ __MSUF (__ieee754_gamma) (_Mdouble_); +extern _Mdouble_ __MSUF_R (__ieee754_gamma) (_Mdouble_, int *); +extern _Mdouble_ __MSUF (__ieee754_hypot) (_Mdouble_, _Mdouble_); +extern _Mdouble_ __MSUF (__ieee754_j0) (_Mdouble_); +extern _Mdouble_ __MSUF (__ieee754_j1) (_Mdouble_); +extern _Mdouble_ __MSUF (__ieee754_jn) (int, _Mdouble_); +extern _Mdouble_ __MSUF (__ieee754_lgamma) (_Mdouble_); +extern _Mdouble_ __MSUF_R (__ieee754_lgamma) (_Mdouble_, int *); +extern _Mdouble_ __MSUF (__ieee754_log) (_Mdouble_); +extern _Mdouble_ __MSUF (__ieee754_log10) (_Mdouble_); +extern _Mdouble_ __MSUF (__ieee754_log2) (_Mdouble_); +extern _Mdouble_ __MSUF (__ieee754_pow) (_Mdouble_, _Mdouble_); +extern _Mdouble_ __MSUF (__ieee754_remainder) (_Mdouble_, _Mdouble_); +extern _Mdouble_ __MSUF (__ieee754_sinh) (_Mdouble_); +extern _Mdouble_ __MSUF (__ieee754_sqrt) (_Mdouble_); +extern _Mdouble_ __MSUF (__ieee754_y0) (_Mdouble_); +extern _Mdouble_ __MSUF (__ieee754_y1) (_Mdouble_); +extern _Mdouble_ __MSUF (__ieee754_yn) (int, _Mdouble_); + +extern _Mdouble_ __MSUF (__ieee754_scalb) (_Mdouble_, _Mdouble_); +extern int __MSUF (__ieee754_ilogb) (_Mdouble_); + +extern int32_t __MSUF (__ieee754_rem_pio2) (_Mdouble_, _Mdouble_ *); + +/* fdlibm kernel functions. */ +extern _Mdouble_ __MSUF (__kernel_sin) (_Mdouble_, _Mdouble_, int); +extern _Mdouble_ __MSUF (__kernel_cos) (_Mdouble_, _Mdouble_); +extern _Mdouble_ __MSUF (__kernel_tan) (_Mdouble_, _Mdouble_, int); + +#if defined __MATH_DECLARING_LONG_DOUBLE || defined __MATH_DECLARING_FLOATN +extern void __MSUF (__kernel_sincos) (_Mdouble_, _Mdouble_, + _Mdouble_ *, _Mdouble_ *, int); +#endif + +#if !defined __MATH_DECLARING_LONG_DOUBLE || defined __MATH_DECLARING_FLOATN +extern int __MSUF (__kernel_rem_pio2) (_Mdouble_ *, _Mdouble_ *, int, + int, int, const int32_t *); +#endif + +/* Internal functions. */ +#if !defined __MATH_DECLARING_LONG_DOUBLE || !defined NO_LONG_DOUBLE +extern _Mdouble_ __MSUF (__copysign) (_Mdouble_ x, _Mdouble_ __y); + +extern inline _Mdouble_ +__MSUF (__copysign) (_Mdouble_ x, _Mdouble_ __y) +{ + return __MSUF (__builtin_copysign) (x, __y); +} +#endif + +/* Return X^2 + Y^2 - 1, computed without large cancellation error. + It is given that 1 > X >= Y >= epsilon / 2, and that X^2 + Y^2 >= + 0.5. */ +extern _Mdouble_ __MSUF (__x2y2m1) (_Mdouble_ x, _Mdouble_ y); + +/* Compute the product of X + X_EPS, X + X_EPS + 1, ..., X + X_EPS + N + - 1, in the form R * (1 + *EPS) where the return value R is an + approximation to the product and *EPS is set to indicate the + approximate error in the return value. X is such that all the + values X + 1, ..., X + N - 1 are exactly representable, and X_EPS / + X is small enough that factors quadratic in it can be + neglected. */ +extern _Mdouble_ __MSUF (__gamma_product) (_Mdouble_ x, _Mdouble_ x_eps, + int n, _Mdouble_ *eps); + +/* Compute lgamma of a negative argument X, if it is in a range + (depending on the floating-point format) for which expansion around + zeros is used, setting *SIGNGAMP accordingly. */ +extern _Mdouble_ __MSUF (__lgamma_neg) (_Mdouble_ x, int *signgamp); + +/* Compute the product of 1 + (T / (X + X_EPS)), 1 + (T / (X + X_EPS + + 1)), ..., 1 + (T / (X + X_EPS + N - 1)), minus 1. X is such that + all the values X + 1, ..., X + N - 1 are exactly representable, and + X_EPS / X is small enough that factors quadratic in it can be + neglected. */ +#if !defined __MATH_DECLARING_FLOAT +extern _Mdouble_ __MSUF (__lgamma_product) (_Mdouble_ t, _Mdouble_ x, + _Mdouble_ x_eps, int n); +#endif + +#undef __MSUF_X +#undef __MSUF_S +#undef __MSUF + +#undef __MSUF_R_X +#undef __MSUF_R_S +#undef __MSUF_R |