Method build_saturation_shim

src/Backends/REFPROP/REFPROPMixtureBackend.cpp:2386–2407

}

shared_ptr<REFPROPMixtureBackend> REFPROPMixtureBackend::build_saturation_shim(int Q) {
    this->check_loaded_fluid();
    if (!(Q == 0 || Q == 1)) {
        throw ValueError(format("build_saturation_shim requires Q==0 or Q==1; got %d", Q));
    }
    const std::vector<double>& x_phase = (Q == 0) ? mole_fractions_liq : mole_fractions_vap;
    // is_cached (operator bool) check first: operator double()/CoolPropDbl() throws if uncached.
    bool have = (Q == 0) ? static_cast<bool>(_rhoLmolar) : static_cast<bool>(_rhoVmolar);
    CoolPropDbl rho_phase = have ? ((Q == 0) ? static_cast<CoolPropDbl>(_rhoLmolar) : static_cast<CoolPropDbl>(_rhoVmolar)) : _HUGE;
    if (!have || !ValidNumber(rho_phase)) {
        throw ValueError("The saturated state has not been set (no two-phase density available).");
    }
    shared_ptr<REFPROPMixtureBackend> shim(new REFPROPMixtureBackend());
    shim->link_to_loaded_fluids(*this);
    // Copy only the first Ncomp entries of the per-phase composition.
    std::vector<CoolPropDbl> x(x_phase.begin(), x_phase.begin() + Ncomp);
    shim->set_mole_fractions(x);
    shim->specify_phase((Q == 0) ? iphase_liquid : iphase_gas);
    shim->update_DmolarT_direct(rho_phase, _T);
    shim->_Q = static_cast<CoolPropDbl>(Q);  // mark the shim as sitting on the envelope
    return shim;
}

double REFPROPMixtureBackend::dpdT_along_saturation_pure(int kph) {
    if (Ncomp != 1) {  // precondition: DPTSATKdll is a single-component saturation routine