github/CoolProp
1
1
Fork 0
mirror of https://github.com/CoolProp/CoolProp.git synced 2026-04-23 03:00:17 -04:00
Code Issues Packages Projects Releases Wiki Activity

Fix logic for Helmholtz evaluation

Browse Source
This commit is contained in:
Ian Bell
2016-10-17 21:45:04 -06:00
parent 1fd55c9eef
commit 0dbd096a56
1 changed files with 1 additions and 1 deletions
Download Patch File Download Diff File Expand all files Collapse all files

2
src/Helmholtz.cpp
View File

@@ -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;