summaryrefslogtreecommitdiff
path: root/newlib/libc/machine/i386/f_frexpf.S
diff options
context:
space:
mode:
Diffstat (limited to 'newlib/libc/machine/i386/f_frexpf.S')
-rw-r--r--newlib/libc/machine/i386/f_frexpf.S43
1 files changed, 43 insertions, 0 deletions
diff --git a/newlib/libc/machine/i386/f_frexpf.S b/newlib/libc/machine/i386/f_frexpf.S
new file mode 100644
index 00000000000..05f0555d69b
--- /dev/null
+++ b/newlib/libc/machine/i386/f_frexpf.S
@@ -0,0 +1,43 @@
+/*
+ * ====================================================
+ * Copyright (C) 1998 by Cygnus Solutions. All rights reserved.
+ *
+ * Permission to use, copy, modify, and distribute this
+ * software is freely granted, provided that this notice
+ * is preserved.
+ * ====================================================
+ */
+
+/*
+Fast version of frexpf using Intel float instructions.
+
+ float _f_frexpf (float x, int *exp);
+
+Function splits x into y * 2 ** z. It then
+returns the value of y and updates *exp with z.
+There is no error checking or setting of errno.
+*/
+
+ #include "i386mach.h"
+
+ .global SYM (_f_frexpf)
+
+SYM (_f_frexpf):
+ pushl ebp
+ movl esp,ebp
+ flds 8(ebp)
+ movl 12(ebp),eax
+
+ fxtract
+ fld1
+ fchs
+ fxch
+ fscale
+ fstp st1
+ fxch
+ fld1
+ faddp
+ fistpl 0(eax)
+
+ leave
+ ret