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

Fix bug in PhaseEnvelopeRoutines::evaluate; closes #971

Browse Source
This commit is contained in:
Ian Bell
2016-02-11 22:42:02 -07:00
parent cfa9dc1d90
commit 4ff587c1df
1 changed files with 5 additions and 3 deletions
Download Patch File Download Diff File Expand all files Collapse all files

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

@@ -360,7 +360,7 @@ void PhaseEnvelopeRoutines::refine(HelmholtzEOSMixtureBackend &HEOS)
}
double PhaseEnvelopeRoutines::evaluate(const PhaseEnvelopeData &env, parameters output, parameters iInput1, double value1, std::size_t &i)
{
int _i = i;
int _i = static_cast<int>(i);
std::vector<double> const *x, *y;
switch (output){
@@ -388,9 +388,11 @@ double PhaseEnvelopeRoutines::evaluate(const PhaseEnvelopeData &env, parameters
}
if ( _i + 2 >= y->size() ){ _i--; }
if ( _i + 1 >= y->size() ){ _i--; }
if ( _i - 1 < 0 ){ i++; }
if ( _i - 1 < 0 ){ _i++; }
return CubicInterp(*x, *y, _i - 1, _i, _i + 1, _i + 2, inval);
double outval = CubicInterp(*x, *y, _i - 1, _i, _i + 1, _i + 2, inval);
i = static_cast<std::size_t>(_i);
return outval;
}
void PhaseEnvelopeRoutines::finalize(HelmholtzEOSMixtureBackend &HEOS)
{