summaryrefslogtreecommitdiff
path: root/newlib/libm/math/s_isinf.c
diff options
context:
space:
mode:
Diffstat (limited to 'newlib/libm/math/s_isinf.c')
-rw-r--r--newlib/libm/math/s_isinf.c26
1 files changed, 26 insertions, 0 deletions
diff --git a/newlib/libm/math/s_isinf.c b/newlib/libm/math/s_isinf.c
new file mode 100644
index 00000000000..87f099566b6
--- /dev/null
+++ b/newlib/libm/math/s_isinf.c
@@ -0,0 +1,26 @@
+/*
+ * isinf(x) returns 1 if x is infinity, else 0;
+ * no branching!
+ * Added by Cygnus Support.
+ */
+
+#include "fdlibm.h"
+
+#ifndef _DOUBLE_IS_32BITS
+
+#ifdef __STDC__
+ int isinf(double x)
+#else
+ int isinf(x)
+ double x;
+#endif
+{
+ __int32_t hx,lx;
+ EXTRACT_WORDS(hx,lx,x);
+ hx &= 0x7fffffff;
+ hx |= (__uint32_t)(lx|(-lx))>>31;
+ hx = 0x7ff00000 - hx;
+ return 1 - (int)((__uint32_t)(hx|(-hx))>>31);
+}
+
+#endif /* _DOUBLE_IS_32BITS */