Blame sysdeps/i386/fpu/s_remquol.S
|
Packit |
6c4009 |
/*
|
|
Packit |
6c4009 |
* Written by Ulrich Drepper <drepper@cygnus.com>.
|
|
Packit |
6c4009 |
* Based on e_remainder by J.T. Conklin <jtc@netbsd.org>.
|
|
Packit |
6c4009 |
* Public domain.
|
|
Packit |
6c4009 |
*/
|
|
Packit |
6c4009 |
|
|
Packit |
6c4009 |
#include <libm-alias-ldouble.h>
|
|
Packit |
6c4009 |
#include <machine/asm.h>
|
|
Packit |
6c4009 |
|
|
Packit |
6c4009 |
#define PARMS 4 /* no space for saved regs */
|
|
Packit |
6c4009 |
#define DVDND PARMS
|
|
Packit |
6c4009 |
#define DVSOR DVDND+12
|
|
Packit |
6c4009 |
#define QUOP DVSOR+12
|
|
Packit |
6c4009 |
|
|
Packit |
6c4009 |
.text
|
|
Packit |
6c4009 |
ENTRY (__remquol)
|
|
Packit |
6c4009 |
|
|
Packit |
6c4009 |
fldt DVSOR(%esp)
|
|
Packit |
6c4009 |
fldt DVDND(%esp)
|
|
Packit |
6c4009 |
1: fprem1
|
|
Packit |
6c4009 |
fstsw %ax
|
|
Packit |
6c4009 |
sahf
|
|
Packit |
6c4009 |
jp 1b
|
|
Packit |
6c4009 |
fstp %st(1)
|
|
Packit |
6c4009 |
/* Compute the congruent of the quotient. */
|
|
Packit |
6c4009 |
movl %eax, %ecx
|
|
Packit |
6c4009 |
shrl $8, %eax
|
|
Packit |
6c4009 |
shrl $12, %ecx
|
|
Packit |
6c4009 |
andl $4, %ecx
|
|
Packit |
6c4009 |
andl $3, %eax
|
|
Packit |
6c4009 |
orl %eax, %ecx
|
|
Packit |
6c4009 |
leal (%ecx,%ecx,2),%ecx
|
|
Packit |
6c4009 |
movl $0xef2a60, %eax
|
|
Packit |
6c4009 |
shrl %cl, %eax
|
|
Packit |
6c4009 |
andl $7, %eax
|
|
Packit |
6c4009 |
movl QUOP(%esp), %ecx
|
|
Packit |
6c4009 |
movl DVDND+8(%esp), %edx
|
|
Packit |
6c4009 |
xorl DVSOR+8(%esp), %edx
|
|
Packit |
6c4009 |
testl $0x8000, %edx
|
|
Packit |
6c4009 |
jz 1f
|
|
Packit |
6c4009 |
negl %eax
|
|
Packit |
6c4009 |
1: movl %eax, (%ecx)
|
|
Packit |
6c4009 |
|
|
Packit |
6c4009 |
ret
|
|
Packit |
6c4009 |
END (__remquol)
|
|
Packit |
6c4009 |
libm_alias_ldouble (__remquo, remquo)
|