assembly optimizations for fmod/remainder functions

diff --git a/src/math/i386/fmod.s b/src/math/i386/fmod.s
new file mode 100644 (file)
index 0000000..069fbfe
--- /dev/null
+++ b/src/math/i386/fmod.s
@@ -0,0 +1,11 @@
+.global fmod
+.type fmod,@function
+fmod:
+       fldl 12(%esp)
+       fldl 4(%esp)
+1:     fprem
+       fstsw %ax
+       sahf
+       jp 1b
+       fstp %st(1)
+       ret

diff --git a/src/math/i386/fmodf.s b/src/math/i386/fmodf.s
new file mode 100644 (file)
index 0000000..d99c80f
--- /dev/null
+++ b/src/math/i386/fmodf.s
@@ -0,0 +1,11 @@
+.global fmodf
+.type fmodf,@function
+fmodf:
+       flds 8(%esp)
+       flds 4(%esp)
+1:     fprem
+       fstsw %ax
+       sahf
+       jp 1b
+       fstp %st(1)
+       ret

diff --git a/src/math/i386/fmodl.s b/src/math/i386/fmodl.s
new file mode 100644 (file)
index 0000000..7e07e7b
--- /dev/null
+++ b/src/math/i386/fmodl.s
@@ -0,0 +1,11 @@
+.global fmodl
+.type fmodl,@function
+fmodl:
+       fldt 16(%esp)
+       fldt 4(%esp)
+1:     fprem
+       fstsw %ax
+       sahf
+       jp 1b
+       fstp %st(1)
+       ret

diff --git a/src/math/i386/remainder.s b/src/math/i386/remainder.s
new file mode 100644 (file)
index 0000000..47ee340
--- /dev/null
+++ b/src/math/i386/remainder.s
@@ -0,0 +1,11 @@
+.global remainder
+.type remainder,@function
+remainder:
+       fldl 12(%esp)
+       fldl 4(%esp)
+1:     fprem1
+       fstsw %ax
+       sahf
+       jp 1b
+       fstp %st(1)
+       ret

diff --git a/src/math/i386/remainderf.s b/src/math/i386/remainderf.s
new file mode 100644 (file)
index 0000000..5b5fc23
--- /dev/null
+++ b/src/math/i386/remainderf.s
@@ -0,0 +1,11 @@
+.global remainderf
+.type remainderf,@function
+remainderf:
+       flds 8(%esp)
+       flds 4(%esp)
+1:     fprem1
+       fstsw %ax
+       sahf
+       jp 1b
+       fstp %st(1)
+       ret

diff --git a/src/math/i386/remainderl.s b/src/math/i386/remainderl.s
new file mode 100644 (file)
index 0000000..0097872
--- /dev/null
+++ b/src/math/i386/remainderl.s
@@ -0,0 +1,11 @@
+.global remainderl
+.type remainderl,@function
+remainderl:
+       fldt 16(%esp)
+       fldt 4(%esp)
+1:     fprem1
+       fstsw %ax
+       sahf
+       jp 1b
+       fstp %st(1)
+       ret

diff --git a/src/math/x86_64/fmodl.s b/src/math/x86_64/fmodl.s
new file mode 100644 (file)
index 0000000..ca81e60
--- /dev/null
+++ b/src/math/x86_64/fmodl.s
@@ -0,0 +1,11 @@
+.global fmodl
+.type fmodl,@function
+fmodl:
+       fldt 24(%rsp)
+       fldt 8(%rsp)
+1:     fprem
+       fstsw %ax
+       sahf
+       jp 1b
+       fstp %st(1)
+       ret

diff --git a/src/math/x86_64/remainderl.s b/src/math/x86_64/remainderl.s
new file mode 100644 (file)
index 0000000..75c1237
--- /dev/null
+++ b/src/math/x86_64/remainderl.s
@@ -0,0 +1,11 @@
+.global remainderl
+.type remainderl,@function
+remainderl:
+       fldt 24(%rsp)
+       fldt 8(%rsp)
+1:     fprem1
+       fstsw %ax
+       sahf
+       jp 1b
+       fstp %st(1)
+       ret