about summary refs log tree commit diff
path: root/src/math
diff options
context:
space:
mode:
authorRich Felker <dalias@aerifal.cx>2012-03-19 05:15:30 -0400
committerRich Felker <dalias@aerifal.cx>2012-03-19 05:15:30 -0400
commit58ff9e8eaf99f2294232be860daa2ca6f3674faf (patch)
treeb2fc48ab1eb651bdb5e023a3547747d558c888df /src/math
parentbc33e617040166e971ec1e6822ac1cc417eb6c9c (diff)
downloadmusl-58ff9e8eaf99f2294232be860daa2ca6f3674faf.tar.gz
musl-58ff9e8eaf99f2294232be860daa2ca6f3674faf.tar.xz
musl-58ff9e8eaf99f2294232be860daa2ca6f3674faf.zip
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.
Diffstat (limited to 'src/math')
-rw-r--r--src/math/i386/ldexp.s1
-rw-r--r--src/math/i386/ldexpf.s1
-rw-r--r--src/math/i386/ldexpl.s1
-rw-r--r--src/math/i386/scalbln.s1
-rw-r--r--src/math/i386/scalblnf.s1
-rw-r--r--src/math/i386/scalblnl.s1
-rw-r--r--src/math/i386/scalbn.s20
-rw-r--r--src/math/i386/scalbnf.s20
-rw-r--r--src/math/i386/scalbnl.s18
9 files changed, 64 insertions, 0 deletions
diff --git a/src/math/i386/ldexp.s b/src/math/i386/ldexp.s
new file mode 100644
index 00000000..c430f784
--- /dev/null
+++ b/src/math/i386/ldexp.s
@@ -0,0 +1 @@
+# see scalbn.s
diff --git a/src/math/i386/ldexpf.s b/src/math/i386/ldexpf.s
new file mode 100644
index 00000000..3f8e4b95
--- /dev/null
+++ b/src/math/i386/ldexpf.s
@@ -0,0 +1 @@
+# see scalbnf.s
diff --git a/src/math/i386/ldexpl.s b/src/math/i386/ldexpl.s
new file mode 100644
index 00000000..86fe5621
--- /dev/null
+++ b/src/math/i386/ldexpl.s
@@ -0,0 +1 @@
+# see scalbnl.s
diff --git a/src/math/i386/scalbln.s b/src/math/i386/scalbln.s
new file mode 100644
index 00000000..c430f784
--- /dev/null
+++ b/src/math/i386/scalbln.s
@@ -0,0 +1 @@
+# see scalbn.s
diff --git a/src/math/i386/scalblnf.s b/src/math/i386/scalblnf.s
new file mode 100644
index 00000000..3f8e4b95
--- /dev/null
+++ b/src/math/i386/scalblnf.s
@@ -0,0 +1 @@
+# see scalbnf.s
diff --git a/src/math/i386/scalblnl.s b/src/math/i386/scalblnl.s
new file mode 100644
index 00000000..86fe5621
--- /dev/null
+++ b/src/math/i386/scalblnl.s
@@ -0,0 +1 @@
+# see scalbnl.s
diff --git a/src/math/i386/scalbn.s b/src/math/i386/scalbn.s
new file mode 100644
index 00000000..e275d14f
--- /dev/null
+++ b/src/math/i386/scalbn.s
@@ -0,0 +1,20 @@
+.global ldexp
+.type ldexp,@function
+ldexp:
+	nop
+
+.global scalbln
+.type scalbln,@function
+scalbln:
+	nop
+
+.global scalbn
+.type scalbn,@function
+scalbn:
+	fildl 12(%esp)
+	fldl 4(%esp)
+	fscale
+	fstp %st(1)
+	fstpl 4(%esp)
+	fldl 4(%esp)
+	ret
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
diff --git a/src/math/i386/scalbnl.s b/src/math/i386/scalbnl.s
new file mode 100644
index 00000000..224b1bef
--- /dev/null
+++ b/src/math/i386/scalbnl.s
@@ -0,0 +1,18 @@
+.global ldexpl
+.type ldexpl,@function
+ldexpl:
+	nop
+
+.global scalblnl
+.type scalblnl,@function
+scalblnl:
+	nop
+
+.global scalbnl
+.type scalbnl,@function
+scalbnl:
+	fildl 16(%esp)
+	fldt 4(%esp)
+	fscale
+	fstp %st(1)
+	ret