double gasinfo::heat_capacity_volume(const state &st) const { double Cv = 0; for (int i = 0; i < nc; i++) { double x = st.rho[i] * Rspecific[i] * beta[i]; Cv += x; } return Cv / st.density(); }
double gasinfo::pressure(const state &st) const { return (gamma_ratio(st) - 1) * st.density() * st.specific_energy(); }
double gasinfo::density(const state &st) const { return st.density(); }