mhd.h

Go to the documentation of this file.

#pragma once
 
// wxm includes
#include "apps/mhd/imhd/imhd.h"
 
real get_p(real gas_gamma, real rho, real E, real V2, real B2);
 
real get_E(real gas_gamma, real rho, real p, real V2, real B2);
 
 
void calculate_grad_pressure(const real* q,
                             const real* aux,
                             real* grad_pressure,
                             const real gas_gamma,
                             const real rho_min = std::numeric_limits<real>::epsilon());
 
real get_temperature(real p, real rho, real Ai, real Zi);
 
void calculate_grad_T(const real* q,
                      const real* aux,
                      real* grad_T,
                      const real gas_gamma,
                      const real Ai = 1.0,
                      const real Zi = 1.0,
                      const real rho_min = std::numeric_limits<real>::epsilon(),
                      const real press_min = std::numeric_limits<real>::epsilon());
 
void ensure_floor_density(const real* q, real num_q, real density_floor, real* mod_q);
 
void ensure_floor_pressure(const real* q,
                           real num_q,
                           real gas_gamma,
                           real pressure_floor,
                           real* mod_q);
 
real getRhoFloored(const real* q,
                   const real rho_min = std::numeric_limits<real>::epsilon());
 
real getPressureFloored(const real* q,
                        const real gas_gamma,
                        const real rho_min = std::numeric_limits<real>::epsilon(),
                        const real press_min = std::numeric_limits<real>::epsilon());
 
real getTemperatureFloored(const real* q,
                           const real gas_gamma,
                           const real Ai = 1.0,
                           const real Zi = 1.0,
                           const real rho_min = std::numeric_limits<real>::epsilon(),
                           const real press_min = std::numeric_limits<real>::epsilon());