- # special cases: 2*x is +-inf, nan or |x| < 0x1p-32
- # check (exponent|0x8000)+2 < 0xbfff+2-32
- movw 12(%esp), %ax
- movw %ax, %dx
- orw $0x8000, %dx
- addw $2, %dx
- cmpw $0xbfff-30, %dx
- jnb 3f
- cmpw $1, %dx
- jbe 1f
- # if |x|<0x1p-32 return 1+x
+ # interesting case: 0x1p-32 <= |x| < 16384
+ # check if (exponent|0x8000) is in [0xbfff-32, 0xbfff+13]
+ mov 12(%esp), %ax
+ or $0x8000, %ax
+ sub $0xbfdf, %ax
+ cmp $45, %ax
+ jbe 2f
+ test %ax, %ax