8 #ifndef COOLPROPFLUID_H_
9 #define COOLPROPFLUID_H_
40 std::vector<double>
cL,
cV;
44 const void get_densities(
double T,
double rho_min,
double rho_crit,
double rho_max,
double& rhoL,
double& rhoV)
const {
45 int Nsoln = -1, Ngood = 0;
46 double rho1 = 0, rho2 = 0, rho3 = 0;
53 if (Nsoln == 1 && rho1 < rho_max && rho1 > rho_crit) {
56 if (rho1 < rho_max && rho1 > rho_crit) {
60 if (rho2 < rho_max && rho2 > rho_crit) {
64 if (Nsoln > 2 && rho3 < rho_max && rho3 > rho_crit) {
69 throw ValueError(
format(
"More than one liquid solution found for critical spline for T=%0.12g", T));
72 throw ValueError(
format(
"No liquid solution found for critical spline for T=%0.12g", T));
82 if (Nsoln == 1 && rho1 > rho_min && rho1 < rho_crit) {
85 if (rho1 > rho_min && rho1 < rho_crit) {
89 if (rho2 > rho_min && rho2 < rho_crit) {
93 if (Nsoln > 2 && rho3 > rho_min && rho3 < rho_crit) {
98 throw ValueError(
format(
"More than one vapor solution found for critical spline for T=%0.12g", T));
101 throw ValueError(
format(
"No vapor solution found for critical spline for T=%0.12g", T));
122 std::vector<CoolPropDbl>
A,
t;
128 std::vector<CoolPropDbl>
A,
B,
n,
m;
160 std::vector<CoolPropDbl>
B,
t,
d;
223 std::vector<CoolPropDbl>
a,
t;
229 std::vector<CoolPropDbl>
a,
t;
233 std::vector<CoolPropDbl>
a,
t;
237 std::vector<CoolPropDbl>
a,
t;
267 std::vector<CoolPropDbl>
b,
t;
271 std::vector<CoolPropDbl>
n,
d,
t;
293 std::vector<CoolPropDbl>
a,
d1,
d2,
t1,
t2,
f,
g,
h,
p,
q,
gamma,
l;
445 assert(R_u < 9 && R_u > 8);
568 #if !defined(NO_FMTLIB) && FMT_VERSION >= 90000
570 return fmt::underlying(type);
574 return fmt::underlying(viscosity);
578 return fmt::underlying(conductivity);