From 0dbd096a56e97c77c4083c0500957296bcbebae2 Mon Sep 17 00:00:00 2001
From: Ian Bell <ian.h.bell@gmail.com>
Date: Mon, 17 Oct 2016 21:45:04 -0600
Subject: [PATCH] Fix logic for Helmholtz evaluation

---
 src/Helmholtz.cpp | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/src/Helmholtz.cpp b/src/Helmholtz.cpp
index f5b13e3b..95a8b23f 100644
--- a/src/Helmholtz.cpp
+++ b/src/Helmholtz.cpp
@@ -161,7 +161,7 @@ void ResidualHelmholtzGeneralizedExponential::all(const CoolPropDbl &tau, const
         
         if (delta_li_in_u){
             CoolPropDbl  ci = el.c, l_double = el.l_double;
-            if (ValidNumber(l_double) && ci > 0){
+            if (ValidNumber(l_double) && l_double > 0 && std::abs(ci) > DBL_EPSILON){
                 const CoolPropDbl u_increment = (el.l_is_int) ? -ci*powInt(delta, el.l_int) : -ci*pow(delta, l_double);
                 const CoolPropDbl du_ddelta_increment = l_double*u_increment*one_over_delta;
                 const CoolPropDbl d2u_ddelta2_increment = (l_double-1)*du_ddelta_increment*one_over_delta;