1 files changed, 34 insertions, 0 deletions
diff --git a/newlib/libm/machine/i386/f_ldexp.S b/newlib/libm/machine/i386/f_ldexp.S
new file mode 100644
index 000000000..af6b05825
--- /dev/null
+++ b/newlib/libm/machine/i386/f_ldexp.S
@@ -0,0 +1,34 @@
+/*
+ * ====================================================
+ * 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 ldexp using Intel float instructions.
+
+   double _f_ldexp (double x, int exp);
+
+Function calculates x * 2 ** exp.
+There is no error checking or setting of errno.
+*/
+
+	#include "i386mach.h"
+
+	.global SYM (_f_ldexp)
+       SOTYPE_FUNCTION(_f_ldexp)
+
+SYM (_f_ldexp):
+	pushl ebp
+	movl esp,ebp
+	fild 16(ebp)
+	fldl 8(ebp)
+	fscale
+	fstp st1
+
+	leave
+	ret