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 Coverity CID 38598

Browse Source 
Commented out the (seemingly unused and dead) parts that responsible for
derivative.
This commit is contained in:
mikekaganski
2015-03-07 12:56:12 +10:00
parent 53fb1bd75c
commit bdeb976f89
1 changed files with 9 additions and 12 deletions
Download Patch File Download Diff File Expand all files Collapse all files

21
src/Backends/Helmholtz/FlashRoutines.h
View File

@@ -125,25 +125,22 @@ public:
class solver_TP_resid : public FuncWrapper1D
{
public:
CoolPropDbl T, p, r, peos, rhomolar, rhor, tau, R_u, delta, dalphar_dDelta;
CoolPropDbl T, p, rhor, tau, R_u;
HelmholtzEOSMixtureBackend *HEOS;
solver_TP_resid(HelmholtzEOSMixtureBackend &HEOS, CoolPropDbl T, CoolPropDbl p){
this->HEOS = &HEOS; this->T = T; this->p = p; this->rhor = HEOS.get_reducing_state().rhomolar;
this->tau = HEOS.get_reducing_state().T/T; this->R_u = HEOS.gas_constant();
};
solver_TP_resid(HelmholtzEOSMixtureBackend &HEOS, CoolPropDbl T, CoolPropDbl p):
HEOS(&HEOS),T(T),p(p),rhor(HEOS.get_reducing_state().rhomolar),
tau(HEOS.get_reducing_state().T/T),R_u(HEOS.gas_constant()) {}
double call(double rhomolar){
this->rhomolar = rhomolar;
delta = rhomolar/rhor; // needed for derivative
// CoolPropDbl delta = rhomolar/rhor; // needed for derivative
HEOS->update_DmolarT_direct(rhomolar, T);
peos = HEOS->p();
r = (peos-p)/p;
return r;
CoolPropDbl peos = HEOS->p();
return (peos-p)/p;
};
double deriv(double rhomolar){
/* double deriv(double rhomolar){
// dp/drho|T / pspecified
return R_u*T*(1+2*delta*HEOS->dalphar_dDelta()+pow(delta, 2)*HEOS->d2alphar_dDelta2())/p;
};
};*/
};
/** A residual function for the f(P, Y) solver