+2015-10-05 Joseph Myers <joseph@codesourcery.com>
+
+ [BZ #887]
+ [BZ #19049]
+ [BZ #19050]
+ * sysdeps/generic/fix-int-fp-convert-zero.h: New file.
+ * sysdeps/ieee754/dbl-64/e_log10.c: Include
+ <fix-int-fp-convert-zero.h>.
+ (__ieee754_log10): Adjust signs as needed if FIX_INT_FP_CONVERT_ZERO.
+ * sysdeps/ieee754/dbl-64/e_log2.c: Include
+ <fix-int-fp-convert-zero.h>.
+ (__ieee754_log2): Adjust signs as needed if FIX_INT_FP_CONVERT_ZERO.
+ * sysdeps/ieee754/dbl-64/s_erf.c: Include
+ <fix-int-fp-convert-zero.h>.
+ (__erfc): Adjust signs as needed if FIX_INT_FP_CONVERT_ZERO.
+ * sysdeps/ieee754/dbl-64/s_logb.c: Include
+ <fix-int-fp-convert-zero.h>.
+ (__logb): Adjust signs as needed if FIX_INT_FP_CONVERT_ZERO.
+ * sysdeps/ieee754/flt-32/e_log10f.c: Include
+ <fix-int-fp-convert-zero.h>.
+ (__ieee754_log10f): Adjust signs as needed if FIX_INT_FP_CONVERT_ZERO.
+ * sysdeps/ieee754/flt-32/e_log2f.c: Include
+ <fix-int-fp-convert-zero.h>.
+ (__ieee754_log2f): Adjust signs as needed if FIX_INT_FP_CONVERT_ZERO.
+ * sysdeps/ieee754/flt-32/s_erff.c: Include
+ <fix-int-fp-convert-zero.h>.
+ (__erfcf): Adjust signs as needed if FIX_INT_FP_CONVERT_ZERO.
+ * sysdeps/ieee754/flt-32/s_logbf.c: Include
+ <fix-int-fp-convert-zero.h>.
+ (__logbf): Adjust signs as needed if FIX_INT_FP_CONVERT_ZERO.
+ * sysdeps/ieee754/ldbl-128ibm/s_erfl.c: Include
+ <fix-int-fp-convert-zero.h>.
+ (__erfcl): Adjust signs as needed if FIX_INT_FP_CONVERT_ZERO.
+ * sysdeps/ieee754/ldbl-128ibm/s_logbl.c: Include
+ <fix-int-fp-convert-zero.h>.
+ (__logbl): Adjust signs as needed if FIX_INT_FP_CONVERT_ZERO.
+ * sysdeps/powerpc/powerpc32/fpu/configure.ac: New file.
+ * sysdeps/powerpc/powerpc32/fpu/configure: New generated file.
+ * sysdeps/powerpc/powerpc32/fpu/fix-int-fp-convert-zero.h: New
+ file.
+ * config.h.in [_LIBC] (HAVE_PPC_FCFID): New macro.
+
2015-10-03 Paul Pluzhnikov <ppluzhnikov@google.com>
* sysdeps/x86_64/fpu/libm-test-ulps: Regenerated.
* The following bugs are resolved with this release:
- 2542, 2543, 2558, 2898, 4404, 6803, 14341, 14912, 15367, 15384, 15786,
- 15918, 16141, 16296, 16347, 16415, 16517, 16519, 16520, 16521, 16620,
- 16734, 16973, 16985, 17118, 17243, 17244, 17250, 17441, 17787, 17886,
- 17887, 17905, 18084, 18086, 18240, 18265, 18370, 18421, 18480, 18525,
- 18595, 18610, 18618, 18647, 18661, 18674, 18675, 18681, 18724, 18757,
- 18778, 18781, 18787, 18789, 18790, 18795, 18796, 18803, 18820, 18823,
- 18824, 18825, 18857, 18863, 18870, 18872, 18873, 18875, 18887, 18921,
- 18951, 18952, 18956, 18961, 18966, 18967, 18969, 18970, 18977, 18980,
- 18981, 18985, 19003, 19016, 19032, 19046, 19059.
+ 887, 2542, 2543, 2558, 2898, 4404, 6803, 14341, 14912, 15367, 15384,
+ 15786, 15918, 16141, 16296, 16347, 16415, 16517, 16519, 16520, 16521,
+ 16620, 16734, 16973, 16985, 17118, 17243, 17244, 17250, 17441, 17787,
+ 17886, 17887, 17905, 18084, 18086, 18240, 18265, 18370, 18421, 18480,
+ 18525, 18595, 18610, 18618, 18647, 18661, 18674, 18675, 18681, 18724,
+ 18757, 18778, 18781, 18787, 18789, 18790, 18795, 18796, 18803, 18820,
+ 18823, 18824, 18825, 18857, 18863, 18870, 18872, 18873, 18875, 18887,
+ 18921, 18951, 18952, 18956, 18961, 18966, 18967, 18969, 18970, 18977,
+ 18980, 18981, 18985, 19003, 19016, 19032, 19046, 19049, 19050, 19059.
* The obsolete header <regexp.h> has been removed. Programs that require
this header must be updated to use <regex.h> instead.
/* The PowerPC64 ELFv2 ABI is being used. */
#undef HAVE_ELFV2_ABI
+/* PowerPC32 uses fcfid for integer to floating point conversions. */
+#define HAVE_PPC_FCFID 0
+
#endif
--- /dev/null
+/* Fix for conversion of integer 0 to floating point. Generic version.
+ 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/>. */
+
+#ifndef FIX_INT_FP_CONVERT_ZERO_H
+#define FIX_INT_FP_CONVERT_ZERO_H 1
+
+/* Define this macro to 1 to work around conversions of integer 0 to
+ floating point returning -0 instead of the correct +0 in some
+ rounding modes. */
+#define FIX_INT_FP_CONVERT_ZERO 0
+
+#endif /* fix-int-fp-convert-zero.h */
#include <math.h>
#include <math_private.h>
+#include <fix-int-fp-convert-zero.h>
static const double two54 = 1.80143985094819840000e+16; /* 0x43500000, 0x00000000 */
static const double ivln10 = 4.34294481903251816668e-01; /* 0x3FDBCB7B, 0x1526E50E */
i = ((u_int32_t) k & 0x80000000) >> 31;
hx = (hx & 0x000fffff) | ((0x3ff - i) << 20);
y = (double) (k + i);
+ if (FIX_INT_FP_CONVERT_ZERO && y == 0.0)
+ y = 0.0;
SET_HIGH_WORD (x, hx);
z = y * log10_2lo + ivln10 * __ieee754_log (x);
return z + y * log10_2hi;
#include <math.h>
#include <math_private.h>
+#include <fix-int-fp-convert-zero.h>
static const double ln2 = 0.69314718055994530942;
static const double two54 = 1.80143985094819840000e+16; /* 43500000 00000000 */
if ((0x000fffff & (2 + hx)) < 3)
{ /* |f| < 2**-20 */
if (f == zero)
- return dk;
+ {
+ if (FIX_INT_FP_CONVERT_ZERO && dk == 0.0)
+ dk = 0.0;
+ return dk;
+ }
R = f * f * (0.5 - 0.33333333333333333 * f);
return dk - (R - f) / ln2;
}
#include <float.h>
#include <math.h>
#include <math_private.h>
+#include <fix-int-fp-convert-zero.h>
static const double
tiny = 1e-300,
ix = hx & 0x7fffffff;
if (ix >= 0x7ff00000) /* erfc(nan)=nan */
{ /* erfc(+-inf)=0,2 */
- return (double) (((u_int32_t) hx >> 31) << 1) + one / x;
+ double ret = (double) (((u_int32_t) hx >> 31) << 1) + one / x;
+ if (FIX_INT_FP_CONVERT_ZERO && ret == 0.0)
+ return 0.0;
+ return ret;
}
if (ix < 0x3feb0000) /* |x|<0.84375 */
#include <math.h>
#include <math_private.h>
+#include <fix-int-fp-convert-zero.h>
double
__logb (double x)
ma = __builtin_clz (ix);
rix -= ma - 12;
}
+ if (FIX_INT_FP_CONVERT_ZERO && rix == 1023)
+ return 0.0;
return (double) (rix - 1023);
}
weak_alias (__logb, logb)
#include <math.h>
#include <math_private.h>
+#include <fix-int-fp-convert-zero.h>
static const float
two25 = 3.3554432000e+07, /* 0x4c000000 */
i = ((u_int32_t)k&0x80000000)>>31;
hx = (hx&0x007fffff)|((0x7f-i)<<23);
y = (float)(k+i);
+ if (FIX_INT_FP_CONVERT_ZERO && y == 0.0f)
+ y = 0.0f;
SET_FLOAT_WORD(x,hx);
z = y*log10_2lo + ivln10*__ieee754_logf(x);
return z+y*log10_2hi;
#include <math.h>
#include <math_private.h>
+#include <fix-int-fp-convert-zero.h>
static const float
ln2 = 0.69314718055994530942,
dk = (float)k;
f = x-(float)1.0;
if((0x007fffff&(15+ix))<16) { /* |f| < 2**-20 */
- if(f==zero) return dk;
+ if(f==zero)
+ {
+ if (FIX_INT_FP_CONVERT_ZERO && dk == 0.0f)
+ dk = 0.0f;
+ return dk;
+ }
R = f*f*((float)0.5-(float)0.33333333333333333*f);
return dk-(R-f)/ln2;
}
#include <float.h>
#include <math.h>
#include <math_private.h>
+#include <fix-int-fp-convert-zero.h>
static const float
tiny = 1e-30,
ix = hx&0x7fffffff;
if(ix>=0x7f800000) { /* erfc(nan)=nan */
/* erfc(+-inf)=0,2 */
- return (float)(((u_int32_t)hx>>31)<<1)+one/x;
+ float ret = (float)(((u_int32_t)hx>>31)<<1)+one/x;
+ if (FIX_INT_FP_CONVERT_ZERO && ret == 0.0f)
+ return 0.0f;
+ return ret;
}
if(ix < 0x3f580000) { /* |x|<0.84375 */
#include <math.h>
#include <math_private.h>
+#include <fix-int-fp-convert-zero.h>
float
__logbf (float x)
though it were normalized. */
rix -= __builtin_clz (ix) - 9;
}
+ if (FIX_INT_FP_CONVERT_ZERO && rix == 127)
+ return 0.0f;
return (float) (rix - 127);
}
weak_alias (__logbf, logbf)
#include <math.h>
#include <math_private.h>
#include <math_ldbl_opt.h>
+#include <fix-int-fp-convert-zero.h>
/* Evaluate P[n] x^n + P[n-1] x^(n-1) + ... + P[0] */
if (ix >= 0x7ff00000)
{ /* erfc(nan)=nan */
/* erfc(+-inf)=0,2 */
- return (long double) ((hx >> 31) << 1) + one / x;
+ long double ret = (long double) ((hx >> 31) << 1) + one / x;
+ if (FIX_INT_FP_CONVERT_ZERO && ret == 0.0L)
+ return 0.0L;
+ return ret;
}
if (ix < 0x3fd00000) /* |x| <1/4 */
#include <math.h>
#include <math_private.h>
#include <math_ldbl_opt.h>
+#include <fix-int-fp-convert-zero.h>
long double
__logbl (long double x)
if ((hxs ^ lx) < 0 && (lx & 0x7fffffffffffffffLL) != 0)
rhx--;
}
+ if (FIX_INT_FP_CONVERT_ZERO && rhx == 1023)
+ return 0.0L;
return (long double) (rhx - 1023);
}
#ifndef __logbl
--- /dev/null
+# This file is generated from configure.ac by Autoconf. DO NOT EDIT!
+ # Local configure fragment for sysdeps/powerpc/powerpc32/fpu.
+
+# Test whether integer to floating point conversions use fcfid.
+{ $as_echo "$as_me:${as_lineno-$LINENO}: checking for fcfid use" >&5
+$as_echo_n "checking for fcfid use... " >&6; }
+if ${libc_cv_ppc_fcfid+:} false; then :
+ $as_echo_n "(cached) " >&6
+else
+ echo 'double foo (int x) { return (double) x; }' > conftest.c
+libc_cv_ppc_fcfid=no
+if { ac_try='${CC-cc} -S $CFLAGS conftest.c -o conftest.s 1>&5'
+ { { eval echo "\"\$as_me\":${as_lineno-$LINENO}: \"$ac_try\""; } >&5
+ (eval $ac_try) 2>&5
+ ac_status=$?
+ $as_echo "$as_me:${as_lineno-$LINENO}: \$? = $ac_status" >&5
+ test $ac_status = 0; }; }; then
+ if grep '[ ]fcfid' conftest.s > /dev/null 2>&1; then
+ libc_cv_ppc_fcfid=yes
+ fi
+fi
+rm -rf conftest*
+fi
+{ $as_echo "$as_me:${as_lineno-$LINENO}: result: $libc_cv_ppc_fcfid" >&5
+$as_echo "$libc_cv_ppc_fcfid" >&6; }
+if test $libc_cv_ppc_fcfid = yes; then
+ $as_echo "#define HAVE_PPC_FCFID 1" >>confdefs.h
+
+fi
--- /dev/null
+GLIBC_PROVIDES dnl See aclocal.m4 in the top level source directory.
+# Local configure fragment for sysdeps/powerpc/powerpc32/fpu.
+
+# Test whether integer to floating point conversions use fcfid.
+AC_CACHE_CHECK([for fcfid use], [libc_cv_ppc_fcfid], [dnl
+echo 'double foo (int x) { return (double) x; }' > conftest.c
+libc_cv_ppc_fcfid=no
+if AC_TRY_COMMAND(${CC-cc} -S $CFLAGS conftest.c -o conftest.s 1>&AS_MESSAGE_LOG_FD); then
+changequote(,)dnl
+ if grep '[ ]fcfid' conftest.s > /dev/null 2>&1; then
+ libc_cv_ppc_fcfid=yes
+ fi
+changequote([,])dnl
+fi
+rm -rf conftest*])
+if test $libc_cv_ppc_fcfid = yes; then
+ AC_DEFINE([HAVE_PPC_FCFID])
+fi
--- /dev/null
+/* Fix for conversion of integer 0 to floating point. PowerPC version.
+ 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/>. */
+
+#ifndef FIX_INT_FP_CONVERT_ZERO_H
+#define FIX_INT_FP_CONVERT_ZERO_H 1
+
+/* The code sequences GCC generates for conversion of integers to
+ floating point result in -0 instead of +0 in FE_DOWNWARD mode when
+ the fcfid instruction is not used, as of GCC 5. See
+ <https://gcc.gnu.org/bugzilla/show_bug.cgi?id=67771>. */
+#define FIX_INT_FP_CONVERT_ZERO (!HAVE_PPC_FCFID)
+
+#endif /* fix-int-fp-convert-zero.h */
projects / thirdparty / glibc.git / commitdiff
