Method PropertyPlot

src/CoolPropPlot.cpp:273–308

}

PropertyPlot::PropertyPlot(const std::string& fluid_name, CoolProp::parameters ykey, CoolProp::parameters xkey, TPLimits tp_limits)
  : xkey_(xkey), ykey_(ykey) {
    this->state_ = Detail::process_fluid_state(fluid_name);
    this->critical_state_ = Detail::get_critical_point(state_);

    xaxis.scale = Detail::default_scale(xkey);
    yaxis.scale = Detail::default_scale(ykey);

    // We are just assuming that all inputs and outputs are in SI units. We
    // take care of any conversions before calling the library and after
    // getting the results.
    int out1 = 0, out2 = 0;
    axis_pair_ = CoolProp::generate_update_pair(xkey, 0, ykey, 1, out1, out2);
    swap_axis_inputs_for_update_ = (out1 == 1);

    const double HI_FACTOR = 2.25;  // Upper default limits: HI_FACTOR*T_crit and HI_FACTOR*p_crit
    const double LO_FACTOR = 1.01;  // Lower default limits: LO_FACTOR*T_triple and LO_FACTOR*p_triple
    switch (tp_limits) {
        case TPLimits::None:
            this->Tp_limits_ = {{Detail::NaN, Detail::NaN}, {Detail::NaN, Detail::NaN}};
            break;
        case TPLimits::Def:
            this->Tp_limits_ = {{LO_FACTOR, HI_FACTOR}, {LO_FACTOR, HI_FACTOR}};
            break;
        case TPLimits::Achp:
            this->Tp_limits_ = {{173.15, 493.15}, {0.25e5, HI_FACTOR}};
            break;
        case TPLimits::Orc:
            this->Tp_limits_ = {{273.15, 673.15}, {0.25e5, HI_FACTOR}};
            break;
    }

    Range2D ranges = get_axis_limits();
    xaxis.range = ranges.x;
    yaxis.range = ranges.y;
}

Range PropertyPlot::isoline_range(CoolProp::parameters key) const {
    if (key == CoolProp::iQ)