From 4420289c7c4b2aa3da74ea5c5f4a89a112da1e3a Mon Sep 17 00:00:00 2001
From: gnzlbg <gonzalobg88@gmail.com>
Date: Wed, 10 Jul 2019 12:59:36 +0200
Subject: [PATCH 1/2] Add FMA tests that cause it to segfault

---
 src/math/fma.rs | 19 +++++++++++++++++++
 tests/unit.rs   | 19 +++++++++++++++++++
 2 files changed, 38 insertions(+)
 create mode 100644 tests/unit.rs

diff --git a/src/math/fma.rs b/src/math/fma.rs
index 07d90f8..6d4471a 100644
--- a/src/math/fma.rs
+++ b/src/math/fma.rs
@@ -205,3 +205,22 @@ pub fn fma(x: f64, y: f64, z: f64) -> f64 {
     }
     scalbn(r, e)
 }
+
+#[cfg(test)]
+mod tests {
+    use super::*;
+    #[test]
+    fn fma_segfault() {
+        // These two inputs cause fma to segfault on release due to overflow:
+        assert_eq!(
+            fma(
+                -0.0000000000000002220446049250313,
+                -0.0000000000000002220446049250313,
+                -0.0000000000000002220446049250313
+            ),
+            -0.00000000000000022204460492503126,
+        );
+
+        assert_eq!(fma(-0.992, -0.992, -0.992), -0.00793599999988632,);
+    }
+}
diff --git a/tests/unit.rs b/tests/unit.rs
new file mode 100644
index 0000000..4e70028
--- /dev/null
+++ b/tests/unit.rs
@@ -0,0 +1,19 @@
+use libm::*;
+
+#[test]
+fn fma_segfault() {
+    // These two inputs cause fma to segfault on release due to overflow:
+    assert_eq!(
+        fma(
+            -0.0000000000000002220446049250313,
+            -0.0000000000000002220446049250313,
+            -0.0000000000000002220446049250313
+        ),
+        -0.00000000000000022204460492503126,
+    );
+
+    assert_eq!(
+        fma(-0.992, -0.992, -0.992),
+        -0.00793599999988632,
+    );
+}

From a52ff87b719d7080203a833478e5effc6084012c Mon Sep 17 00:00:00 2001
From: gnzlbg <gonzalobg88@gmail.com>
Date: Wed, 10 Jul 2019 13:01:02 +0200
Subject: [PATCH 2/2] Fix overflow bug in fma

---
 src/math/fma.rs |  4 ++--
 tests/unit.rs   | 19 -------------------
 2 files changed, 2 insertions(+), 21 deletions(-)
 delete mode 100644 tests/unit.rs

diff --git a/src/math/fma.rs b/src/math/fma.rs
index 6d4471a..99f77bc 100644
--- a/src/math/fma.rs
+++ b/src/math/fma.rs
@@ -126,8 +126,8 @@ pub fn fma(x: f64, y: f64, z: f64) -> f64 {
     } else {
         /* r -= z */
         let t = rlo;
-        rlo -= zlo;
-        rhi = rhi - zhi - (t < rlo) as u64;
+        rlo = rlo.wrapping_sub(zlo);
+        rhi = rhi.wrapping_sub(zhi.wrapping_sub((t < rlo) as u64));
         if (rhi >> 63) != 0 {
             rlo = (-(rlo as i64)) as u64;
             rhi = (-(rhi as i64)) as u64 - (rlo != 0) as u64;
diff --git a/tests/unit.rs b/tests/unit.rs
deleted file mode 100644
index 4e70028..0000000
--- a/tests/unit.rs
+++ /dev/null
@@ -1,19 +0,0 @@
-use libm::*;
-
-#[test]
-fn fma_segfault() {
-    // These two inputs cause fma to segfault on release due to overflow:
-    assert_eq!(
-        fma(
-            -0.0000000000000002220446049250313,
-            -0.0000000000000002220446049250313,
-            -0.0000000000000002220446049250313
-        ),
-        -0.00000000000000022204460492503126,
-    );
-
-    assert_eq!(
-        fma(-0.992, -0.992, -0.992),
-        -0.00793599999988632,
-    );
-}