2 * Written by J.T. Conklin <jtc@netbsd.org>.
3 * Adapted for exp2 by Ulrich Drepper <drepper@cygnus.com>.
4 * Adapted for x86-64 by Andreas Jaeger <aj@suse.de>.
8 #include <machine/asm.h>
10 ENTRY(__ieee754_exp2l)
12 /* I added the following ugly construct because exp(+-Inf) resulted
13 in NaN. The ugliness results from the bright minds at Intel.
14 For the i686 the code can be written better.
15 -- drepper@cygnus.com. */
16 fxam /* Is NaN or +-Inf? */
21 je 1f /* Is +-Inf, jump. */
24 fsubr %st,%st(1) /* fract(x) */
26 f2xm1 /* 2^(fract(x)) - 1 */
28 faddp /* 2^(fract(x)) */
33 1: testl $0x200, %eax /* Test sign. */
34 jz 2f /* If positive, jump. */
36 fldz /* Set result to 0. */