45 lines
720 B
ArmAsm
45 lines
720 B
ArmAsm
|
/*
|
||
|
|
* Public domain.
|
||
|
|
*/
|
||
|
|
|
||
|
|
#include <libm-alias-ldouble.h>
|
||
|
|
#include <machine/asm.h>
|
||
|
|
|
||
|
|
#define PARMS 4 /* no space for saved regs */
|
||
|
|
#define DVDND PARMS
|
||
|
|
#define DVSOR DVDND+12
|
||
|
|
#define QUOP DVSOR+12
|
||
|
|
|
||
|
|
.text
|
||
|
|
ENTRY (__remquol)
|
||
|
|
|
||
|
|
fldt DVSOR(%esp)
|
||
|
|
fldt DVDND(%esp)
|
||
|
|
1: fprem1
|
||
|
|
fstsw %ax
|
||
|
|
sahf
|
||
|
|
jp 1b
|
||
|
|
fstp %st(1)
|
||
|
|
/* Compute the congruent of the quotient. */
|
||
|
|
movl %eax, %ecx
|
||
|
|
shrl $8, %eax
|
||
|
|
shrl $12, %ecx
|
||
|
|
andl $4, %ecx
|
||
|
|
andl $3, %eax
|
||
|
|
orl %eax, %ecx
|
||
|
|
leal (%ecx,%ecx,2),%ecx
|
||
|
|
movl $0xef2a60, %eax
|
||
|
|
shrl %cl, %eax
|
||
|
|
andl $7, %eax
|
||
|
|
movl QUOP(%esp), %ecx
|
||
|
|
movl DVDND+8(%esp), %edx
|
||
|
|
xorl DVSOR+8(%esp), %edx
|
||
|
|
testl $0x8000, %edx
|
||
|
|
jz 1f
|
||
|
|
negl %eax
|
||
|
|
1: movl %eax, (%ecx)
|
||
|
|
|
||
|
|
ret
|
||
|
|
END (__remquol)
|
||
|
|
libm_alias_ldouble (__remquo, remquo)
|