equation_of_state.h

Go to the documentation of this file.

/*******************************************************************************
 * This file is part of SWIFT.
 * Copyright (c) 2016   Matthieu Schaller (schaller@strw.leidenuniv.nl).
 *
 * This program is free software: you can redistribute it and/or modify
 * it under the terms of the GNU Lesser General Public License as published
 * by the Free Software Foundation, either version 3 of the License, or
 * (at your option) any later version.
 *
 * This program is distributed in the hope that it will be useful,
 * but WITHOUT ANY WARRANTY; without even the implied warranty of
 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
 * GNU General Public License for more details.
 *
 * You should have received a copy of the GNU Lesser General Public License
 * along with this program.  If not, see <http://www.gnu.org/licenses/>.
 *
 ******************************************************************************/
#ifndef SWIFT_IDEAL_GAS_EQUATION_OF_STATE_H
#define SWIFT_IDEAL_GAS_EQUATION_OF_STATE_H
 
/* Some standard headers. */
#include <math.h>
 
/* Local headers. */
#include "adiabatic_index.h"
//#include "common_io.h"
#include "inline.h"
//#include "physical_constants.h"
 
extern struct eos_parameters eos;
 
struct eos_parameters {};
 
__attribute__((always_inline, const)) INLINE static float
gas_internal_energy_from_entropy(float density, float entropy) {
 
  return entropy * pow_gamma_minus_one(density) *
         hydro_one_over_gamma_minus_one;
}
 
__attribute__((always_inline, const)) INLINE static float
gas_pressure_from_entropy(float density, float entropy) {
 
  return entropy * pow_gamma(density);
}
 
__attribute__((always_inline, const)) INLINE static float
gas_entropy_from_pressure(float density, float pressure) {
 
  return pressure * pow_minus_gamma(density);
}
 
__attribute__((always_inline, const)) INLINE static float
gas_soundspeed_from_entropy(float density, float entropy) {
 
  return sqrtf(hydro_gamma * pow_gamma_minus_one(density) * entropy);
}
 
__attribute__((always_inline, const)) INLINE static float
gas_entropy_from_internal_energy(float density, float u) {
 
  return hydro_gamma_minus_one * u * pow_minus_gamma_minus_one(density);
}
 
__attribute__((always_inline, const)) INLINE static float
gas_pressure_from_internal_energy(float density, float u) {
 
  return hydro_gamma_minus_one * u * density;
}
 
__attribute__((always_inline, const)) INLINE static float
gas_internal_energy_from_pressure(float density, float pressure) {
 
  return hydro_one_over_gamma_minus_one * pressure / density;
}
 
__attribute__((always_inline, const)) INLINE static float
gas_soundspeed_from_internal_energy(float density, float u) {
 
  return sqrtf(u * hydro_gamma * hydro_gamma_minus_one);
}
 
__attribute__((always_inline, const)) INLINE static float
gas_soundspeed_from_pressure(float density, float P) {
 
  return sqrtf(hydro_gamma * P / density);
}
 
INLINE static void eos_init(struct eos_parameters *e,
                            const struct phys_const *phys_const,
                            const struct unit_system *us,
                            struct swift_params *params) {}
INLINE static void eos_print(const struct eos_parameters *e) {
 
//  message("Equation of state: Ideal gas.");
 
//  message("Adiabatic index gamma: %f.", hydro_gamma);
}
 
#if defined(HAVE_HDF5)
INLINE static void eos_print_snapshot(hid_t h_grpsph,
                                      const struct eos_parameters *e) {
 
  io_write_attribute_f(h_grpsph, "Adiabatic index", hydro_gamma);
 
  io_write_attribute_s(h_grpsph, "Equation of state", "Ideal gas");
}
#endif
 
#endif /* SWIFT_IDEAL_GAS_EQUATION_OF_STATE_H */