resolve warnings

update

update

Author: tgross35

crates/libm-test/benches/random.rs

@@ -126,6 +126,7 @@ libm_macros::for_each_function! {
         | fdimf128
         | fdimf16
         | floorf128
+        | fmaf16
         | floorf16
         | rintf128
         | rintf16

crates/util/src/main.rs

@@ -96,6 +96,7 @@ fn do_eval(basis: &str, op: &str, inputs: &[&str]) {
             | fdimf16
             | floorf128
             | floorf16
+            | fmaf16
             | rintf128
             | rintf16
             | sqrtf128

src/math/generic/fma.rs

@@ -1,39 +1,34 @@
-use super::super::fenv::{
-    FE_INEXACT, FE_TONEAREST, FE_UNDERFLOW, feclearexcept, fegetround, feraiseexcept, fetestexcept,
-};
-use super::super::{CastFrom, CastInto, DFloat, Float, HFloat, Int, IntTy, MinInt};
+use super::super::fenv::{FE_TONEAREST, fegetround};
+use super::super::{CastFrom, CastInto, DFloat, Float, HFloat, IntTy, MinInt};
 
 /// FMA implementation when a hardware-backed larger float type is available.
 pub fn fma_big<F, B>(x: F, y: F, z: F) -> F
 where
     F: Float + HFloat<D = B>,
     B: Float + DFloat<H = F>,
-    // F: Float + CastInto<B>,
-    // B: Float + CastInto<F> + CastFrom<F>,
     B::Int: CastInto<i32>,
     i32: CastFrom<i32>,
 {
     let one = IntTy::<B>::ONE;
 
-    let xy: B;
-    let mut result: B;
-    let mut ui: B::Int;
-    let e: i32;
-
-    xy = x.widen() * y.widen();
-    result = xy + z.widen();
-    ui = result.to_bits();
-    e = result.exp().signed();
+    let xy: B = x.widen() * y.widen();
+    let result: B = xy + z.widen();
+    let mut ui: B::Int = result.to_bits();
+    let re = result.exp();
     let zb: B = z.widen();
 
     let prec_diff = B::SIG_BITS - F::SIG_BITS;
     let excess_prec = ui & ((one << prec_diff) - one);
-    let x = one << (prec_diff - 1);
-
-    // Common case: the larger precision is fine
-    if excess_prec != x
-        || e == i32::cast_from(F::EXP_MAX)
+    let halfway = one << (prec_diff - 1);
+
+    // Common case: the larger precision is fine if...
+    // This is not a halfway case
+    if excess_prec != halfway
+        // Or the result is NaN
+        || re == B::EXP_MAX
+        // Or the result is exact
         || (result - xy == zb && result - zb == xy)
+        // Or the mode is something other than round to nearest
         || fegetround() != FE_TONEAREST
     {
         // TODO: feclearexcept