github/CoolProp
1
1
Fork 0
mirror of https://github.com/CoolProp/CoolProp.git synced 2026-01-22 12:28:04 -05:00
Code Issues Packages Projects Releases Wiki Activity

Fix first criticality condition line solver; closes #768

Browse Source
This commit is contained in:
Ian Bell
2015-08-06 10:08:41 -04:00
parent 66153844cc
commit b585bc744b
1 changed files with 6 additions and 2 deletions
Download Patch File Download Diff File Expand all files Collapse all files

8
src/Backends/Helmholtz/HelmholtzEOSMixtureBackend.cpp
View File

@@ -3088,8 +3088,12 @@ std::vector<CoolProp::SimpleState> HelmholtzEOSMixtureBackend::find_all_critical
{
std::string errstr;
OneDimObjective resid_L0(*this, 0.5);
double tau_L0 = Halley(resid_L0, 0.76, 1e-10, 100, errstr);
double tau_L01 = Newton(resid_L0, 0.76, 1e-10, 100, errstr);
double tau0 = 0.66;
resid_L0.call(tau0);
if (resid_L0.deriv(tau0) > 0){
tau0 += 0.1;
}
double tau_L0 = Halley(resid_L0, tau0, 1e-10, 100, errstr);
L0CurveTracer tracer(*this, tau_L0, 0.5);
tracer.trace();