1 files changed, 58 insertions, 0 deletions
diff --git a/newlib/libm/mathfp/sf_frexp.c b/newlib/libm/mathfp/sf_frexp.c
new file mode 100644
index 00000000000..c2751f65c5c
--- /dev/null
+++ b/newlib/libm/mathfp/sf_frexp.c
@@ -0,0 +1,58 @@
+
+/* @(#)z_frexpf.c 1.0 98/08/13 */
+/******************************************************************
+ * frexp
+ *
+ * Input:
+ *   d   - floating point value
+ *   exp - exponent value
+ *
+ * Output:
+ *   A floating point value in the range [0.5, 1).
+ *
+ * Description:
+ *   This routine breaks a floating point value into a number f and
+ *   an exponent exp such that d = f * 2 ^ exp.
+ *
+ *****************************************************************/
+
+#include "fdlibm.h"
+#include "zmath.h"
+
+float frexpf (float d, int *exp)
+{
+  float f;
+  __int32_t wf, wd;
+
+  GET_FLOAT_WORD (wd, d);
+
+  /* Get the exponent. */
+  *exp = ((wd & 0x7f800000) >> 23) - 126;
+
+  /* Get the mantissa. */ 
+  wf = wd & 0x7fffff;  
+  wf |= 0x3f000000;
+
+  SET_FLOAT_WORD (f, wf);
+
+  /* Check for special values. */
+  switch (numtestf (f))
+    {
+      case NAN:
+      case INF:
+        errno = EDOM;
+        *exp = 0;
+        return (f);
+    }
+
+  return (f);
+}
+
+#ifdef _DOUBLE_IS_32BITS
+
+double frexp (double x, int *exp)
+{
+  return (double) frexpf ((float) x, exp);
+}
+
+#endif /* defined(_DOUBLE_IS_32BITS) */