+2015-04-13 Stefan Liebler <stli@linux.vnet.ibm.com>
+
+ [BZ #6792]
+ * math/w_log1p.c: New file.
+ * math/w_log1pf.c: Likewise.
+ * math/w_log1pl.c: Likewise.
+ * math/Makefile (libm-calls): Add w_log1p.
+ * math/s_log1pl.c (log1pl): Remove weak_alias.
+ * sysdeps/i386/fpu/s_log1p.S (log1p): Likewise.
+ * sysdeps/i386/fpu/s_log1pf.S (log1pf): Likewise.
+ * sysdeps/i386/fpu/s_log1pl.S (log1pl): Likewise.
+ * sysdeps/x86_64/fpu/s_log1pl.S (log1pl): Likewise.
+ * sysdeps/ieee754/dbl-64/s_log1p.c (log1p): Likewise.
+ [NO_LONG_DOUBLE] (log1pl): Likewise.
+ * sysdeps/ieee754/flt-32/s_log1pf.c (log1pf): Likewise.
+ * sysdeps/ieee754/ldbl-128/s_log1pl.c (log1pl): Likewise.
+ * sysdeps/ieee754/ldbl-64-128/s_log1pl.c
+ (log1p): Remove long_double_symbol.
+ * sysdeps/ieee754/ldbl-128ibm/s_log1pl.c (log1pl): Likewise.
+ * sysdeps/ieee754/ldbl-64-128/w_log1pl.c: New file.
+ * sysdeps/ieee754/ldbl-128ibm/w_log1pl.c: Likewise.
+ * sysdeps/m68k/m680x0/fpu/s_log1p.c: Define empty weak_alias to
+ remove weak_alias for corresponding log1p function.
+ * sysdeps/m68k/m680x0/fpu/s_log1pf.c: Likewise.
+ * sysdeps/m68k/m680x0/fpu/s_log1pl.c: Likewise.
+ * sysdeps/ia64/fpu/w_log1p.c: New file.
+ * sysdeps/ia64/fpu/w_log1pf.c: Likewise.
+ * sysdeps/ia64/fpu/w_log1pl.c: Likewise.
+ * math/libm-test.inc (log1p_test_data): Add errno expectations.
+
2015-04-10 Joseph Myers <joseph@codesourcery.com>
[BZ #18247]
* The following bugs are resolved with this release:
- 4719, 13064, 14094, 14841, 14906, 15319, 15467, 15790, 15969, 16351,
+ 4719, 6792, 13064, 14094, 14841, 14906, 15319, 15467, 15790, 15969, 16351,
16512, 16560, 16783, 16850, 17090, 17195, 17269, 17523, 17542, 17569,
17588, 17596, 17620, 17621, 17628, 17631, 17711, 17776, 17779, 17792,
17836, 17912, 17916, 17930, 17932, 17944, 17949, 17964, 17965, 17967,
e_ilogb \
k_cos k_rem_pio2 k_sin k_tan s_asinh s_atan s_cbrt \
s_ceil s_cos s_erf s_expm1 s_fabs \
- s_floor s_log1p s_logb \
+ s_floor s_log1p w_log1p s_logb \
s_nextafter s_nexttoward s_rint s_scalbln w_scalbln \
s_significand s_sin s_tan s_tanh w_acos w_acosh w_asin \
w_atan2 w_atanh w_cosh w_drem w_exp w_exp2 w_exp10 w_fmod \
static const struct test_f_f_data log1p_test_data[] =
{
- TEST_f_f (log1p, -1, minus_infty, DIVIDE_BY_ZERO_EXCEPTION),
- TEST_f_f (log1p, -2, qnan_value, INVALID_EXCEPTION),
- TEST_f_f (log1p, -max_value, qnan_value, INVALID_EXCEPTION),
- TEST_f_f (log1p, minus_infty, qnan_value, INVALID_EXCEPTION),
+ TEST_f_f (log1p, -1, minus_infty, DIVIDE_BY_ZERO_EXCEPTION|ERRNO_ERANGE),
+ TEST_f_f (log1p, -2, qnan_value, INVALID_EXCEPTION|ERRNO_EDOM),
+ TEST_f_f (log1p, -max_value, qnan_value, INVALID_EXCEPTION|ERRNO_EDOM),
+ TEST_f_f (log1p, minus_infty, qnan_value, INVALID_EXCEPTION|ERRNO_EDOM),
TEST_f_f (log1p, plus_infty, plus_infty),
TEST_f_f (log1p, qnan_value, qnan_value, NO_INEXACT_EXCEPTION),
__set_errno (ENOSYS);
return 0.0;
}
-weak_alias (__log1pl, log1pl)
stub_warning (log1pl)
--- /dev/null
+/* Wrapper for __log1p that handles setting errno.
+ Copyright (C) 2015 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/>. */
+
+#include <errno.h>
+#include <math.h>
+#include <math_private.h>
+
+double
+__w_log1p (double x)
+{
+ if (__glibc_unlikely (islessequal (x, -1.0)))
+ {
+ if (x == -1.0)
+ __set_errno (ERANGE);
+ else
+ __set_errno (EDOM);
+ }
+
+ return __log1p (x);
+}
+weak_alias (__w_log1p, log1p)
+
+#ifdef NO_LONG_DOUBLE
+strong_alias (__w_log1p, __log1pl)
+weak_alias (__w_log1p, log1pl)
+#endif
--- /dev/null
+/* Wrapper for __log1pf that handles setting errno.
+ Copyright (C) 2015 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/>. */
+
+#include <errno.h>
+#include <math.h>
+#include <math_private.h>
+
+float
+__w_log1pf (float x)
+{
+ if (__glibc_unlikely (islessequal (x, -1.0f)))
+ {
+ if (x == -1.0f)
+ __set_errno (ERANGE);
+ else
+ __set_errno (EDOM);
+ }
+
+ return __log1pf (x);
+}
+weak_alias (__w_log1pf, log1pf)
--- /dev/null
+/* Wrapper for __log1pl that handles setting errno.
+ Copyright (C) 2015 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/>. */
+
+#include <errno.h>
+#include <math.h>
+#include <math_private.h>
+
+long double
+__w_log1pl (long double x)
+{
+ if (__glibc_unlikely (islessequal (x, -1.0L)))
+ {
+ if (x == -1.0L)
+ __set_errno (ERANGE);
+ else
+ __set_errno (EDOM);
+ }
+
+ return __log1pl (x);
+}
+weak_alias (__w_log1pl, log1pl)
ret
END (__log1p)
-weak_alias (__log1p, log1p)
ret
END (__log1pf)
-weak_alias (__log1pf, log1pf)
ret
END (__log1pl)
-weak_alias (__log1pl, log1pl)
--- /dev/null
+/* Wrapper for __log1p that handles setting errno.
+ Copyright (C) 2015 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/>. */
+
+/* Nothing to do.
+ errno is set in sysdeps/ia64/fpu/s_log1p.S. */
--- /dev/null
+/* Wrapper for __log1pf that handles setting errno.
+ Copyright (C) 2015 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/>. */
+
+/* Nothing to do.
+ errno is set in sysdeps/ia64/fpu/s_log1pf.S. */
--- /dev/null
+/* Wrapper for __log1pl that handles setting errno.
+ Copyright (C) 2015 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/>. */
+
+/* Nothing to do.
+ errno is set in sysdeps/ia64/fpu/s_log1pl.S. */
else
return k * ln2_hi - ((hfsq - (s * (hfsq + R) + (k * ln2_lo + c))) - f);
}
-weak_alias (__log1p, log1p)
-#ifdef NO_LONG_DOUBLE
-strong_alias (__log1p, __log1pl)
-weak_alias (__log1p, log1pl)
-#endif
if(k==0) return f-(hfsq-s*(hfsq+R)); else
return k*ln2_hi-((hfsq-(s*(hfsq+R)+(k*ln2_lo+c)))-f);
}
-weak_alias (__log1pf, log1pf)
z = z + e * C1;
return (z);
}
-
-weak_alias (__log1pl, log1pl)
z = z + e * C1;
return (z);
}
-
-long_double_symbol (libm, __log1pl, log1pl);
--- /dev/null
+/* Wrapper for __log1pl that handles setting errno.
+ Copyright (C) 2015 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/>. */
+
+#include <math_ldbl_opt.h>
+#undef weak_alias
+#define weak_alias(n,a)
+#include <math/w_log1pl.c>
+long_double_symbol (libm, __w_log1pl, log1pl);
#include <math_ldbl_opt.h>
-#undef weak_alias
-#define weak_alias(n,a)
#include <sysdeps/ieee754/ldbl-128/s_log1pl.c>
-long_double_symbol (libm, __log1pl, log1pl);
--- /dev/null
+/* Wrapper for __log1pl that handles setting errno.
+ Copyright (C) 2015 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/>. */
+
+#include <math_ldbl_opt.h>
+#undef weak_alias
+#define weak_alias(n,a)
+#include <math/w_log1pl.c>
+long_double_symbol (libm, __w_log1pl, log1pl);
#define FUNC log1p
+#undef weak_alias
+#define weak_alias(a,b)
#include <s_atan.c>
#define FUNC log1pf
+#undef weak_alias
+#define weak_alias(a,b)
#include <s_atanf.c>
#define FUNC log1pl
+#undef weak_alias
+#define weak_alias(a,b)
#include <s_atanl.c>
ret
END (__log1pl)
-weak_alias (__log1pl, log1pl)