Avoid one more overflow

CalebBell · May 16, 2021 · 3a39cd7 · 3a39cd7
1 parent 87378cc
commit 3a39cd7
Show file tree

Hide file tree

Showing 2 changed files with 6 additions and 1 deletion.
diff --git a/chemicals/dippr.py b/chemicals/dippr.py
@@ -884,7 +884,7 @@ def EQ115(T, A, B, C=0, D=0, E=0, order=0):
        DIPPR/AIChE
     '''
     if order == 0:
-        return exp(A+B/T+C*log(T)+D*T**2 + E/T**2)
+        return trunc_exp(A+B/T+C*log(T)+D*T**2 + E/T**2)
     elif order == 1:
         x0 = T**2
         x1 = 1/x0

diff --git a/tests/test_dippr.py b/tests/test_dippr.py
@@ -163,6 +163,11 @@ def test_EQ115_more():
     d3_numerical = derivative(lambda T: EQ115(T, *args, order=2), 300.0, dx=1e-2)
     assert_close(d3_analytical, 0.0007252940633608988, rtol=1e-13)
     assert_close(d3_analytical, d3_numerical)
+
+    # Check case avoid overflow
+    kwargs = {'T': 400.05, 'A': 1.0, 'B': 1.0, 'C': 1.0, 'D': 1.0, 'E': 1.0, 'order': 0}
+    EQ115(**kwargs)
+
 
 
 def test_EQ116_more():