summaryrefslogtreecommitdiff
diff options
context:
space:
mode:
authorRich Felker <dalias@aerifal.cx>2012-03-19 06:22:54 -0400
committerRich Felker <dalias@aerifal.cx>2012-03-19 06:22:54 -0400
commitbe5b01f855b9a6460f8ed575bea50dd01030cb50 (patch)
tree085cb865c1ee70c5549080e949254e26f3e79178
parent1bf4dad32727ee17e6067caac6600ca0c5ab9f2d (diff)
downloadmusl-be5b01f855b9a6460f8ed575bea50dd01030cb50.tar.gz
fix broken modf family functions
-rw-r--r--src/math/i386/modf.s31
-rw-r--r--src/math/i386/modff.s31
-rw-r--r--src/math/i386/modfl.s31
3 files changed, 66 insertions, 27 deletions
diff --git a/src/math/i386/modf.s b/src/math/i386/modf.s
index b88e7841..8e3ea50d 100644
--- a/src/math/i386/modf.s
+++ b/src/math/i386/modf.s
@@ -1,15 +1,28 @@
.global modf
.type modf,@function
modf:
- mov 12(%esp),%eax
+ mov 12(%esp),%ecx
fldl 4(%esp)
- fld1
- fld %st(1)
-1: fprem
- fstsw %ax
+ fxam
+ fnstsw %ax
sahf
- jp 1b
- fstp %st(1)
- fsubr %st(1)
- fstpl (%eax)
+ jnp 1f
+ jc 2f
+1: fld1
+ fld %st(1)
+1: fprem
+ fnstsw %ax
+ sahf
+ jp 1b
+ fstp %st(1)
+ fsubr %st(0),%st(1)
+ fxch %st(1)
+ fstpl (%ecx)
+ ret
+
+2: fstpl (%ecx)
+ mov 4(%ecx),%eax
+ and $0x80000000,%eax
+ mov %eax,4(%esp)
+ flds 4(%esp)
ret
diff --git a/src/math/i386/modff.s b/src/math/i386/modff.s
index d29b4b60..35a08ca9 100644
--- a/src/math/i386/modff.s
+++ b/src/math/i386/modff.s
@@ -1,15 +1,28 @@
.global modff
.type modff,@function
modff:
- mov 8(%esp),%eax
+ mov 8(%esp),%ecx
flds 4(%esp)
- fld1
- fld %st(1)
-1: fprem
- fstsw %ax
+ fxam
+ fnstsw %ax
sahf
- jp 1b
- fstp %st(1)
- fsubr %st(1)
- fstps (%eax)
+ jnp 1f
+ jc 2f
+1: fld1
+ fld %st(1)
+1: fprem
+ fnstsw %ax
+ sahf
+ jp 1b
+ fstp %st(1)
+ fsubr %st(0),%st(1)
+ fxch %st(1)
+ fstps (%ecx)
+ ret
+
+2: fstps (%ecx)
+ mov (%ecx),%eax
+ and $0x80000000,%eax
+ mov %eax,4(%esp)
+ flds 4(%esp)
ret
diff --git a/src/math/i386/modfl.s b/src/math/i386/modfl.s
index f9380083..d650debb 100644
--- a/src/math/i386/modfl.s
+++ b/src/math/i386/modfl.s
@@ -1,15 +1,28 @@
.global modfl
.type modfl,@function
modfl:
- mov 16(%esp),%eax
+ mov 16(%esp),%ecx
fldt 4(%esp)
- fld1
- fld %st(1)
-1: fprem
- fstsw %ax
+ fxam
+ fnstsw %ax
sahf
- jp 1b
- fstp %st(1)
- fsubr %st(1)
- fstpt (%eax)
+ jnp 1f
+ jc 2f
+1: fld1
+ fld %st(1)
+1: fprem
+ fnstsw %ax
+ sahf
+ jp 1b
+ fstp %st(1)
+ fsubr %st(0),%st(1)
+ fxch %st(1)
+ fstpt (%ecx)
+ ret
+
+2: fstpt (%ecx)
+ mov 6(%ecx),%eax
+ and $0x80000000,%eax
+ mov %eax,4(%esp)
+ flds 4(%esp)
ret