From cb20cb9f25c9f18cf9891dcfef8795a2a68b6006 Mon Sep 17 00:00:00 2001
From: =?UTF-8?q?Sebastian=20H=C3=BCtter?= <sebastian.huetter@ovgu.de>
Date: Tue, 17 Dec 2019 22:23:19 +0100
Subject: [PATCH] Add range checks for MathSpecial::fm_exp

---
 src/math_special.cpp | 2 ++
 1 file changed, 2 insertions(+)

diff --git a/src/math_special.cpp b/src/math_special.cpp
index bf11a1ad45..d4abc36f25 100644
--- a/src/math_special.cpp
+++ b/src/math_special.cpp
@@ -538,6 +538,8 @@ double MathSpecial::exp2_x86(double x)
 double MathSpecial::fm_exp(double x)
 {
 #if defined(__BYTE_ORDER__) && (__BYTE_ORDER__ == __ORDER_LITTLE_ENDIAN__)
+    if (x < -1022.0/FM_DOUBLE_LOG2OFE) return 0;
+    if (x > 1023.0/FM_DOUBLE_LOG2OFE) return INFINITY;
     return exp2_x86(FM_DOUBLE_LOG2OFE * x);
 #else
     return ::exp(x);