Skip to content
Snippets Groups Projects
Commit 80b371bd authored by Matthieu Schaller's avatar Matthieu Schaller
Browse files

Defined additional functions to get hydrodynamical quantities

parent 79575985
Branches
Tags
1 merge request!208Test 125
...@@ -283,3 +283,39 @@ __attribute__((always_inline)) INLINE static float hydro_get_internal_energy( ...@@ -283,3 +283,39 @@ __attribute__((always_inline)) INLINE static float hydro_get_internal_energy(
return entropy * pow_gamma_minus_one(p->rho) * hydro_one_over_gamma_minus_one; return entropy * pow_gamma_minus_one(p->rho) * hydro_one_over_gamma_minus_one;
} }
/**
* @brief Returns the pressure of a particle
*
* @param p The particle of interest
* @param dt Time since the last kick
*/
__attribute__((always_inline)) INLINE static float hydro_get_pressure(
const struct part *restrict p, float dt) {
return p->force.P_over_rho2 * p->rho * p->rho / p->rho_dh;
}
/**
* @brief Returns the entropy of a particle
*
* @param p The particle of interest
* @param dt Time since the last kick
*/
__attribute__((always_inline)) INLINE static float hydro_get_entropy(
const struct part *restrict p, float dt) {
return p->entropy + p->force.entropy_dt * dt;
}
/**
* @brief Returns the sound speed of a particle
*
* @param p The particle of interest
* @param dt Time since the last kick
*/
__attribute__((always_inline)) INLINE static float hydro_get_soundspeed(
const struct part *restrict p, float dt) {
return p->force.soundspeed;
}
...@@ -29,8 +29,8 @@ __attribute__((always_inline)) INLINE static void hydro_debug_particle( ...@@ -29,8 +29,8 @@ __attribute__((always_inline)) INLINE static void hydro_debug_particle(
p->x[0], p->x[1], p->x[2], p->v[0], p->v[1], p->v[2], xp->v_full[0], p->x[0], p->x[1], p->x[2], p->v[0], p->v[1], p->v[2], xp->v_full[0],
xp->v_full[1], xp->v_full[2], p->a_hydro[0], p->a_hydro[1], p->a_hydro[2], xp->v_full[1], xp->v_full[2], p->a_hydro[0], p->a_hydro[1], p->a_hydro[2],
p->h, p->density.wcount, p->density.wcount_dh, p->mass, p->rho_dh, p->rho, p->h, p->density.wcount, p->density.wcount_dh, p->mass, p->rho_dh, p->rho,
p->force.P_over_rho2 * p->rho * p->rho / p->rho_dh, p->force.P_over_rho2, hydro_get_pressure(p, 0.), p->force.P_over_rho2, p->entropy,
p->entropy, p->force.entropy_dt, p->force.soundspeed, p->density.div_v, p->force.entropy_dt, p->force.soundspeed, p->density.div_v,
p->density.rot_v[0], p->density.rot_v[1], p->density.rot_v[2], p->density.rot_v[0], p->density.rot_v[1], p->density.rot_v[2],
p->force.balsara, p->force.v_sig, p->force.h_dt, p->ti_begin, p->ti_end); p->force.balsara, p->force.v_sig, p->force.h_dt, p->ti_begin, p->ti_end);
} }
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment