summaryrefslogtreecommitdiff
path: root/src/math/i386
diff options
context:
space:
mode:
authorSzabolcs Nagy <nsz@port70.net>2012-12-16 17:30:29 +0100
committerSzabolcs Nagy <nsz@port70.net>2012-12-16 17:30:29 +0100
commit58bba42d1bd14e1ab01f3249ffc98afdbf841a6a (patch)
tree8790505fc53cfcd7a3bd6ecd4b1f7877018e2e17 /src/math/i386
parent525ad96e0ef6e46b6ac28438ef2ffa1d5317fb64 (diff)
downloadmusl-58bba42d1bd14e1ab01f3249ffc98afdbf841a6a.tar.gz
math: x86_64 version of expl, fixed some comments in the i386 version
Diffstat (limited to 'src/math/i386')
-rw-r--r--src/math/i386/expl.s4
-rw-r--r--src/math/i386/expm1f.s2
-rw-r--r--src/math/i386/expm1l.s2
3 files changed, 4 insertions, 4 deletions
diff --git a/src/math/i386/expl.s b/src/math/i386/expl.s
index 3f2f707d..8ceb40d2 100644
--- a/src/math/i386/expl.s
+++ b/src/math/i386/expl.s
@@ -1,6 +1,7 @@
# exp(x) = 2^hi + 2^hi (2^lo - 1)
# where hi+lo = log2e*x with 128bit precision
# exact log2e*x calculation depends on nearest rounding mode
+# using the exact multiplication method of Dekker and Veltkamp
.global expl
.type expl,@function
@@ -35,7 +36,6 @@ expl:
# should be 0x1.71547652b82fe178p0 == 0x3fff b8aa3b29 5c17f0bc
# it will be wrong on non-nearest rounding mode
3: fldl2e
-# subl $32, %esp
subl $44, %esp
# hi = log2e_hi*x
# 2^hi = exp2l(hi)
@@ -54,7 +54,7 @@ expl:
fldt 16(%esp)
# fpu stack: 2^hi x hi
# exact mult: x*log2e
- fld %st(1) # x
+ fld %st(1)
# c = 0x1p32+1
pushl $0x41f00000
pushl $0x00100000
diff --git a/src/math/i386/expm1f.s b/src/math/i386/expm1f.s
index 89c2d290..f335a3e5 100644
--- a/src/math/i386/expm1f.s
+++ b/src/math/i386/expm1f.s
@@ -1 +1 @@
-# see expm1.s
+# see exp.s
diff --git a/src/math/i386/expm1l.s b/src/math/i386/expm1l.s
index 89c2d290..f335a3e5 100644
--- a/src/math/i386/expm1l.s
+++ b/src/math/i386/expm1l.s
@@ -1 +1 @@
-# see expm1.s
+# see exp.s