Print this page
*** 26,36 ****
/*
* Copyright 2011, Richard Lowe
*/
! /* Functions in this file are duplicated in libm.m4. Keep them in sync */
#ifndef _LIBM_INLINES_H
#define _LIBM_INLINES_H
#ifdef __GNUC__
--- 26,36 ----
/*
* Copyright 2011, Richard Lowe
*/
! /* Functions in this file are duplicated in locallibm.il. Keep them in sync */
#ifndef _LIBM_INLINES_H
#define _LIBM_INLINES_H
#ifdef __GNUC__
*** 45,113 ****
#define _LO_WORD(x) ((uint32_t *)&x)[0]
#define _HI_WORD(x) ((uint32_t *)&x)[1]
#define _HIER_WORD(x) ((uint32_t *)&x)[2]
extern __inline__ double
- __ieee754_sqrt(double a)
- {
- double ret;
-
- __asm__ __volatile__("fsqrt\n\t" : "=t" (ret) : "0" (a));
- return (ret);
- }
-
- extern __inline__ double
__inline_sqrt(double a)
{
double ret;
! __asm__ __volatile__("fsqrt\n\t" : "=t" (ret) : "0" (a));
return (ret);
}
extern __inline__ double
! __d_sqrt_(double *a)
{
! double ret;
!
! __asm__ __volatile__("fsqrt\n\t" : "=t" (ret) : "0" (*a));
! return (ret);
}
extern __inline__ float
__inline_sqrtf(float a)
{
float ret;
! __asm__ __volatile__("fsqrt\n\t" : "=t" (ret) : "0" (a));
return (ret);
}
extern __inline__ double
__inline_rint(double a)
{
- double ret;
-
__asm__ __volatile__(
! "andl $0x7fffffff,%2\n\t"
! "cmpl $0x43300000,%2\n\t"
"jae 1f\n\t"
"frndint\n\t"
"1: fwait\n\t"
! : "=t" (ret)
! : "0" (a), "r" (_HI_WORD(a)));
! return (ret);
! }
!
! extern __inline__ short
! __inline_fstsw(void)
! {
! short ret;
!
! __asm__ __volatile__("fstsw %0\n\t" : "=r" (ret));
! return (ret);
}
/*
* 00 - 24 bits
* 01 - reserved
--- 45,91 ----
#define _LO_WORD(x) ((uint32_t *)&x)[0]
#define _HI_WORD(x) ((uint32_t *)&x)[1]
#define _HIER_WORD(x) ((uint32_t *)&x)[2]
extern __inline__ double
__inline_sqrt(double a)
{
double ret;
! __asm__ __volatile__("fsqrt\n\t" : "=t" (ret) : "0" (a) : "cc");
return (ret);
}
extern __inline__ double
! __ieee754_sqrt(double a)
{
! return (__inline_sqrt(a));
}
extern __inline__ float
__inline_sqrtf(float a)
{
float ret;
! __asm__ __volatile__("fsqrt\n\t" : "=t" (ret) : "0" (a) : "cc");
return (ret);
}
extern __inline__ double
__inline_rint(double a)
{
__asm__ __volatile__(
! "andl $0x7fffffff,%1\n\t"
! "cmpl $0x43300000,%1\n\t"
"jae 1f\n\t"
"frndint\n\t"
"1: fwait\n\t"
! : "+t" (a), "+&r" (_HI_WORD(a))
! :
! : "cc");
! return (a);
}
/*
* 00 - 24 bits
* 01 - reserved
*** 153,331 ****
}
extern __inline__ double
ceil(double d)
{
- double ret;
short rd = __swap87RD(fp_positive);
! __asm__ __volatile__("frndint" : "=t" (ret) : "0" (d));
__swap87RD(rd);
! return (ret);
}
extern __inline__ double
copysign(double d1, double d2)
{
__asm__ __volatile__(
"andl $0x7fffffff,%0\n\t" /* %0 <-- hi_32(abs(d)) */
"andl $0x80000000,%1\n\t" /* %1[31] <-- sign_bit(d2) */
"orl %1,%0\n\t" /* %0 <-- hi_32(copysign(x,y)) */
! : "+r" (_HI_WORD(d1))
! : "r" (_HI_WORD(d2)));
return (d1);
}
extern __inline__ double
fabs(double d)
{
! double ret;
!
! __asm__ __volatile__("fabs\n\t" : "=t" (ret) : "0" (d));
! return (ret);
}
extern __inline__ float
fabsf(float d)
{
! float ret;
!
! __asm__ __volatile__("fabs\n\t" : "=t" (ret) : "0" (d));
! return (ret);
}
extern __inline__ long double
fabsl(long double d)
{
! long double ret;
!
! __asm__ __volatile__("fabs\n\t" : "=t" (ret) : "0" (d));
! return (ret);
}
extern __inline__ int
finite(double d)
{
! int ret;
__asm__ __volatile__(
! "notl %1\n\t"
! "andl $0x7ff00000,%1\n\t"
! "negl %1\n\t"
! "shrl $31,%1\n\t"
! : "=r" (ret)
! : "0" (_HI_WORD(d)));
return (ret);
}
extern __inline__ double
floor(double d)
{
- double ret;
short rd = __swap87RD(fp_negative);
! __asm__ __volatile__("frndint" : "=t" (ret) : "0" (d));
__swap87RD(rd);
! return (ret);
! }
!
! /*
! * branchless __isnan
! * ((0x7ff00000-[((lx|-lx)>>31)&1]|ahx)>>31)&1 = 1 iff x is NaN
! */
! extern __inline__ int
! isnan(double d)
! {
! int ret;
!
! __asm__ __volatile__(
! "movl %1,%%ecx\n\t"
! "negl %%ecx\n\t" /* ecx <-- -lo_32(x) */
! "orl %%ecx,%1\n\t"
! "shrl $31,%1\n\t" /* 1 iff lx != 0 */
! "andl $0x7fffffff,%2\n\t" /* ecx <-- hi_32(abs(x)) */
! "orl %2,%1\n\t"
! "subl $0x7ff00000,%1\n\t"
! "negl %1\n\t"
! "shrl $31,%1\n\t"
! : "=r" (ret)
! : "0" (_HI_WORD(d)), "r" (_LO_WORD(d))
! : "ecx");
!
! return (ret);
}
extern __inline__ int
isnanf(float f)
{
- int ret;
-
__asm__ __volatile__(
"andl $0x7fffffff,%0\n\t"
"negl %0\n\t"
"addl $0x7f800000,%0\n\t"
"shrl $31,%0\n\t"
! : "=r" (ret)
! : "0" (f));
!
! return (ret);
! }
!
! extern __inline__ int
! isinf(double d)
! {
! int ret;
! __asm__ __volatile__(
! "andl $0x7fffffff,%1\n\t" /* set first bit to 0 */
! "cmpl $0x7ff00000,%1\n\t"
! "pushfl\n\t"
! "popl %0\n\t"
! "cmpl $0,%2\n\t" /* is lo_32(x) = 0? */
! "pushfl\n\t"
! "popl %2\n\t" /* bit 6 of ecx <-- lo_32(x) == 0 */
! "andl %2,%0\n\t"
! "andl $0x40,%0\n\t"
! "shrl $6,%0\n\t"
! : "=r" (ret)
! : "0" (_HI_WORD(d)), "r" (_LO_WORD(d)));
!
! return (ret);
}
extern __inline__ double
rint(double a) {
! double ret;
!
! __asm__ __volatile__(
! "andl $0x7fffffff,%2\n\t"
! "cmpl $0x43300000,%2\n\t"
! "jae 1f\n\t"
! "frndint\n\t"
! "1: fwait\n\t"
! : "=t" (ret)
! : "0" (a), "r" (_HI_WORD(a)));
!
! return (ret);
}
extern __inline__ double
scalbn(double d, int n)
{
! double ret, dummy;
__asm__ __volatile__(
! "fildl %3\n\t" /* Convert N to extended */
"fxch\n\t"
"fscale\n\t"
! : "=t" (ret), "=u" (dummy)
! : "0" (d), "m" (n));
! return (ret);
}
extern __inline__ int
signbit(double d)
{
--- 131,244 ----
}
extern __inline__ double
ceil(double d)
{
short rd = __swap87RD(fp_positive);
! __asm__ __volatile__("frndint" : "+t" (d), "+r" (rd) : : "cc");
__swap87RD(rd);
! return (d);
}
extern __inline__ double
copysign(double d1, double d2)
{
__asm__ __volatile__(
"andl $0x7fffffff,%0\n\t" /* %0 <-- hi_32(abs(d)) */
"andl $0x80000000,%1\n\t" /* %1[31] <-- sign_bit(d2) */
"orl %1,%0\n\t" /* %0 <-- hi_32(copysign(x,y)) */
! : "+&r" (_HI_WORD(d1)), "+r" (_HI_WORD(d2))
! :
! : "cc");
return (d1);
}
extern __inline__ double
fabs(double d)
{
! __asm__ __volatile__("fabs\n\t" : "+t" (d) : : "cc");
! return (d);
}
extern __inline__ float
fabsf(float d)
{
! __asm__ __volatile__("fabs\n\t" : "+t" (d) : : "cc");
! return (d);
}
extern __inline__ long double
fabsl(long double d)
{
! __asm__ __volatile__("fabs\n\t" : "+t" (d) : : "cc");
! return (d);
}
extern __inline__ int
finite(double d)
{
! int ret = _HI_WORD(d);
__asm__ __volatile__(
! "notl %0\n\t"
! "andl $0x7ff00000,%0\n\t"
! "negl %0\n\t"
! "shrl $31,%0\n\t"
! : "+r" (ret)
! :
! : "cc");
return (ret);
}
extern __inline__ double
floor(double d)
{
short rd = __swap87RD(fp_negative);
! __asm__ __volatile__("frndint" : "+t" (d), "+r" (rd) : : "cc");
__swap87RD(rd);
! return (d);
}
extern __inline__ int
isnanf(float f)
{
__asm__ __volatile__(
"andl $0x7fffffff,%0\n\t"
"negl %0\n\t"
"addl $0x7f800000,%0\n\t"
"shrl $31,%0\n\t"
! : "+r" (f)
! :
! : "cc");
! return (f);
}
extern __inline__ double
rint(double a) {
! return (__inline_rint(a));
}
extern __inline__ double
scalbn(double d, int n)
{
! double dummy;
__asm__ __volatile__(
! "fildl %2\n\t" /* Convert N to extended */
"fxch\n\t"
"fscale\n\t"
! : "+t" (d), "=u" (dummy)
! : "m" (n)
! : "cc");
! return (d);
}
extern __inline__ int
signbit(double d)
{
*** 339,394 ****
}
extern __inline__ double
sqrt(double d)
{
! double ret;
! __asm__ __volatile__("fsqrt" : "=t" (ret) : "0" (d));
! return (ret);
}
extern __inline__ float
sqrtf(float f)
{
! float ret;
! __asm__ __volatile__("fsqrt" : "=t" (ret) : "0" (f));
! return (ret);
}
extern __inline__ long double
sqrtl(long double ld)
{
! long double ret;
! __asm__ __volatile__("fsqrt" : "=t" (ret) : "0" (ld));
! return (ret);
}
extern __inline__ int
isnanl(long double ld)
{
! int ret;
__asm__ __volatile__(
! "andl $0x00007fff,%1\n\t"
! "jz 1f\n\t" /* jump if __exp is all 0 */
! "xorl $0x00007fff,%1\n\t"
! "jz 2f\n\t" /* jump if __exp is all 1 */
! "testl $0x80000000,%2\n\t"
"jz 3f\n\t" /* jump if leading bit is 0 */
! "movl $0,%1\n\t"
"jmp 1f\n\t"
! "2:\n\t" /* note that %eax = 0 from before */
! "cmpl $0x80000000,%2\n\t" /* what is first half of __significand? */
"jnz 3f\n\t" /* jump if not equal to 0x80000000 */
! "testl $0xffffffff,%3\n\t" /* is second half of __significand 0? */
"jnz 3f\n\t" /* jump if not equal to 0 */
"jmp 1f\n\t"
"3:\n\t"
! "movl $1,%1\n\t"
"1:\n\t"
! : "=r" (ret)
! : "0" (_HIER_WORD(ld)), "r" (_HI_WORD(ld)), "r" (_LO_WORD(ld)));
return (ret);
}
#ifdef __cplusplus
--- 252,303 ----
}
extern __inline__ double
sqrt(double d)
{
! return (__inline_sqrt(d));
}
extern __inline__ float
sqrtf(float f)
{
! return (__inline_sqrtf(f));
}
extern __inline__ long double
sqrtl(long double ld)
{
! __asm__ __volatile__("fsqrt" : "+t" (ld) : : "cc");
! return (ld);
}
extern __inline__ int
isnanl(long double ld)
{
! int ret = _HIER_WORD(ld);
__asm__ __volatile__(
! "andl $0x00007fff,%0\n\t"
! "jz 1f\n\t" /* jump if exp is all 0 */
! "xorl $0x00007fff,%0\n\t"
! "jz 2f\n\t" /* jump if exp is all 1 */
! "testl $0x80000000,%1\n\t"
"jz 3f\n\t" /* jump if leading bit is 0 */
! "movl $0,%0\n\t"
"jmp 1f\n\t"
! "2:\n\t" /* note that %0 = 0 from before */
! "cmpl $0x80000000,%1\n\t" /* what is first half of significand? */
"jnz 3f\n\t" /* jump if not equal to 0x80000000 */
! "testl $0xffffffff,%2\n\t" /* is second half of significand 0? */
"jnz 3f\n\t" /* jump if not equal to 0 */
"jmp 1f\n\t"
"3:\n\t"
! "movl $1,%0\n\t"
"1:\n\t"
! : "+&r" (ret)
! : "r" (_HI_WORD(ld)), "r" (_LO_WORD(ld))
! : "cc");
return (ret);
}
#ifdef __cplusplus