1 files changed, 48 insertions, 0 deletions
diff --git a/newlib/libm/mathfp/sf_isinf.c b/newlib/libm/mathfp/sf_isinf.c
new file mode 100644
index 000000000..84ab3f71c
--- /dev/null
+++ b/newlib/libm/mathfp/sf_isinf.c
@@ -0,0 +1,48 @@
+
+/* @(#)z_isinff.c 1.0 98/08/13 */
+/******************************************************************
+ * isinff
+ *
+ * Input:
+ *   x - pointer to a floating point value
+ *
+ * Output:
+ *   An integer that indicates if the number is infinite.
+ *
+ * Description:
+ *   This routine returns an integer that indicates if the number
+ *   passed in is infinite (1) or is finite (0).
+ *
+ *****************************************************************/
+
+#include "fdlibm.h"
+#include "zmath.h"
+
+int 
+_DEFUN (isinff, (float),
+	float x)
+{
+  __uint32_t wx;
+  int exp;
+
+  GET_FLOAT_WORD (wx, x);
+  exp = (wx & 0x7f800000) >> 23;
+
+  if ((exp == 0x7f8) && !(wx & 0xf0000))
+    return (1);
+  else
+    return (0);
+}
+
+#ifdef _DOUBLE_IS_32BITS
+
+int
+_DEFUN (isinf, (double),
+        double x)
+{
+  return isinff ((float) x);
+}
+
+#endif /* defined(_DOUBLE_IS_32BITS) */
+
+