about summary refs log tree commit diff
path: root/sysdeps/x86_64/fpu/s_ceill.S
blob: 16dbecd56dd331200e381fe96724060701a735e5 (plain) (blame)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
/*
 * Public domain.
 */

#include <libm-alias-ldouble.h>
#include <machine/asm.h>


ENTRY(__ceill)
	fldt	8(%rsp)

	fnstenv	-28(%rsp)		/* store fpu environment */

	/* 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	-28(%rsp),%edx
	andl	$0xfbff,%edx
	movl	%edx,-32(%rsp)
	fldcw	-32(%rsp)		/* load modified control word */

	frndint				/* round */

	/* Preserve "invalid" exceptions from sNaN input.  */
	fnstsw
	andl	$0x1, %eax
	orl	%eax, -24(%rsp)

	fldenv	-28(%rsp)		/* restore original environment */

	ret
END (__ceill)
libm_alias_ldouble (__ceil, ceil)