exp(inf), exp(-inf), exp(nan) used to raise wrong flags
exp2:
fldl 4(%esp)
1: pushl $0x467ff000
exp2:
fldl 4(%esp)
1: pushl $0x467ff000
xorl %eax,%eax
pushl $0x80000000
push %eax
fld %st(1)
fabs
xorl %eax,%eax
pushl $0x80000000
push %eax
fld %st(1)
fabs
+ ja 3f # |x| > 16380
+ jp 2f # x is nan (avoid invalid except in fistp)
fld %st(0)
fistpl 8(%esp)
fildl 8(%esp)
fld %st(0)
fistpl 8(%esp)
fildl 8(%esp)
add %eax,8(%esp)
f2xm1
fld1
add %eax,8(%esp)
f2xm1
fld1
+ faddp # 2^(x-rint(x))
+ fldt (%esp) # 2^rint(x)
- mov 9(%esp),%ah
- and $0x7f,%ah
- cmp $0x7f,%ah
- jne 1f
- decb 9(%esp)
- fstp %st(0)
- fldt (%esp)
-1: fld %st(0)
+ fld1
+ mov 8(%esp),%ax
+ and $0x7fff,%ax
+ cmp $0x7fff,%ax
+ je 1f # x = +-inf
+ fld %st(1)
- fxch %st(1)
- fsub %st(1)
+ fxch %st(2)
+ fsub %st(2) # st(0)=x-rint(x), st(1)=1, st(2)=rint(x)
+ faddp # 2^(x-rint(x))
+1: fscale
fstp %st(1)
add $12,%esp
ret
fstp %st(1)
add $12,%esp
ret
mov %eax,-20(%rsp)
xor %eax,%eax
mov %eax,-24(%rsp)
mov %eax,-20(%rsp)
xor %eax,%eax
mov %eax,-24(%rsp)
fld %st(1)
fabs
fucom %st(1)
fld %st(1)
fabs
fucom %st(1)
fstp %st(0)
fstp %st(0)
sahf
fstp %st(0)
fstp %st(0)
sahf
+ ja 3f # |x| > 16380
+ jp 2f # x is nan (avoid invalid except in fistp)
fld %st(0)
fistpl -16(%rsp)
fildl -16(%rsp)
fld %st(0)
fistpl -16(%rsp)
fildl -16(%rsp)
add %eax,-16(%rsp)
f2xm1
fld1
add %eax,-16(%rsp)
f2xm1
fld1
+ faddp # 2^(x-rint(x))
+ fldt -24(%rsp) # 2^rint(x)
- mov -15(%rsp),%ah
- and $0x7f,%ah
- cmp $0x7f,%ah
- jne 1f
- decb -15(%rsp)
- fstp %st(0)
- fldt -24(%rsp)
-1: fld %st(0)
+ fld1
+ mov -15(%rsp),%ax
+ and $0x7fff,%ax
+ cmp $0x7fff,%ax
+ je 1f # x = +-inf
+ fld %st(1)
- fxch %st(1)
- fsub %st(1)
+ fxch %st(2)
+ fsub %st(2) # st(0)=x-rint(x), st(1)=1, st(2)=rint(x)
+ faddp # 2^(x-rint(x))
+1: fscale