asm for scalbn family

unlike some implementations, these functions perform the equivalent of
gcc's -ffloat-store on the result before returning. this is necessary
to raise underflow/overflow/inexact exceptions, perform the correct
rounding with denormals, etc.

1 files changed, 20 insertions, 0 deletions

diff --git a/src/math/i386/scalbnf.s b/src/math/i386/scalbnf.s
new file mode 100644
index 00000000..40232b6a
--- /dev/null
+++ b/src/math/i386/scalbnf.s
@@ -0,0 +1,20 @@
+.global ldexpf
+.type ldexpf,@function
+ldexpf:
+	nop
+
+.global scalblnf
+.type scalblnf,@function
+scalblnf:
+	nop
+
+.global scalbnf
+.type scalbnf,@function
+scalbnf:
+	fildl 8(%esp)
+	flds 4(%esp)
+	fscale
+	fstp %st(1)
+	fstps 4(%esp)
+	flds 4(%esp)
+	ret