about summary refs log tree commit diff
path: root/sysdeps/i386/fpu/s_ceil.S
diff options
context:
space:
mode:
Diffstat (limited to 'sysdeps/i386/fpu/s_ceil.S')
-rw-r--r--sysdeps/i386/fpu/s_ceil.S32
1 files changed, 32 insertions, 0 deletions
diff --git a/sysdeps/i386/fpu/s_ceil.S b/sysdeps/i386/fpu/s_ceil.S
new file mode 100644
index 0000000000..b0159128aa
--- /dev/null
+++ b/sysdeps/i386/fpu/s_ceil.S
@@ -0,0 +1,32 @@
+/*
+ * Written by J.T. Conklin <jtc@netbsd.org>.
+ * Public domain.
+ */
+
+#include <machine/asm.h>
+
+RCSID("$NetBSD: s_ceil.S,v 1.4 1995/05/08 23:52:13 jtc Exp $")
+
+ENTRY(__ceil)
+	fldl	4(%esp)
+	subl	$8,%esp
+
+	fstcw	4(%esp)			/* 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(%esp),%edx
+	andl	$0xfbff,%edx
+	movl	%edx,(%esp)
+	fldcw	(%esp)			/* load modified control word */
+
+	frndint				/* round */
+
+	fldcw	4(%esp)			/* restore original control word */
+
+	addl	$8,%esp
+	ret
+END (__ceil)
+weak_alias (__ceil, ceil)