@@ -1497,6 +1497,29 @@ exp2 1023
exp2 -1074
exp2 16383
exp2 -16400
+exp2 -126.125
+exp2 -126.25
+exp2 -126.375
+exp2 -126.5
+exp2 -126.625
+exp2 -126.75
+exp2 -126.875
+exp2 -1022.125
+exp2 -1022.25
+exp2 -1022.375
+exp2 -1022.5
+exp2 -1022.625
+exp2 -1022.75
+exp2 -1022.875
+exp2 -0x3.fe4e8p+8
+exp2 -0x3.fe513p+8
+exp2 -16382.125
+exp2 -16382.25
+exp2 -16382.375
+exp2 -16382.5
+exp2 -16382.625
+exp2 -16382.75
+exp2 -16382.875
exp2 0x1p-10
exp2 -0x1p-10
exp2 0x1p-20
@@ -38,9 +38,17 @@ __ieee754_exp2l (long double x)
{
int intx = (int) x;
long double fractx = x - intx;
+ long double result;
if (fabsl (fractx) < LDBL_EPSILON / 4.0L)
- return __scalbnl (1.0L + fractx, intx);
- return __scalbnl (__ieee754_expl (M_LN2l * fractx), intx);
+ result = __scalbnl (1.0L + fractx, intx);
+ else
+ result = __scalbnl (__ieee754_expl (M_LN2l * fractx), intx);
+ if (result < LDBL_MIN)
+ {
+ long double force_underflow = result * result;
+ math_force_eval (force_underflow);
+ }
+ return result;
}
else
{
@@ -6,7 +6,24 @@
#include <machine/asm.h>
+ .section .rodata.cst8,"aM",@progbits,8
+
+ .p2align 3
+ .type dbl_min,@object
+dbl_min: .byte 0, 0, 0, 0, 0, 0, 0x10, 0
+ ASM_SIZE_DIRECTIVE(dbl_min)
+
+#ifdef PIC
+# define MO(op) op##@GOTOFF(%ecx)
+#else
+# define MO(op) op
+#endif
+
+ .text
ENTRY(__ieee754_exp2)
+#ifdef PIC
+ LOAD_PIC_REG (cx)
+#endif
fldl 4(%esp)
/* I added the following ugly construct because exp(+-Inf) resulted
in NaN. The ugliness results from the bright minds at Intel.
@@ -27,7 +44,22 @@ ENTRY(__ieee754_exp2)
faddp /* 2^(fract(x)) */
fscale /* e^x */
fstp %st(1)
- ret
+ fldl MO(dbl_min)
+ fld %st(1)
+ fucompp
+ fnstsw
+ sahf
+ jnc 3f
+ subl $8, %esp
+ cfi_adjust_cfa_offset (8)
+ fld %st(0)
+ fmul %st(0)
+ fstpl (%esp)
+ fstpl (%esp)
+ fldl (%esp)
+ addl $8, %esp
+ cfi_adjust_cfa_offset (-8)
+3: ret
1: testl $0x200, %eax /* Test sign. */
jz 2f /* If positive, jump. */
@@ -6,7 +6,24 @@
#include <machine/asm.h>
+ .section .rodata.cst4,"aM",@progbits,4
+
+ .p2align 2
+ .type flt_min,@object
+flt_min: .byte 0, 0, 0x80, 0
+ ASM_SIZE_DIRECTIVE(flt_min)
+
+#ifdef PIC
+# define MO(op) op##@GOTOFF(%ecx)
+#else
+# define MO(op) op
+#endif
+
+ .text
ENTRY(__ieee754_exp2f)
+#ifdef PIC
+ LOAD_PIC_REG (cx)
+#endif
flds 4(%esp)
/* I added the following ugly construct because exp(+-Inf) resulted
in NaN. The ugliness results from the bright minds at Intel.
@@ -27,7 +44,22 @@ ENTRY(__ieee754_exp2f)
faddp /* 2^(fract(x)) */
fscale /* e^x */
fstp %st(1)
- ret
+ flds MO(flt_min)
+ fld %st(1)
+ fucompp
+ fnstsw
+ sahf
+ jnc 3f
+ subl $4, %esp
+ cfi_adjust_cfa_offset (4)
+ fld %st(0)
+ fmul %st(0)
+ fstps (%esp)
+ fstps (%esp)
+ flds (%esp)
+ addl $4, %esp
+ cfi_adjust_cfa_offset (-4)
+3: ret
1: testl $0x200, %eax /* Test sign. */
jz 2f /* If positive, jump. */
@@ -6,7 +6,24 @@
#include <machine/asm.h>
+ .section .rodata.cst16,"aM",@progbits,16
+ .p2align 4
+ .type ldbl_min,@object
+ldbl_min: .byte 0, 0, 0, 0, 0, 0, 0, 0x80, 0x1, 0
+ .byte 0, 0, 0, 0, 0, 0
+ ASM_SIZE_DIRECTIVE(ldbl_min)
+
+#ifdef PIC
+# define MO(op) op##@GOTOFF(%ecx)
+#else
+# define MO(op) op
+#endif
+
+ .text
ENTRY(__ieee754_exp2l)
+#ifdef PIC
+ LOAD_PIC_REG (cx)
+#endif
fldt 4(%esp)
/* I added the following ugly construct because exp(+-Inf) resulted
in NaN. The ugliness results from the bright minds at Intel.
@@ -35,7 +52,17 @@ ENTRY(__ieee754_exp2l)
faddp /* 2^(fract(x)) */
fscale /* e^x */
fstp %st(1)
- ret
+ /* Ensure underflow for tiny result. */
+ fldt MO(ldbl_min)
+ fld %st(1)
+ fucompp
+ fnstsw
+ sahf
+ jnc 4f
+ fld %st
+ fmul %st
+ fstp %st
+4: ret
1: testl $0x200, %eax /* Test sign. */
jz 2f /* If positive, jump. */
@@ -120,7 +120,15 @@ __ieee754_exp2 (double x)
if (!unsafe)
return result;
else
- return result * scale_u.d;
+ {
+ result *= scale_u.d;
+ if (result < DBL_MIN)
+ {
+ double force_underflow = result * result;
+ math_force_eval (force_underflow);
+ }
+ return result;
+ }
}
else
/* Return x, if x is a NaN or Inf; or overflow, otherwise. */
@@ -109,7 +109,15 @@ __ieee754_exp2f (float x)
if (!unsafe)
return result;
else
- return result * scale_u.f;
+ {
+ result *= scale_u.f;
+ if (result < FLT_MIN)
+ {
+ float force_underflow = result * result;
+ math_force_eval (force_underflow);
+ }
+ return result;
+ }
}
/* Exceptional cases: */
else if (isless (x, himark))
@@ -7,6 +7,20 @@
#include <machine/asm.h>
+ .section .rodata.cst16,"aM",@progbits,16
+ .p2align 4
+ .type ldbl_min,@object
+ldbl_min: .byte 0, 0, 0, 0, 0, 0, 0, 0x80, 0x1, 0
+ .byte 0, 0, 0, 0, 0, 0
+ ASM_SIZE_DIRECTIVE(ldbl_min)
+
+#ifdef PIC
+# define MO(op) op##(%rip)
+#else
+# define MO(op) op
+#endif
+
+ .text
ENTRY(__ieee754_exp2l)
fldt 8(%rsp)
/* I added the following ugly construct because exp(+-Inf) resulted
@@ -36,7 +50,16 @@ ENTRY(__ieee754_exp2l)
faddp /* 2^(fract(x)) */
fscale /* e^x */
fstp %st(1)
- ret
+ /* Ensure underflow for tiny result. */
+ fldt MO(ldbl_min)
+ fld %st(1)
+ fucomip %st(1), %st
+ fstp %st
+ jnc 4f
+ fld %st
+ fmul %st
+ fstp %st
+4: ret
1: testl $0x200, %eax /* Test sign. */
jz 2f /* If positive, jump. */