94 lines
1.4 KiB
ArmAsm
94 lines
1.4 KiB
ArmAsm
.global expm1l
|
|
.type expm1l,@function
|
|
expm1l:
|
|
fldt 4(%esp)
|
|
fldl2e
|
|
fmulp
|
|
mov $0xc2820000,%eax
|
|
push %eax
|
|
flds (%esp)
|
|
pop %eax
|
|
fucomp %st(1)
|
|
fnstsw %ax
|
|
sahf
|
|
fld1
|
|
jb 1f
|
|
# x*log2e < -65, return -1 without underflow
|
|
fstp %st(1)
|
|
fchs
|
|
ret
|
|
1: fld %st(1)
|
|
fabs
|
|
fucom %st(1)
|
|
fnstsw %ax
|
|
fstp %st(0)
|
|
fstp %st(0)
|
|
sahf
|
|
ja 1f
|
|
f2xm1
|
|
ret
|
|
1: call 1f
|
|
fld1
|
|
fsubrp
|
|
ret
|
|
|
|
.global exp2l
|
|
.global __exp2l
|
|
.hidden __exp2l
|
|
.type exp2l,@function
|
|
exp2l:
|
|
__exp2l:
|
|
fldt 4(%esp)
|
|
1: sub $12,%esp
|
|
fld %st(0)
|
|
fstpt (%esp)
|
|
mov 8(%esp),%ax
|
|
and $0x7fff,%ax
|
|
cmp $0x3fff+13,%ax
|
|
jb 4f # |x| < 8192
|
|
cmp $0x3fff+15,%ax
|
|
jae 3f # |x| >= 32768
|
|
fsts (%esp)
|
|
cmpl $0xc67ff800,(%esp)
|
|
jb 2f # x > -16382
|
|
movl $0x5f000000,(%esp)
|
|
flds (%esp) # 0x1p63
|
|
fld %st(1)
|
|
fsub %st(1)
|
|
faddp
|
|
fucomp %st(1)
|
|
fnstsw
|
|
sahf
|
|
je 2f # x - 0x1p63 + 0x1p63 == x
|
|
movl $1,(%esp)
|
|
flds (%esp) # 0x1p-149
|
|
fdiv %st(1)
|
|
fstps (%esp) # raise underflow
|
|
2: fld1
|
|
fld %st(1)
|
|
frndint
|
|
fxch %st(2)
|
|
fsub %st(2) # st(0)=x-rint(x), st(1)=1, st(2)=rint(x)
|
|
f2xm1
|
|
faddp # 2^(x-rint(x))
|
|
1: fscale
|
|
fstp %st(1)
|
|
add $12,%esp
|
|
ret
|
|
3: xor %eax,%eax
|
|
4: cmp $0x3fff-64,%ax
|
|
fld1
|
|
jb 1b # |x| < 0x1p-64
|
|
fstpt (%esp)
|
|
fistl 8(%esp)
|
|
fildl 8(%esp)
|
|
fsubrp %st(1)
|
|
addl $0x3fff,8(%esp)
|
|
f2xm1
|
|
fld1
|
|
faddp # 2^(x-rint(x))
|
|
fldt (%esp) # 2^rint(x)
|
|
fmulp
|
|
add $12,%esp
|
|
ret
|