From be5b01f855b9a6460f8ed575bea50dd01030cb50 Mon Sep 17 00:00:00 2001
From: Rich Felker <dalias@aerifal.cx>
Date: Mon, 19 Mar 2012 06:22:54 -0400
Subject: [PATCH] fix broken modf family functions

---
 src/math/i386/modf.s  | 31 ++++++++++++++++++++++---------
 src/math/i386/modff.s | 31 ++++++++++++++++++++++---------
 src/math/i386/modfl.s | 31 ++++++++++++++++++++++---------
 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
-- 
2.20.1