assembly optimizations for fmod/remainder functions

[musl] / src / math / x86_64 / fmodl.s

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