From a288d5c6f7a85f0e6c0c86bb9d5417c9e7a3eac5 Mon Sep 17 00:00:00 2001
From: Ray Speth <speth@mit.edu>
Date: Thu, 14 Apr 2016 18:55:33 -0400
Subject: [PATCH] [1D] Make Jacobian atol independent of compiler optimizations

If compiling with -funsafe-math-optimizations, this method of attempting to
calculate the floating point epsilon value would fail.
---
 src/oneD/MultiJac.cpp | 6 +-----
 1 file changed, 1 insertion(+), 5 deletions(-)

diff --git a/src/oneD/MultiJac.cpp b/src/oneD/MultiJac.cpp
index dd2ad2317..9775bdd50 100644
--- a/src/oneD/MultiJac.cpp
+++ b/src/oneD/MultiJac.cpp
@@ -24,11 +24,7 @@ MultiJac::MultiJac(OneDim& r)
     m_elapsed = 0.0;
     m_nevals = 0;
     m_age = 100000;
-    doublereal ff = 1.0;
-    while (1.0 + ff != 1.0) {
-        ff *= 0.5;
-    }
-    m_atol = sqrt(ff);
+    m_atol = sqrt(std::numeric_limits<double>::epsilon());
     m_rtol = 1.0e-5;
 }