about summary refs log tree commit diff
diff options
context:
space:
mode:
-rw-r--r--ChangeLog2
-rw-r--r--sysdeps/ieee754/dbl-64/wordsize-64/s_finite.c32
2 files changed, 34 insertions, 0 deletions
diff --git a/ChangeLog b/ChangeLog
index b7a245b8b7..c4ea055c9d 100644
--- a/ChangeLog
+++ b/ChangeLog
@@ -1,5 +1,7 @@
 2011-09-15  Ulrich Drepper  <drepper@gmail.com>
 
+	* sysdeps/ieee754/dbl-64/wordsize-64/s_finite.c: New file.
+
 	* sysdeps/unix/sysv/linux/bits/in.h (IP_MULTICAST_ALL): Define.
 	Patch mostly by Neil Horman <nhorman@tuxdriver.com>.
 
diff --git a/sysdeps/ieee754/dbl-64/wordsize-64/s_finite.c b/sysdeps/ieee754/dbl-64/wordsize-64/s_finite.c
new file mode 100644
index 0000000000..585f63e571
--- /dev/null
+++ b/sysdeps/ieee754/dbl-64/wordsize-64/s_finite.c
@@ -0,0 +1,32 @@
+/*
+ * ====================================================
+ * 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.
+ * ====================================================
+ */
+
+/*
+ * finite(x) returns 1 is x is finite, else 0;
+ * no branching!
+ */
+
+#include "math.h"
+#include "math_private.h"
+
+int
+__finite(double x)
+{
+  int64_t lx;
+  EXTRACT_WORDS64(lx,x);
+  return (int)((uint64_t)((lx&INT64_C(0x7fffffffffffffff))-INT64_C(0x7ff0000000000000))>>63);
+}
+hidden_def (__finite)
+weak_alias (__finite, finite)
+#ifdef NO_LONG_DOUBLE
+strong_alias (__finite, __finitel)
+weak_alias (__finite, finitel)
+#endif