aboutsummaryrefslogtreecommitdiff
path: root/sysdeps/x86_64/fpu/s_ceill.S
diff options
context:
space:
mode:
Diffstat (limited to 'sysdeps/x86_64/fpu/s_ceill.S')
-rw-r--r--sysdeps/x86_64/fpu/s_ceill.S31
1 files changed, 31 insertions, 0 deletions
diff --git a/sysdeps/x86_64/fpu/s_ceill.S b/sysdeps/x86_64/fpu/s_ceill.S
new file mode 100644
index 0000000000..efc8dd995e
--- /dev/null
+++ b/sysdeps/x86_64/fpu/s_ceill.S
@@ -0,0 +1,31 @@
+/*
+ * Written by J.T. Conklin <jtc@netbsd.org>.
+ * Changes for long double by Ulrich Drepper <drepper@cygnus.com>
+ * Changes for x86-64 by Andreas Jaeger <aj@suse.de>
+ * Public domain.
+ */
+
+#include <machine/asm.h>
+
+
+ENTRY(__ceill)
+ fldt 8(%rsp)
+
+ fstcw -4(%rsp) /* store fpu control word */
+
+ /* We use here %edx although only the low 1 bits are defined.
+ But none of the operations should care and they are faster
+ than the 16 bit operations. */
+ movl $0x0800,%edx /* round towards +oo */
+ orl -4(%rsp),%edx
+ andl $0xfbff,%edx
+ movl %edx,-8(%rsp)
+ fldcw -8(%rsp) /* load modified control word */
+
+ frndint /* round */
+
+ fldcw -4(%rsp) /* restore original control word */
+
+ ret
+END (__ceill)
+weak_alias (__ceill, ceill)