about summary refs log tree commit diff
path: root/src/math
diff options
context:
space:
mode:
Diffstat (limited to 'src/math')
-rw-r--r--src/math/i386/modf.s28
-rw-r--r--src/math/i386/modff.s28
-rw-r--r--src/math/i386/modfl.s28
3 files changed, 84 insertions, 0 deletions
diff --git a/src/math/i386/modf.s b/src/math/i386/modf.s
new file mode 100644
index 00000000..8e3ea50d
--- /dev/null
+++ b/src/math/i386/modf.s
@@ -0,0 +1,28 @@
+.global modf
+.type modf,@function
+modf:
+	mov 12(%esp),%ecx
+	fldl 4(%esp)
+	fxam
+	fnstsw %ax
+	sahf
+        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
new file mode 100644
index 00000000..35a08ca9
--- /dev/null
+++ b/src/math/i386/modff.s
@@ -0,0 +1,28 @@
+.global modff
+.type modff,@function
+modff:
+	mov 8(%esp),%ecx
+	flds 4(%esp)
+	fxam
+	fnstsw %ax
+	sahf
+        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
new file mode 100644
index 00000000..d650debb
--- /dev/null
+++ b/src/math/i386/modfl.s
@@ -0,0 +1,28 @@
+.global modfl
+.type modfl,@function
+modfl:
+	mov 16(%esp),%ecx
+	fldt 4(%esp)
+	fxam
+	fnstsw %ax
+	sahf
+        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