#pragma once namespace teqp { const double N_A = 6.02214076e23; ///< Avogadro's number ///< Gas constant, according to CODATA 2019, in the given number type template<typename NumType> const auto get_R_gas() { const NumType k_B = 1.380649e-23; ///< Boltzmann constant const NumType N_A_ = 6.02214076e23; ///< Avogadro's number return static_cast<NumType>(N_A_ * k_B); }; }; // namespace teqp