From 324ff738ec5f8547a8e866290c39ad2262960e98 Mon Sep 17 00:00:00 2001
From: Ben Kimock <kimockb@gmail.com>
Date: Fri, 19 Apr 2024 18:24:06 -0400
Subject: [PATCH] Remove force_eval uses for raising fp exceptions

---
 src/math/asinh.rs      | 4 ----
 src/math/asinhf.rs     | 4 ----
 src/math/atan.rs       | 5 -----
 src/math/atanf.rs      | 4 ----
 src/math/atanh.rs      | 7 +------
 src/math/atanhf.rs     | 7 +------
 src/math/cosf.rs       | 5 -----
 src/math/cosh.rs       | 2 --
 src/math/coshf.rs      | 3 ---
 src/math/exp.rs        | 6 ------
 src/math/exp2.rs       | 5 -----
 src/math/exp2f.rs      | 3 ---
 src/math/expf.rs       | 5 -----
 src/math/expm1.rs      | 3 ---
 src/math/expm1f.rs     | 3 ---
 src/math/floorf.rs     | 2 --
 src/math/fmaf.rs       | 3 +--
 src/math/ilogb.rs      | 2 --
 src/math/ilogbf.rs     | 2 --
 src/math/log1p.rs      | 4 ----
 src/math/log1pf.rs     | 4 ----
 src/math/mod.rs        | 2 +-
 src/math/nextafter.rs  | 9 ---------
 src/math/nextafterf.rs | 9 ---------
 src/math/sin.rs        | 8 --------
 src/math/sincos.rs     | 7 -------
 src/math/sincosf.rs    | 8 --------
 src/math/sinf.rs       | 8 --------
 src/math/tan.rs        | 8 --------
 src/math/tanf.rs       | 8 --------
 src/math/tanh.rs       | 2 --
 src/math/tanhf.rs      | 2 --
 src/math/tgamma.rs     | 2 --
 src/math/trunc.rs      | 3 ---
 src/math/truncf.rs     | 2 --
 35 files changed, 4 insertions(+), 157 deletions(-)

diff --git a/src/math/asinh.rs b/src/math/asinh.rs
index 0abd80c2..d7f9dd91 100644
--- a/src/math/asinh.rs
+++ b/src/math/asinh.rs
@@ -26,10 +26,6 @@ pub fn asinh(mut x: f64) -> f64 {
     } else if e >= 0x3ff - 26 {
         /* |x| >= 0x1p-26, up to 1.6ulp error in [0.125,0.5] */
         x = log1p(x + x * x / (sqrt(x * x + 1.0) + 1.0));
-    } else {
-        /* |x| < 0x1p-26, raise inexact if x != 0 */
-        let x1p120 = f64::from_bits(0x4770000000000000);
-        force_eval!(x + x1p120);
     }
 
     if sign {
diff --git a/src/math/asinhf.rs b/src/math/asinhf.rs
index 09c77823..8dd6692d 100644
--- a/src/math/asinhf.rs
+++ b/src/math/asinhf.rs
@@ -25,10 +25,6 @@ pub fn asinhf(mut x: f32) -> f32 {
     } else if i >= 0x3f800000 - (12 << 23) {
         /* |x| >= 0x1p-12, up to 1.6ulp error in [0.125,0.5] */
         x = log1pf(x + x * x / (sqrtf(x * x + 1.0) + 1.0));
-    } else {
-        /* |x| < 0x1p-12, raise inexact if x!=0 */
-        let x1p120 = f32::from_bits(0x7b800000);
-        force_eval!(x + x1p120);
     }
 
     if sign {
diff --git a/src/math/atan.rs b/src/math/atan.rs
index 4259dc71..bb4413bc 100644
--- a/src/math/atan.rs
+++ b/src/math/atan.rs
@@ -83,11 +83,6 @@ pub fn atan(x: f64) -> f64 {
         /* |x| < 0.4375 */
         if ix < 0x3e40_0000 {
             /* |x| < 2^-27 */
-            if ix < 0x0010_0000 {
-                /* raise underflow for subnormal x */
-                force_eval!(x as f32);
-            }
-
             return x;
         }
 
diff --git a/src/math/atanf.rs b/src/math/atanf.rs
index d042b3bc..321d39dd 100644
--- a/src/math/atanf.rs
+++ b/src/math/atanf.rs
@@ -63,10 +63,6 @@ pub fn atanf(mut x: f32) -> f32 {
         /* |x| < 0.4375 */
         if ix < 0x39800000 {
             /* |x| < 2**-12 */
-            if ix < 0x00800000 {
-                /* raise underflow for subnormal x */
-                force_eval!(x * x);
-            }
             return x;
         }
         -1
diff --git a/src/math/atanh.rs b/src/math/atanh.rs
index b984c4ac..d3a85415 100644
--- a/src/math/atanh.rs
+++ b/src/math/atanh.rs
@@ -15,12 +15,7 @@ pub fn atanh(x: f64) -> f64 {
     let mut y = f64::from_bits(u & 0x7fff_ffff_ffff_ffff);
 
     if e < 0x3ff - 1 {
-        if e < 0x3ff - 32 {
-            /* handle underflow */
-            if e == 0 {
-                force_eval!(y as f32);
-            }
-        } else {
+        if e >= 0x3ff - 32 {
             /* |x| < 0.5, up to 1.7ulp error */
             y = 0.5 * log1p(2.0 * y + 2.0 * y * y / (1.0 - y));
         }
diff --git a/src/math/atanhf.rs b/src/math/atanhf.rs
index a1aa314a..46d18b94 100644
--- a/src/math/atanhf.rs
+++ b/src/math/atanhf.rs
@@ -15,12 +15,7 @@ pub fn atanhf(mut x: f32) -> f32 {
     x = f32::from_bits(u);
 
     if u < 0x3f800000 - (1 << 23) {
-        if u < 0x3f800000 - (32 << 23) {
-            /* handle underflow */
-            if u < (1 << 23) {
-                force_eval!((x * x) as f32);
-            }
-        } else {
+        if u >= 0x3f800000 - (32 << 23) {
             /* |x| < 0.5, up to 1.7ulp error */
             x = 0.5 * log1pf(2.0 * x + 2.0 * x * x / (1.0 - x));
         }
diff --git a/src/math/cosf.rs b/src/math/cosf.rs
index 424fa42e..778f8749 100644
--- a/src/math/cosf.rs
+++ b/src/math/cosf.rs
@@ -28,8 +28,6 @@ const C4_PIO2: f64 = 4. * FRAC_PI_2; /* 0x401921FB, 0x54442D18 */
 pub fn cosf(x: f32) -> f32 {
     let x64 = x as f64;
 
-    let x1p120 = f32::from_bits(0x7b800000); // 0x1p120f === 2 ^ 120
-
     let mut ix = x.to_bits();
     let sign = (ix >> 31) != 0;
     ix &= 0x7fffffff;
@@ -37,9 +35,6 @@ pub fn cosf(x: f32) -> f32 {
     if ix <= 0x3f490fda {
         /* |x| ~<= pi/4 */
         if ix < 0x39800000 {
-            /* |x| < 2**-12 */
-            /* raise inexact if x != 0 */
-            force_eval!(x + x1p120);
             return 1.;
         }
         return k_cosf(x64);
diff --git a/src/math/cosh.rs b/src/math/cosh.rs
index 2fb568ab..2b0c2a28 100644
--- a/src/math/cosh.rs
+++ b/src/math/cosh.rs
@@ -18,8 +18,6 @@ pub fn cosh(mut x: f64) -> f64 {
     /* |x| < log(2) */
     if w < 0x3fe62e42 {
         if w < 0x3ff00000 - (26 << 20) {
-            let x1p120 = f64::from_bits(0x4770000000000000);
-            force_eval!(x + x1p120);
             return 1.;
         }
         let t = expm1(x); // exponential minus 1
diff --git a/src/math/coshf.rs b/src/math/coshf.rs
index e7b68458..c45fcfe7 100644
--- a/src/math/coshf.rs
+++ b/src/math/coshf.rs
@@ -9,8 +9,6 @@ use super::k_expo2f;
 /// Angles are specified in radians.
 #[cfg_attr(all(test, assert_no_panic), no_panic::no_panic)]
 pub fn coshf(mut x: f32) -> f32 {
-    let x1p120 = f32::from_bits(0x7b800000); // 0x1p120f === 2 ^ 120
-
     /* |x| */
     let mut ix = x.to_bits();
     ix &= 0x7fffffff;
@@ -20,7 +18,6 @@ pub fn coshf(mut x: f32) -> f32 {
     /* |x| < log(2) */
     if w < 0x3f317217 {
         if w < (0x3f800000 - (12 << 23)) {
-            force_eval!(x + x1p120);
             return 1.;
         }
         let t = expm1f(x);
diff --git a/src/math/exp.rs b/src/math/exp.rs
index d4994277..f45f257c 100644
--- a/src/math/exp.rs
+++ b/src/math/exp.rs
@@ -84,7 +84,6 @@ const P5: f64 = 4.13813679705723846039e-08; /* 0x3E663769, 0x72BEA4D0 */
 #[cfg_attr(all(test, assert_no_panic), no_panic::no_panic)]
 pub fn exp(mut x: f64) -> f64 {
     let x1p1023 = f64::from_bits(0x7fe0000000000000); // 0x1p1023 === 2 ^ 1023
-    let x1p_149 = f64::from_bits(0x36a0000000000000); // 0x1p-149 === 2 ^ -149
 
     let hi: f64;
     let lo: f64;
@@ -106,13 +105,10 @@ pub fn exp(mut x: f64) -> f64 {
             return x;
         }
         if x > 709.782712893383973096 {
-            /* overflow if x!=inf */
             x *= x1p1023;
             return x;
         }
         if x < -708.39641853226410622 {
-            /* underflow if x!=-inf */
-            force_eval!((-x1p_149 / x) as f32);
             if x < -745.13321910194110842 {
                 return 0.;
             }
@@ -137,8 +133,6 @@ pub fn exp(mut x: f64) -> f64 {
         hi = x;
         lo = 0.;
     } else {
-        /* inexact if x!=0 */
-        force_eval!(x1p1023 + x);
         return 1. + x;
     }
 
diff --git a/src/math/exp2.rs b/src/math/exp2.rs
index e0e385df..83198562 100644
--- a/src/math/exp2.rs
+++ b/src/math/exp2.rs
@@ -336,7 +336,6 @@ pub fn exp2(mut x: f64) -> f64 {
     // union {double f; uint64_t i;} u = {x};
     // union {uint32_t u; int32_t i;} k;
     let x1p1023 = f64::from_bits(0x7fe0000000000000);
-    let x1p52 = f64::from_bits(0x4330000000000000);
     let _0x1p_149 = f64::from_bits(0xb6a0000000000000);
 
     /* Filter out exceptional cases. */
@@ -356,10 +355,6 @@ pub fn exp2(mut x: f64) -> f64 {
         }
         if ui >> 63 != 0 {
             /* x <= -1022 */
-            /* underflow */
-            if x <= -1075.0 || x - x1p52 + x1p52 != x {
-                force_eval!((_0x1p_149 / x) as f32);
-            }
             if x <= -1075.0 {
                 return 0.0;
             }
diff --git a/src/math/exp2f.rs b/src/math/exp2f.rs
index f4867b80..7ef2b909 100644
--- a/src/math/exp2f.rs
+++ b/src/math/exp2f.rs
@@ -102,9 +102,6 @@ pub fn exp2f(mut x: f32) -> f32 {
         }
         if ui >= 0x80000000 {
             /* x < -126 */
-            if ui >= 0xc3160000 || (ui & 0x0000ffff != 0) {
-                force_eval!(f32::from_bits(0x80000001) / x);
-            }
             if ui >= 0xc3160000 {
                 /* x <= -150 */
                 return 0.0;
diff --git a/src/math/expf.rs b/src/math/expf.rs
index a53aa90a..9a5fc8b8 100644
--- a/src/math/expf.rs
+++ b/src/math/expf.rs
@@ -33,7 +33,6 @@ const P2: f32 = -2.7667332906e-3; /* -0xb55215.0p-32 */
 #[cfg_attr(all(test, assert_no_panic), no_panic::no_panic)]
 pub fn expf(mut x: f32) -> f32 {
     let x1p127 = f32::from_bits(0x7f000000); // 0x1p127f === 2 ^ 127
-    let x1p_126 = f32::from_bits(0x800000); // 0x1p-126f === 2 ^ -126  /*original 0x1p-149f    ??????????? */
     let mut hx = x.to_bits();
     let sign = (hx >> 31) as i32; /* sign bit of x */
     let signb: bool = sign != 0;
@@ -53,8 +52,6 @@ pub fn expf(mut x: f32) -> f32 {
             return x;
         }
         if signb {
-            /* underflow */
-            force_eval!(-x1p_126 / x);
             if hx >= 0x42cff1b5 {
                 /* x <= -103.972084f */
                 return 0.;
@@ -84,8 +81,6 @@ pub fn expf(mut x: f32) -> f32 {
         hi = x;
         lo = 0.;
     } else {
-        /* raise inexact */
-        force_eval!(x1p127 + x);
         return 1. + x;
     }
 
diff --git a/src/math/expm1.rs b/src/math/expm1.rs
index 42608509..09f049f1 100644
--- a/src/math/expm1.rs
+++ b/src/math/expm1.rs
@@ -82,9 +82,6 @@ pub fn expm1(mut x: f64) -> f64 {
         c = (hi - x) - lo;
     } else if hx < 0x3c900000 {
         /* |x| < 2**-54, return x */
-        if hx < 0x00100000 {
-            force_eval!(x);
-        }
         return x;
     } else {
         c = 0.0;
diff --git a/src/math/expm1f.rs b/src/math/expm1f.rs
index 3fc2a247..0d6d18af 100644
--- a/src/math/expm1f.rs
+++ b/src/math/expm1f.rs
@@ -84,9 +84,6 @@ pub fn expm1f(mut x: f32) -> f32 {
         c = (hi - x) - lo;
     } else if hx < 0x33000000 {
         /* when |x|<2**-25, return x */
-        if hx < 0x00800000 {
-            force_eval!(x * x);
-        }
         return x;
     } else {
         k = 0;
diff --git a/src/math/floorf.rs b/src/math/floorf.rs
index dfdab91a..421d75dc 100644
--- a/src/math/floorf.rs
+++ b/src/math/floorf.rs
@@ -24,13 +24,11 @@ pub fn floorf(x: f32) -> f32 {
         if (ui & m) == 0 {
             return x;
         }
-        force_eval!(x + f32::from_bits(0x7b800000));
         if ui >> 31 != 0 {
             ui += m;
         }
         ui &= !m;
     } else {
-        force_eval!(x + f32::from_bits(0x7b800000));
         if ui >> 31 == 0 {
             ui = 0;
         } else if ui << 1 != 0 {
diff --git a/src/math/fmaf.rs b/src/math/fmaf.rs
index 2848f2ae..0ef8cf75 100644
--- a/src/math/fmaf.rs
+++ b/src/math/fmaf.rs
@@ -26,7 +26,6 @@
  */
 
 use core::f32;
-use core::ptr::read_volatile;
 
 use super::fenv::{
     feclearexcept, fegetround, feraiseexcept, fetestexcept, FE_INEXACT, FE_TONEAREST, FE_UNDERFLOW,
@@ -74,7 +73,7 @@ pub fn fmaf(x: f32, y: f32, mut z: f32) -> f32 {
         if e < 0x3ff - 126 && e >= 0x3ff - 149 && fetestexcept(FE_INEXACT) != 0 {
             feclearexcept(FE_INEXACT);
             // prevent `xy + vz` from being CSE'd with `xy + z` above
-            let vz: f32 = unsafe { read_volatile(&z) };
+            let vz = force_eval!(z);
             result = xy + vz as f64;
             if fetestexcept(FE_INEXACT) != 0 {
                 feraiseexcept(FE_UNDERFLOW);
diff --git a/src/math/ilogb.rs b/src/math/ilogb.rs
index 7d74dcfb..23bc39ab 100644
--- a/src/math/ilogb.rs
+++ b/src/math/ilogb.rs
@@ -9,7 +9,6 @@ pub fn ilogb(x: f64) -> i32 {
     if e == 0 {
         i <<= 12;
         if i == 0 {
-            force_eval!(0.0 / 0.0);
             return FP_ILOGB0;
         }
         /* subnormal x */
@@ -20,7 +19,6 @@ pub fn ilogb(x: f64) -> i32 {
         }
         e
     } else if e == 0x7ff {
-        force_eval!(0.0 / 0.0);
         if (i << 12) != 0 {
             FP_ILOGBNAN
         } else {
diff --git a/src/math/ilogbf.rs b/src/math/ilogbf.rs
index 0fa58748..244256f3 100644
--- a/src/math/ilogbf.rs
+++ b/src/math/ilogbf.rs
@@ -9,7 +9,6 @@ pub fn ilogbf(x: f32) -> i32 {
     if e == 0 {
         i <<= 9;
         if i == 0 {
-            force_eval!(0.0 / 0.0);
             return FP_ILOGB0;
         }
         /* subnormal x */
@@ -20,7 +19,6 @@ pub fn ilogbf(x: f32) -> i32 {
         }
         e
     } else if e == 0xff {
-        force_eval!(0.0 / 0.0);
         if (i << 9) != 0 {
             FP_ILOGBNAN
         } else {
diff --git a/src/math/log1p.rs b/src/math/log1p.rs
index 4fd1c73e..9051933e 100644
--- a/src/math/log1p.rs
+++ b/src/math/log1p.rs
@@ -95,10 +95,6 @@ pub fn log1p(x: f64) -> f64 {
         }
         if hx << 1 < 0x3ca00000 << 1 {
             /* |x| < 2**-53 */
-            /* underflow if subnormal */
-            if (hx & 0x7ff00000) == 0 {
-                force_eval!(x as f32);
-            }
             return x;
         }
         if hx <= 0xbfd2bec4 {
diff --git a/src/math/log1pf.rs b/src/math/log1pf.rs
index 500e8eea..85da5977 100644
--- a/src/math/log1pf.rs
+++ b/src/math/log1pf.rs
@@ -50,10 +50,6 @@ pub fn log1pf(x: f32) -> f32 {
         }
         if ix << 1 < 0x33800000 << 1 {
             /* |x| < 2**-24 */
-            /* underflow if subnormal */
-            if (ix & 0x7f800000) == 0 {
-                force_eval!(x * x);
-            }
             return x;
         }
         if ix <= 0xbe95f619 {
diff --git a/src/math/mod.rs b/src/math/mod.rs
index 05ebb708..c82c3ca0 100644
--- a/src/math/mod.rs
+++ b/src/math/mod.rs
@@ -1,6 +1,6 @@
 macro_rules! force_eval {
     ($e:expr) => {
-        unsafe { ::core::ptr::read_volatile(&$e) }
+        core::hint::black_box($e)
     };
 }
 
diff --git a/src/math/nextafter.rs b/src/math/nextafter.rs
index 05762619..5351ea1e 100644
--- a/src/math/nextafter.rs
+++ b/src/math/nextafter.rs
@@ -23,15 +23,6 @@ pub fn nextafter(x: f64, y: f64) -> f64 {
         ux_i += 1;
     }
 
-    let e = ux_i >> 52 & 0x7ff;
-    // raise overflow if ux.f is infinite and x is finite
-    if e == 0x7ff {
-        force_eval!(x + x);
-    }
     let ux_f = f64::from_bits(ux_i);
-    // raise underflow if ux.f is subnormal or zero
-    if e == 0 {
-        force_eval!(x * x + ux_f * ux_f);
-    }
     ux_f
 }
diff --git a/src/math/nextafterf.rs b/src/math/nextafterf.rs
index 8ba38335..2853390b 100644
--- a/src/math/nextafterf.rs
+++ b/src/math/nextafterf.rs
@@ -23,15 +23,6 @@ pub fn nextafterf(x: f32, y: f32) -> f32 {
         ux_i += 1;
     }
 
-    let e = ux_i & 0x7f80_0000_u32;
-    // raise overflow if ux_f is infinite and x is finite
-    if e == 0x7f80_0000_u32 {
-        force_eval!(x + x);
-    }
     let ux_f = f32::from_bits(ux_i);
-    // raise underflow if ux_f is subnormal or zero
-    if e == 0 {
-        force_eval!(x * x + ux_f * ux_f);
-    }
     ux_f
 }
diff --git a/src/math/sin.rs b/src/math/sin.rs
index a53843dc..f2efb900 100644
--- a/src/math/sin.rs
+++ b/src/math/sin.rs
@@ -42,8 +42,6 @@ use super::{k_cos, k_sin, rem_pio2};
 //      TRIG(x) returns trig(x) nearly rounded
 #[cfg_attr(all(test, assert_no_panic), no_panic::no_panic)]
 pub fn sin(x: f64) -> f64 {
-    let x1p120 = f64::from_bits(0x4770000000000000); // 0x1p120f === 2 ^ 120
-
     /* High word of x. */
     let ix = (f64::to_bits(x) >> 32) as u32 & 0x7fffffff;
 
@@ -51,12 +49,6 @@ pub fn sin(x: f64) -> f64 {
     if ix <= 0x3fe921fb {
         if ix < 0x3e500000 {
             /* |x| < 2**-26 */
-            /* raise inexact if x != 0 and underflow if subnormal*/
-            if ix < 0x00100000 {
-                force_eval!(x / x1p120);
-            } else {
-                force_eval!(x + x1p120);
-            }
             return x;
         }
         return k_sin(x, 0.0, 0);
diff --git a/src/math/sincos.rs b/src/math/sincos.rs
index ff5d87a1..72df4f09 100644
--- a/src/math/sincos.rs
+++ b/src/math/sincos.rs
@@ -25,13 +25,6 @@ pub fn sincos(x: f64) -> (f64, f64) {
     if ix <= 0x3fe921fb {
         /* if |x| < 2**-27 * sqrt(2) */
         if ix < 0x3e46a09e {
-            /* raise inexact if x!=0 and underflow if subnormal */
-            let x1p120 = f64::from_bits(0x4770000000000000); // 0x1p120 == 2^120
-            if ix < 0x00100000 {
-                force_eval!(x / x1p120);
-            } else {
-                force_eval!(x + x1p120);
-            }
             return (x, 1.0);
         }
         return (k_sin(x, 0.0, 0), k_cos(x, 0.0));
diff --git a/src/math/sincosf.rs b/src/math/sincosf.rs
index 9a4c3610..a2e5089c 100644
--- a/src/math/sincosf.rs
+++ b/src/math/sincosf.rs
@@ -38,14 +38,6 @@ pub fn sincosf(x: f32) -> (f32, f32) {
     if ix <= 0x3f490fda {
         /* |x| < 2**-12 */
         if ix < 0x39800000 {
-            /* raise inexact if x!=0 and underflow if subnormal */
-
-            let x1p120 = f32::from_bits(0x7b800000); // 0x1p120 == 2^120
-            if ix < 0x00100000 {
-                force_eval!(x / x1p120);
-            } else {
-                force_eval!(x + x1p120);
-            }
             return (x, 1.0);
         }
         return (k_sinf(x as f64), k_cosf(x as f64));
diff --git a/src/math/sinf.rs b/src/math/sinf.rs
index 6e20be2a..7a14c523 100644
--- a/src/math/sinf.rs
+++ b/src/math/sinf.rs
@@ -28,8 +28,6 @@ const S4_PIO2: f64 = 4. * FRAC_PI_2; /* 0x401921FB, 0x54442D18 */
 pub fn sinf(x: f32) -> f32 {
     let x64 = x as f64;
 
-    let x1p120 = f32::from_bits(0x7b800000); // 0x1p120f === 2 ^ 120
-
     let mut ix = x.to_bits();
     let sign = (ix >> 31) != 0;
     ix &= 0x7fffffff;
@@ -38,12 +36,6 @@ pub fn sinf(x: f32) -> f32 {
         /* |x| ~<= pi/4 */
         if ix < 0x39800000 {
             /* |x| < 2**-12 */
-            /* raise inexact if x!=0 and underflow if subnormal */
-            force_eval!(if ix < 0x00800000 {
-                x / x1p120
-            } else {
-                x + x1p120
-            });
             return x;
         }
         return k_sinf(x64);
diff --git a/src/math/tan.rs b/src/math/tan.rs
index 5a72f680..d9b52ce4 100644
--- a/src/math/tan.rs
+++ b/src/math/tan.rs
@@ -41,19 +41,11 @@ use super::{k_tan, rem_pio2};
 //      TRIG(x) returns trig(x) nearly rounded
 #[cfg_attr(all(test, assert_no_panic), no_panic::no_panic)]
 pub fn tan(x: f64) -> f64 {
-    let x1p120 = f32::from_bits(0x7b800000); // 0x1p120f === 2 ^ 120
-
     let ix = (f64::to_bits(x) >> 32) as u32 & 0x7fffffff;
     /* |x| ~< pi/4 */
     if ix <= 0x3fe921fb {
         if ix < 0x3e400000 {
             /* |x| < 2**-27 */
-            /* raise inexact if x!=0 and underflow if subnormal */
-            force_eval!(if ix < 0x00100000 {
-                x / x1p120 as f64
-            } else {
-                x + x1p120 as f64
-            });
             return x;
         }
         return k_tan(x, 0.0, 0);
diff --git a/src/math/tanf.rs b/src/math/tanf.rs
index 10de59c3..426bea4d 100644
--- a/src/math/tanf.rs
+++ b/src/math/tanf.rs
@@ -28,8 +28,6 @@ const T4_PIO2: f64 = 4. * FRAC_PI_2; /* 0x401921FB, 0x54442D18 */
 pub fn tanf(x: f32) -> f32 {
     let x64 = x as f64;
 
-    let x1p120 = f32::from_bits(0x7b800000); // 0x1p120f === 2 ^ 120
-
     let mut ix = x.to_bits();
     let sign = (ix >> 31) != 0;
     ix &= 0x7fffffff;
@@ -38,12 +36,6 @@ pub fn tanf(x: f32) -> f32 {
         /* |x| ~<= pi/4 */
         if ix < 0x39800000 {
             /* |x| < 2**-12 */
-            /* raise inexact if x!=0 and underflow if subnormal */
-            force_eval!(if ix < 0x00800000 {
-                x / x1p120
-            } else {
-                x + x1p120
-            });
             return x;
         }
         return k_tanf(x64, false);
diff --git a/src/math/tanh.rs b/src/math/tanh.rs
index 980c6855..7ca52a83 100644
--- a/src/math/tanh.rs
+++ b/src/math/tanh.rs
@@ -40,8 +40,6 @@ pub fn tanh(mut x: f64) -> f64 {
         t = -t / (t + 2.0);
     } else {
         /* |x| is subnormal */
-        /* note: the branch above would not raise underflow in [0x1p-1023,0x1p-1022) */
-        force_eval!(x as f32);
         t = x;
     }
 
diff --git a/src/math/tanhf.rs b/src/math/tanhf.rs
index fc94e3dd..b6df5452 100644
--- a/src/math/tanhf.rs
+++ b/src/math/tanhf.rs
@@ -27,8 +27,6 @@ pub fn tanhf(mut x: f32) -> f32 {
         let t = expm1f(-2. * x);
         -t / (t + 2.)
     } else {
-        /* |x| is subnormal */
-        force_eval!(x * x);
         x
     };
     if sign {
diff --git a/src/math/tgamma.rs b/src/math/tgamma.rs
index 3f38c0b1..f0354f28 100644
--- a/src/math/tgamma.rs
+++ b/src/math/tgamma.rs
@@ -166,8 +166,6 @@ pub fn tgamma(mut x: f64) -> f64 {
     if ix >= 0x40670000 {
         /* |x| >= 184 */
         if sign {
-            let x1p_126 = f64::from_bits(0x3810000000000000); // 0x1p-126 == 2^-126
-            force_eval!((x1p_126 / x) as f32);
             if floor(x) * 0.5 == floor(x * 0.5) {
                 return 0.0;
             } else {
diff --git a/src/math/trunc.rs b/src/math/trunc.rs
index f7892a2c..7e8a1e2c 100644
--- a/src/math/trunc.rs
+++ b/src/math/trunc.rs
@@ -10,8 +10,6 @@ pub fn trunc(x: f64) -> f64 {
             return unsafe { ::core::intrinsics::truncf64(x) }
         }
     }
-    let x1p120 = f64::from_bits(0x4770000000000000); // 0x1p120f === 2 ^ 120
-
     let mut i: u64 = x.to_bits();
     let mut e: i64 = (i >> 52 & 0x7ff) as i64 - 0x3ff + 12;
     let m: u64;
@@ -26,7 +24,6 @@ pub fn trunc(x: f64) -> f64 {
     if (i & m) == 0 {
         return x;
     }
-    force_eval!(x + x1p120);
     i &= !m;
     f64::from_bits(i)
 }
diff --git a/src/math/truncf.rs b/src/math/truncf.rs
index 20d5b73b..d72b1a71 100644
--- a/src/math/truncf.rs
+++ b/src/math/truncf.rs
@@ -10,7 +10,6 @@ pub fn truncf(x: f32) -> f32 {
             return unsafe { ::core::intrinsics::truncf32(x) }
         }
     }
-    let x1p120 = f32::from_bits(0x7b800000); // 0x1p120f === 2 ^ 120
 
     let mut i: u32 = x.to_bits();
     let mut e: i32 = (i >> 23 & 0xff) as i32 - 0x7f + 9;
@@ -26,7 +25,6 @@ pub fn truncf(x: f32) -> f32 {
     if (i & m) == 0 {
         return x;
     }
-    force_eval!(x + x1p120);
     i &= !m;
     f32::from_bits(i)
 }