| /** | |
| * This file has no copyright assigned and is placed in the Public Domain. | |
| * This file is part of the w64 mingw-runtime package. | |
| * No warranty is given; refer to the file DISCLAIMER within this package. | |
| */ | |
| #ifdef _WIN64 | |
| .file "exp2.S" | |
| .text | |
| .align 4 | |
| .globl _exp2 | |
| .def _exp2; .scl 2; .type 32; .endef | |
| _exp2: | |
| movsd %xmm0,-16(%rsp) | |
| fldl -16(%rsp) | |
| fxam /* Is NaN or +-Inf? */ | |
| fstsw %ax | |
| movb $0x45, %dh | |
| andb %ah, %dh | |
| cmpb $0x05, %dh | |
| je 1f /* Is +-Inf, jump. */ | |
| fld %st | |
| frndint /* int(x) */ | |
| fsubr %st,%st(1) /* fract(x) */ | |
| fxch | |
| f2xm1 /* 2^(fract(x)) - 1 */ | |
| fld1 | |
| faddp /* 2^(fract(x)) */ | |
| fscale /* e^x */ | |
| fstp %st(1) | |
| fstpl -16(%rsp) | |
| movsd -16(%rsp),%xmm0 | |
| ret | |
| 1: testl $0x200, %eax /* Test sign. */ | |
| jz 2f /* If positive, jump. */ | |
| fstp %st | |
| fldz /* Set result to 0. */ | |
| 2: fstpl -16(%rsp) | |
| movsd -16(%rsp),%xmm0 | |
| ret | |
| #else | |
| .file "exp2.S" | |
| .text | |
| .align 4 | |
| .globl _exp2 | |
| .def _exp2; .scl 2; .type 32; .endef | |
| _exp2: | |
| fldl 4(%esp) | |
| /* I added the following ugly construct because exp(+-Inf) resulted | |
| in NaN. The ugliness results from the bright minds at Intel. | |
| For the i686 the code can be written better. | |
| -- drepper@cygnus.com. */ | |
| fxam /* Is NaN or +-Inf? */ | |
| fstsw %ax | |
| movb $0x45, %dh | |
| andb %ah, %dh | |
| cmpb $0x05, %dh | |
| je 1f /* Is +-Inf, jump. */ | |
| fld %st | |
| frndint /* int(x) */ | |
| fsubr %st,%st(1) /* fract(x) */ | |
| fxch | |
| f2xm1 /* 2^(fract(x)) - 1 */ | |
| fld1 | |
| faddp /* 2^(fract(x)) */ | |
| fscale /* e^x */ | |
| fstp %st(1) | |
| ret | |
| 1: testl $0x200, %eax /* Test sign. */ | |
| jz 2f /* If positive, jump. */ | |
| fstp %st | |
| fldz /* Set result to 0. */ | |
| 2: ret | |
| #endif |