diff options
Diffstat (limited to 'math')
| -rw-r--r-- | math/w_exp_compat.c | 41 | ||||
| -rw-r--r-- | math/w_expf_compat.c | 37 | ||||
| -rw-r--r-- | math/w_expl_compat.c | 48 |
3 files changed, 118 insertions, 8 deletions
diff --git a/math/w_exp_compat.c b/math/w_exp_compat.c new file mode 100644 index 0000000000..98aa5b76a0 --- /dev/null +++ b/math/w_exp_compat.c @@ -0,0 +1,41 @@ +/* Copyright (C) 2011-2017 Free Software Foundation, Inc. + This file is part of the GNU C Library. + Contributed by Ulrich Drepper <drepper@gmail.com>, 2011. + + 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 <math.h> +#include <math_private.h> +#include <math-svid-compat.h> + +#if LIBM_SVID_COMPAT +/* wrapper exp */ +double +__exp (double x) +{ + double z = __ieee754_exp (x); + if (__builtin_expect (!isfinite (z) || z == 0, 0) + && isfinite (x) && _LIB_VERSION != _IEEE_) + return __kernel_standard (x, x, 6 + !!signbit (x)); + + return z; +} +hidden_def (__exp) +weak_alias (__exp, exp) +# ifdef NO_LONG_DOUBLE +strong_alias (__exp, __expl) +weak_alias (__exp, expl) +# endif +#endif diff --git a/math/w_expf_compat.c b/math/w_expf_compat.c new file mode 100644 index 0000000000..a38ff4082d --- /dev/null +++ b/math/w_expf_compat.c @@ -0,0 +1,37 @@ +/* Copyright (C) 2011-2017 Free Software Foundation, Inc. + This file is part of the GNU C Library. + Contributed by Ulrich Drepper <drepper@gmail.com>, 2011. + + 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 <math.h> +#include <math_private.h> +#include <math-svid-compat.h> + +#if LIBM_SVID_COMPAT +/* wrapper expf */ +float +__expf (float x) +{ + float z = __ieee754_expf (x); + if (__builtin_expect (!isfinite (z) || z == 0, 0) + && isfinite (x) && _LIB_VERSION != _IEEE_) + return __kernel_standard_f (x, x, 106 + !!signbit (x)); + + return z; +} +hidden_def (__expf) +weak_alias (__expf, expf) +#endif diff --git a/math/w_expl_compat.c b/math/w_expl_compat.c index 70096a820c..ec076deaad 100644 --- a/math/w_expl_compat.c +++ b/math/w_expl_compat.c @@ -1,13 +1,45 @@ +/* w_expl.c -- long double version of w_exp.c. + * Conversion to long double by Ulrich Drepper, + * Cygnus Support, drepper@cygnus.com. + */ + +/* + * ==================================================== + * Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved. + * + * Developed at SunPro, a Sun Microsystems, Inc. business. + * Permission to use, copy, modify, and distribute this + * software is freely granted, provided that this notice + * is preserved. + * ==================================================== + */ + +#if defined(LIBM_SCCS) && !defined(lint) +static char rcsid[] = "$NetBSD: $"; +#endif + +/* + * wrapper expl(x) + */ + #include <math.h> -#include <stdio.h> -#include <errno.h> +#include <math_private.h> +#include <math-svid-compat.h> -long double -__expl(long double x) +#if LIBM_SVID_COMPAT +long double __expl(long double x) /* wrapper exp */ { - fputs ("__expl not implemented\n", stderr); - __set_errno (ENOSYS); - return 0.0; -} +# ifdef _IEEE_LIBM + return __ieee754_expl(x); +# else + long double z = __ieee754_expl (x); + if (__glibc_unlikely (!isfinite (z) || z == 0) + && isfinite (x) && _LIB_VERSION != _IEEE_) + return __kernel_standard_l (x, x, 206 + !!signbit (x)); + return z; +# endif +} +hidden_def (__expl) weak_alias (__expl, expl) +#endif |