MeltingCaloric.h

Go to the documentation of this file.

#ifndef MELTING_CALORIC_H
#define MELTING_CALORIC_H
 
#include <cmath>
#include <cstddef>
#include <limits>
#include <memory>
#include <optional>
#include <utility>
#include <vector>
#include "CoolProp/superancillary/superancillary.h"
#include <Eigen/Dense>
 
namespace CoolProp {
 
class HelmholtzEOSMixtureBackend;  // forward decl
 
class MeltingCaloric
{
   public:
    MeltingCaloric() = default;
 
    void sample(HelmholtzEOSMixtureBackend& H, std::size_t n);
 
    [[nodiscard]] std::size_t n_samples() const {
        return m_lnp.size();
    }
    [[nodiscard]] double sample_lnp(std::size_t i) const {
        return m_lnp[i];
    }
    [[nodiscard]] double sample_T(std::size_t i) const {
        return m_T[i];
    }
    [[nodiscard]] double sample_rho(std::size_t i) const {
        return m_rho[i];
    }
    [[nodiscard]] double sample_h(std::size_t i) const {
        return m_h[i];
    }
    [[nodiscard]] double sample_s(std::size_t i) const {
        return m_s[i];
    }
 
    void build(HelmholtzEOSMixtureBackend& H);
    [[nodiscard]] bool built() const {
        return m_built;
    }
 
    [[nodiscard]] double lnp_min() const {
        return m_T_approx ? m_T_approx->xmin() : 0.0;
    }
    [[nodiscard]] double lnp_max() const {
        return m_T_approx ? m_T_approx->xmax() : 0.0;
    }
 
    [[nodiscard]] double eval_T(double lnp) const {
        return m_T_approx ? m_T_approx->eval(lnp) : std::numeric_limits<double>::quiet_NaN();
    }
    [[nodiscard]] double eval_rho(double lnp) const {
        return m_rho_approx ? m_rho_approx->eval(lnp) : std::numeric_limits<double>::quiet_NaN();
    }
    [[nodiscard]] double eval_h(double lnp) const {
        return m_h_approx ? m_h_approx->eval(lnp) : std::numeric_limits<double>::quiet_NaN();
    }
    [[nodiscard]] double eval_s(double lnp) const {
        return m_s_approx ? m_s_approx->eval(lnp) : std::numeric_limits<double>::quiet_NaN();
    }
 
    bool seed_for_hs(double s_cache, double h_cache, double& T0, double& rho0) const;
 
    [[nodiscard]] std::optional<std::pair<double, double>> stamp() const {
        return m_stamp;
    }
 
    [[nodiscard]] double curve_Tmin() const {
        return m_curve_Tmin;
    }
 
   protected:
    std::vector<double> m_lnp, m_T, m_rho, m_h, m_s;
 
    using Approx = CoolProp::superancillary::ChebyshevApproximation1D<Eigen::ArrayXd>;
    bool m_built = false;
    std::optional<Approx> m_T_approx, m_rho_approx, m_h_approx, m_s_approx;
    std::optional<std::pair<double, double>> m_stamp;
    double m_curve_Tmin = 0.0;
};
 
std::shared_ptr<MeltingCaloric> get_melting_caloric_cached(HelmholtzEOSMixtureBackend& H);
 
}  // namespace CoolProp
#endif