euler1d_riemann_problem_analytic_solution.h

Go to the documentation of this file.

#ifndef WXM_APPS_5MOMENT_EULER1D_RIEMANN_PROBLEM_ANALYTIC_SOLUTION_H
#define WXM_APPS_5MOMENT_EULER1D_RIEMANN_PROBLEM_ANALYTIC_SOLUTION_H
 
// General includes
#include "apps/wmapplication.h"
 
// Specific includes
// #include "apps/vlasovmaxwell/vlasovmaxwell.h"
 
namespace wxm
{
namespace apps
{
namespace five_moment
{
 
class Euler1D_Riemann_Problem_Analytic_Solution : public WmApplication
{
public:
    Euler1D_Riemann_Problem_Analytic_Solution();
 
    ~Euler1D_Riemann_Problem_Analytic_Solution() override;
 
    void setup(const WxCryptSet& wxc) override;
 
    // const std::vector<int>& getInputVariableIndexes(int flag) const override
    // {
    //  return _variables;
    // }
 
    const std::vector<int>& getOutputVariableIndexes(int flag) const override
    {
        return _variables;
    }
 
    // const std::vector<int>& getAuxiliaryVariableIndexes(int flag) const override
    // {
    //     return _auxVariables;
    // }
 
    void evaluate_function(const real* q,
                           const real* aux,
                           const solverVariables_t* pSV,
                           real* result) const override;
 
protected:
    real _gas_gamma; // gamma
 
    real _xc; // x center location (where jump is)
 
    std::vector<real> _w_left; // input left state primitive variables, (rho_left, u_left,
                               // v_left=0, w_left=0, press_left)
    std::vector<real> _w_right; // input right state primitive variables, (rho_right,
                                // u_right, v_right=0, w_right=0, press_right)
 
    std::vector<real> _q_left;  // input left conserved variables, (rho_left, px_left,
                                // py_left=0, pz_left=0,  ener_left)
    std::vector<real> _q_right; // input right conserved variables, rho_right, px_right,
                                // py_right=0, pz_right=0, ener_right)
 
    std::vector<real> _wave_speeds;           // calculated wave speeds
    std::vector<std::vector<real>> _q_states; // calculated conserved variable states
    std::vector<std::vector<real>> _w_states; // calcualtes primitive variable states
 
    std::vector<int> _variables; // conserved variable vector to output to
 
private:
    Euler1D_Riemann_Problem_Analytic_Solution&
    operator=(const Euler1D_Riemann_Problem_Analytic_Solution& var);
    Euler1D_Riemann_Problem_Analytic_Solution(
        const Euler1D_Riemann_Problem_Analytic_Solution& var);
};
 
} // namespace five_moment
} // namespace apps
} // namespace wxm
 
#endif // WXM_APPS_5MOMENT_EULER1D_RIEMANN_PROBLEM_ANALYTIC_SOLUTION_H