update

Author: tgross35

crates/libm-test/src/gen/case_list.rs

@@ -432,11 +432,17 @@ fn hypot_cases() -> Vec<TestCase<op::hypot::Routine>> {
     let mut v = vec![];
     TestCase::append_pairs(
         &mut v,
-        &[(
-            // Case that can overflow exponent if wrapping arithmetic is not used
-            (hf64!("-0x1.800f800f80100p+1023"), hf64!("0x1.8354835473720p+996")),
-            Some(hf64!("0x1.800f800f80100p+1023")),
-        )],
+        &[
+            // Cases that can overflow exponent if wrapping arithmetic is not used
+            (
+                (hf64!("-0x1.800f800f80100p+1023"), hf64!("0x1.8354835473720p+996")),
+                Some(hf64!("0x1.800f800f80100p+1023")),
+            ),
+            (
+                (hf64!("0x1.201b201b201c0p+0"), hf64!("0x1.b028b028b02a0p-1")),
+                Some(hf64!("0x1.6821e821e8230p+0")),
+            ),
+        ],
     );
     v
 }

src/math/hypot.rs

@@ -17,7 +17,7 @@
 use super::sqrt;
 #[allow(unused_imports)] // msrv compat
 use super::support::Float;
-use super::support::cold_path;
+use super::support::{Hexf, cold_path};
 
 #[cfg_attr(all(test, assert_no_panic), no_panic::no_panic)]
 pub fn hypot(x: f64, y: f64) -> f64 {
@@ -185,7 +185,7 @@ fn as_hypot_hard(x: f64, y: f64, flag: FExcept) -> f64 {
 
     if ls >= 0 {
         lm <<= ls;
-        m2 += lm.wrapping_mul(lm);
+        m2 = m2.wrapping_add(lm.wrapping_mul(lm));
     } else {
         cold_path();
         let lm2: u128 = (lm as u128) * (lm as u128);