This is the mail archive of the
glibc-cvs@sourceware.org
mailing list for the glibc project.
GNU C Library master sources branch, master, updated. glibc-2.14-385-g1004d18
- From: drepper at sourceware dot org
- To: glibc-cvs at sourceware dot org
- Date: 18 Oct 2011 02:31:04 -0000
- Subject: GNU C Library master sources branch, master, updated. glibc-2.14-385-g1004d18
This is an automated email from the git hooks/post-receive script. It was
generated because a ref change was pushed to the repository containing
the project "GNU C Library master sources".
The branch, master has been updated
via 1004d1821ea9093f6b75e68047dc0b3d093578f7 (commit)
from 228a984d543439308632ade9a77a171c88c11afc (commit)
Those revisions listed above that are new to this repository have
not appeared on any other notification email; so we list those
revisions in full, below.
- Log -----------------------------------------------------------------
http://sources.redhat.com/git/gitweb.cgi?p=glibc.git;a=commitdiff;h=1004d1821ea9093f6b75e68047dc0b3d093578f7
commit 1004d1821ea9093f6b75e68047dc0b3d093578f7
Author: Ulrich Drepper <drepper@gmail.com>
Date: Mon Oct 17 22:30:44 2011 -0400
Add optimized 64-bit frexp
diff --git a/ChangeLog b/ChangeLog
index 081b414..c47826c 100644
--- a/ChangeLog
+++ b/ChangeLog
@@ -1,5 +1,7 @@
2011-10-17 Ulrich Drepper <drepper@gmail.com>
+ * sysdeps/ieee754/dbl-64/wordsize-64/s_frexp.c: New file.
+
* sysdeps/x86_64/fpu/math_private.h: Relax asm requirements for
recently added interfaces.
* sysdeps/x86_64/fpu/bits/mathinline.h: Likewise.
diff --git a/sysdeps/ieee754/dbl-64/wordsize-64/s_frexp.c b/sysdeps/ieee754/dbl-64/wordsize-64/s_frexp.c
new file mode 100644
index 0000000..76f7610
--- /dev/null
+++ b/sysdeps/ieee754/dbl-64/wordsize-64/s_frexp.c
@@ -0,0 +1,67 @@
+/* Copyright (C) 2011 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, write to the Free
+ Software Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA
+ 02111-1307 USA. */
+
+#include <inttypes.h>
+#include "math.h"
+#include "math_private.h"
+
+/*
+ * for non-zero, finite x
+ * x = frexp(arg,&exp);
+ * return a double fp quantity x such that 0.5 <= |x| <1.0
+ * and the corresponding binary exponent "exp". That is
+ * arg = x*2^exp.
+ * If arg is inf, 0.0, or NaN, then frexp(arg,&exp) returns arg
+ * with *exp=0.
+ */
+
+
+double
+__frexp (double x, int *eptr)
+{
+ int64_t ix;
+ EXTRACT_WORDS64 (ix, x);
+ int32_t ex = 0x7ff & (ix >> 52);
+ int e = 0;
+
+ if (__builtin_expect (ex != 0x7ff && x != 0.0, 1))
+ {
+ /* Not zero and finite. */
+ e = ex - 1022;
+ if (__builtin_expect (ex == 0, 0))
+ {
+ /* Subnormal. */
+ x *= 0x1p54;
+ EXTRACT_WORDS64 (ix, x);
+ ex = 0x7ff & (ix >> 52);
+ e = ex - 1022 - 54;
+ }
+
+ ix = (ix & INT64_C (0x800fffffffffffff)) | INT64_C (0x3fe0000000000000);
+ INSERT_WORDS64 (x, ix);
+ }
+
+ *eptr = e;
+ return x;
+}
+weak_alias (__frexp, frexp)
+#ifdef NO_LONG_DOUBLE
+strong_alias (__frexp, __frexpl)
+weak_alias (__frexp, frexpl)
+#endif
-----------------------------------------------------------------------
Summary of changes:
ChangeLog | 2 +
.../dbl-64/{e_atanh.c => wordsize-64/s_frexp.c} | 73 +++++++++----------
2 files changed, 37 insertions(+), 38 deletions(-)
copy sysdeps/ieee754/dbl-64/{e_atanh.c => wordsize-64/s_frexp.c} (50%)
hooks/post-receive
--
GNU C Library master sources